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.

TTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATC  >  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)
        cGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTTGCACCACAGCAGGACTGAGGCTGCCAAACGGGGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGgtgt                                                                                                                       <  2:4264074/137‑1 (MQ=38)
           tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGTGTGTCCGCGAAACAGGCGCCGCGGTGCATGGCCCCGGCTTTGCCGGTGttcg                                                                                                                   >  2:10847414/1‑135 (MQ=255)
           tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGCCCGCGAAACAAGAGCCACGGTGCATGGCCCCGGCTTTTGCAGTGtttt                                                                                                                  >  1:10837755/1‑136 (MQ=255)
           tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGATGACCACGAAACAAGCACCACAGTGCATGGCACCGACTTTGTCCGTTTTTc                                                                                                                 >  1:1138912/1‑137 (MQ=255)
             aaGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGCCCGCGAAACAAGCGCCACGGTGCACGGCACCAGCTTTGCCGGTGttcgcg                                                                                                                 >  1:11774809/1‑133 (MQ=255)
             aaGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGCCCACGAAACAGGCGCCGCGGTGTATGGCACCGGCTTTGCCGGTGttcgcg                                                                                                                 >  2:15529424/1‑133 (MQ=38)
             aaGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGGGAAACAGGGGCCGCGGTGTATGGCACCGGCTTTTCCGGTGtttgcg                                                                                                                >  1:13578118/1‑134 (MQ=38)
                     gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCACGGTGCCCGGCACCGGCTTTGCCGGTGTTCGCAGGCAcgcc                                                                                                      >  2:8178206/1‑137 (MQ=25)
                     gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCACGGTGCATGGCACCGGCTTTGTCGGTGTTCGCGGGCAcgcc                                                                                                      >  2:254030/1‑137 (MQ=25)
                     gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCGAAACAGGCACCGCAGTGTATGGCACCGGCTTTGCCGGTGTTCGCGGGCAcgcc                                                                                                      >  1:12830544/1‑137 (MQ=25)
                       tGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGTGTGCCCCCGAAACAAGAGCCACAGTGCATGGCCCCCGCTTTTTCTGTGTTCTCGGGAAcgcccc                                                                                                    >  1:1591178/1‑135 (MQ=255)
                       tGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGTCCGCGAAACAAACGCCACGGTGCCCGGCACCCGCTTTTTCTGTGTTTTCGGGCAcacccc                                                                                                    >  1:7210195/1‑135 (MQ=255)
                       tGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAcgcccg                                                                                                    >  2:8600903/1‑135 (MQ=25)
                           gCCATTCAGAGCTTGAAAGAAGCTTTGGCACCACAGCGGGATGGGGGCTGACAAGCGGGGGGAGTGGGCGTTCCCGAGAAAAGGGCGCCGCGGTGCATGGCACCGGCTTTGCTGGTGTTCGCGGGCACTCCCGCTcc                                                                                                 <  1:9136802/137‑1 (MQ=12)
                            ccATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTccc                                                                                                >  2:2146854/1‑137 (MQ=21)
                            ccATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCGCGGTTCACGGCACCGGCTTTGCCAGTGTTCGCAGGCACGCCCGCTccc                                                                                                >  2:10757387/1‑137 (MQ=18)
                              aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCcac                                                                                              >  2:8386896/1‑137 (MQ=21)
                              aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCATGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCcac                                                                                              >  2:11755414/1‑137 (MQ=21)
                                 cAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCAAGAAACAGGAGCCGCGGGGCATGGCACCGGATTTGCCAGTGTTCGCGGGCACGCCCGCTCCCACAGa                                                                                           >  1:2601955/1‑137 (MQ=21)
                                       cTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGCCCCGGGTTTGTCGGTGTTCGCGGGCCCGTCCGCTCCCACAGACTGGtt                                                                                     >  2:7634847/1‑137 (MQ=21)
                                         gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGACCGATTGc                                                                                   >  2:5207217/1‑137 (MQ=18)
                                         gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCACGGCACCGGCTTTGCCAGTGTTCGCGGGCACGCCCACTCCCACAGAATGATTGc                                                                                   >  2:16211417/1‑137 (MQ=18)
                                         gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCAAAACAGGAGCCGCAGTGTATGGAACCGGATTTGCCAGTGTTCGCAGGCACGCCCCCTCCCACAGAATGGTTga                                                                                   >  1:14015697/1‑136 (MQ=16)
                                         gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAAAAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCACTCCCACAGACTGATTGc                                                                                   >  2:15649237/1‑137 (MQ=21)
                                         gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGACCGCGAAAAAGGCACCGCGGTGCATGGCACCAGCTTTGGCGGTTTTCACGGGCCCGCCCGATCCCACAGGCTGATTGc                                                                                   >  2:11411147/1‑137 (MQ=255)
                                             gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGTGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGCACCAGCTTTTCCGGTGTTCTCGGGCACACCCACTCCCCCAAACCGATTTCCAtt                                                                               >  1:8251449/1‑137 (MQ=255)
                                             gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGTGTGCCCACGAAACAAGAGCCGCGGTGTATGGCACCAGCCTTTCCTGTGTTTTCGGGCACACCCACTCCCAAAAACTGATTGCCAtt                                                                               >  1:3842728/1‑137 (MQ=255)
                                             gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGGCCGCGAAACAAGAGCCGCGGTGTATTGCACCAGCTTTGCCGGTGTTTGCGGGACCGCCCACTCCCACAGAATGAGTGACAtt                                                                               >  2:14493285/1‑137 (MQ=255)
                                             gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAAACTGATTGCCAtt                                                                               >  2:17476067/1‑137 (MQ=21)
                                             gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCACGGTGCACTGCACCAGCTTTTCCGGTGTTCGCGGGCACGCCCCCTCCCACAAAATGATTGCCAtt                                                                               >  2:8098455/1‑137 (MQ=18)
                                             gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGAGCCGAGGTGCACGGCAACAGCTTTGCCGGTGTTCGTGGGCACGCCCGCTCCCACAGAATGAGTGCCAtt                                                                               >  2:14394106/1‑137 (MQ=18)
                                             gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCAAAAAAAGCGCCGCGGTGCATGGCACCGGCTTTTCCGGTGTTTGCGGGCACACCCGCTCCCACAGAATGATTGCCAtt                                                                               >  2:12625641/1‑137 (MQ=18)
                                             gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGTCCGCAAAACAGGAGCCACGGTGCATGGCACCGGCTTTTCCGGTGTTTTCGGGCACGCCCACTCCCACCGACGGATTGCCCtt                                                                               >  1:17232763/1‑137 (MQ=21)
                                             gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCGAAACAAGAGCCACAGTGCACGGCACCAGATTTGCCAGTGTTTGCGGGCACGCCCACTCCCACAAAATGATTGCCact                                                                               >  2:6169870/1‑135 (MQ=255)
                                             gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCACGAAACAAGCACCACGGTGTCTGGCACCAACTTTTCCGGTGTTTTTGGGCACACCCACTCCCCCAAAATGATTGCCAtt                                                                               >  1:5430246/1‑137 (MQ=255)
                                                 ccTTGACCCCACAACAGGCCTGAGGCTGCCAAACGGTGGGGAGGGGCGTACCCGCAAAAAGTACCCCGCGGTACGTAACACCGGCTTGGCTGGGGTCAGTGGGCACGCCCGCTCCCACAGACTGATTACCATTAAGc                                                                           <  2:3107489/137‑1 (MQ=255)
                                                   ttGGCACCACAGCGGGACTGAGGCTGCAAAAAGGTGGGAGCGGGTGTGCCAGAGAAACGGGCGCTGTGGTGCGGGGCACCGGTTTTGCTGGTGTTCGCGGGAACGCCCGCTCCCACAGACTGATTGCCATTAAGCgg                                                                         <  2:10918848/137‑1 (MQ=255)
                                                   ttGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGATTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGAATGATTGCCATTAAGCgg                                                                         >  2:1846857/1‑137 (MQ=21)
                                                   ttGGAACCACAGCAGGACTGAGGTTGCCAACCGGTGGGAGTGGGCGTACCCGCAAAACAGGCGCCGTGGTGCATGGCACTGGTTTTGCCGGTGTTTGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCgg                                                                         <  2:17361721/137‑1 (MQ=21)
                                                   tgggCACCACAGCAGGACTGAGGCTGCCAAATTGTGGGAGCGGGCGTGCCTGCGAAACAGTCGCCGCGGTGCGTGGCACCGGTTTTGCAGTTGTTCGTGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCgg                                                                         <  2:3475749/135‑1 (MQ=18)
                                                    tGGCACCACAGCGGGATTGAGGTTGCCAAATGGTGGGAGCGGTTGTTCCCGCGAAACGGGCGCCGTGGGGCAGGGAACTGGCTTTGCTGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCggg                                                                        <  2:1833063/137‑1 (MQ=16)
                                                      gaaccacAGAAGGATGGGGGCTACAAAATGGGGGGAGCGGGTGTTCCCGAGAAACAGGCGCCGTGGGGCGTGGCACTGGTTTTTCCGGTTTTCGCGGGCACCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                      <  1:9341229/135‑1 (MQ=255)
                                                      gCACCACAGCAGTACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACGGGCGCCGTGGTACATGGCACTGGCTTTTCCGGTGTTCGCGGGCACTCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                      <  2:11886665/137‑1 (MQ=18)
                                                      gCACCACAGCAGGATTGAGGCTGCAAAATGGTGGGGGCGGGTGTGCCCGCGAAAAGGGCGCCGCGGTGCATGGCACCGGTTTTTCTGTTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                      <  2:8899140/137‑1 (MQ=16)
                                                      gCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGAGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                      <  2:16713815/137‑1 (MQ=21)
                                                      gCACCACAGAAGGGTGGGTGCGGCAAAATGGGGGGAGTGGGAGTTCCCGAGAAACAGGCGCCGGGGGGTGTGTCACTGGTTTTGTTGTTGTTCGTGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                      <  1:2070509/137‑1 (MQ=255)
                                                      gCACCACAGAAGGATTGAGGCTGCAAACCGGTGGGAGGGGGCGTTCCCGAGAAAAGGGCGCCACGGGGCGGGGCACTGGTTTTGCTGGGGTTTGCGGGCACTCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                      <  1:12384849/137‑1 (MQ=255)
                                                      acaccccAGCAGGACTGGGTTTGCCAAACGTGGGGATCGGTTTTCCCTGCGAAACGGGCGCCGTGGTGCGTGCCACCGTTTTTACCGGTGTTCGCGGGCACACCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                      <  1:6619821/136‑1 (MQ=255)
                                                      acaccacAGAAGGGTGGAGGTTGCAAAATGGGGGGAGGGGGTGTTCACGAAAAAAGTGCGCCGTGGGGCAGGGCACTGGTTTTGCGGGTGTTCGTGGGCACCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                      <  2:5584688/136‑1 (MQ=255)
                                                      aaacaacaGAAGGATGGGGGCGGCAAAATGGGGGGGGGGGGGGGGCCCGCGAAAAGGGCGCGGTGGTGAGGGGCATTGGTTTTGCGGGTGTTCGTGGACACCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                      <  1:3337782/132‑1 (MQ=255)
                                                              gcagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCACGGTGCCTGGCACCGGCTTTGTCGGTGTTCGAAGGCACACCCGCTCCCACAGGATGATTGACATTAAGCGGGACTATCGCTa                                                              >  2:1312640/1‑137 (MQ=18)
                                                              gcagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAAGCGCCACGGTGCATGGCACCGGCTTTGTCGGTGTTTACGGGCACGCCCGCTCCCACAGACCGATTGCCATTAAACGGGACTATCGccc                                                              >  2:1020820/1‑135 (MQ=18)
                                                                   aCTGAGGCTGCCAAACGGTGGGAGCGGGGGTGCCCACAAAACAAGAGCCACGGTGTATGGCACCAGATTTTTCGGTGTTCACGGGCACACCCACCCCCCCAAACTGATTGCCATTAAACGGGCCCATCGCTGCAGTg                                                         >  1:4112236/1‑137 (MQ=255)
                                                                     tGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCGAAACAGGCGCCACAGTGCACGGCACCAGCTTTGCCGGTGTTTGCGGGCACAACCGCTTCCACAGACTGATTGCCACTCAGCAGGCCTACCACTACAATGAc                                                       >  2:1346109/1‑137 (MQ=255)
                                                                     ggAGGCTGCAAAATGGTGGGGGGGGGTGTGCCCGAAAAACGGGCGCCGCGGTGCATGGCATTGGTTTTTCTGGTGTTCGAGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAc                                                       <  2:14789684/136‑1 (MQ=255)
                                                                      gAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCAGTTCACGGCACCGGCTTTGCCGGTTTTCGCGGGCCCGCCCGCTCCCACCGACTGATTGCCATTAAGCGGGCCTATCGCCACAGTGAca                                                      >  2:9880844/1‑136 (MQ=18)
                                                                          cTGCCAAACGGTGGGAGCGGGCGTTCCCGCGAAACAGACGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGTGGGCACGCCCGCTCCCAAAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGcac                                                  >  2:7252958/1‑135 (MQ=21)
                                                                          cTGCCAAACGGTGGGAGCGGGCGTGTCCGCGAAACAGGCGCCGCGGTGCCTGGCCCCGGCTTTGCCGGTGTTTGCGGGGACGCCCGCTTCCACAAACTGATTGCCATTAAGCGGGCCTATTGCTACAGCGACGGccc                                                  >  2:6722324/1‑137 (MQ=37)
                                                                          cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGAATGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGccc                                                  >  2:17104684/1‑137 (MQ=18)
                                                                          cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCACGGCACCGGATTTGCCAGTGTTCGCGGGCACGCCCGCTCCCACAAAATGATTGCAATTCAACGGGCCTATCGCCACAATGACGGccc                                                  >  2:12508491/1‑137 (MQ=17)
                                                                          cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAGGCGCCACGGTGCATGGCACCGGATTTGTCGGTGTTTGCGGGCACACCCACTCCCACAGACTGGTTTCCCTTCAGCGGAACTGTTGCTGCAGTGAAGGgcc                                                  >  1:18470448/1‑137 (MQ=9)
                                                                               aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCCCAGACTGATTGCCCTTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttt                                             >  2:18772385/1‑137 (MQ=21)
                                                                               aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCACGtttt                                             >  2:9333810/1‑137 (MQ=21)
                                                                               aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTCAGCGGGCCTCTCGCCACAGTGACGGCCCGtttt                                             >  2:13888989/1‑137 (MQ=18)
                                                                               aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCGCGGTGCATGGCCCCGGCTTTGCCGGTGTTCGCGGGCACACCCGCTCCCCCAGAATGATTGCCATTCAGCGGGCCTATCGCTACCGTGACGGACCGtttt                                             >  1:4223472/1‑137 (MQ=18)
                                                                               aaaCGGTGGGAGCGGGCGTGCCCGCAAAACAGGCGCCGCGGTGCCCGGCACCGGCTTTGCCGGTGTTCGCGGGCACACCCGCTCCCAAAGAATGATTGCCATTAAGTGGGCCTATCGCTACAATGACGGCCCGtttt                                             >  2:2726313/1‑137 (MQ=18)
                                                                                 aCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACCGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt                                           >  2:15937897/1‑137 (MQ=18)
                                                                                   tttGGGAGCGGGCGTCCCCGCGAAATAGGCGTCGTGGGGCAGGGCACTGGTTTTGCTGGTGTTCGCGGGCACTCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgt                                         <  1:5784409/135‑1 (MQ=255)
                                                                                        gAGCGGGCGTGCCCGCGAAACAAGCGTCCCAGTGCCTGGCCCCAGCTTTGTCCGTGTTCGCCGGCACACCCGCTCCCACCAACTGATTGTCATTAAGCGGGCCTCCCTCTCCAGTGGCGGACCGTTTTTTGTTGCCt                                    >  2:400945/1‑137 (MQ=9)
                                                                                          gCGGGCGTGCCCACGAAAAAAGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTTTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTCAGCGGGCCTCTCGCTACAGTTGCGGCCCGTTTTTTGTTGCCTCg                                  >  2:18213356/1‑137 (MQ=18)
                                                                                            gggTGTGCCTGAGAAACAGGCGCCGCGGTGCGTGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCt                                <  2:16479936/137‑1 (MQ=25)
                                                                                                    ccGCGAAAAGGGCGCCGTGGTGCGTGGCATTGTTTTTGCGGGTGTTCGTGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATc                        <  2:15865049/137‑1 (MQ=255)
                                                                                                    ccGCGAAAAAGGCGCCGGGGTGCGTGGAACCGGTTTTCCTGGTGTTCGTGGGAACACCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATc                        <  2:1423918/137‑1 (MQ=255)
                                                                                                     tgcgAAACAGGTTCCGTGGGGCAGGGCACTGTTTTTGTGGGTGTTTGGGGGCACCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  2:12828264/136‑1 (MQ=255)
                                                                                                     cgcgAAACGGGCGCCGCGGTGCATGGCACCGGTTTTTCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  2:6310582/137‑1 (MQ=25)
                                                                                                     cgcgAAACAGGTGCCGCGGTGCATGGCACCGTCTTTACCGTTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  2:1977617/137‑1 (MQ=37)
                                                                                                     cgcgAAACAGGCGCCGCGTTGCATGGCACCGGTTTTGCTGGTGTTCGAGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  2:15171002/137‑1 (MQ=25)
                                                                                                     cgcgAAACAGGCGCCGCGGTGCATGACACCGGCTTTTCCGTTGTTCGCGGGCACCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  2:15655868/137‑1 (MQ=255)
                                                                                                     cgcaAAACGGGCGCCGCGGTGCATGCCACCGTTTTTTCCGTTGTTCTCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  1:18116023/137‑1 (MQ=255)
                                                                                                     cgagAAACGGGCTCCGTGGTGCATGGAACTGGTTTTGCTGGTGTTCGGGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  2:1234948/134‑1 (MQ=25)
                                                                                                     cgagAAACAGGTGCCGCGGTGAGGGGCACCGGTTTTGTTGGGGTTTGTGGGCATGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  1:10226358/134‑1 (MQ=255)
                                                                                                     cacgAAACAGTCGCCGCGGTGCATGACACCGGCTTTGTCGTTTTCCGCGGGCACTCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  1:8593520/135‑1 (MQ=255)
                                                                                                     agcgAAACAGGCGCCGTGGTGCATGGCATCGGTTTTGCTGGTGTTTGCGGGCACCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  1:12434153/136‑1 (MQ=255)
                                                                                                     agagAAACGGGCGCTGGGGTGCAGGGCACTGTTTTTGCTGGTGTTCGCGGGCACCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  1:14902472/134‑1 (MQ=255)
                                                                                                      gcgAAACGGGCGCTGTGGGGCAGGGCACTGGCTTTGCTGGTTTTCGCGGGCACCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                      <  1:15761641/137‑1 (MQ=255)
                                                                                                      gcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                      <  2:17732192/137‑1 (MQ=25)
                                                                                                                cttcgcgGGGCAGGGCACTGGTTTTGCTGGTGTTCGCGGACACCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa            <  1:10923403/134‑1 (MQ=255)
                                                                                                                ctccgcGGTGCATGGCACCGGCTTTTCCGGTGTTCGTGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa            <  2:6490289/135‑1 (MQ=37)
                                                                                                                ctccgcGGTACGTGGCACTGGCTCTGCCGGTTTTCGGGGGCACTCCCGCTCACACAGACTGAGTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa            <  1:5500727/135‑1 (MQ=255)
                                                                                                                cgccgtGGTGCGTGACACCGTCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa            <  2:7882999/137‑1 (MQ=37)
                                                                                                                cgccgcGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa            >  2:10681648/1‑137 (MQ=37)
                                                                                                                caccgcGGTGCATGGCACCGGCTTGTCCGGGGTTCGTGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa            <  2:7414084/135‑1 (MQ=25)
                                                                                                                     cgGTGCATGGCACCTGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACg       <  2:13705471/137‑1 (MQ=37)
                                                                                                                     cgGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACg       <  2:17891435/137‑1 (MQ=37)
                                                                                                                                ||                                                                                                                                     
TTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATC  >  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: