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,4800CA95.2% 449.9 / NA 207intergenic (+57/‑107)trpS CDS (Tryptophan‑‑tRNA ligase)/zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain)–/–
Reads supporting (aligned to +/- strand):  ref base C (10/0);  new base A (101/96);  total (111/96)
Fisher's exact test for biased strand distribution p-value = 1.97e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*AE015451_phaM11,499,4810AC97.5% 480.2 / NA 207intergenic (+58/‑106)trpS CDS (Tryptophan‑‑tRNA ligase)/zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain)–/–
Reads supporting (aligned to +/- strand):  ref base A (4/0);  major base C (107/91);  minor base T (0/5);  total (111/96)
Fisher's exact test for biased strand distribution p-value = 2.24e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 207 total aligned reads displayed.

TTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTG‑T‑T‑T‑CGGGCAC‑G‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATC  >  AE015451_phaM1/1499352‑1499608
                                                                                                                                ||                                                                                                                                     
ttCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCttt                                                                                                                                <  2:5230731/137‑1 (MQ=255)
 tCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTg                                                                                                                               <  2:10421838/137‑1 (MQ=255)
 tCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTg                                                                                                                               <  1:9372552/137‑1 (MQ=255)
   gACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGcc                                                                                                                             <  2:6937/137‑1 (MQ=255)
     cAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCgg                                                                                                                           <  2:4894192/137‑1 (MQ=255)
       aCGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGgtg                                                                                                                         <  2:1069215/137‑1 (MQ=255)
        cGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTTGCACCACAGCAGGACTGAGGCTGCCAAACGGGGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGgtg‑t                                                                                                                       <  2:4264074/137‑1 (MQ=38)
           tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGTGTGTCCGCGAAACAGGCGCCGCGGTGCATGGCCCCGGCTTTGCCGGTG‑t‑tcg                                                                                                                   >  2:10847414/1‑135 (MQ=255)
           tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGCCCGCGAAACAAGAGCCACGGTGCATGGCCCCGGCTTTTGCAGTG‑t‑t‑tt                                                                                                                  >  1:10837755/1‑136 (MQ=255)
           tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGATGACCACGAAACAAGCACCACAGTGCATGGCACCGACTTTGTCCGTT‑T‑T‑T‑c                                                                                                                 >  1:1138912/1‑137 (MQ=255)
             aaGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGCCCGCGAAACAAGCGCCACGGTGCACGGCACCAGCTTTGCCGGTG‑t‑tcgcg                                                                                                                 >  1:11774809/1‑133 (MQ=255)
             aaGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGCCCACGAAACAGGCGCCGCGGTGTATGGCACCGGCTTTGCCGGTG‑t‑tcgcg                                                                                                                 >  2:15529424/1‑133 (MQ=38)
             aaGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGGGAAACAGGGGCCGCGGTGTATGGCACCGGCTTTTCCGGTG‑t‑t‑tgcg                                                                                                                >  1:13578118/1‑134 (MQ=38)
                     gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCACGGTGCCCGGCACCGGCTTTGCCGGTG‑T‑TCG‑CAGGCAc‑g‑cc                                                                                                      >  2:8178206/1‑137 (MQ=25)
                     gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCACGGTGCATGGCACCGGCTTTGTCGGTG‑T‑TCG‑CGGGCAc‑g‑cc                                                                                                      >  2:254030/1‑137 (MQ=25)
                     gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCGAAACAGGCACCGCAGTGTATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGCAc‑g‑cc                                                                                                      >  1:12830544/1‑137 (MQ=25)
                       tGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGTGTGCCCCCGAAACAAGAGCCACAGTGCATGGCCCCCGCTTTTTCTGTG‑T‑TCT‑CGGGAAc‑g‑cccc                                                                                                    >  1:1591178/1‑135 (MQ=255)
                       tGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGTCCGCGAAACAAACGCCACGGTGCCCGGCACCCGCTTTTTCTGTGTT‑T‑T‑CGGGCAc‑a‑cccc                                                                                                    >  1:7210195/1‑135 (MQ=255)
                       tGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGCAc‑g‑cccg                                                                                                    >  2:8600903/1‑135 (MQ=25)
                           gCCATTCAGAGCTTGAAAGAAGCTTTGGCACCACAGCGGGATGGGGGCTGACAAGCGGGGGGAGTGGGCGTTCCCGAGAAAAGGGCGCCGCGGTGCATGGCACCGGCTTTGCTGGTG‑T‑TCG‑CGGGCAC‑TCCCGCTcc                                                                                                 <  1:9136802/137‑1 (MQ=12)
                            ccATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGCAC‑GCCCGCTccc                                                                                                >  2:2146854/1‑137 (MQ=21)
                            ccATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCGCGGTTCACGGCACCGGCTTTGCCAGTG‑T‑TCG‑CAGGCAC‑GCCCGCTccc                                                                                                >  2:10757387/1‑137 (MQ=18)
                              aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGCAC‑GCCCGCTCCcac                                                                                              >  2:8386896/1‑137 (MQ=21)
                              aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCATGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGCAC‑GCCCGCTCCcac                                                                                              >  2:11755414/1‑137 (MQ=21)
                                 cAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCAAGAAACAGGAGCCGCGGGGCATGGCACCGGATTTGCCAGTG‑T‑TCG‑CGGGCAC‑GCCCGCTCCCACAGa                                                                                           >  1:2601955/1‑137 (MQ=21)
                                       cTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGCCCCGGGTTTGTCGGTG‑T‑TCG‑CGGGCCC‑GTCCGCTCCCACAGACTGGtt                                                                                     >  2:7634847/1‑137 (MQ=21)
                                         gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGCAA‑GCCCGCTCCCACAGACCGATTGc                                                                                   >  2:5207217/1‑137 (MQ=18)
                                         gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCACGGCACCGGCTTTGCCAGTG‑T‑TCG‑CGGGCAC‑GCCCACTCCCACAGAATGATTGc                                                                                   >  2:16211417/1‑137 (MQ=18)
                                         gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCAAAACAGGAGCCGCAGTGTATGGAACCGGATTTGCCAGTG‑T‑TCG‑CAGGCAC‑GCCCCCTCCCACAGAATGGTTga                                                                                   >  1:14015697/1‑136 (MQ=16)
                                         gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAAAAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGCAC‑GCCCACTCCCACAGACTGATTGc                                                                                   >  2:15649237/1‑137 (MQ=21)
                                         gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGACCGCGAAAAAGGCACCGCGGTGCATGGCACCAGCTTTGGCGGTT‑T‑TCA‑CGGGCCC‑GCCCGATCCCACAGGCTGATTGc                                                                                   >  2:11411147/1‑137 (MQ=255)
                                             gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGTGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGCACCAGCTTTTCCGGTG‑T‑TCT‑CGGGCACAC‑CCACTCCCCCAAACCGATTTCCAtt                                                                               >  1:8251449/1‑137 (MQ=255)
                                             gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGTGTGCCCACGAAACAAGAGCCGCGGTGTATGGCACCAGCCTTTCCTGTGTT‑T‑T‑CGGGCACAC‑CCACTCCCAAAAACTGATTGCCAtt                                                                               >  1:3842728/1‑137 (MQ=255)
                                             gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGGCCGCGAAACAAGAGCCGCGGTGTATTGCACCAGCTTTGCCGGTG‑T‑T‑TGCGGGACC‑GCCCACTCCCACAGAATGAGTGACAtt                                                                               >  2:14493285/1‑137 (MQ=255)
                                             gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑T‑CGCGGGCAC‑GCCCGCTCCCACAAACTGATTGCCAtt                                                                               >  2:17476067/1‑137 (MQ=21)
                                             gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCACGGTGCACTGCACCAGCTTTTCCGGTG‑T‑TCG‑CGGGCAC‑GCCCCCTCCCACAAAATGATTGCCAtt                                                                               >  2:8098455/1‑137 (MQ=18)
                                             gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGAGCCGAGGTGCACGGCAACAGCTTTGCCGGTG‑T‑TCG‑TGGGCAC‑GCCCGCTCCCACAGAATGAGTGCCAtt                                                                               >  2:14394106/1‑137 (MQ=18)
                                             gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCAAAAAAAGCGCCGCGGTGCATGGCACCGGCTTTTCCGGTG‑T‑T‑TGCGGGCACAC‑CCGCTCCCACAGAATGATTGCCAtt                                                                               >  2:12625641/1‑137 (MQ=18)
                                             gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGTCCGCAAAACAGGAGCCACGGTGCATGGCACCGGCTTTTCCGGTGTT‑T‑T‑CGGGCAC‑GCCCACTCCCACCGACGGATTGCCCtt                                                                               >  1:17232763/1‑137 (MQ=21)
                                             gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCGAAACAAGAGCCACAGTGCACGGCACCAGATTTGCCAGTG‑T‑T‑TGCGGGCAC‑GCCCACTCCCACAAAATGATTGCCact                                                                               >  2:6169870/1‑135 (MQ=255)
                                             gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCACGAAACAAGCACCACGGTGTCTGGCACCAACTTTTCCGGTGTT‑T‑T‑TGGGCACAC‑CCACTCCCCCAAAATGATTGCCAtt                                                                               >  1:5430246/1‑137 (MQ=255)
                                                 ccTTGACCCCACAACAGGCCTGAGGCTGCCAAACGGTGGGGAGGGGCGTACCCGCAAAAAGTACCCCGCGGTACGTAACACCGGCTTGGCTGGGG‑TCA‑G‑TGGGCAC‑GCCCGCTCCCACAGACTGATTACCATTAAGc                                                                           <  2:3107489/137‑1 (MQ=255)
                                                   ttGGCACCACAGCGGGACTGAGGCTGCAAAAAGGTGGGAGCGGGTGTGCCAGAGAAACGGGCGCTGTGGTGCGGGGCACCGGTTTTGCTGGTG‑T‑T‑CGCGGGAAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCgg                                                                         <  2:10918848/137‑1 (MQ=255)
                                                   ttGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGATTTGCCGGTG‑T‑TCG‑CGGGCAC‑GCCCGCTCCCACAGAATGATTGCCATTAAGCgg                                                                         >  2:1846857/1‑137 (MQ=21)
                                                   ttGGAACCACAGCAGGACTGAGGTTGCCAACCGGTGGGAGTGGGCGTACCCGCAAAACAGGCGCCGTGGTGCATGGCACTGGTTTTGCCGGTG‑T‑T‑TGCGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCgg                                                                         <  2:17361721/137‑1 (MQ=21)
                                                   tgggCACCACAGCAGGACTGAGGCTGCCAAATTGTGGGAGCGGGCGTGCCTGCGAAACAGTCGCCGCGGTGCGTGGCACCGGTTTTGCAGTTG‑T‑TCG‑TGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCgg                                                                         <  2:3475749/135‑1 (MQ=18)
                                                    tGGCACCACAGCGGGATTGAGGTTGCCAAATGGTGGGAGCGGTTGTTCCCGCGAAACGGGCGCCGTGGGGCAGGGAACTGGCTTTGCTGGTG‑T‑TCG‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCggg                                                                        <  2:1833063/137‑1 (MQ=16)
                                                      gaaccacAGAAGGATGGGGGCTACAAAATGGGGGGAGCGGGTGTTCCCGAGAAACAGGCGCCGTGGGGCGTGGCACTGGTTTTTCCGGTT‑T‑TCG‑CGGGCAC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                      <  1:9341229/135‑1 (MQ=255)
                                                      gCACCACAGCAGTACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACGGGCGCCGTGGTACATGGCACTGGCTTTTCCGGTG‑T‑TCG‑CGGGCAC‑TCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                      <  2:11886665/137‑1 (MQ=18)
                                                      gCACCACAGCAGGATTGAGGCTGCAAAATGGTGGGGGCGGGTGTGCCCGCGAAAAGGGCGCCGCGGTGCATGGCACCGGTTTTTCTGTTG‑T‑T‑CGCGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                      <  2:8899140/137‑1 (MQ=16)
                                                      gCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑T‑CGAGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                      <  2:16713815/137‑1 (MQ=21)
                                                      gCACCACAGAAGGGTGGGTGCGGCAAAATGGGGGGAGTGGGAGTTCCCGAGAAACAGGCGCCGGGGGGTGTGTCACTGGTTTTGTTGTTG‑T‑T‑CGTGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                      <  1:2070509/137‑1 (MQ=255)
                                                      gCACCACAGAAGGATTGAGGCTGCAAACCGGTGGGAGGGGGCGTTCCCGAGAAAAGGGCGCCACGGGGCGGGGCACTGGTTTTGCTGGGG‑T‑T‑TGCGGGCAC‑TCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                      <  1:12384849/137‑1 (MQ=255)
                                                      acaccccAGCAGGACTGGGTTTGCCAAACGTGGGGATCGGTTTTCCCTGCGAAACGGGCGCCGTGGTGCGTGCCACCGTTTTTACCGGTG‑T‑TCG‑CGGGCAC‑ACCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                      <  1:6619821/136‑1 (MQ=255)
                                                      acaccacAGAAGGGTGGAGGTTGCAAAATGGGGGGAGGGGGTGTTCACGAAAAAAGTGCGCCGTGGGGCAGGGCACTGGTTTTGCGGGTG‑T‑TCG‑TGGGCAC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                      <  2:5584688/136‑1 (MQ=255)
                                                      aaacaacaGAAGGATGGGGGCGGCAAAATGGGGGGGGGGGGGGGGCCCGCGAAAAGGGCGCGGTGGTGAGGGGCATTGGTTTTGCGGGTG‑T‑TCG‑TGGACAC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                      <  1:3337782/132‑1 (MQ=255)
                                                              gcagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCACGGTGCCTGGCACCGGCTTTGTCGGTG‑T‑TCG‑AAGGCACAC‑CCGCTCCCACAGGATGATTGACATTAAGCGGGACTATCGCTa                                                              >  2:1312640/1‑137 (MQ=18)
                                                              gcagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAAGCGCCACGGTGCATGGCACCGGCTTTGTCGGTG‑T‑T‑TACGGGCAC‑GCCCGCTCCCACAGACCGATTGCCATTAAACGGGACTATCGccc                                                              >  2:1020820/1‑135 (MQ=18)
                                                                   aCTGAGGCTGCCAAACGGTGGGAGCGGGGGTGCCCACAAAACAAGAGCCACGGTGTATGGCACCAGATTTTTCGGTG‑T‑TCA‑CGGGCACAC‑CCACCCCCCCAAACTGATTGCCATTAAACGGGCCCATCGCTGCAGTg                                                         >  1:4112236/1‑137 (MQ=255)
                                                                     tGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCGAAACAGGCGCCACAGTGCACGGCACCAGCTTTGCCGGTG‑T‑T‑TGCGGGCACAA‑CCGCTTCCACAGACTGATTGCCACTCAGCAGGCCTACCACTACAATGAc                                                       >  2:1346109/1‑137 (MQ=255)
                                                                     ggAGGCTGCAAAATGGTGGGGGGGGGTGTGCCCGAAAAACGGGCGCCGCGGTGCATGGCATTGGTTTTTCTGGTG‑T‑T‑CGAGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAc                                                       <  2:14789684/136‑1 (MQ=255)
                                                                      gAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCAGTTCACGGCACCGGCTTTGCCGGTT‑T‑TCG‑CGGGCCC‑GCCCGCTCCCACCGACTGATTGCCATTAAGCGGGCCTATCGCCACAGTGAca                                                      >  2:9880844/1‑136 (MQ=18)
                                                                          cTGCCAAACGGTGGGAGCGGGCGTTCCCGCGAAACAGACGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑TGGGCAC‑GCCCGCTCCCAAAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGcac                                                  >  2:7252958/1‑135 (MQ=21)
                                                                          cTGCCAAACGGTGGGAGCGGGCGTGTCCGCGAAACAGGCGCCGCGGTGCCTGGCCCCGGCTTTGCCGGTG‑T‑T‑TGCGGGGAC‑GCCCGCTTCCACAAACTGATTGCCATTAAGCGGGCCTATTGCTACAGCGACGGccc                                                  >  2:6722324/1‑137 (MQ=37)
                                                                          cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGCAC‑GCCCGCTCCCACAGAATGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGccc                                                  >  2:17104684/1‑137 (MQ=18)
                                                                          cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCACGGCACCGGATTTGCCAGTG‑T‑TCG‑CGGGCAC‑GCCCGCTCCCACAAAATGATTGCAATTCAACGGGCCTATCGCCACAATGACGGccc                                                  >  2:12508491/1‑137 (MQ=17)
                                                                          cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAGGCGCCACGGTGCATGGCACCGGATTTGTCGGTG‑T‑T‑TGCGGGCACAC‑CCACTCCCACAGACTGGTTTCCCTTCAGCGGAACTGTTGCTGCAGTGAAGGgcc                                                  >  1:18470448/1‑137 (MQ=9)
                                                                               aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGCAC‑GCCCGCTCCCCCAGACTGATTGCCCTTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttt                                             >  2:18772385/1‑137 (MQ=21)
                                                                               aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCACGtttt                                             >  2:9333810/1‑137 (MQ=21)
                                                                               aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑T‑CGCGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTCAGCGGGCCTCTCGCCACAGTGACGGCCCGtttt                                             >  2:13888989/1‑137 (MQ=18)
                                                                               aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCGCGGTGCATGGCCCCGGCTTTGCCGGTG‑T‑TCG‑CGGGCAC‑ACCCGCTCCCCCAGAATGATTGCCATTCAGCGGGCCTATCGCTACCGTGACGGACCGtttt                                             >  1:4223472/1‑137 (MQ=18)
                                                                               aaaCGGTGGGAGCGGGCGTGCCCGCAAAACAGGCGCCGCGGTGCCCGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGCACAC‑CCGCTCCCAAAGAATGATTGCCATTAAGTGGGCCTATCGCTACAATGACGGCCCGtttt                                             >  2:2726313/1‑137 (MQ=18)
                                                                                 aCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGCAC‑GCCCGCTCCCACAGACCGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt                                           >  2:15937897/1‑137 (MQ=18)
                                                                                   tttGGGAGCGGGCGTCCCCGCGAAATAGGCGTCGTGGGGCAGGGCACTGGTTTTGCTGGTG‑T‑TCG‑CGGGCAC‑TCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgt                                         <  1:5784409/135‑1 (MQ=255)
                                                                                        gAGCGGGCGTGCCCGCGAAACAAGCGTCCCAGTGCCTGGCCCCAGCTTTGTCCGTG‑T‑TCG‑CCGGCACAC‑CCGCTCCCACCAACTGATTGTCATTAAGCGGGCCTCCCTCTCCAGTGGCGGACCGTTTTTTGTTGCCt                                    >  2:400945/1‑137 (MQ=9)
                                                                                          gCGGGCGTGCCCACGAAAAAAGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTT‑T‑TCG‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTCAGCGGGCCTCTCGCTACAGTTGCGGCCCGTTTTTTGTTGCCTCg                                  >  2:18213356/1‑137 (MQ=18)
                                                                                            gggTGTGCCTGAGAAACAGGCGCCGCGGTGCGTGGCACCGGCTTTGCCGGTG‑T‑T‑CGCGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCt                                <  2:16479936/137‑1 (MQ=25)
                                                                                                    ccGCGAAAAGGGCGCCGTGGTGCGTGGCATTGTTTTTGCGGGTG‑T‑TCG‑TGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATc                        <  2:15865049/137‑1 (MQ=255)
                                                                                                    ccGCGAAAAAGGCGCCGGGGTGCGTGGAACCGGTTTTCCTGGTG‑T‑TCG‑TGGGAAC‑ACCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATc                        <  2:1423918/137‑1 (MQ=255)
                                                                                                     tgcgAAACAGGTTCCGTGGGGCAGGGCACTGTTTTTGTGGGTG‑T‑T‑TGGGGGCAC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  2:12828264/136‑1 (MQ=255)
                                                                                                     cgcgAAACGGGCGCCGCGGTGCATGGCACCGGTTTTTCCGGTG‑T‑T‑CGCGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  2:6310582/137‑1 (MQ=25)
                                                                                                     cgcgAAACAGGTGCCGCGGTGCATGGCACCGTCTTTACCGTTG‑T‑TCG‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  2:1977617/137‑1 (MQ=37)
                                                                                                     cgcgAAACAGGCGCCGCGTTGCATGGCACCGGTTTTGCTGGTG‑T‑TCG‑AGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  2:15171002/137‑1 (MQ=25)
                                                                                                     cgcgAAACAGGCGCCGCGGTGCATGACACCGGCTTTTCCGTTG‑T‑TCG‑CGGGCAC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  2:15655868/137‑1 (MQ=255)
                                                                                                     cgcaAAACGGGCGCCGCGGTGCATGCCACCGTTTTTTCCGTTG‑T‑TCT‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  1:18116023/137‑1 (MQ=255)
                                                                                                     cgagAAACGGGCTCCGTGGTGCATGGAACTGGTTTTGCTGGTG‑T‑T‑CGGGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  2:1234948/134‑1 (MQ=25)
                                                                                                     cgagAAACAGGTGCCGCGGTGAGGGGCACCGGTTTTGTTGGGG‑T‑T‑TGTGGGCAT‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  1:10226358/134‑1 (MQ=255)
                                                                                                     cacgAAACAGTCGCCGCGGTGCATGACACCGGCTTTGTCGTTT‑T‑CCG‑CGGGCAC‑TCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  1:8593520/135‑1 (MQ=255)
                                                                                                     agcgAAACAGGCGCCGTGGTGCATGGCATCGGTTTTGCTGGTG‑T‑T‑TGCGGGCAC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  1:12434153/136‑1 (MQ=255)
                                                                                                     agagAAACGGGCGCTGGGGTGCAGGGCACTGTTTTTGCTGGTG‑T‑TCG‑CGGGCAC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  1:14902472/134‑1 (MQ=255)
                                                                                                      gcgAAACGGGCGCTGTGGGGCAGGGCACTGGCTTTGCTGGTT‑T‑TCG‑CGGGCAC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                      <  1:15761641/137‑1 (MQ=255)
                                                                                                      gcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                      <  2:17732192/137‑1 (MQ=25)
                                                                                                                cttcgcgGGGCAGGGCACTGGTTTTGCTGGTG‑T‑TCG‑CGGACAC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa            <  1:10923403/134‑1 (MQ=255)
                                                                                                                ctccgcGGTGCATGGCACCGGCTTTTCCGGTG‑T‑T‑CGTGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa            <  2:6490289/135‑1 (MQ=37)
                                                                                                                ctccgcGGTACGTGGCACTGGCTCTGCCGGTT‑T‑TCG‑GGGGCACTC‑CCGCTCACACAGACTGAGTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa            <  1:5500727/135‑1 (MQ=255)
                                                                                                                cgccgtGGTGCGTGACACCGTCTTTGCCGGTG‑T‑T‑CGCGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa            <  2:7882999/137‑1 (MQ=37)
                                                                                                                cgccgcGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa            >  2:10681648/1‑137 (MQ=37)
                                                                                                                caccgcGGTGCATGGCACCGGCTTGTCCGGGG‑T‑TCG‑TGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa            <  2:7414084/135‑1 (MQ=25)
                                                                                                                     cgGTGCATGGCACCTGCTTTGCCGGTG‑T‑T‑CGCGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACg       <  2:13705471/137‑1 (MQ=37)
                                                                                                                     cgGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACg       <  2:17891435/137‑1 (MQ=37)
                                                                                                                                ||                                                                                                                                     
TTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTG‑T‑T‑T‑CGGGCAC‑G‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATC  >  AE015451_phaM1/1499352‑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: ‑