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.C94.7% 743.4 / 40.3 229intergenic (+83/‑81)trpS CDS (Tryptophan‑‑tRNA ligase)/zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain)–/–
Reads supporting (aligned to +/- strand):  ref base . (11/1);  new base C (147/69);  total (159/70)
Fisher's exact test for biased strand distribution p-value = 1.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 243 total aligned reads displayed.

GATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTT‑T‑CGGG‑C‑A‑C‑G‑C‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTC  >  AE015451_phaM1/1499373‑1499637
                                                                                                                                            |                                                                                                                                    
gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGCCCCCGAAACAAGAGCCACTGTGCACGGCACCAGACTTTTCGGTGTT‑TGCGGG‑A‑A‑c‑g‑c‑c                                                                                                                                   >  1:10365934/1‑137 (MQ=255)
gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGTCTGGCACCAGCTTTTTCGGTGTTCG‑CGGG‑C‑A‑c‑g‑c‑c                                                                                                                                   >  2:2199249/1‑137 (MQ=25)
gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCAAAACAGGGGCCGCGGTGCATGGCACCGGCTTTGTCGGTGTTCG‑CAGG‑C‑A‑c‑g‑c‑c                                                                                                                                   >  2:5692553/1‑137 (MQ=25)
gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCAAAACAAGCGCCAAGGTGCCTGGCACCGGATTTGCCGGTGTTCG‑CGGG‑C‑A‑c‑g‑c‑c                                                                                                                                   >  1:8788219/1‑137 (MQ=25)
  tGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGAGCCACGGTGCACGGCCCCAGACTTTTCGGTGTT‑TGCGGG‑G‑A‑c‑g‑c‑ccg                                                                                                                                 >  1:2876421/1‑135 (MQ=255)
         aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGG‑C‑A‑A‑GCC‑CGCTCCcac                                                                                                                           >  2:16161304/1‑137 (MQ=18)
         aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGAGAAAAAGGCGCCGAGGTGCATGGGAACGGCTTTGCCGGTGTT‑TGCGGG‑G‑A‑C‑GAC‑CGCTCCcac                                                                                                                           >  1:13416872/1‑137 (MQ=25)
              gagCCTGAACGAAGCCTTGGCACCACAACAGGATTGAGGCTGCCAAATGGTGGGAGCGGGCGTGCCCGAGAAACGGGCGCCGCGGTGCATGGCACTGGTTTTGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACt                                                                                                                      <  2:11572861/137‑1 (MQ=18)
                    gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCTCGAAACAAGCGCCGCGGTGCATGGCACCGACTTTTCCGGTGTTCG‑CGGG‑C‑ACC‑C‑C‑CACTCCCACAAACTCATTGc                                                                                                                >  2:16403143/1‑137 (MQ=255)
                    gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCACGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGG‑C‑A‑C‑GCC‑CGCTCCCACAGAATGATTGc                                                                                                                >  2:8638097/1‑137 (MQ=21)
                       cGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCAAAACAGGTGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGG‑C‑A‑C‑ACC‑CGCTCCCAAAGACTGATTgtcct                                                                                                             >  2:15640509/1‑133 (MQ=18)
                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGAATGATTGCCAtt                                                                                                            >  2:9156560/1‑137 (MQ=21)
                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGAGCCGCGGTGTATGGCACCGGCTTTTCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CACTCCCACAAACTGATTGCCact                                                                                                            >  1:12922201/1‑135 (MQ=18)
                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACCGGCGCCGCGGTGCATGGCACCGGATTTGCCGGTGTT‑TGCAGG‑G‑A‑C‑GCC‑CGCTCCCACAGACTGAGTGGCagt                                                                                                            >  2:7477864/1‑135 (MQ=18)
                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACAAGCGCCGAGGTGTATGGCACCTGATTTTCCGGTGTTCG‑CAGG‑CAA‑A‑C‑C‑CGCTCCCACAGAATGATTGCCCtt                                                                                                            >  1:3845641/1‑137 (MQ=12)
                            cattGGCACAAAAGCAGGAATGAGGCTGCAAAACGGTGGGAGCGGGCATGCCCGCGAAACAGGCGCCGCGGTGCATGGACCCGGGTTTGACGGTGTTCT‑CGGG‑C‑A‑C‑GCC‑CAGTCCCACAGACCTATTGCCATTAAgt                                                                                                        >  2:12078630/3‑136 (MQ=16)
                              ttGACACCACAGAAGAACTGAGGCTGCCAAACGGTGGGAGCGGGCTTTCCTAAGAAAAAGTCGCCGTGGGACAGGACACCGGTTTTCCCGGTGTT‑CGCGGA‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCgg                                                                                                      <  2:77658/137‑1 (MQ=255)
                                 gCACCACAGCAGGGCTGAGGCTGCCAAACGTTGGGAGCGGGCGGCCCCGCGAAACAGGCCCCAAGGGGCAGGGCACTGGCTTTGCCGGTGTT‑TGCGGG‑C‑A‑C‑ACC‑CGCTCCCGCAGACTGATTGCCATTAAGCGGGcc                                                                                                   <  2:3931265/137‑1 (MQ=255)
                                 gCACCACAGCAGGACGGAGGCTGCCAAATGGTGGGAGCGGGTGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGTTTTGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                                   <  2:8857660/137‑1 (MQ=21)
                                 gCACCACAGAAGGATGGGGGTTGCAAAATGGGGGGAGGGGGTGTCCCCGAAAAACAGGCGCCGCGGTGAAGGGCACTGGTTTTGCTGGTGTT‑CGCGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                                   <  2:1483340/137‑1 (MQ=12)
                                 gCACCACAGAAGGATGGGGGCTGCAAAATGGGGGGAGGGGGTGTGCCCGCGAAACGGGCGCCGTGGGGCGTGACACCGGTTTTGCGGGTGTT‑CGTGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                                   <  2:5991782/137‑1 (MQ=12)
                                 acaccacATCAGGACTGAGGCGGCAAAATGGGGGGAGCGGGTGTGCCAGAGAAACAGGCACAGCGGGACATGGCACCGGTTTTGCCGGTTTTCG‑CGGG‑A‑A‑C‑CCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                                   <  2:9717769/136‑1 (MQ=255)
                                      acaGCAGGACTGAGGCTGCCAAATGGTGGGAGCGGGTGTGCCCGCGAAACAGGCGCCGTGGGGCATGGCACTGGCTTTGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                                              <  2:16224513/137‑1 (MQ=21)
                                         gcagGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCGAAACAAGCACCACGGTGCCCGGCACCCGCTTTTCCGGTGTTCT‑CGGG‑C‑A‑CAC‑C‑CCCCCCCCCAAAATGATTGCCACTCCACGGGACTATCACtc                                                                                           >  1:11976457/1‑136 (MQ=255)
                                             gACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGCCCACGAAACAAGTACCACGGTGCATGGCACCAGGTTTGCCGGTGTTCG‑CGGG‑C‑A‑CAC‑C‑CACCCCCCCAGACCGAGTGCCCTTCAGCGGGACTATCGCTaaaat                                                                                       >  1:1816872/1‑133 (MQ=9)
                                             gACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCAAAACAAGAGCCACGGTGCATGGCACCAGGTTTGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CTCTCCCCCAGAATGATTGCCATTAAGAGGGCCTCTCGCTAAAgg                                                                                       >  1:7944955/1‑136 (MQ=12)
                                                     cTGCCAAACGGTGGGAGCGGGCGTTCCCACGAAACAGGAGCCGCAGTGCATGGCACCCGATTTGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CCCCCCCACAAAATGATTGCAATTAAGCGGGCCTCTTGCTACAATGAAGGacc                                                                               >  2:15661874/1‑137 (MQ=12)
                                                     cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑TGCGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCCTTCAGCGGGCCTATCGCTACAGTGACGGccc                                                                               >  2:7075035/1‑137 (MQ=21)
                                                     cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGtcc                                                                               >  2:12217738/1‑137 (MQ=21)
                                                     cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAGGGGCCGCGGGGCATGGCACCGGGTTTGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGGCTGATTGCCATTAAGAGGGCCTATCGTTACAATGACGGcac                                                                               >  2:15941076/1‑135 (MQ=17)
                                                     cTGCCAAACGGTGGGAGCGGGCATGCCCGCAAAACAAGCACCACAGTGCCCGGCACCAGCCTTGCCCGTGTTCC‑CGGG‑C‑A‑C‑GCC‑CACCCCCCCAAACCGATTGTCACTAAACAGGCCCATCACTCCAATGACAGccc                                                                               >  2:685039/1‑137 (MQ=255)
                                                     cTGCCAAACGGTGGGAGCGGGAGTGCCCACGAAACAAGCGCCACGGTGCATGGCACCGGATTTGTCGGTGTTCT‑CGGG‑C‑A‑C‑GCC‑CACTCCCACAAAACGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGccc                                                                               >  2:2668583/1‑137 (MQ=255)
                                                          gaaCGGTGGGAGTGGGTGTGGCCGAGGAAAAGGCGCCGCGGTGCATGGAACCGGATTTGAAGGTGTTCG‑CGGG‑C‑A‑A‑G‑CACGCTCCCACAAACTTATTGCCATTAAGTGGGCCTATCACTACAGTGACAGCCCGtttt                                                                          >  2:15439402/2‑137 (MQ=255)
                                                          aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCACGGCACCAGCTTTGCCGGTGTTCG‑CGGG‑C‑A‑CAC‑C‑CGCTCCCACAGACTGATTGCCACTAAGCGGGAATATCGATACAGTGACGGAACGtttt                                                                          >  2:8547538/1‑137 (MQ=17)
                                                                  ggAACGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGTCTTTGCCGGTGTT‑CGTGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGcc                                                                  <  2:16713039/137‑1 (MQ=21)
                                                                     gCGGGCGTGCCCACGAAACAAGCGCCGCGGTGCACGGCACCGGATTTGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAAACTGATTGCCATTCAGCGGGCCTACCGCTACAGTGACGGCCCGTTTTTTTTTGCCTCg                                                               >  2:4021771/1‑137 (MQ=18)
                                                                             gCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGaa                                                       >  1:18377998/1‑137 (MQ=25)
                                                                                cgtgAAACAGGCGCCGTGGGGCGTGGCACCGGCTTTGCTGGTGTT‑TGCGGG‑C‑A‑C‑TCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    <  1:16280599/134‑1 (MQ=255)
                                                                                cgcgAAATGGGCACCGCGGTGCATGGCACCGATTTTCCCGGTGTTCG‑GGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    <  2:2797099/137‑1 (MQ=25)
                                                                                cgcgAAACGGGCTCCACGGTGCATGGCACCGGCTTTACTGTTTTT‑CGTGGG‑C‑A‑C‑TCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    <  1:3391724/137‑1 (MQ=255)
                                                                                cgcgAAACGGGCGCCGTGGGGCAGGACACCGGTTTTGTAGGTGTTCG‑CGGA‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    <  2:1746102/137‑1 (MQ=255)
                                                                                cgcgAAACAGGCGCCGTGGTGCGTGGCACCGTTTTTGCCGGTGTTCG‑TGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    <  1:15377013/137‑1 (MQ=37)
                                                                                cgcgAAACAGCCCCCGTGGGGAATGGCACTGGCTTTGCTGGTGTTCA‑GGGG‑C‑A‑T‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    <  1:8479008/137‑1 (MQ=255)
                                                                                cgagAAACGGGCGCGGTGGTGCATGGCACCGGTTTGGCGGGTTTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    <  2:4806553/134‑1 (MQ=25)
                                                                                cgaaaaaCAGGAGCTGTGGTGCGTGGCACCGGTTTTGCCGGTGTT‑CGTGGG‑C‑A‑C‑ACC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    <  1:982101/133‑1 (MQ=255)
                                                                                 gcgAAACAGGCGCCACGGTGCATGGCACCGGTTTTGCTGGTGTT‑CGCGGA‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                                   <  2:8053469/137‑1 (MQ=255)
                                                                                 acaaaaCAGGCGCCCCGATGCAGGGCACCGCCTTTGCAGGTGTCCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTTACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                                   <  2:3343609/134‑1 (MQ=25)
                                                                                   gAAAAAGGCGCCGCGGTGCGTGGCACCGGCTTTGCGGGGGTT‑CGCGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGa                                                 <  2:17958142/137‑1 (MQ=37)
                                                                                     aataGGCGCCGCAGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGATACAGTGACGGCCAGTTGTTTGTTGCCTCGCTAACGAATCGTGCCt                                               >  2:10389638/4‑137 (MQ=21)
                                                                                           cgctgtggTGCGTGGCACCGGCTTTGCCGTTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                         <  1:13832102/131‑1 (MQ=25)
                                                                                           cgccgtGGTGCATGGCACCGGTTTGGCTGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                         <  2:7873477/137‑1 (MQ=37)
                                                                                              cgcgGTGCATGGCACCGGCTTTGCCGGTGTT‑CGGGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGaa                                      <  2:17178491/137‑1 (MQ=37)
                                                                                                     cATGGCACCGGCTTTGTCGGTGTTAG‑CGGG‑C‑A‑C‑TCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                               <  1:6980392/137‑1 (MQ=255)
                                                                                                     cATGGCACCGGCTTTGCCGGTGTC‑CGCGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGAATCCCTTAGAACGCTATc                               <  2:12989992/137‑1 (MQ=37)
                                                                                                       tggcaccGGCTTTTCCGTTGTT‑CGTGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCaa                             <  2:14593273/131‑1 (MQ=37)
                                                                                                       gggcaccGGCTTTTCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCaa                             <  2:15795445/131‑1 (MQ=37)
                                                                                                        ggcaccGGCTTTGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg                            >  2:5953570/6‑137 (MQ=37)
                                                                                                        ggcaccGGCTTTGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg                            >  1:10193723/6‑137 (MQ=37)
                                                                                                             cGGCTTTGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGAt                       >  1:14354154/1‑137 (MQ=38)
                                                                                                             cGGCTTTGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGAt                       >  2:14298372/1‑137 (MQ=38)
                                                                                                                cTTTGTCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  2:10904530/1‑137 (MQ=255)
                                                                                                                cTTTGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    >  1:14493433/1‑105 (MQ=37)
                                                                                                                cTTTGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                         >  2:3397235/1‑116 (MQ=38)
                                                                                                                cTTTGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  1:6084168/1‑137 (MQ=39)
                                                                                                                cTTTGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  2:11385035/1‑137 (MQ=39)
                                                                                                                cTTTGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  2:9324748/1‑137 (MQ=39)
                                                                                                                cTTTGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  1:11478442/1‑137 (MQ=39)
                                                                                                                cTTTGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  1:187898/1‑137 (MQ=39)
                                                                                                                cTTTGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  2:11638010/1‑137 (MQ=39)
                                                                                                                cTTTGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  1:13795327/1‑137 (MQ=39)
                                                                                                                cTTTGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  2:4498712/1‑137 (MQ=39)
                                                                                                                cTTTGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  1:9341786/1‑137 (MQ=39)
                                                                                                                cTTTGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  1:15655455/1‑137 (MQ=39)
                                                                                                                cTTTGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  2:17072888/1‑137 (MQ=39)
                                                                                                                cTTTGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  1:7690791/1‑137 (MQ=39)
                                                                                                                cTTTGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCACTTAGAACGCTATCAAGCAGATCTg                    >  2:10725640/1‑137 (MQ=38)
                                                                                                                cTTTGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCATATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  2:17811088/1‑137 (MQ=38)
                                                                                                                cTTTGCCGGTGTT‑CGCGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGGTTGCCATCAAGCGGGCCTATTGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGCATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  2:17507578/1‑137 (MQ=38)
                                                                                                                 ttggCCGGGGTT‑T‑GGGGAC‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATGAAGCGGGCTTATCGCTACAGTGACGGCCCGTTTGTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGTTATCAAGCAGATCTGa                   <  2:5378544/134‑1 (MQ=255)
                                                                                                                  ttGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa                  >  2:18175226/1‑137 (MQ=39)
                                                                                                                  ttGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa                  >  2:14149559/1‑137 (MQ=39)
                                                                                                                  ttGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa                  >  2:3842495/1‑137 (MQ=39)
                                                                                                                  ttGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa                  >  1:10395739/1‑137 (MQ=39)
                                                                                                                  ttGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa                  >  2:12329714/1‑137 (MQ=39)
                                                                                                                  ttGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa                  >  1:17836665/1‑137 (MQ=39)
                                                                                                                  ttGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa                  >  2:11940172/1‑137 (MQ=39)
                                                                                                                  ttGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa                  >  1:5625481/1‑137 (MQ=39)
                                                                                                                  ttGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa                  >  1:1695341/1‑137 (MQ=39)
                                                                                                                  ttGCCGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa                  >  1:12668476/1‑137 (MQ=39)
                                                                                                                  ttGCCGGTGTT‑CGCGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa                  >  2:6449007/1‑137 (MQ=39)
                                                                                                                     ccGGTGTTCG‑CGGG‑C‑A‑C‑GCC‑CGCTCACACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCGCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATATGAAACg               >  2:2830242/1‑137 (MQ=255)
                                                                                                                          tgttcgcGGG‑C‑A‑C‑TCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTccc           <  1:18011307/131‑1 (MQ=255)
                                                                                                                          tgttcgcGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTccc           <  1:4155272/131‑1 (MQ=255)
                                                                                                                            ttcgcGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGa         >  1:16642225/5‑137 (MQ=255)
                                                                                                                             tcgcGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAc        >  1:3463661/4‑137 (MQ=255)
                                                                                                                              cgcGGG‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAct       >  2:13312789/3‑137 (MQ=255)
                                                                                                                                 ggg‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActtct    >  2:472074/1‑137 (MQ=255)
                                                                                                                                 ggg‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActtct    >  2:2549475/1‑137 (MQ=255)
                                                                                                                                 ggg‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActtct    >  1:6778454/1‑137 (MQ=255)
                                                                                                                                   g‑C‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACttcttc  >  2:17698262/1‑137 (MQ=255)
                                                                                                                                            |                                                                                                                                    
GATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTT‑T‑CGGG‑C‑A‑C‑G‑C‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTC  >  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 < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑