Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,584,9280T.60.8% 32.2 / 28.9 51coding (1559/2280 nt)ydePputative oxidoreductase YdeP
Reads supporting (aligned to +/- strand):  ref base T (10/10);  new base . (13/18);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 7.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.60e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AGCCTTACCTGACGGCGTCATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAACCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCA  >  NC_000913/1584621‑1585289
                                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                                                                                                    
aGCCTTACCTGACGGCGTCCATCCAGCGCCTTTTCAGCAGGCTGCATTTATCAGGTGAAAACCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGAt                                                                                                                                                                                                                                                                                                              <  1:320552/383‑1 (MQ=255)
                          cgcCTTTCAGCAGCTGCATTTATCAGGTGAAAACCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGc                                                                                                                                                                                                                                                                                                                               >  1:320559/1‑334 (MQ=255)
                                               ttATCAGGTGAAAACCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCt                                                                                                                                                                                                                                                                                               <  1:320558/348‑1 (MQ=255)
                                                    aGGTGAAAACCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGCCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTTGCATGTCAAtct                                                                                                                                                                                                                                                                                                                                                                   >  1:1137456/1‑275 (MQ=255)
                                                                      tGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGaa                                                                                                                                                                                                                                                                                                                                           <  1:1137454/281‑1 (MQ=255)
                                                                                gcTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATGAATCCTCAACGGTTACCGCCTGCGCACC‑‑Gttttt                                                                                                                                                                                                                                                                                                                                                                                  >  1:1137461/1‑234 (MQ=255)
                                                                                                 aaCTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTTg                                                                                                                                                                                                                                                                                                                                                                               <  1:1137460/220‑1 (MQ=255)
                                                                                                                   tCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTGCATGTc                                                                                                                                                                                                                                                                                                                                                                        <  1:1137464/207‑1 (MQ=255)
                                                                                                                           aTTGCGAAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCTGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTACAACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTGTCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGa                                                                                                                                                                                                                                >  1:320566/1‑332 (MQ=255)
                                                                                                                             tGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATGAATCCTCAACGGTTCCGCCTGCGCACC‑‑GTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGagag                                                                                                                                                                                                                                                                                                         >  1:1137467/1‑258 (MQ=255)
                                                                                                                                          cATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGTTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGcgctcgct                                                                                                                                                                                                                                       >  1:320570/1‑310 (MQ=255)
                                                                                                                                               tCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGCCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCa                                                                                                                                                                                                                                                                                                                    >  1:320573/1‑229 (MQ=255)
                                                                                                                                                   ccaccaGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑G‑‑TTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACa                                                                                                                                                                <  1:320565/366‑1 (MQ=255)
                                                                                                                                                     accaGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCACCGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑G‑‑‑TTTTGCATGTCAATCTCGCTACGTccc                                                                                                                                                                                                                                                                                                                                                        <  1:320567/189‑1 (MQ=255)
                                                                                                                                                        aGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACgg                                                                                                                                                                                                                                                      <  1:1137468/282‑1 (MQ=255)
                                                                                                                                                            aCTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACg                                                                                                                                                                                                                                                       >  1:320575/1‑279 (MQ=255)
                                                                                                                                                            aCTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCa                                                                                                                                                                                                                                                                          >  1:1137476/1‑259 (MQ=255)
                                                                                                                                                               cccAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGAt                                                                                                                                                                                                                                                                                                              <  1:1137469/221‑1 (MQ=255)
                                                                                                                                                                   ggCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAgcg                                                                                                                                                                                                                                             <  1:1137471/280‑1 (MQ=255)
                                                                                                                                                                   ggCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACgtgt                                                                                                                                                                                                                                                                                 <  1:1137470/245‑1 (MQ=255)
                                                                                                                                                                     cTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑G‑‑TTTTTGCATGTCAACTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCCGATCTGGCATTGCCCAGCGCAAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGa                                                                                                                    <  1:320571/394‑1 (MQ=255)
                                                                                                                                                                       aCCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGaa                                                                                                                                                                                                                                                                                                                                           <  1:1137473/185‑1 (MQ=255)
                                                                                                                                                                        ccACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑G‑‑TTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGtt                                                                                                                                                                                                                                                           <  1:1137474/260‑1 (MQ=255)
                                                                                                                                                                         cACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTGTTTTTTTGCATGTCAATCTCGCTACGTccc                                                                                                                                                                                                                                                                                                                                                        <  1:320574/174‑1 (MQ=255)
                                                                                                                                                                          aCGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACgt                                                                                                                                                                                                                                                                                   >  1:320576/1‑238 (MQ=255)
                                                                                                                                                                               ctGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTCCCGATCTGGCATTGCCAGCGCAAA‑‑TGTCCCCACTGCAGATCAATGCt                                                                                                                                                                                 >  1:320577/1‑321 (MQ=255)
                                                                                                                                                                               ctGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTAcc                                                                                                                                                                                                                                                                              >  1:1137479/1‑237 (MQ=255)
                                                                                                                                                                               ctGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAgcg                                                                                                                                                                                                                                             >  1:1137480/1‑266 (MQ=255)
                                                                                                                                                                                                    gATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTc                                                                                                                                                                                                                                    <  1:1137478/254‑1 (MQ=255)
                                                                                                                                                                                                                          ctctGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTGCATGTTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACgt                                                                                                                                                                                                                                                                                   >  1:1137483/1‑190 (MQ=255)
                                                                                                                                                                                                                          ctctGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACgt                                                                                                                                                                                                                                                                                   >  1:1137484/1‑189 (MQ=255)
                                                                                                                                                                                                                          ctctGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACgt                                                                                                                                                                                                                                                                                   >  1:1137482/1‑189 (MQ=255)
                                                                                                                                                                                                                           ctgattttCAGCATTACACCGGCGGGTTTTAACAACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAAATTGCGTTAACGGTAACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCACTGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATg                                                                                                                                 >  1:320582/6‑332 (MQ=255)
                                                                                                                                                                                                                                 ttCAGCATTACACCGGCGGGTTTTAACACGCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACGCCTGCGCACC‑‑‑‑TGTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTAc                                                                                                                                                                 <  1:320580/287‑1 (MQ=255)
                                                                                                                                                                                                                                 ttCAGCATTACACGGCGGGTTTTTAACACGCACGCGAGGCATGAATCATCGACATTGAATCCTCAAGGTTACCGCCTGCGCACC‑‑G‑‑TTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATAttgcttg                                                                                                                                                       <  1:320579/298‑1 (MQ=255)
                                                                                                                                                                                                                                  tCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACa                                                                                                                                                                                                                                                                                                                 >  1:1137487/1‑151 (MQ=255)
                                                                                                                                                                                                                                              cggcggGTTTTAACACGCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑TGTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACgg                                                                                                                                                                                                                                                      <  1:320581/197‑1 (MQ=255)
                                                                                                                                                                                                                                               gcggggTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAATCTAATTGCGTTAACGGTACAGCGCTCGCTCCCGATCTGGCATTGCCAGCGCAAATTGCCCCCACTGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGt                                         >  1:320586/4‑392 (MQ=255)
                                                                                                                                                                                                                                                      ttttAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTGCCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCaaa                                                                                                                                                                                                          <  1:1137485/231‑1 (MQ=255)
                                                                                                                                                                                                                                                      ttttAACACGCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTGCATGTCAATCTCGCTAGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATAttgcttg                                                                                                                                                       <  1:1137486/279‑1 (MQ=255)
                                                                                                                                                                                                                                                            cacgccacgcGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑G‑‑‑TTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCaaa                                                                                                                                                                                                          <  1:320583/222‑1 (MQ=255)
                                                                                                                                                                                                                                                              acgccacgcGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTGTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGt                                                                                                                          <  1:320584/302‑1 (MQ=255)
                                                                                                                                                                                                                                                                  cacgcGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGc                                                                                                                                             <  1:1137488/279‑1 (MQ=255)
                                                                                                                                                                                                                                                                   acgcGAGGCATGAATCATCGACATGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAAtt                                                                                                                                                                                                                                                                >  1:1137493/1‑168 (MQ=255)
                                                                                                                                                                                                                                                                     gcgAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCACGCAAAATTGCCCCCCATGCAGATc                                                                                                                                                                                       <  1:1137491/234‑1 (MQ=255)
                                                                                                                                                                                                                                                                     gcgAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑G‑‑TTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATg                                                                                                                                 <  1:1137490/287‑1 (MQ=255)
                                                                                                                                                                                                                                                                               aaTCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATGCCCCCACTGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGaa                                                                           >  1:320587/1‑331 (MQ=255)
                                                                                                                                                                                                                                                                                    tcGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑G‑‑TTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAAGATAGGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATc                                                                                                                                                                                       <  1:1137492/218‑1 (MQ=255)
                                                                                                                                                                                                                                                                                            gAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAATTGCCCCACTGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCGAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCTGTACATTAGAGTGTCCACGTAGTGGGCa  >  1:320589/1‑390 (MQ=255)
                                                                                                                                                                                                                                                                                              aTCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTATTGCGTTAACGGTACAGcgctcgc                                                                                                                                                                                                                                        >  1:1137494/1‑163 (MQ=255)
                                                                                                                                                                                                                                                                                                    aaCGGTTACCGCCTGCGCACC‑‑GTTTTTTTGCATGTTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAATCTAATTGCGTTAACGGTTACAgcgc                                                                                                                                                                                                                                            >  1:320590/1‑157 (MQ=255)
                                                                                                                                                                                                                                                                                                         ttACCGCCTGCGCACC‑‑TGTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCaa                                                                                                                                                                                                           <  1:320588/178‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                                            tttttttGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCt                                                                                              >  1:1137498/1‑271 (MQ=255)
                                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                                                                                                    
AGCCTTACCTGACGGCGTCATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAACCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACC‑‑GTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCA  >  NC_000913/1584621‑1585289

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: