Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451_phaM1 1,499,508 +C intergenic (+85/‑79) trpS CDS (Tryptophan‑‑tRNA ligase) → / → zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain) –/–

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM11,499,5061.C92.8% 903.6 / 77.2 308intergenic (+83/‑81)trpS CDS (Tryptophan‑‑tRNA ligase)/zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain)–/–
Reads supporting (aligned to +/- strand):  ref base . (18/4);  new base C (171/111);  total (193/115)
Fisher's exact test for biased strand distribution p-value = 6.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 325 total aligned reads displayed.

GATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCAC‑GGCTTTGCCGGTG‑TT‑T‑CGGG‑CAC‑G‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTC  >  AE015451_phaM1/1499373‑1499637
                                                                                                                                            |                                                                                                                                   
gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGTCCGCGAAACAAGAGCCGCGGTGCATGGCACC‑GGTTTTTCCGGTG‑TT‑TACGGG‑CAc‑g‑cc                                                                                                                                   >  1:6044929/1‑137 (MQ=255)
gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTTCCCACTAAACAAGCACCGCATTGCATAGCACC‑GACCTTGCCGGTG‑TTCG‑CGGG‑CAc‑g‑cc                                                                                                                                   >  1:15383535/1‑137 (MQ=25)
gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTG‑TTCG‑CGGG‑CAc‑a‑cc                                                                                                                                   >  2:19005410/1‑137 (MQ=25)
gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTACCGGTG‑TTCG‑CGGG‑CAc‑g‑cc                                                                                                                                   >  2:19349402/1‑137 (MQ=25)
gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCACGGTGCATGGCACC‑GGCTTTGCCGGTG‑TTCG‑CGGG‑CAc‑a‑cc                                                                                                                                   >  2:13168356/1‑137 (MQ=25)
gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACAGGCGCCACGGTGCACGGCCCC‑GGCTTTTCCAGTG‑TTCT‑CGGG‑TAc‑g‑cc                                                                                                                                   >  2:2472108/1‑137 (MQ=38)
gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGGCCGCGAAAAAAGAGCCACAGTGCCCGGCACC‑GAATTTTCCAGTGTTT‑T‑CGGG‑CAc‑a‑cc                                                                                                                                   >  1:7548722/1‑137 (MQ=255)
 aTGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCGCGGTGCATGGCCCC‑GGCTTTGCCGGTG‑TTCG‑CGGG‑CAc‑g‑ccc                                                                                                                                  >  2:3069346/1‑136 (MQ=25)
       ccATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCAAAAAAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTccc                                                                                                                             >  2:13985222/1‑137 (MQ=21)
         aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGTCTGGCACC‑GGGTTTGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCACTCCcac                                                                                                                           >  2:1792298/1‑137 (MQ=21)
         aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGAGGTGCATGGCACC‑GGCTTTGCCAGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCcac                                                                                                                           >  2:5956804/1‑137 (MQ=21)
         aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCAAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTG‑TTCG‑CGGG‑CAC‑GTCCGCTCCcac                                                                                                                           >  2:6933881/1‑137 (MQ=21)
             agagCCGAAACGAAGCTTTGGCACCACAGCAGGACGAAGGCTGCCAAACGGGGGGAGCGGTTGTACCCCCGAAAAAGGCGCCGCGGGGCATGACACC‑GTCTTTGCTGGGG‑TTCG‑AGGA‑CAC‑ACCCGCTCCCACAGAc                                                                                                                       <  2:18096541/137‑1 (MQ=255)
                   tGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTG‑TTCG‑CGGG‑CAA‑GCCCGCTCCCACAGACTGATTg                                                                                                                 >  2:16317676/1‑137 (MQ=18)
                    gAACGACGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACAGGCGCCACGGTGCACGGCACC‑GGATTTGCCGGTG‑TTCG‑CGGG‑AAC‑GCCCGCTCCCACAGAATGATTGc                                                                                                                >  2:12003050/1‑137 (MQ=21)
                    gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCTCTCCCACAGACTGATTGc                                                                                                                >  2:12420983/1‑137 (MQ=21)
                       cGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCCC‑GGCTTTGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGccct                                                                                                             >  2:15877932/1‑135 (MQ=21)
                       cGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGGCTTGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCac                                                                                                             >  2:1661530/1‑136 (MQ=21)
                       cGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAAACTGATTGCCAt                                                                                                             >  2:13814346/1‑137 (MQ=18)
                       cGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAAGCGCCACGGTGCCTGGCACC‑GGCTTTGTCGGTG‑TTCT‑CGGG‑CCC‑GCCCGCTCCCACAGACTGATTGCCAt                                                                                                             >  2:10207892/1‑137 (MQ=21)
                        gAAGCTTTGGCACCACAGCAGGATTGGGGCGGCAAAACGGTGGGAACGGGCGTTCCCGCGAAGCAGTCGCCGTGGTGCGTGGCACC‑GGCTTTGCTGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCAtt                                                                                                            <  2:9231633/137‑1 (MQ=12)
                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGTCCGCGAAACAAGAGCCGCAGTGCCTGGCACC‑GGCTTTTCCGGTGTTT‑T‑TGGG‑CAC‑GCCCGCTCCCCCAAAACGAGTGCCact                                                                                                            >  1:13031249/1‑135 (MQ=255)
                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGTCCGCGAAAAAAGCGCCGTGGTGCCCGGCACC‑TGCTTTGCCGGTG‑TT‑TGCGGG‑CAC‑GCCCGCTCCCACCAACCGAGTGCCAtt                                                                                                            >  1:17826374/1‑137 (MQ=19)
                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGTCCGCGAAAAAGGAGCCGCGGTGCATGGTACC‑GGCTTTGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGAATGGTTGCCact                                                                                                            >  1:8670328/1‑135 (MQ=21)
                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCAtt                                                                                                            >  2:9291196/1‑137 (MQ=21)
                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCACGGTGCATGGCCCC‑GGCTTTGCCAGTG‑TTCG‑CAGG‑CAC‑GCCCTCCCCCACAGACCGATTGCCact                                                                                                            >  2:5698429/1‑135 (MQ=21)
                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCACGGTGCACGGCACC‑GGCTTTTCCGGTG‑TTCG‑CGGG‑CAC‑GCCCACTCCCACAAACTGATTGCCCtt                                                                                                            >  2:16740160/1‑137 (MQ=18)
                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAAAGGTGGGAGCGGGCGTGTCCGCGAAACAGACGCTGAGATGTCTGTCACC‑GGCTTTGCCGGTG‑TTCA‑CGGG‑TAC‑GCCCGCTCCCAAAGACCGATTGCCAtt                                                                                                            >  2:10588537/1‑137 (MQ=19)
                              ttGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGTCGCGGTGCATGGCACC‑GGCTTTGCCAGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGAATGATTGCCACTAAGCgg                                                                                                      >  2:15292133/1‑137 (MQ=18)
                              ttGGCACCACAGAAGGATTGGGGTGGCAAAATGGTGGGAGTGGGTGTTCCCGAGAAACAGGCGCTG‑GGGGGAGGGGCACGGGTTTTGCTGGTG‑TT‑TGGGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCgg                                                                                                      <  1:1105792/137‑1 (MQ=255)
                                 gCACCACAGCAGGACGGAGGCTGCAAAACGGGGGGAGCGGGTGTGCCCGCGAAACGGGCGCCGTGTTGCGTGGCACC‑GGTTTTTCCGGTG‑TT‑TGCGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                                   <  2:8721340/137‑1 (MQ=21)
                                 gCACCACAGAAGGATGGAGGCTGCAAAATGGGGGGAGCGGGCGTGCCCGAGAAAAAGGCACAGTGGGGCGGGAAACC‑GGTTTTGCCGGTG‑TT‑C‑GGGGACAC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                                   <  1:8793533/137‑1 (MQ=255)
                                 gCACCAAAGCAGGACTGAGGCTGCCAAACGGTGGGGACGGGCGTGCCCGCGAAACATGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTA‑TTCG‑TGGA‑CAC‑GCCCGCTCCCACAGACTGATTACCATTAAGCGAGcc                                                                                                   <  2:18855286/137‑1 (MQ=17)
                                  caccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGTCCGCGAAACAGGAGACACGGTGCATGGCACC‑AGCTTTGCCGGTG‑TT‑TGCGGG‑TACAC‑CCTCTCCCCCAGACTGACTGCCCTTCCGCGGGCCt                                                                                                  >  1:7326042/1‑137 (MQ=255)
                                  caccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGTATGGCACC‑GGCTTTGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGAACGATTGCCATTAAGCGGGCCt                                                                                                  >  2:16908586/1‑137 (MQ=18)
                                  aaccacAGAAGGACGGAGGCTACAAAATGGGGGGAGCGGGTGTGCCTGCGAAACGGGCGCCGCGGTGCATGGAACC‑GGTTTTGCTGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCt                                                                                                  <  2:7794666/136‑1 (MQ=16)
                                      acaGCAGGACTGAGGCTACTAAATGGTGGGAGCGGGCGTTCCCTCGAAACAGGCCCCGCGGTACGTGGCACC‑GGCTTTGCCGGTG‑TTCG‑TGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                                              <  2:4428047/137‑1 (MQ=18)
                                      acaGCAGGACTGAGGCTACCAAACGGTGGGAGTGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGTTTGGCCGGGT‑TT‑AGCGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                                              <  2:17594403/137‑1 (MQ=18)
                                      acaGCAGGACTAAGGCTGCCAAACGGTGGGAATGGGCGTTCCCGCAAAACATGCGCCGCGGTACATGACACC‑GGCTTTGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                                              <  2:7616976/137‑1 (MQ=21)
                                      acaGCAGAACGGAGGTGGCAAAAGGGGGGGGGTGGGTGTGCCCGAAAAACAGGCACCGGGGGGCATGGCACC‑GGTTTTGCTGGTG‑TC‑CGAGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                                              <  1:18358297/137‑1 (MQ=12)
                                        agcagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCACGGTGCATGGCACC‑GGCCTTGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCACTCCCACAGACTGAGTGCCATTCAAAGGGACTCTCGCt                                                                                            >  2:3263343/1‑137 (MQ=17)
                                             gACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCACGAAAAAAGCGCCGCGGTGCATGGCCCC‑GGCCTTTCCGGTG‑TTCG‑CGGG‑CAC‑ACCCACCCCCACAAACCGATTTCCCTTCAACAGGCCTCTCGCTAccat                                                                                       >  1:9433310/1‑134 (MQ=255)
                                                ggAGGCTGCCAAACGGGGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGACACC‑GGCTTTGCCGTTG‑TT‑CGCGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAc                                                                                    <  2:14410429/136‑1 (MQ=18)
                                                     cTGCCAAACGGTGGGAGCGGGCGTTCCCGCGAAACAGGGGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTG‑TTCG‑CGGG‑CACAC‑CCGATCCCACAGGATGATTGTCATTAAGTGGGCCTATCGCTGAAGTGGCGGtcc                                                                               >  2:5326698/1‑137 (MQ=9)
                                                     cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACCAGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTCAGCGGGCCTATCGCTCCAGTGACGGccc                                                                               >  2:1500370/1‑137 (MQ=17)
                                                     cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCGTTGCCGGTG‑TTCA‑CGGG‑CAA‑GCCCGCTCCCACAAACTGATTGCCATTCAGCGGGCCTATCGCTCAAGTGACGGccc                                                                               >  2:9024319/1‑137 (MQ=18)
                                                       gCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGCACC‑GGCTTTGCCGGTG‑TT‑CGCGGG‑CAC‑GACCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTCCAGTGACGGCCCGt                                                                             >  2:18547803/1‑137 (MQ=21)
                                                          aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCAGTGCATGGCCCC‑GGCTTTGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGACTGCCATTAAGAAGGACTATTGCTAAAGTGACGGTCCGtttt                                                                          >  2:9729314/1‑137 (MQ=17)
                                                           aaCGGTGGGAGCGGGCGTGTCCGCGGAACAGGCGCCGCGGTGCATGGCACC‑GGATTTGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGAGGGCCTATCGCTACAATGACGGCCCGctttt                                                                         >  2:11547479/1‑137 (MQ=18)
                                                           aaCGGTGGGAGCGGGCGTGTCCACGAAACAGGCGCCGCGGTGCATGGCACC‑GGCCTTGCCGGTG‑TTCG‑CAGG‑CACAC‑CCACTCCCACAAACCGATTGCCACTCAACAGGCCCATCACTCCAGTGACAGCCCAttttt                                                                         >  2:6589226/1‑137 (MQ=9)
                                                               gTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTG‑TT‑CGCGGG‑CAC‑GCCCGCTCCCACAGACTGATTTCCATTAAGCGGGCCTCTCGCTACAATGACGGCCCGGTTTttgtt                                                                     >  2:16147671/1‑137 (MQ=18)
                                                                ggggAGCGGGCGTGCCCGCGAAACAGGCGCCGTGGGGAATGGAACC‑GTCCTTGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgttg                                                                    <  2:9127012/136‑1 (MQ=21)
                                                                  ggAGCGGGCGGGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTTCCGGTG‑TTCG‑CGGG‑CAC‑TCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGcc                                                                  <  2:9485577/137‑1 (MQ=21)
                                                                     gCGGGTGTGTCCGCGAAATAGGCGCCGCGGTGCATGGAACC‑GGCTTTGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGAATGATTGCCATTAAGCGGGACTATCGCTACAGTGACGGCCCGTTTTTTGTTGCATCg                                                               >  2:19213887/1‑137 (MQ=21)
                                                                     gCGGGCGTGCCCACGAAAAAAGCGCCACGGTGCATGGCACC‑AGCTTTTCCGGTG‑TTCG‑CGGG‑CAC‑ACCCGCTCCCACAGAATGAGTGCCATTAAGCGGGCCTATCGCTACAATGACGGACCATTTTTTGTTGCCCCg                                                               >  2:16086950/1‑137 (MQ=255)
                                                                      tgtttttGCCCGTGAAACAGGCCCCGTGGGGCGCGACACC‑GGCTTTGCCGGTG‑TT‑CGAGGA‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTAGCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGc                                                              <  2:10692027/131‑1 (MQ=19)
                                                                      cGGGCTTGCCCGCGAAACAGGCGCCGCGGGGCATGACACC‑GGCTTGGCTGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGc                                                              <  1:14500790/137‑1 (MQ=21)
                                                                         gtgTGCCTTCGAAACGTGCGCCGTGGTGCGTGGCACT‑GTCTTTGCTGGTG‑TT‑CGTGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTaa                                                           <  1:2380521/135‑1 (MQ=19)
                                                                                cgtgAAACAGGCGCTGTGGTGCATGGCACC‑GTCTTTGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    <  1:14066113/134‑1 (MQ=25)
                                                                                cgcgAAACAGGCGCTACGGTGCATGGCACC‑GGCTTTGCCGGTG‑TTCG‑TGGA‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACNGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    <  1:18785082/137‑1 (MQ=25)
                                                                                cgcgAAACAGGCGCCGCGGTGCATGGCACT‑GGTTTTGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    <  2:14559781/137‑1 (MQ=25)
                                                                                cgcaAAACAGGCGCCACGGTGCATGGCACC‑GGCTTTGCCGGTG‑TTCG‑TGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    <  1:13625716/137‑1 (MQ=25)
                                                                                cgagAAAAAGGCGTCGTGGTGCAGGGCACT‑GGTTTTGCTGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    <  1:16714405/134‑1 (MQ=25)
                                                                                agagAAACAGTCGCCGCGGTGCATGGCACT‑GGCTTTGCAGGTG‑TT‑CGCGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    <  2:7872206/134‑1 (MQ=25)
                                                                                 tcgAAAAAGGCTCCGCGGTGCGGGGAACC‑GGTTTTTCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                                   <  2:4996916/136‑1 (MQ=25)
                                                                                 gggAAACGGTCGCCGTGGTGCATGGCACT‑GTCTTTGCTGTTT‑TTCG‑CGGA‑CAC‑TCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                                   <  1:13475424/135‑1 (MQ=255)
                                                                                 gcgAAACAGGCGCCGCGGTGCAGGGCACC‑GGCTTTTCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                                   <  2:11731394/137‑1 (MQ=37)
                                                                                 gcaaaaCAGGCCCCGCGGTACAGGGCACC‑GGCTTTCCCGGTT‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                                   <  2:5654137/134‑1 (MQ=25)
                                                                                 gagAAACGGGCGCCGTGGGGCGTGGCACC‑GGTTTGGCTGTTG‑TT‑TGGGGG‑CAC‑ACCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                                   <  1:5203727/135‑1 (MQ=255)
                                                                                   gAAACAGGCGCCGCGGGGCATGGCACC‑GGTTTTGCTGGTG‑TTCG‑TGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGa                                                 <  1:17650694/137‑1 (MQ=37)
                                                                                     aaCAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    >  1:540212/1‑132 (MQ=37)
                                                                                     aaCAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACt                                               >  2:13618278/1‑137 (MQ=37)
                                                                                      acggtcgccgcGGTGCATTGCACC‑GGTTTTGCTGGGG‑TTCG‑TGGG‑CAC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTc                                              <  1:19440212/132‑1 (MQ=255)
                                                                                           cgccgtGGTACGTGGCACC‑GGCTTTTCCGGTT‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                         <  2:13136217/137‑1 (MQ=25)
                                                                                           cgccgtGGGGCATGGCACT‑GGTTTTGCTGGTG‑TT‑C‑GGGGACAC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                         <  1:1905471/137‑1 (MQ=255)
                                                                                             ctgtggTGCGTGGCCCT‑GGTTTTTCCGGTG‑TT‑TGCGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGa                                       <  1:3389339/133‑1 (MQ=255)
                                                                                             ctgcgGGGCAGGGCACC‑GGTTTTGCCGGTG‑TTCG‑AGGG‑CAC‑TCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGa                                       <  1:16125544/135‑1 (MQ=255)
                                                                                             ccgcGGTGCATGGCACC‑GGCTTTGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGa                                       >  1:10129997/1‑137 (MQ=37)
                                                                                                     cATGGCACC‑GGCTTTGCCGGTG‑TTCG‑TGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                               <  2:12333394/137‑1 (MQ=37)
                                                                                                     cATGGCACC‑GGCTTTGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                               <  2:12693782/137‑1 (MQ=37)
                                                                                                     cATGGCACC‑GGCTTTGCCGGTG‑TT‑CGCGGG‑CAC‑TCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                               <  2:11666021/137‑1 (MQ=37)
                                                                                                        ggcacc‑GGCTTTGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg                            >  1:3699995/6‑137 (MQ=37)
                                                                                                        ggcacc‑GGCTTTGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg                            >  2:11200355/6‑137 (MQ=37)
                                                                                                           acc‑GGCTTTGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAg                         >  2:4722300/3‑137 (MQ=38)
                                                                                                           acc‑GGCTTTGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAg                         >  2:17053346/3‑137 (MQ=38)
                                                                                                             c‑GGCTTTGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGAt                       >  1:7617038/1‑137 (MQ=38)
                                                                                                               ggCTTTGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATc                      >  1:12129885/1‑137 (MQ=38)
                                                                                                                 cTTTGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAATGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  2:8436614/1‑137 (MQ=38)
                                                                                                                 cTTTGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCGTGACTCCCTTAGGACGCTATCAAGCAGATCTg                    >  2:14887565/1‑137 (MQ=38)
                                                                                                                 cTTTGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCCAGCAGATCTg                    >  1:14072510/1‑137 (MQ=38)
                                                                                                                 cTTTGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  1:1463526/1‑137 (MQ=39)
                                                                                                                 cTTTGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  1:14975006/1‑137 (MQ=39)
                                                                                                                 cTTTGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  1:12739307/1‑137 (MQ=39)
                                                                                                                   ttGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa                  >  2:1342042/1‑137 (MQ=39)
                                                                                                                   ttGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa                  >  1:6894022/1‑137 (MQ=39)
                                                                                                                   ttGCCGGTG‑TTCG‑CGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa                  >  1:18545355/1‑137 (MQ=39)
                                                                                                                               tcgcGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAc        >  2:1711713/4‑137 (MQ=255)
                                                                                                                               tcgcGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAc        >  1:10605428/4‑137 (MQ=255)
                                                                                                                                cgcGGG‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAct       >  2:2107025/3‑137 (MQ=255)
                                                                                                                                   ggg‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActtct    >  1:3929480/1‑137 (MQ=255)
                                                                                                                                            |                                                                                                                                   
GATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCAC‑GGCTTTGCCGGTG‑TT‑T‑CGGG‑CAC‑G‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTC  >  AE015451_phaM1/1499373‑1499637

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑