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.4% 825.7 / 71.6 267intergenic (+83/‑81)trpS CDS (Tryptophan‑‑tRNA ligase)/zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain)–/–
Reads supporting (aligned to +/- strand):  ref base . (19/1);  new base C (155/89);  total (177/90)
Fisher's exact test for biased strand distribution p-value = 2.75e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 282 total aligned reads displayed.

GATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTG‑TT‑‑T‑CGGGCAC‑G‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTT  >  AE015451_phaM1/1499373‑1499636
                                                                                                                                           |                                                                                                                                  
gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTC‑G‑CGGGGAc‑g‑cc                                                                                                                                  >  1:12830507/1‑137 (MQ=37)
gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCGAAACAAGCGCCGCGGTGCACGGCACCAGATTTTCCAGTG‑TTC‑T‑CGGGCAc‑g‑cc                                                                                                                                  >  2:14381565/1‑137 (MQ=25)
 aTGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCAAAACAGGCGCCGAGGTGCATGGCACCGGCTTTGCCGGTG‑TTC‑G‑CGGGCAc‑g‑ccc                                                                                                                                 >  2:8668969/1‑136 (MQ=25)
     ggCCATTAAGATCCGGAACGAAGCTTTGGCACCACAACAGGACTGAGGCTTCCAAACGGTGGGAGCGGGCGTGCCCACGAAACGGGCACTGCGGTGCATGGCACTGTTTTTGCTGGTG‑TT‑‑TGGGGGCAc‑g‑cccgctc                                                                                                                             <  2:17977249/137‑6 (MQ=12)
         aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCcac                                                                                                                          >  2:8386896/1‑137 (MQ=21)
            cAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGACCGCGAAACAGGCGCCGCGGTGTCTGGCACCAGCTTTTCCGGTG‑TTC‑G‑CGGGCAC‑GTCCGCTCCcacaaa                                                                                                                       >  1:554812/1‑135 (MQ=21)
                    gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCACGGCACCGGCTTTGCCAGTG‑TTC‑G‑CGGGCAC‑GCCCACTCCCACAGAATGATTGc                                                                                                               >  2:16211417/1‑137 (MQ=18)
                    gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGAGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGTCGGTG‑TTC‑G‑CGGGCACAC‑CCGCTCCCACAGACTGATTGc                                                                                                               >  2:5387511/1‑137 (MQ=21)
                    gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGTCGGTG‑TTC‑G‑CAGGCACAC‑CCGCTCCCACAGCCTGACTGc                                                                                                               >  2:7622155/1‑137 (MQ=18)
                    gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGACCGCGAAAAAGGCACCGCGGTGCATGGCACCAGCTTTGGCGGTT‑TTC‑A‑CGGGCCC‑GCCCGATCCCACAGGCTGATTGc                                                                                                               >  2:11411147/1‑137 (MQ=255)
                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGGCCGCGAAACAAGAGCCGCGGTGTATTGCACCAGCTTTGCCGGTG‑TT‑‑TGCGGGACC‑GCCCACTCCCACAGAATGAGTGACAtt                                                                                                           >  2:14493285/1‑137 (MQ=255)
                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGCCCACAAAAAAAGCGCCACGGGGGACGGCAACAGCTTTGTCAGTT‑TTC‑G‑CGGGCACAC‑CCACTCCCAAAAACTGATTTCCAtt                                                                                                           >  1:4370863/1‑137 (MQ=255)
                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCACGGTGCACTGCACCAGCTTTTCCGGTG‑TTC‑G‑CGGGCAC‑GCCCCCTCCCACAAAATGATTGCCAtt                                                                                                           >  2:8098455/1‑137 (MQ=18)
                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCAAAAAAAGCGCCGCGGTGCATGGCACCGGCTTTTCCGGTG‑TT‑‑TGCGGGCACAC‑CCGCTCCCACAGAATGATTGCCAtt                                                                                                           >  2:12625641/1‑137 (MQ=18)
                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TT‑‑CGCGGGCAC‑GCCCGCTCCCACAGACTGATTGCCact                                                                                                           >  2:15842965/1‑135 (MQ=21)
                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGTCCGCAAAACAGGAGCCACGGTGCATGGCACCGGCTTTTCCGGTGTTT‑‑T‑CGGGCAC‑GCCCACTCCCACCGACGGATTGCCCtt                                                                                                           >  1:17232763/1‑137 (MQ=21)
                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCGAAACAAGAGCCACAGTGCACGGCACCAGATTTGCCAGTG‑TT‑‑TGCGGGCAC‑GCCCACTCCCACAAAATGATTGCCact                                                                                                           >  2:6169870/1‑135 (MQ=255)
                              ttGGCACCACAGCGGGACTGAGGCTGCAAAAAGGTGGGAGCGGGTGTGCCAGAGAAACGGGCGCTGTGGTGCGGGGCACCGGTTTTGCTGGTG‑TT‑‑CGCGGGAAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCgg                                                                                                     <  2:10918848/137‑1 (MQ=255)
                              tgggaaccacAGAAGGAGTGGGGCGACAAAAGGGGGGGAGAGGGTGTTCCCGCGAAACAGGCTCCGTGGGGCATGGCACCGGTTTTGCTGGTG‑TTC‑G‑TGGGCAC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCgg                                                                                                     <  1:9022621/132‑1 (MQ=255)
                               tGGCACCACAGCGGGATTGAGGTTGCCAAATGGTGGGAGCGGTTGTTCCCGCGAAACGGGCGCCGTGGGGCAGGGAACTGGCTTTGCTGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCggg                                                                                                    <  2:1833063/137‑1 (MQ=16)
                               tGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGTGGGCTTGCCCGCGAAACGGGCGCCGCGGTGCATGGCACTTGCTTTGCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCggg                                                                                                    <  2:17134834/137‑1 (MQ=21)
                                 gaaccacAGAAGGATGGGGGCTACAAAATGGGGGGAGCGGGTGTTCCCGAGAAACAGGCGCCGTGGGGCGTGGCACTGGTTTTTCCGGTT‑TTC‑G‑CGGGCAC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                                  <  1:9341229/135‑1 (MQ=255)
                                 gCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TT‑‑CGAGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                                  <  2:16713815/137‑1 (MQ=21)
                                 gCACCACAGAAGGGTGGGTGCGGCAAAATGGGGGGAGTGGGAGTTCCCGAGAAACAGGCGCCGGGGGGTGTGTCACTGGTTTTGTTGTTG‑TT‑‑CGTGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                                  <  1:2070509/137‑1 (MQ=255)
                                 gCACCACAGAAGGATTGAGGCTGCAAACCGGTGGGAGGGGGCGTTCCCGAGAAAAGGGCGCCACGGGGCGGGGCACTGGTTTTGCTGGGG‑TT‑‑TGCGGGCAC‑TCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                                  <  1:12384849/137‑1 (MQ=255)
                                 gCACAACAGCAGGACGGAGGCTGCCAAACGTTGGGAGCGGGTGTGCCCGCGAAACGGGCTCCGCGGTGCATGGCACCGGCTTTGCTGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                                  <  2:14017001/137‑1 (MQ=18)
                                 acaccccAGCAGGACTGGGTTTGCCAAACGTGGGGATCGGTTTTCCCTGCGAAACGGGCGCCGTGGTGCGTGCCACCGTTTTTACCGGTG‑TTC‑G‑CGGGCAC‑ACCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                                  <  1:6619821/136‑1 (MQ=255)
                                 aaacaacaGAAGGATGGGGGCGGCAAAATGGGGGGGGGGGGGGGGCCCGCGAAAAGGGCGCGGTGGTGAGGGGCATTGGTTTTGCGGGTG‑TTC‑G‑TGGACAC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                                  <  1:3337782/132‑1 (MQ=255)
                                   accacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGTCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGAATGATTGCCATTAAGCGGGCCTa                                                                                                >  2:15524719/1‑137 (MQ=18)
                                    ccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCGAAACAGGCGCCGCGGTGCATGGTACCGGCTTTGCCGGTG‑TTCGT‑GGGGAAG‑C‑CCGCTCCCACAGACCGATTGCCATCAAAAGGGCCTAt                                                                                               >  1:14536208/1‑137 (MQ=18)
                                         gcagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTCAGCGGGCCCATCGCTa                                                                                          >  2:13691235/1‑137 (MQ=18)
                                              aCTGAGGCTGCCAAACGGTGGGAGCGGGGGTGCCCACAAAACAAGAGCCACGGTGTATGGCACCAGATTTTTCGGTG‑TTC‑A‑CGGGCACAC‑CCACCCCCCCAAACTGATTGCCATTAAACGGGCCCATCGCTGCAGTg                                                                                     >  1:4112236/1‑137 (MQ=255)
                                                 gAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCAGTTCACGGCACCGGCTTTGCCGGTT‑TTC‑G‑CGGGCCC‑GCCCGCTCCCACCGACTGATTGCCATTAAGCGGGCCTATCGCCACAGTGAca                                                                                  >  2:9880844/1‑136 (MQ=18)
                                                     cTGCCAAACGGTGGGAGCGGGCGTTCCCGCGAAACAGACGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTC‑G‑TGGGCAC‑GCCCGCTCCCAAAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGcac                                                                              >  2:7252958/1‑135 (MQ=21)
                                                     cTGCCAAACGGTGGGAGCGGGCGTGTCCGCGAAACAGGCGCCGCGGTGTATGGCAACGGATTTGTCGGTG‑TT‑‑TGCGGGCAC‑GCCCGCTCCCACAGACCGATTGCCATTCAGCGGGCCTGTCGATACAGTGACGGccc                                                                              >  2:2099695/1‑137 (MQ=21)
                                                     cTGCCAAACGGTGGGAGCGGGCGTGTCCGCGAAACAGGCGCCGCGGTGCCTGGCCCCGGCTTTGCCGGTG‑TT‑‑TGCGGGGAC‑GCCCGCTTCCACAAACTGATTGCCATTAAGCGGGCCTATTGCTACAGCGACGGccc                                                                              >  2:6722324/1‑137 (MQ=37)
                                                     cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCACGGCACCGGATTTGCCAGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAAAATGATTGCAATTCAACGGGCCTATCGCCACAATGACGGccc                                                                              >  2:12508491/1‑137 (MQ=17)
                                                     cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCGCGGTGCATGGCACCGGCTTTTCCGGTT‑TTC‑G‑CGGGCCC‑GCCCGCTCCCACAGACTGATTGCCATTCAGCGGGCCTCTCGCTCCAGTGACGGccc                                                                              >  2:14764576/1‑137 (MQ=17)
                                                          aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGTAGTGCATGGCACCAGCTTTGCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTTCCACAGAATGATTGCCATTAAGCGGGACCACCGCTACAGTGACGGCCCGtttt                                                                         >  2:11839400/1‑137 (MQ=17)
                                                          aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCCCAGACTGATTGCCCTTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttt                                                                         >  2:18772385/1‑137 (MQ=21)
                                                          aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCACGtttt                                                                         >  2:9333810/1‑137 (MQ=21)
                                                          aaaCGGTGGGAGCGGGCGTGCCCCCGAAAACGGCGCCACGGTGCCTGGCACCGGCTTTGCCGGTG‑TTC‑T‑CGGGCAC‑GCCCACCCCCACAGAACGATTGCCATTCAACGGGCCTCTCGCTACAGTGGCGGTCCAtttt                                                                         >  2:10754704/1‑137 (MQ=12)
                                                           aaCGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCGCGGTGCACGGCACCGGCCTTTCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAAAACGATTGACCTTAAGCGGGCCTATCGCTCCAGTGACGGCCCGttttt                                                                        >  2:6925553/1‑137 (MQ=17)
                                                            aCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGACCGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt                                                                       >  2:15937897/1‑137 (MQ=18)
                                                              tttGGGAGCGGGCGTCCCCGCGAAATAGGCGTCGTGGGGCAGGGCACTGGTTTTGCTGGTG‑TTC‑G‑CGGGCAC‑TCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgt                                                                     <  1:5784409/135‑1 (MQ=255)
                                                                   gAGCGGGCGTGCCCGCGAAACAAGCGTCCCAGTGCCTGGCCCCAGCTTTGTCCGTG‑TTC‑G‑CCGGCACAC‑CCGCTCCCACCAACTGATTGTCATTAAGCGGGCCTCCCTCTCCAGTGGCGGACCGTTTTTTGTTGCCt                                                                >  2:400945/1‑137 (MQ=9)
                                                                     gCGGGCGTGCCCACGAAACAGGCGCCGCGGTGCCTGGCACCGGCTTTGCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGACCCCACAAAATGGTTGCCATTAAGCGGGGCTATCGCTACAGTGAAGGCCCGGTTTTTGTTGCCTCg                                                              >  2:10089338/1‑137 (MQ=18)
                                                                       gggTGTGCCTGAGAAACAGGCGCCGCGGTGCGTGGCACCGGCTTTGCCGGTG‑TT‑‑CGCGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCt                                                            <  2:16479936/137‑1 (MQ=25)
                                                                              cccGCGAAACAGGCACCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TT‑‑TGCGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAAt                                                     <  2:10593367/137‑1 (MQ=25)
                                                                                tgcgAAACAGGTTCCGTGGGGCAGGGCACTGTTTTTGTGGGTG‑TT‑‑TGGGGGCAC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                   <  2:12828264/136‑1 (MQ=255)
                                                                                cgcgAAACAGGTGCCGCGGTGCATGGCACCGTCTTTACCGTTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                   <  2:1977617/137‑1 (MQ=37)
                                                                                cgcgAAACAGGCGCCGCGGTGCATGACACCGGCTTTTCCGTTG‑TTC‑G‑CGGGCAC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                   <  2:15655868/137‑1 (MQ=255)
                                                                                cgagAAACGGGCTCCGTGGTGCATGGAACTGGTTTTGCTGGTG‑TT‑‑CGGGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                   <  2:1234948/134‑1 (MQ=25)
                                                                                cacgAAACAGTCGCCGCGGTGCATGACACCGGCTTTGTCGTTT‑TCC‑G‑CGGGCAC‑TCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                   <  1:8593520/135‑1 (MQ=255)
                                                                                 gcgAAACGGGCGCTGTGGGGCAGGGCACTGGCTTTGCTGGTT‑TTC‑G‑CGGGCAC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                                  <  1:15761641/137‑1 (MQ=255)
                                                                                     aaCAGGCGCCGAGGGGCGTGGCACCGTCTTTGCCGGTT‑TCC‑G‑CGGACAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACt                                              <  1:16757533/137‑1 (MQ=255)
                                                                                           cttcgcgGGGCAGGGCACTGGTTTTGCTGGTG‑TTC‑G‑CGGACAC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                        <  1:10923403/134‑1 (MQ=255)
                                                                                           ctccgcGGTGCATGGCACCGGCTTTTCCGGTG‑TT‑‑CGTGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                        <  2:6490289/135‑1 (MQ=37)
                                                                                           ctccgcGGTACGTGGCACTGGCTCTGCCGGTT‑TTC‑G‑GGGGCACTC‑CCGCTCACACAGACTGAGTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                        <  1:5500727/135‑1 (MQ=255)
                                                                                           cgctgcGGTGCATGGCACCGGCTTTGCTGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACATACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                        <  2:17446809/137‑1 (MQ=25)
                                                                                           cgccgtGGTGCGTGACACCGTCTTTGCCGGTG‑TT‑‑CGCGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                        <  2:7882999/137‑1 (MQ=37)
                                                                                           cgccgtGGTGCATGGCACTGTCTTTTCTGGTT‑TTC‑G‑CGGACAC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                        <  1:2432918/137‑1 (MQ=255)
                                                                                           cgccgcGGTGCGTGGCACGGGTTTTGCAGGTG‑TT‑‑CGTGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                        <  1:10425024/137‑1 (MQ=37)
                                                                                                cgGTGCATGGCACCGGCTTTGCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACg                                   <  2:17891395/137‑1 (MQ=37)
                                                                                                     cATGGCACCGGTTTTCCGGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                              <  1:7686484/137‑1 (MQ=37)
                                                                                                     cATGGCACCGGCTTTGCCGGTG‑TT‑‑CGCGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                              <  2:6417235/137‑1 (MQ=37)
                                                                                                          caccGGCTTTGCCGGGG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCa                         <  1:11420721/134‑1 (MQ=38)
                                                                                                           accGGCTTTGCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTACCGAATCATGACTCCCTTAGAACGCTATCAAGCAg                        >  1:14210515/3‑137 (MQ=38)
                                                                                                           accGGCTTTGCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAg                        >  1:13585468/3‑137 (MQ=38)
                                                                                                             cGGCTTTGCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGAt                      >  2:48268/1‑137 (MQ=38)
                                                                                                                cTTTGCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCGTGACTACCTTAGAACGCTATCAAGCAGATCTg                   >  2:15637092/1‑137 (MQ=38)
                                                                                                                cTTTGCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                   >  1:10052060/1‑137 (MQ=39)
                                                                                                                cTTTGCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                   >  2:2191682/1‑137 (MQ=39)
                                                                                                                cTTTGCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                   >  2:278876/1‑137 (MQ=39)
                                                                                                                cTTTGCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                   >  2:3548856/1‑137 (MQ=39)
                                                                                                                cTTTGCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                   >  1:12067149/1‑137 (MQ=39)
                                                                                                                cTTTGCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                   >  1:366654/1‑137 (MQ=39)
                                                                                                                cTTTGCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                   >  1:14339335/1‑137 (MQ=39)
                                                                                                                cTTTGCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                   >  1:15105480/1‑137 (MQ=39)
                                                                                                                cTTTGCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                   >  1:2748043/1‑137 (MQ=39)
                                                                                                                cTTTGCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                   >  2:10823602/1‑137 (MQ=39)
                                                                                                                cTTTGCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                   >  1:2707465/1‑137 (MQ=39)
                                                                                                                cTTTGCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                   >  2:13319116/1‑137 (MQ=39)
                                                                                                                cTTTGCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                   >  1:18466759/1‑137 (MQ=39)
                                                                                                                cTTTGCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCATTAGGACGCTATCAAGCAGATCTg                   >  2:13724964/1‑137 (MQ=38)
                                                                                                                cTTTGCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCATTAGAACGCTATCAAGCAGATCTg                   >  1:4452831/1‑137 (MQ=38)
                                                                                                                cTTTGCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGAATCCCTTAGAACGCTATCAAGCAGATCTg                   >  1:7335606/1‑137 (MQ=38)
                                                                                                                  ttGCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa                 >  2:15151939/1‑137 (MQ=39)
                                                                                                                  ttGCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa                 >  1:15624613/1‑137 (MQ=39)
                                                                                                                  ttGCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa                 >  1:10546976/1‑137 (MQ=39)
                                                                                                                  ttGCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACG‑ATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaaa                >  1:18858225/1‑137 (MQ=38)
                                                                                                                  ttGCCGGTG‑TTC‑G‑CGGGCAC‑GCCCGCTCCCACAGAATGATTGCAATTGAGCGGGCATATCGCTACAGTGACGGCAAGGTTTTTGTTGCCTCGCTAACGAATCATGACTACCTTAGGACGCTATCAAGTAGATCTGaa                 >  2:4386115/1‑137 (MQ=25)
                                                                                                                            tgttcgcGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTccc          <  2:13460306/131‑1 (MQ=255)
                                                                                                                              ttagcGGACAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTGTGTTGCCTCGCTAACGAATCATGACTNCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGa        <  1:18194793/133‑1 (MQ=255)
                                                                                                                               tcgcGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAc       >  1:18035305/4‑137 (MQ=255)
                                                                                                                               tcgcGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAc       <  2:2667833/134‑1 (MQ=255)
                                                                                                                               tcgcGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAc       >  2:17720065/4‑137 (MQ=255)
                                                                                                                               tcgcGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCNATCAAGCAGATCTGAAACGTCCCGAc       >  1:18617284/4‑137 (MQ=255)
                                                                                                                                   gggCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActtct   >  2:18612934/1‑137 (MQ=255)
                                                                                                                                    ggCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActtctt  <  1:12504200/137‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                  
GATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTG‑TT‑‑T‑CGGGCAC‑G‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTT  >  AE015451_phaM1/1499373‑1499636

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: ‑