Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451_phaM1 1,499,498 1 bp→CG intergenic (+75/‑89) 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,4971.C77.6% 467.0 / 195.9 216intergenic (+74/‑90)trpS CDS (Tryptophan‑‑tRNA ligase)/zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain)–/–
Reads supporting (aligned to +/- strand):  ref base . (21/27);  new base C (130/36);  total (151/65)
Fisher's exact test for biased strand distribution p-value = 9.86e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*AE015451_phaM11,499,4980TG85.0% 440.0 / 27.6 216intergenic (+75/‑89)trpS CDS (Tryptophan‑‑tRNA ligase)/zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain)–/–
Reads supporting (aligned to +/- strand):  ref base T (20/8);  new base G (119/39);  total (151/65)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 232 total aligned reads displayed.

TGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑CACGGCTTTGCCGGTG‑TT‑T‑CG‑GGCAC‑G‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTT  >  AE015451_phaM1/1499363‑1499633
                                                                                                                                          |||                                                                                                                                         
tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGTCCCCGAAACAAGCGCCGCGGTGCCTTGC‑CCCGGCTCTTCCGGTG‑tt‑tg                                                                                                                                            >  1:11151487/1‑136 (MQ=255)
          gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGCCCCCGAAACAAGAGCCACTGTGCACGGC‑ACCAGACTTTTCGGTG‑TT‑TGCG‑GGAAc‑g‑cc                                                                                                                               >  1:10365934/1‑137 (MQ=255)
          gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACAAAACAAGCGCCGCTGTGCATGGC‑ACCAGCTTTTCCCGTG‑TTCG‑CG‑GGCAc‑g‑cc                                                                                                                               >  1:14698481/1‑137 (MQ=25)
          gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGTCCACGAAACAAGCGCCACGGTGCCTTGCACCCGCCTTTCCCGTTT‑TC‑T‑CG‑GGcccacc                                                                                                                                 >  1:7926179/1‑132 (MQ=255)
            tGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCATGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTTGCCGGTG‑TTCG‑CG‑GGCAc‑g‑cccg                                                                                                                             >  1:13546675/1‑135 (MQ=25)
            tGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGTCCGCGAAACAGGCGCCGCGGTGCCTGGC‑ACCGGCTTTGCCGGTG‑TTCG‑CG‑GGCAc‑g‑ccca                                                                                                                             >  1:18477335/1‑135 (MQ=37)
                 ccATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGTCCGCGAAAAAGGCGCCGCGGTGCATGGC‑AACGGCTTTGCCGGTG‑TTCG‑CG‑GGCACAC‑CCGCTccc                                                                                                                         >  2:506186/1‑137 (MQ=25)
                       agggCCTGAACGAAGCCTTGGCACCACAACATGATTGAGGCTGCCAAACGGTGGGAGCGGGTGGGCCCGCGAAACGGGCGCCGCGGTGCATGGC‑ACTGGTTTTGCCGGTG‑TTCG‑TG‑GGAAC‑GCCCGCTCCCACAGAc                                                                                                                   <  2:12209045/134‑1 (MQ=255)
                        gagCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAAAAGGCGCCGCGGTGCATGGC‑AACGGCTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACt                                                                                                                  >  2:12727240/1‑137 (MQ=25)
                        gagCCTGAACGAAGCCTTGGCACCACAACAGGATTGAGGCTGCCAAATGGTGGGAGCGGGCGTGCCCGAGAAACGGGCGCCGCGGTGCATGGC‑ACTGGTTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACt                                                                                                                  <  2:11572861/137‑1 (MQ=18)
                        aagCCTGAACAAAACCTTGGCACCAAGGAAGGACGGAGGCTGCAAAACGGGGGGAGCGGGCGTGCCCGCGAAACAGACACCGCGGTGCATGGC‑ACCGGCTTTGCCGGTG‑TTCG‑CG‑GACAC‑GCCCGCTCCCACAGACt                                                                                                                  <  2:17107008/136‑1 (MQ=9)
                                 cGAAGCTTTGTCACCACAGCAGGGTTGAGGCTACCAAACGGTGGGGGCGGGCGTGCGCGCGAAACGGGCGCCGTTGTGCGTGGC‑ACCGGCTTTGCTGGTG‑TT‑TGCG‑GGCAT‑GCCCGCTCCCACAGACTGATTGCCAt                                                                                                         <  2:15808892/137‑1 (MQ=16)
                                 cGAAGCATTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAATGGATGTCGCGGTGCATGGC‑ACCGGCTTTGCCGGTG‑TTCG‑CG‑GGCGC‑GCCCGCTCCCACGGACTGATTGCCAt                                                                                                         >  2:10516894/1‑137 (MQ=18)
                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACCGGCGCCACGGTTCATGGC‑ACCAGCTTTGCCGGTG‑TTCT‑CA‑GGTAC‑ACCCCCTTCCACAGACTGATTGGCAtt                                                                                                        >  1:13167018/1‑137 (MQ=21)
                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGC‑ACCAGCTTTGCCGGTG‑TTCT‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCAtt                                                                                                        >  2:16580361/1‑137 (MQ=21)
                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTTGCCGGTG‑TTCT‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCAtt                                                                                                        >  2:15614547/1‑137 (MQ=25)
                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACCGGCGCCGCGGTGCATGGC‑ACCGGATTTGCCGGTG‑TT‑TGCA‑GGGAC‑GCCCGCTCCCACAGACTGAGTGGCagt                                                                                                        >  2:7477864/1‑135 (MQ=18)
                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCAAAAAAAGAACCTCAGTGCATGGC‑ACCAGCTTTTCCAGTGTTT‑T‑CG‑GGCACAC‑CCCCTCCCACCAACCTATTGcccct                                                                                                        >  2:5433856/1‑134 (MQ=255)
                                  aaaGCTTTGGCACCACAGCAGGATGGGGGCTGCAAAATGGTGGGAGGGGTCGGCCCTGCAAAACGGGCGCCGTGGGGAATGGC‑ACTGGTTTTGCTGGTG‑TTCG‑TG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCAtt                                                                                                        <  1:12792947/136‑1 (MQ=12)
                                      cattGGCACAAAAGCAGGAATGAGGCTGCAAAACGGTGGGAGCGGGCATGCCCGCGAAACAGGCGCCGCGGTGCATGGA‑CCCGGGTTTGACGGTG‑TTCT‑CG‑GGCAC‑GCCCAGTCCCACAGACCTATTGCCATTAAgt                                                                                                    >  2:12078630/3‑136 (MQ=16)
                                        ttGACACCACAGAAGAACTGAGGCTGCCAAACGGTGGGAGCGGGCTTTCCTAAGAAAAAGTCGCCGTGGGACAGGAC‑ACCGGTTTTCCCGGTG‑TT‑CGCG‑GACAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCgg                                                                                                  <  2:77658/137‑1 (MQ=255)
                                           gCACCACGGCGGGATTGAGGTGACCAAATGGTGGGAGGGGTTGTGCCCGCGAAACAGGAGCCAAGTGGCATGGA‑ACTGGTTTTTCCGGTT‑TTCG‑CG‑GGCAC‑ACCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                               <  2:14700943/137‑1 (MQ=255)
                                           gCACCACAGCAGGGCTGAGGCTGCCAAACGTTGGGAGCGGGCGGCCCCGCGAAACAGGCCCCAAGGGGCAGGGC‑ACTGGCTTTGCCGGTG‑TT‑TGCG‑GGCAC‑ACCCGCTCCCGCAGACTGATTGCCATTAAGCGGGcc                                                                                               <  2:3931265/137‑1 (MQ=255)
                                           gCACCACAGCAGGACTGAGTCTCCCAAACGGTGGGAGCGGGTGTGCCCGCGAAACGGGCGCCGCGGTGCAGGAA‑ACCGGCTTGGCCGTTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                               <  2:2851360/137‑1 (MQ=18)
                                           gCACCACAGCAGGACGGAGGCTGCCAAATGGTGGGAGCGGGTGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGTTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                               <  2:8857660/137‑1 (MQ=21)
                                           acaccacATCAGGACTGAGGCGGCAAAATGGGGGGAGCGGGTGTGCCAGAGAAACAGGCACAGCGGGACATGGC‑ACCGGTTTTGCCGGTT‑TTCG‑CG‑GGAAC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                               <  2:9717769/136‑1 (MQ=255)
                                            caccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGC‑ACCGGCTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGAATGATTGCCATTAAGCGGGCCt                                                                                              >  2:15651896/1‑137 (MQ=21)
                                                   gcagGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCGAAACAAGCACCACGGTGCCCGGC‑ACCCGCTTTTCCGGTG‑TTCT‑CG‑GGCACAC‑CCCCCCCCCCAAAATGATTGCCACTCCACGGGACTATCACtc                                                                                       >  1:11976457/1‑136 (MQ=255)
                                                           gAGGCTGCAAAACGGGGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGGGCATGAC‑ACCGGCTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACg                                                                               <  2:16166847/137‑1 (MQ=18)
                                                               cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTTGCCGGTG‑TT‑TGCG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCCTTCAGCGGGCCTATCGCTACAGTGACGGccc                                                                           >  2:7075035/1‑137 (MQ=21)
                                                                    aaaCGGTGGGAGCGGGCGTTCCCGCGAAACAGGCGCCACAGTGCACGGC‑ACCGGCTTTGCCTGTG‑TTCA‑CG‑GGCAC‑GCCCGCCCCCCCCGACTGATTGCCCTTCAACAGGCCTATCGCTACAGTGACGGCCCGtttt                                                                      >  2:8547585/1‑137 (MQ=255)
                                                                              aGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTTGCCGGTG‑TTCG‑CG‑GGCACAC‑CCGCTCCCACAGACCGATTGCCATTAAGCGGGCCTATCGCTAAAATGACGGCCCGTTTTTTGTTGCCTc                                                            >  2:18046116/1‑137 (MQ=18)
                                                                               gCGGGCGTGCCCACGAAACAAGCGCCGCGGTGCACGGC‑ACCGGATTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAAACTGATTGCCATTCAGCGGGCCTACCGCTACAGTGACGGCCCGTTTTTTTTTGCCTCg                                                           >  2:4021771/1‑137 (MQ=18)
                                                                                cGGGCATTCCCGCAAAACAGGCGCCGCGGCGCATGGC‑ACCTACTTTACCGGTG‑TTCG‑CG‑GGCACTC‑CCACTCCCACAGACTGATTACCATTAAACGGGCCTATCGCTACACTGACGGCCCGTTTTTTGTTGCCTCGc                                                          <  2:15756961/137‑1 (MQ=255)
                                                                                       gCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCACGTTTTTTGTTGCCTCGCTAACGaa                                                   >  1:18377949/1‑137 (MQ=25)
                                                                                          ttcgAAACAGGCGCTACGGGGCATGGC‑ACTGGCTTTGCTGGTG‑TT‑CGCG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                <  2:17489204/135‑1 (MQ=25)
                                                                                          cgtgAAACAGGCGCCGAGGTGCCTTGC‑ACCGGCTGTGCTGGTG‑TT‑C‑GGTGGCAC‑GCCCGCTCTCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                <  2:4228083/134‑1 (MQ=255)
                                                                                          cgcgAAATGGGCACCGCGGTGCATGGC‑ACCGATTTTCCCGGTG‑TTCG‑GG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                <  2:2797099/137‑1 (MQ=25)
                                                                                          cgcgAAACGGGTGCCGTGGGAGATGGC‑ACCGGTTTGGTTGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                <  1:1644419/137‑1 (MQ=21)
                                                                                          cgcgAAACAGCCCCCGTGGGGAATGGC‑ACTGGCTTTGCTGGTG‑TTCA‑GG‑GGCAT‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                <  1:8479008/137‑1 (MQ=255)
                                                                                          cgaaaaaCAGGAGCTGTGGTGCGTGGC‑ACCGGTTTTGCCGGTG‑TT‑CGTG‑GGCAC‑ACCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                <  1:982101/133‑1 (MQ=255)
                                                                                           gcgAAACAGGCGCCGTGTTCCGTGTC‑ACTGGCTTGGCTGTGG‑TT‑CGCG‑GGCAC‑TCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGCGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                               <  2:14250331/137‑1 (MQ=255)
                                                                                           gcgAAACAGGCGCCACGGTGCATGGC‑ACCGGTTTTGCTGGTG‑TT‑CGCG‑GACAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                               <  2:8053469/137‑1 (MQ=255)
                                                                                           acaaaaCAGGCGCCCCGATGCAGGGC‑ACCGCCTTTGCAGGTG‑TCCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTTACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                               <  2:3343609/134‑1 (MQ=25)
                                                                                             gAAAAAGGCGCCGCGGTACAGGGC‑ACTGGTTTTTCTGTTG‑TTCG‑CG‑GACAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGa                                             <  2:8738887/137‑1 (MQ=255)
                                                                                                aCAGGCGCCGCGGTGCATGGC‑ACCGGGTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAAACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGAATc                                          >  2:9215777/1‑137 (MQ=25)
                                                                                                aCAGGCGCCGCGGTGCATGGC‑ACCGGCTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTc                                          >  1:3223687/1‑137 (MQ=37)
                                                                                                     cgctgtggTGCGTGGC‑ACCGGCTTTGCCGTTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                     <  1:13832102/131‑1 (MQ=25)
                                                                                                     cgccgtGGTGCATGGC‑ACCGGTTTGGCTGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                     <  2:7873477/137‑1 (MQ=37)
                                                                                                     cgccgcGGTGCATGGC‑ACCGTCTTTGCTGGTG‑TTCG‑CG‑GACAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                     <  2:17925980/137‑1 (MQ=255)
                                                                                                        cgcgGTGCATGGC‑ACCGGCTTTGCCGGTG‑TT‑CGGG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGaa                                  <  2:17178491/137‑1 (MQ=37)
                                                                                                               cATGGC‑ACCGGCTTTGTCGGTG‑TTAG‑CG‑GGCAC‑TCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                           <  1:6980392/137‑1 (MQ=255)
                                                                                                               cATGGC‑ACCGGCTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                           <  2:13012918/137‑1 (MQ=37)
                                                                                                                   ggcaccGGCTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg                        >  2:5953570/6‑137 (MQ=37)
                                                                                                                   ggcaccGGCTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg                        >  1:10193723/6‑137 (MQ=37)
                                                                                                                      accGGCTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAg                     >  1:9082585/3‑137 (MQ=38)
                                                                                                                      accGGCTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAg                     >  1:10824106/3‑137 (MQ=38)
                                                                                                                        cGGCTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGAt                   >  1:14354154/1‑137 (MQ=38)
                                                                                                                        cGGCTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGAt                   >  1:8282516/1‑137 (MQ=38)
                                                                                                                        cGGCTTTGCCGGTG‑TT‑CGCG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCATATCGCTACAGTGACGGCCAGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGAt                   >  2:16763282/1‑137 (MQ=38)
                                                                                                                           cTTTGTCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  2:10904530/1‑137 (MQ=255)
                                                                                                                           cTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCCCAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTCGAACGCTATCAAGCAGATCTg                >  2:17811087/1‑137 (MQ=38)
                                                                                                                           cTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                >  1:14493433/1‑105 (MQ=37)
                                                                                                                           cTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                <  2:14493433/105‑1 (MQ=37)
                                                                                                                           cTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                     >  2:3397235/1‑116 (MQ=38)
                                                                                                                           cTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                     <  1:3397235/116‑1 (MQ=38)
                                                                                                                           cTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  2:2753269/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  2:11385035/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  2:11624348/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  2:6853634/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  2:16009809/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  2:3791724/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  1:18014491/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  1:15655455/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  1:12581058/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  2:17550587/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  1:9341786/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  1:1620795/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  2:9324748/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  2:163945/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCACTTAGAACGCTATCAAGCAGATCTg                >  2:10725640/1‑137 (MQ=38)
                                                                                                                           cTTTGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCAACAGACTGATTGCCATTACGCGGGCCTATCGATACAGTGACGACCCGTTTTTTGTTGCCTCGCTAACGAATCATGAATCCCTTAGAACGCTATCAAGCATATCTg                >  1:2451806/1‑137 (MQ=37)
                                                                                                                           cTTTGCCGGTG‑TT‑CGCG‑GGCAC‑GCCCGCTCCCACAGACTGATCGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCCTGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  2:4777418/1‑137 (MQ=38)
                                                                                                                             ttGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa              >  2:12329714/1‑137 (MQ=39)
                                                                                                                             ttGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa              >  1:6604616/1‑137 (MQ=39)
                                                                                                                             ttGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa              >  1:5625481/1‑137 (MQ=39)
                                                                                                                             ttGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa              >  2:6019695/1‑137 (MQ=39)
                                                                                                                             ttGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa              >  2:11940172/1‑137 (MQ=39)
                                                                                                                             ttGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa              >  2:2205429/1‑137 (MQ=39)
                                                                                                                             ttGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa              >  2:15412386/1‑137 (MQ=39)
                                                                                                                             ttGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa              >  1:17836626/1‑137 (MQ=39)
                                                                                                                             ttGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa              >  1:391340/1‑137 (MQ=39)
                                                                                                                             ttGCCGGTG‑TTCG‑CG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa              >  1:1695341/1‑137 (MQ=39)
                                                                                                                             ttGCCGGTG‑TT‑CGCG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa              >  2:6449007/1‑137 (MQ=39)
                                                                                                                             ttGCCGGTG‑TT‑CGCG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa              >  2:18322513/1‑137 (MQ=39)
                                                                                                                                   gtg‑tT‑TGCG‑GGCAC‑TCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTcc        <  1:3800781/137‑1 (MQ=255)
                                                                                                                                      tgtttgcG‑GGCAT‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTccc       <  1:18011277/131‑1 (MQ=255)
                                                                                                                                      tgttcgcG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTccc       <  1:4155272/131‑1 (MQ=255)
                                                                                                                                         tcgcG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAc    >  1:7600147/4‑137 (MQ=255)
                                                                                                                                          cgcG‑GGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAct   >  2:14138265/3‑137 (MQ=255)
                                                                                                                                          |||                                                                                                                                         
TGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑CACGGCTTTGCCGGTG‑TT‑T‑CG‑GGCAC‑G‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTT  >  AE015451_phaM1/1499363‑1499633

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