Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451_phaM1 1,499,480 2 bp→AC intergenic (+57/‑106) 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,4800CA91.9% 590.9 / ‑0.4 286intergenic (+57/‑107)trpS CDS (Tryptophan‑‑tRNA ligase)/zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain)–/–
Reads supporting (aligned to +/- strand):  ref base C (17/6);  new base A (131/131);  total (148/138)
Fisher's exact test for biased strand distribution p-value = 3.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*AE015451_phaM11,499,4810AC95.5% 701.4 / ‑3.0 286intergenic (+58/‑106)trpS CDS (Tryptophan‑‑tRNA ligase)/zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain)–/–
Reads supporting (aligned to +/- strand):  ref base A (7/6);  new base C (141/132);  total (148/138)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 286 total aligned reads displayed.

CGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCAC‑GGCTTTG‑CCGG‑TG‑T‑T‑T‑CGGG‑CA‑C‑G‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATC  >  AE015451_phaM1/1499347‑1499608
                                                                                                                                     ||                                                                                                                                          
  cAGTTCGACAACGCTGAAGCACGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCATGCCCGCGGAAAAGGCGCCGTGGTGTATGGCAAC‑GGc                                                                                                                                       >  2:11059863/1‑137 (MQ=255)
    gTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCtt                                                                                                                                     <  1:16075181/137‑1 (MQ=255)
    gTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCtt                                                                                                                                     >  2:15866356/1‑137 (MQ=255)
    gTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCtt                                                                                                                                     <  1:13893984/137‑1 (MQ=255)
    gTTCGACAACGCTGAAGCCCGTGATGCGACCATTCAGATCCTGAACGAAGCCGTGGCACCACAGCAGGACGGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCtt                                                                                                                                     <  2:8855977/137‑1 (MQ=38)
      tCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTg                                                                                                                                   <  1:4708741/137‑1 (MQ=255)
      tCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTg                                                                                                                                   <  1:8371441/137‑1 (MQ=255)
      tCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTg                                                                                                                                   <  1:8767752/137‑1 (MQ=255)
      tCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTg                                                                                                                                   <  1:1542127/137‑1 (MQ=255)
      tCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTg                                                                                                                                   <  2:18362769/137‑1 (MQ=255)
      tCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTg                                                                                                                                   <  2:3848216/137‑1 (MQ=255)
      tCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTg                                                                                                                                   <  2:4222761/137‑1 (MQ=255)
      tCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTg                                                                                                                                   <  1:10382599/137‑1 (MQ=255)
            aCGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTG‑CCGg‑tg                                                                                                                           <  1:18042769/137‑1 (MQ=255)
                tGAAGCCCGTGATGCGGCCATTCAGAGCCTGCACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTG‑CCGG‑TG‑t‑tcg                                                                                                                     >  2:13896775/1‑135 (MQ=255)
                tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTG‑CCGG‑TG‑t‑tcg                                                                                                                     >  2:6507713/1‑135 (MQ=255)
                          gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGTCCGCGAAACAAGAGCCGCGGTGCATGGCACC‑GGTTTTT‑CCGG‑TG‑T‑T‑TACGGG‑CA‑c‑g‑cc                                                                                                      >  1:6044929/1‑137 (MQ=255)
                          gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGCCCGCGAAACAGGAGCCGGGGGGCACGGCACC‑CGCTTTT‑TCGG‑TG‑T‑T‑TGCGGG‑GA‑c‑g‑cc                                                                                                      >  1:128771/1‑137 (MQ=255)
                          gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTG‑CCGG‑TG‑T‑TCG‑CGGG‑CA‑c‑a‑cc                                                                                                      >  2:19005410/1‑137 (MQ=25)
                          gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCACGGTGCATGGCACC‑GGCTTTG‑CCGG‑TG‑T‑TCG‑CGGG‑CA‑c‑a‑cc                                                                                                      >  2:13168356/1‑137 (MQ=25)
                          gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACAAAAAAAGCGCCCCGGTGCATGGCCCC‑AGCTTTT‑TCAG‑TT‑T‑TCG‑CGGG‑CA‑c‑g‑cc                                                                                                      >  1:7585657/1‑137 (MQ=25)
                                 ctaTTCAGAGCTTGAACGAAGCCTGGGCACCACAGCAGGACTGAGGCTGCCAAACGGGGGGAGTGGGTGTGCCCGCGAAACGGGCGCTGCTGGGCCTGGCACT‑GGTTTTG‑CTGG‑TG‑T‑TAG‑CGGG‑CA‑C‑GCCCGCTccc                                                                                                <  2:12141575/135‑1 (MQ=255)
                                 ccATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCAGTGGCTGGCACC‑GGCTTTG‑CCTG‑TT‑T‑TCT‑CGGG‑CA‑C‑GCCCGCTccc                                                                                                >  2:15223691/1‑137 (MQ=25)
                                 ccATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCAAAAAAGGCGCCGCGGTGCATGGCACC‑GGCTTTG‑CCGG‑TG‑T‑TCG‑CGGG‑CA‑C‑GCCCGCTccc                                                                                                >  2:13985222/1‑137 (MQ=21)
                                   aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTG‑CCAG‑TG‑T‑TCG‑CGGG‑CA‑C‑GCCCGCTCCcaa                                                                                              >  2:14395035/1‑136 (MQ=21)
                                   aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGAGGTGCATGGCACC‑GGCTTTG‑CCAG‑TG‑T‑TCG‑CGGG‑CA‑C‑GCCCGCTCCcac                                                                                              >  2:5956804/1‑137 (MQ=21)
                                             tGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTG‑CCGG‑TG‑T‑TCG‑CGGG‑CA‑A‑GCCCGCTCCCACAGACTGATTg                                                                                    >  2:16317676/1‑137 (MQ=18)
                                              gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGTGAAACAGGAGCAGCGGTGTATGGCACC‑GGCTTTG‑CCGG‑TG‑T‑T‑TGCGGG‑CA‑C‑GCCCGCTCCCACGGGATGATTGc                                                                                   >  2:3079256/1‑137 (MQ=21)
                                              gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTG‑CCGG‑TG‑T‑TCG‑CGGG‑CA‑C‑GCCCTCTCCCACAGACTGATTGc                                                                                   >  2:12420983/1‑137 (MQ=21)
                                              gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTG‑CCGG‑TG‑T‑TCG‑CGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGc                                                                                   >  2:9157170/1‑137 (MQ=21)
                                              gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCACGGCACC‑GGCTTTG‑CCGG‑TG‑T‑TCG‑CGGG‑CA‑C‑GCCCGCTCCCACAGAATGATTGc                                                                                   >  2:14989420/1‑137 (MQ=21)
                                              gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCAAAAAAAGCGCCACGGTGCACGGCCCC‑TGATCTG‑CCGG‑TT‑T‑TCG‑CGGG‑CA‑C‑GCCCGCTCACCCAGAATGATTGc                                                                                   >  1:14574739/1‑137 (MQ=19)
                                                 cGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACAAGCGCCACGGTGCCTTGCCCC‑AGCTTTG‑CCGG‑TT‑T‑TCG‑CTGG‑CC‑CAC‑CCGCTCCAACAGAATGATTGccct                                                                                >  2:1661529/1‑135 (MQ=255)
                                                 cGAAGCCTTGGCACCAAAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGCAACAAGCGCCTCGGTGCATGGCAAC‑GGTTTTG‑CCAG‑TA‑T‑TCG‑CGGG‑CA‑CAA‑CCGCTCCCACAGACCGGTTTCCAt                                                                                >  2:12339615/1‑137 (MQ=18)
                                                  ggcaccTTGTCACCACAATAATACTGAGTTGACCAAATGTTGGAGGCGGGAGTGCCCGCGAAAAAGGCGCCGTGGTGCAGGGCCCT‑GGTTGTT‑TAGG‑TG‑T‑T‑C‑GGGGACA‑A‑GCCCGCTCCCACAGACTGATTGCCAtt                                                                               <  2:5112977/133‑1 (MQ=255)
                                                  gcaGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTAGGAGCGGGCGTGCCCGCGAAAAAAGCGCCGCCGTGCATGGCAAC‑GGCCTTG‑CCGG‑TG‑T‑TCT‑CGGG‑CA‑C‑GCCCCCTCCCACAGACTGATCGCCAtt                                                                               >  2:13459259/3‑137 (MQ=21)
                                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGTCCGCGAAAAAAGCGCCACGGTGCATGGCACC‑AGCTTTT‑CCGG‑TGTT‑T‑T‑TGGG‑CA‑CAC‑CCGCTCCCCCAGACTGACTTCCAtt                                                                               >  1:14189225/1‑137 (MQ=255)
                                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGTATGGCACC‑GGCTTTG‑CCGG‑TG‑T‑TCG‑CGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCAtt                                                                               >  2:17841217/1‑137 (MQ=21)
                                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGAGCCACAGTGTCTGGCCCC‑AGGTTTT‑CCGG‑TG‑T‑TCG‑CGGG‑AA‑C‑GCCCGCTCCCACAGACTGATTGTCAtt                                                                               >  1:13518501/1‑137 (MQ=18)
                                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACAGGCGCCGCGGTGCCTGGCACC‑GGCTTTG‑CCGG‑TG‑T‑TCG‑CGGG‑CA‑C‑GCCCGCTCCCACAGAATGATTGCCAtt                                                                               >  2:17929630/1‑137 (MQ=21)
                                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACAGGCGCCGCGGTGCATGGCACC‑AGCTTTG‑CCGG‑TG‑T‑T‑TGCGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCAtt                                                                               >  2:10497540/1‑137 (MQ=25)
                                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAAAGGTGGGAGCGGGCGTGTCCGCGAAACAGACGCTGAGATGTCTGTCACC‑GGCTTTG‑CCGG‑TG‑T‑TCA‑CGGG‑TA‑C‑GCCCGCTCCCAAAGACCGATTGCCAtt                                                                               >  2:10588537/1‑137 (MQ=19)
                                                  gAAGCCTTGCCACCACAACAGGACTGAGGCTGCAAACGGGGGGGAGCGGGCGTGCCCGCGAAACGGGCGCCGCGGGGCAGGGCACC‑GGCTTTG‑CAGG‑TG‑TCC‑G‑CGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCAtt                                                                               <  2:9231687/137‑1 (MQ=12)
                                                      ccTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAGGCGCCGCGGTGCATGGCAAC‑GGATTTG‑CCAG‑TG‑T‑TCG‑CGGG‑CA‑C‑GCCCGCTCCCACAGACCGATTGCCATTCAGc                                                                           >  2:17165786/1‑137 (MQ=21)
                                                      ccTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGAGAAACAGGCGCCGCGGTGTATGGCACC‑GGCTTTG‑CCAG‑TG‑T‑TCG‑CGGG‑CC‑C‑GCCCGCTCCCACAGGCTGATTGTCATTAAGc                                                                           >  2:10096653/1‑137 (MQ=18)
                                                      ccTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAAAAAGCGCCACGGTGCACGGCACC‑GGATTTT‑TCGG‑TG‑T‑TCT‑CGGG‑CA‑C‑GCCCACTCCCACAAACTGATTGCCCTTAAga                                                                           >  1:9116470/1‑136 (MQ=255)
                                                       cTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGACCACGAAACAAGCGCCGCGGTGCCTGGTACC‑GGCTTTT‑TCGGTTT‑T‑T‑T‑TGGG‑TA‑C‑GCCCCCCCCCACAAAATAATTGCCCTTCAGCg                                                                          >  1:3447035/1‑137 (MQ=255)
                                                        ttGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGCAAC‑GGATTTG‑CCGG‑TG‑T‑TCG‑CGGG‑AA‑C‑GCCCACTCCCACAGAATGATTGCCCTTAAGCgg                                                                         >  2:2554355/1‑137 (MQ=21)
                                                        ttGGCACCACAGAAGGATTGGGGTGGCAAAATGGTGGGAGTGGGTGTTCCCGAGAAACAGGCGCTG‑GGGGGAGGGGCACGGGTTTTG‑CTGG‑TG‑T‑T‑TGGGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCgg                                                                         <  1:1105792/137‑1 (MQ=255)
                                                          ggCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACAGGCGCCGCGGTGCATGGCACC‑AGCCTTG‑CCGG‑TG‑T‑TCG‑CGGG‑CA‑C‑GCCCGCCCCCACAAACTGAGTGCCAATAAGCgggg                                                                       >  2:785359/1‑136 (MQ=21)
                                                          ggCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGTCCGCGAAAAAGGCGCCGCGGTGCATGGCACC‑GGGTTTG‑CCGG‑TG‑T‑TCG‑CGGG‑CA‑C‑GCCCGCTTCCACAGACTGATTGCCATTCAGCgggg                                                                       >  2:7964594/1‑136 (MQ=21)
                                                           gaaccacGGAAGGATTGGGGTTGCCAAACGGTGGGAGTGGGCGGGCCCGCGAAAAGGGCGCTGTGGTGAGTGGCACT‑GGCTTTG‑CTGG‑TG‑T‑TCG‑GGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                      <  1:12141205/135‑1 (MQ=12)
                                                           gCACCGCAGCGGGGCGGGGGTTGCCATCTGGGGGGAGTGGTTGTGCCTGCGAAAAAGGCGCTTTGGGGCGTGGCACC‑GGTTTTG‑CTGG‑TG‑T‑TCT‑CGGG‑CA‑C‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                      <  1:19251508/137‑1 (MQ=255)
                                                           gCACCAAAGCAGGACTGAGGTTGCCAAACGGTGGGAGCGGGCGTGCCTGCGAAAAGGGCGCCGCGGTGCATGGCACC‑GGCTTTG‑CCGG‑TG‑T‑T‑CGCGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                      <  2:4859200/137‑1 (MQ=18)
                                                            caccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGTCCGCGAAACAGGAGACACGGTGCATGGCACC‑AGCTTTG‑CCGG‑TG‑T‑T‑TGCGGG‑TA‑CAC‑CCTCTCCCCCAGACTGACTGCCCTTCCGCGGGCCt                                                                     >  1:7326042/1‑137 (MQ=255)
                                                            caccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAGGCGCCGCGGTGCCTGGCACC‑GGCTTTG‑CCGG‑TG‑T‑T‑CGCGGG‑CA‑C‑GCCCGATCCCACAGACTGATTGCCCTTAAGCGGGctt                                                                     >  2:8626533/1‑135 (MQ=18)
                                                             accacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACAAGAGCCGCGGTGCATGGCACC‑AGCTTTG‑CCGG‑TG‑T‑TCG‑CGGG‑CA‑A‑GCCCACTCCCACAGAACGATTGCCACTCAGCGGGCCTa                                                                    >  1:18645930/1‑137 (MQ=14)
                                                                acaGCAGGACTGAGGCTACTAAATGGTGGGAGCGGGCGTTCCCTCGAAACAGGCCCCGCGGTACGTGGCACC‑GGCTTTG‑CCGG‑TG‑T‑TCG‑TGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                 <  2:4428047/137‑1 (MQ=18)
                                                                   gcagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTG‑CCGG‑TG‑T‑TCG‑CGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTa                                                              >  2:9120134/1‑137 (MQ=21)
                                                                   gcagGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCGAAACAAGCGCCGCGGTGTACGGCACC‑AGCTTTT‑CCGG‑TG‑T‑TCG‑CGGG‑CA‑CAC‑CCGCTCCCCCAAACCGACTTCCACTCAGCGGGCCTATCGCtc                                                              >  2:1613244/1‑136 (MQ=255)
                                                                         cTGAGGCTGCCAAACGGTGGGAGCGGGCGTGTCCGCGAAACAGGCGCCACGGTGCATGGCACC‑GGCTTTG‑CCAG‑TG‑T‑TCG‑CGGG‑CA‑C‑GCCCGCTCCCACAGACTGACTGCCACTCAGCGGGACTATCACTACAGTGa                                                        >  2:16831169/1‑137 (MQ=17)
                                                                          tGAGTCTGCCAAACGTGGGAAGCGGGTTTGCCCGCGAAACGGGCGCCGCGGTGCGTGGCACC‑GGCTTTT‑CCGT‑TG‑T‑TCG‑CGGG‑CA‑C‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAc                                                       <  2:3502440/137‑1 (MQ=255)
                                                                            aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGTACGGCACC‑GGCTTTG‑CCGG‑TG‑T‑TCG‑CGGG‑CA‑C‑GCCCTCTCCCCCAAAACTACTTCCATTAAAAAGGCCTATCACCACAGTGACgg                                                     >  2:6797797/1‑137 (MQ=9)
                                                                            aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCACCGCGGTGCATGGCACC‑GGATTTG‑CCGG‑TG‑T‑TCG‑CTGG‑CA‑C‑GCCCGCTCCCCCCAAATGACTTACATTAAGTGGGCCTACCGCTACAGTGACgg                                                     >  2:17111026/1‑137 (MQ=9)
                                                                               ctacaaaaCGGGGGGGGGGGGGGTGCCCGAGAAACAGGCGCTGCGGGGAATGGCACT‑GTTTTTT‑CCGG‑TG‑T‑TCG‑TGGG‑CA‑C‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGccc                                                  <  1:17479890/132‑1 (MQ=255)
                                                                               cTGCCAAACGGTGGGAGCGGGTGTGTCCTTGAAACAGGTGTCTCGGTGTCTGGTACC‑GGCTTTT‑TCGG‑TG‑T‑T‑TGCGGG‑TA‑C‑GCCCGCTCCCCCAGACCGGTTGCCATTCAGCGGGTCTCTCTCTACAGTGACGGccc                                                  >  2:15960531/1‑137 (MQ=255)
                                                                                 gCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGCACC‑GGCTTTG‑CCGG‑TG‑T‑T‑CGCGGG‑CA‑C‑GACCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTCCAGTGACGGCCCGt                                                >  2:18547803/1‑137 (MQ=21)
                                                                                    aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCAGTGCATGGCCCC‑GGCTTTG‑CCGG‑TG‑T‑TCG‑CGGG‑CA‑C‑GCCCGCTCCCACAGACTGACTGCCATTAAGAAGGACTATTGCTAAAGTGACGGTCCGtttt                                             >  2:9729314/1‑137 (MQ=17)
                                                                                     aaCGGTGGGAGCGGGCGTGTCCGCGGAACAGGCGCCGCGGTGCATGGCACC‑GGATTTG‑CCGG‑TG‑T‑TCG‑CGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGAGGGCCTATCGCTACAATGACGGCCCGctttt                                            >  2:11547479/1‑137 (MQ=18)
                                                                                         gTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTG‑CCGG‑TG‑T‑T‑CGCGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTTCCATTAAGCGGGCCTCTCGCTACAATGACGGCCCGGTTTttgtt                                        >  2:16147671/1‑137 (MQ=18)
                                                                                          tGGGAGCGGGTGTGCCCGAGAAAAGGGCGCCGCGGGGCGTGACACT‑GGCTTTT‑CCGG‑TT‑T‑TCG‑CGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgttg                                       <  2:16683899/137‑1 (MQ=21)
                                                                                          tGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGTGGTGCATGGCACC‑GGCTTTG‑CCGG‑TG‑T‑TCG‑CGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgttg                                       <  2:19275366/137‑1 (MQ=21)
                                                                                            ggAGTGGGCGTTCCTGTAAAACAGGCTCTGCGGTGCAGGACCCT‑GTTTTTT‑TTGG‑TG‑T‑T‑T‑GGGGACACC‑C‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGcc                                     <  1:3802861/137‑1 (MQ=255)
                                                                                               gCGGGTGTGTCCGCGAAATAGGCGCCGCGGTGCATGGAACC‑GGCTTTG‑CCGG‑TG‑T‑TCG‑CGGG‑CA‑C‑GCCCGCTCCCACAGAATGATTGCCATTAAGCGGGACTATCGCTACAGTGACGGCCCGTTTTTTGTTGCATCg                                  >  2:19213887/1‑137 (MQ=21)
                                                                                                tgtttttGCCCGTGAAACAGGCCCCGTGGGGCGCGACACC‑GGCTTTG‑CCGG‑TG‑T‑T‑CGAGGA‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTAGCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGc                                 <  2:10692027/131‑1 (MQ=19)
                                                                                                ggggTGTGCCGGCGAAAAGGGCGCCGTGGTACAGGGCACT‑GGTTTTG‑CCGG‑TG‑T‑T‑C‑GGGGACA‑C‑TCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGc                                 <  1:13431228/136‑1 (MQ=255)
                                                                                                   gtgTGCCTTCGAAACGTGCGCCGTGGTGCGTGGCACT‑GTCTTTG‑CTGG‑TG‑T‑T‑CGTGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTaa                              <  1:2380521/135‑1 (MQ=19)
                                                                                                        cccGCGAAACGTGCGCCGCGGTGCATGGCACC‑GTTTTTG‑CTGG‑TG‑T‑TCG‑TGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGTCTCGCTAACGAAt                         <  2:9419672/137‑1 (MQ=25)
                                                                                                        cccGAGAAAAGGCCGCCGGGGGGCATGGAACC‑GGTTGTG‑CCGG‑TG‑T‑TAG‑GGGG‑CA‑C‑TCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAAt                         <  1:17764803/137‑1 (MQ=255)
                                                                                                          cgtgAAACAGGCGCCATGGTGCATGGCACC‑GGTTTTT‑CTGG‑TG‑T‑TCG‑TGGG‑CA‑T‑TCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  2:8974696/134‑1 (MQ=255)
                                                                                                          cgcgAAACGGTCGCCGTGGTGCATGGCACC‑GTCTTTG‑CCGT‑TT‑T‑TCG‑CGGA‑CA‑C‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACNGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  1:18811845/137‑1 (MQ=255)
                                                                                                          cgcgAAACAGGTGCCGAGGTGAGTGTCACA‑GGCTTTG‑CCGG‑TG‑T‑TCG‑TGGA‑CA‑C‑TCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  2:7664033/137‑1 (MQ=255)
                                                                                                          cgcgAAACAGGCGCCGCGGTGCATGGCACT‑GGTTTTG‑CCGG‑TG‑T‑TCG‑CGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  2:14559781/137‑1 (MQ=25)
                                                                                                          cgcgAAAAAGGCGCCGCGGTACATGGCACC‑GGCTTTT‑CTGG‑TG‑T‑TCG‑CGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  1:18364348/137‑1 (MQ=25)
                                                                                                          cgcaaaaCAGACCCCGCGTGGAATGGCACC‑GGTTTTG‑CCGG‑TG‑T‑TCG‑TGCG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  2:702107/133‑1 (MQ=21)
                                                                                                          cgcaAAACAGGCGCCACGGTGCATGGCACC‑GGCTTTG‑CCGG‑TG‑T‑TCG‑TGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  1:13625716/137‑1 (MQ=25)
                                                                                                          cgcaAAAAAGGCGCCGCGGTGAATGACACC‑GGTTTTG‑CTGG‑TT‑T‑T‑TGAGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  2:6432181/137‑1 (MQ=255)
                                                                                                          cgagAAAAAGGCGTCGTGGTGCAGGGCACT‑GGTTTTG‑CTGG‑TG‑T‑TCG‑CGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  1:16714405/134‑1 (MQ=25)
                                                                                                           gcgAAACGGGCGCCGCGGTGCATGGCACC‑GGCTTTG‑CCGG‑TG‑T‑TCG‑CGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                      <  2:13636591/137‑1 (MQ=25)
                                                                                                           gcgAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTG‑CCGG‑TG‑T‑TCG‑CGGG‑CA‑C‑ACCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                      <  2:18740373/137‑1 (MQ=25)
                                                                                                           gcgAAACAGGCGCCGCGGTGCAGGGCACC‑GGCTTTT‑CCGG‑TG‑T‑TCG‑CGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                      <  2:11731394/137‑1 (MQ=37)
                                                                                                           gcgAAACAGGCGCCGCGGGGCATGGCACC‑GGCTTTG‑CCGG‑TG‑T‑T‑CGCGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                      <  2:19466132/137‑1 (MQ=37)
                                                                                                           gcaaaaCAGGCCCCGCGGTACAGGGCACC‑GGCTTTC‑CCGG‑TT‑T‑TCG‑CGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                      <  2:5654137/134‑1 (MQ=25)
                                                                                                               aaCAGGCGCCGCGGTGCATGGCACC‑GGCTTTG‑CCGG‑TG‑T‑TCG‑CGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       >  1:540212/1‑132 (MQ=37)
                                                                                                               aaCAGGCGCCGCGGTGCATGGCACC‑GGCTTTG‑CCGG‑TG‑T‑T‑CGCGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  2:540212/132‑1 (MQ=25)
                                                                                                                     cgccgtGGTACGTGGCACC‑GGCTTTT‑CCGG‑TT‑T‑TCG‑CGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa            <  2:13136217/137‑1 (MQ=25)
                                                                                                                     ccccctggTGCATGG‑CAC‑GGGTTTGACCGG‑TG‑TCC‑G‑CGGG‑CA‑C‑ACCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa            <  1:7484363/131‑1 (MQ=255)
                                                                                                                       ccgcGGTGCATGGCACC‑GGCTTTG‑CCGG‑TG‑T‑TCG‑CGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGa          >  1:10129997/1‑137 (MQ=37)
                                                                                                                               cATGGCACC‑GGTTTTG‑CCGG‑TG‑T‑TCG‑CGGA‑CA‑C‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc  <  1:6737480/137‑1 (MQ=255)
                                                                                                                               cATGGCACC‑GGCTTTG‑CCGG‑TG‑T‑TCG‑CGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc  <  2:12693782/137‑1 (MQ=37)
                                                                                                                                     ||                                                                                                                                          
CGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCAC‑GGCTTTG‑CCGG‑TG‑T‑T‑T‑CGGG‑CA‑C‑G‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATC  >  AE015451_phaM1/1499347‑1499608

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