Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM14,889,6230AG30.8% 250.7 / 29.0 162intergenic (+69/‑84)hyi CDS (hydroxypyruvate isomerase)/glxR CDS (tartronate semialdehyde reductase)–/–
Reads supporting (aligned to +/- strand):  ref base A (62/48);  new base G (0/49);  total (62/100)
Fisher's exact test for biased strand distribution p-value = 5.09e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 162 total aligned reads displayed.

GGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGC  >  AE015451_phaM1/4889488‑4889756
                                                                                                                                       |                                                                                                                                      
ggAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAg                                                                                                                                       <  1:10043436/137‑1 (MQ=255)
  aGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGtt                                                                                                                                     <  2:16681114/137‑1 (MQ=255)
  aGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGtt                                                                                                                                     <  1:5158875/137‑1 (MQ=255)
    tACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGGGTCGCCCTCGAACACCGAAAAAAACAGTGCCACCCCCCCAGCTGt                                                                                                                                   >  2:8282693/1‑137 (MQ=255)
    tACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGt                                                                                                                                   <  1:3873067/137‑1 (MQ=255)
        aGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGc                                                                                                                               <  2:14336107/137‑1 (MQ=255)
        aGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGc                                                                                                                               <  1:15966237/137‑1 (MQ=255)
         gCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCt                                                                                                                              <  1:7757492/137‑1 (MQ=255)
         gCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCt                                                                                                                              <  1:1573883/137‑1 (MQ=255)
              tGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTTGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGAAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTcgcg                                                                                                                        <  2:17610208/137‑1 (MQ=255)
               gACCACCACAGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGGGTCGTCCACAGACACCAGAAAAACCGGTGTCACCCACCTAGTTTTcttcttttcgg                                                                                                                        >  2:13038681/1‑132 (MQ=255)
                 ccacccccGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGGGTCGCCCAAGAAAACCCGCAAAGCAGGTGCCACCCAACAAGTTGTCTTTTTCTCGGGc                                                                                                                     >  2:17934227/1‑137 (MQ=255)
                 ccaccactGAAGTGGTCTTGGGCGGGCGGAAAACCCACAACGCAATCTGAGAAACAACTAGGGACCCTGTGGGAGCGGGCTCGCCCGAGAACACCGGCAAACCCGGTGCCACCCACCGGGTTGTCTGCTTCGCGGGc                                                                                                                     <  2:13366931/137‑1 (MQ=25)
                 ccaccaccGAAGCGGGTCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGTGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGc                                                                                                                     <  2:14360751/137‑1 (MQ=255)
                 ccaccaccGAAGCGGGCCTGTGTTGGCTAAAACCCCACAACACACCCTAAGGAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGc                                                                                                                     <  2:19402334/137‑1 (MQ=25)
                 ccaccaccGAAGCGGGCCTGGGCTGGCCGAAAACCCACAACACAATCTGAGGAACACCTCGGGACCCTGGGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGGGTTGTCTGCTTCGCGGGc                                                                                                                     <  2:12363470/137‑1 (MQ=37)
                  caccacgGAGGCGGGCTGGGGTTGGCTAAAACCCCACAACGCACTCGGAGAAACAACTTGGGACTTTGTGGGAGCGGGCTCGCCCGAGAACACTGGAAAAGCTGGGGCCATCCACTGGGTTGTTTGTTTCGGGGGCg                                                                                                                    <  2:6214420/137‑1 (MQ=255)
                    ccacaGAAGAGGTCCTGGTTTGGCGGAAAACCCACAACCAAATCTGAGGAACACCTCGTGACTTTGTGGGAGCGGGCTCGCCCGCGAAAACCGGCAAGGCAGGGGCCACCCACCGAGTTGTTGGCTTCGTGGGCGAg                                                                                                                  <  2:19352718/137‑1 (MQ=255)
                       ccGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGccc                                                                                                               <  2:5511208/137‑1 (MQ=255)
                        cGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGGGTCTTCCCCGAACACCAGCCAAACCGGTGGCCTCCCCCGAGTTGTCTGCTTCTCGGGTGGGcccc                                                                                                              >  2:13730483/1‑136 (MQ=255)
                          aaGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGGGTCGGCCCCAAACCCCAGCAAAACCGGTGGCCCCCACCCAAGTTTCTGATTTACGGGAGAGCCCCCt                                                                                                            >  2:13301228/1‑137 (MQ=255)
                             cGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGGGTCGCCCGCGAACACCGGAAAAACCAGTGCACTCCCCCAAGTTGGTTTCTTCTCAGGAGAGCCCGATccc                                                                                                         >  2:14422798/1‑137 (MQ=255)
                             cGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGGATCGACCCCGAACACCGGAAAAACCAGTGTCACCCACCGAGGTGTCTTTTTTGCGGGAGAGCCCACcccc                                                                                                         >  2:14558101/1‑137 (MQ=255)
                                  cTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGGGGCACCCCCCGACACCCGAAAAACCCGTTCCCCCCCCAAAACTGTCCTCTTCTCGGGCGAGACCCCTCCcccaag                                                                                                    >  2:9041259/1‑135 (MQ=255)
                                         ggttaaaacCCCACAACGCAATCTGGGGAACAACTCGGGACCCTTGGGGAGCGGGCTCGCCCACGAACACTGACAAAGCTGGTTCCATCCACCGGTTTGTCTGCTTCGCGGGAGAGCCCGCTCCCACAGGGGGCCaa                                                                                             <  2:8322994/132‑1 (MQ=255)
                                         ggtgaaaacCCCAAAACACAATCTGAGGAACAACCCGGACCCCTGGGGGGACTGGCTTACCCCAAAACACATGCAAAACTCGCCCCACCCCCCGATTTGCCTTCTTCTCGGCCAACCCCGCTCCCACAGGGGGCCaa                                                                                             <  2:16075084/132‑1 (MQ=255)
                                         ggtgaaaaaCCCACAACCAAATCGGAGAAACAATGGGGAACCCTTGGGGACTGGGTTCGCCGGAGAACACCGGAAAAGCTGGTGCCCTCCACCGGGTTGTTTTCTTCGTGGGCGATCCCGCTCCCACAGGGGGCCaa                                                                                             <  2:6455222/132‑1 (MQ=255)
                                              aaaaCCCACACCGCAGTCTGAGAAACACTTTGGGACCTTGTGGGAGTGGGTTCCCCCGAAAACACCGAAAAACCCGGGGCCATCCACGGGGTTGCCTTCTTCGTGGGCGAGCCCGCTCCCACAGGGGGCCAAAgctt                                                                                        <  2:17667729/137‑1 (MQ=255)
                                              aaaaCCCACAACGCAGCCTGAGAACCAACTTGGGACCTTTTGGGAGCGGTCTCGCCCGAGAACACCGGCAAAGCCGTTGCCACCCACCGGGTTTTCGGCTTCGAGGGCGACCCCGCTCCCACAGGGGGCCAAAgctt                                                                                        <  2:19156732/137‑1 (MQ=255)
                                                    cacaACGCACTCTAAGAAACAATTGGTGCCCTTGGGGAAGTGGACTTTCCCACGAACACTGACAAACCCGGCGCCCTCAACGGATTGTTCTGCTTTGTGGGCGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTga                                                                                  <  2:5745267/137‑1 (MQ=255)
                                                      caACACACTTTGAGGAACACTTCGGGACCTTGGGGGGGCGGTCTTGCCCGCGAACACCGGAAAAGCCGGTGCCATCCACCGGGTTGTCTGCTTTGTGGGCGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTgaga                                                                                <  2:19387049/137‑1 (MQ=25)
                                                                ggagaAACAATTGGGGACCTTGTGGGAGGGGGCTCCCCCAAAAACACCGGAAAAGCCGGTGCCACCCACTGGGTTGTTTTTTTCGTGGGCGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAAtt                                                                      <  2:2425626/136‑1 (MQ=255)
                                                                  agaaacagcTCGGGACCTTTGGGGGGTGGGCCCTCCAAAAAACACCGAAAACGCCGGTCCAATCCCCTGATTTTTCTGCTGCGTGGCCGCTCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTAc                                                                    <  1:13255173/134‑1 (MQ=255)
                                                                  agaaacaatTTGGGACTTTGTGGGAGGGGTCTCGCCCGGAAACACTGGCAAAGTCGGGCCCGTCCACTGGGTTGTCTGCTTTGTGGGAGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTAc                                                                    <  1:18518320/134‑1 (MQ=255)
                                                                              gaaCCTTGTGGGAGCGGGCTCGCCCGAGAACACCGGAAAATCTGGTGCCATCCCCCGGTTTGTTTGTTTCGGGGGAAACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAAtt                                                        <  2:3292147/135‑1 (MQ=255)
                                                                               gACCTTGTGGGAGTGGGCTCTCCCGCGAACACCGGCAAATCTGGTGCCACCCCCTGGTTTGTTTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAAttt                                                       <  2:4626617/137‑1 (MQ=255)
                                                                               gACCCTGTGGGAGCTGTCTCACCTTCAAAAACCGGCAAAGTTGGTGCCACCCCCCGGGTTTTCTGTGTCGTGGGCGAACCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAAttt                                                       <  2:18359137/137‑1 (MQ=255)
                                                                               gACCCTGGGGGAGTGGCCTCCCCCGAAAACACCGAAAAATTCGGTACCACCCACCGGTTTTTCTGTTCGGCGGGCAACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAAttt                                                       <  2:3968826/137‑1 (MQ=255)
                                                                               aaCCCGGTGGGGGGGGGCTCGCCCGCGAACACAGGCAAAGCTGTTGCCACCCACCGGGTGGTTTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAAttt                                                       <  2:8587411/136‑1 (MQ=255)
                                                                                   ttgtgGGGGTGGTCTCTCCCGAAAACACTGGCAAATCTGGTCCCACCCCCCGGTTTTTTTTTTTCGCGGACAACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCa                                                   <  1:11932446/136‑1 (MQ=255)
                                                                                          ggtggtcTCCCCCGAGAACACTGGAAAATCTGTTGCCATCCACGGAGTTGTTTGTTCCGCGGGCAATCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTaa                                            <  1:17695548/131‑1 (MQ=255)
                                                                                                      tgtgAACACCGGAAAGGCTGTTGCCGTCCACCGGGTTGTTTGCTTCGCGGGCAAGCCCGCTCCCACAGTGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGTTTTCAt                                <  2:11964565/134‑1 (MQ=37)
                                                                                                      ggatggcacCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACCAAGAGGTAATTTCTCATGGCTAAAATAGGTTTCAt                                >  2:283972/7‑137 (MQ=255)
                                                                                                       caaaaaacCGACAAGACCGGTACCATCCCCTGGTTTTTTGCTTCCGCGGCCAACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATc                               <  1:15213081/131‑1 (MQ=255)
                                                                                                           acacCGGCAAATCGGTTGCTATCCCCCGAGTTTCCTGCTTCGCGGGCGAGCCCGCTCCCACAGTGGGCCAAGGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCTTTTTCATCGGCa                           <  2:7038591/137‑1 (MQ=255)
                                                                                                           acacCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCa                           >  1:19428625/1‑137 (MQ=255)
                                                                                                            cacTGGCAAATCGGGTGCCGTCCACCGGGTTGTTGGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCAc                          <  2:9168110/137‑1 (MQ=255)
                                                                                                             acCGGCAAAGCCGGTGCCATCCACAGAGGTGTCTGCTTCGCGAGCGAGCCCGCTCCCAAAGGGGGCCAAAGCTTGCTTGAGATATCACAAGTACAAAGAGGTAAGTTCTAATGGATAAGATCAGTTTCATCGGCAcc                         >  2:9658617/1‑137 (MQ=38)
                                                                                                             acCAGCCAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCACGGGCGAGCCCGATCCCCAAGGGGGCCAAAGCTTGCTTGAGATATCACAATCCCAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCAcc                         >  1:4097029/1‑137 (MQ=255)
                                                                                                                ggCAAATCCGGTGCCATCCACTGGGTGGTCTGTTTCGTGGGGGGTCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGc                      <  1:9670059/137‑1 (MQ=255)
                                                                                                                ggCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGc                      >  1:11732722/1‑137 (MQ=255)
                                                                                                                ggCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGc                      >  1:3936743/1‑137 (MQ=255)
                                                                                                                ggCAAACCTGGTGCCGTCCACTGGGTTTTTGGTTTTGGGGGCGATCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGc                      <  1:9670058/137‑1 (MQ=255)
                                                                                                                gcaaaaGCTGGTGCAATCCACAAAGTTGTCTGCTCCGCGGGAAACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGc                      <  2:7374837/134‑1 (MQ=255)
                                                                                                                 gaaaaGCTGTTGCCATCCCTCGGGTTGTTGGCTTCGGGGGCGCCCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGca                     <  2:4211983/135‑1 (MQ=255)
                                                                                                                  cAAAGCCGGTGCCATCCACCGAGTTGTCTTCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcat                    >  1:7099186/1‑137 (MQ=255)
                                                                                                                  cAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACACAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcat                    >  1:7162193/1‑137 (MQ=255)
                                                                                                                  cAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcat                    >  2:17969221/1‑137 (MQ=255)
                                                                                                                  cAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCAACGAcat                    >  2:2988538/1‑137 (MQ=255)
                                                                                                                  cAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCAAAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCCTCGGCACCGGcat                    >  2:4529554/1‑137 (MQ=255)
                                                                                                                   aaaGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcatc                   >  1:5719790/1‑137 (MQ=255)
                                                                                                                    aaGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcatca                  >  2:7129529/1‑137 (MQ=255)
                                                                                                                     atccGGTGCCCCCCCCCGGTTTTTCTGTTTCGCGGGCGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcatcat                 <  1:8859069/135‑1 (MQ=255)
                                                                                                                     aGCTGGTGCCATCCACCGCGTTGCCTGCTTCCCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAACTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcatcat                 <  2:9152790/137‑1 (MQ=255)
                                                                                                                     aGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcatcat                 >  2:4509044/1‑137 (MQ=255)
                                                                                                                      tctggTGCCGCCCACCGGGTTGTCTGCTCTGCGGCCAACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATg                <  2:3417265/134‑1 (MQ=255)
                                                                                                                      tccGGTGCCCTCCATGGTGTTGTCTGCTTCGGGGGAGATCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATg                <  2:12276579/136‑1 (MQ=255)
                                                                                                                      gctggggCCACCCACGGATTGGTCTGCTTCGGGGGCGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATg                <  2:14912650/131‑1 (MQ=255)
                                                                                                                      gctggTGCCCTCCACCGGTTTGTTGTTTTTGAGGACGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATg                <  1:15184493/134‑1 (MQ=255)
                                                                                                                      gctggTGCCATCCACCGGGTTGTCTGTTTGGGGGGCAACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATg                <  2:1204803/134‑1 (MQ=255)
                                                                                                                      gCCGGTGCCATCCACCGAGTTGTCGGTTTTGCGGAAAAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATg                <  2:12591770/137‑1 (MQ=255)
                                                                                                                      cccGTTGCCATCCACCGAGTTTTCTGCTTCGCGGGAGATCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATg                <  2:8932838/136‑1 (MQ=255)
                                                                                                                      cccGGTGCCACCCACCGAGTTGTCTGCTTCGCGGTCGGGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATg                <  2:18947128/136‑1 (MQ=255)
                                                                                                                        cGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATggg              >  2:17319378/1‑137 (MQ=255)
                                                                                                                         tgTGCCATCCACCGGTTGGTTTGTTTCTCGGGAGAACCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGc             <  2:15186245/136‑1 (MQ=255)
                                                                                                                         ggTGCCGTCCCCTGGGTTGTTTGTTTGGTGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGc             <  1:6819283/137‑1 (MQ=255)
                                                                                                                          gTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGTATCATGGGCa            >  2:2491801/1‑137 (MQ=255)
                                                                                                                            gCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAg          >  2:12296311/1‑137 (MQ=255)
                                                                                                                            gCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAg          >  1:14536715/1‑137 (MQ=255)
                                                                                                                            gCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAg          >  1:14710480/1‑137 (MQ=255)
                                                                                                                            gCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGACATCATGGGCAAg          >  1:247805/1‑137 (MQ=255)
                                                                                                                            accATCCACCGGGTTTTTTGTTTGCTGGGCGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAg          <  2:2698583/136‑1 (MQ=255)
                                                                                                                             ccccccACCGGTTTGTTTTTTTCGGGGGCGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc         <  2:5094080/133‑1 (MQ=255)
                                                                                                                             ccacccacTGGGTTGTCTGCTTCGAGGACGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc         <  1:12720341/137‑1 (MQ=255)
                                                                                                                             ccATCCCCCGAGTTGTCTGCTTCGTGGTCGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc         <  2:14810344/137‑1 (MQ=255)
                                                                                                                             ccATCCACTGATTTTTTTGTTTCGCGGGCGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc         <  2:10762908/137‑1 (MQ=255)
                                                                                                                              cgtcaaCCGATTTGTCGGTTTCGAGGGTGGCCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGcc        <  2:13970464/132‑1 (MQ=255)
                                                                                                                              cATCCCCGGGGTTGTCTGCTTCGCGGGCGAACCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGcc        <  2:13407310/137‑1 (MQ=255)
                                                                                                                              cATCCACTGAGTTTTTGGTTTCGCGGGCGATCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGcc        <  2:14789000/137‑1 (MQ=255)
                                                                                                                              cATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGcc        >  2:7036355/1‑137 (MQ=255)
                                                                                                                              cATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGcc        <  2:7691002/137‑1 (MQ=255)
                                                                                                                                   aCCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACCATTCCAAATAGGTAATTTCTCATGGCTACAATCGGTTTCATCGGCACCGGCCTCATGGGTAAGCCCATgg   >  1:63228/1‑137 (MQ=255)
                                                                                                                                   aCCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATgg   >  2:18023705/1‑137 (MQ=255)
                                                                                                                                   aCCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATgg   >  1:7837355/1‑137 (MQ=255)
                                                                                                                                   aCCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATgg   >  1:63267/1‑137 (MQ=255)
                                                                                                                                   aCCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATgg   >  1:5856484/1‑137 (MQ=255)
                                                                                                                                   aCCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATgg   >  1:14419868/1‑137 (MQ=255)
                                                                                                                                    ccGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGc  >  2:7340219/1‑137 (MQ=255)
                                                                                                                                    ccGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGc  >  2:3871175/1‑137 (MQ=255)
                                                                                                                                    ccGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGc  >  1:10519656/1‑137 (MQ=255)
                                                                                                                                       |                                                                                                                                      
GGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGC  >  AE015451_phaM1/4889488‑4889756

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: