Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,848,3000CT79.7% 116.8 / 28.5 59G126E (GGG→GAG) ‡ydjAputative oxidoreductase YdjA
Reads supporting (aligned to +/- strand):  ref base C (4/8);  new base T (25/22);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.94e-02
Rejected as consensus: Frequency below/above cutoff threshold.

GATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAAAAATGC  >  NC_000913/1848061‑1848542
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   
gagaaGTGAGAGTGTCTGAATTCCTGCGCCTTTGCGCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAAct                                                                                                                                                                                                                                          <  2:171670/248‑1 (MQ=255)
           gtgtCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTg                                                                                                                                                                                                                              >  1:8707/1‑251 (MQ=255)
           gtgtCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTg                                                                                                                                                                                                                              <  1:296643/251‑1 (MQ=255)
                 gAATTCCTGCGCCTTTGCGCACAATCCATACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTTCTGCCATTTGCATCGc                                                                                                                                                                                                                        <  2:54401/251‑1 (MQ=255)
                  aaTTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGcc                                                                                                                                                                                                                       >  2:226743/1‑251 (MQ=255)
                   aTTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGcc                                                                                                                                                                                                                       >  1:59704/1‑250 (MQ=255)
                     tCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCAt                                                                                                                                                                                                                     <  2:8707/250‑1 (MQ=255)
                           gcCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGGGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACcgcg                                                                                                                                                                                                              >  1:102210/1‑251 (MQ=255)
                            ccTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCgcgc                                                                                                                                                                                                             >  1:59406/1‑251 (MQ=255)
                                    cacaATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGcc                                                                                                                                                                                                     >  2:152536/1‑250 (MQ=255)
                                    cacaATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCa                                                                                                                                                                                                     <  2:20388/251‑1 (MQ=255)
                                      caagccAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGa                                                                                                                                                                                                   <  2:102210/247‑1 (MQ=255)
                                          ccAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACAtt                                                                                                                                                                                               <  1:181530/251‑1 (MQ=255)
                                            aGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTc                                                                                                                                                                                             <  1:405356/251‑1 (MQ=255)
                                                  tttCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTc                                                                                                                                                                                       >  1:27625/1‑251 (MQ=255)
                                                     cgcgACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCa                                                                                                                                                                                    <  1:56113/251‑1 (MQ=255)
                                                            aTTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCggg                                                                                                                                                                             >  1:66178/1‑251 (MQ=255)
                                                               aTCAGAAATAAGTTACAAACGGCGTCGGGTCCTGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAAGTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTCCCGGACTTTCAGTTAATGCTCCACTGCACCCGATGCCGCCAAACTGCTGGGCATCTGCTGCCATTTGCATCGCCATGACCCCGCATCCGGCAGGCATTTCCTGTACCCCGCGCGGGAc                                                                                                                                                                           <  2:213698/250‑1 (MQ=255)
                                                               aTCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGAc                                                                                                                                                                           <  2:113828/250‑1 (MQ=255)
                                                                 cAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACttt                                                                                                                                                                        >  2:45429/1‑251 (MQ=255)
                                                                 cAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACttt                                                                                                                                                                        <  2:27625/251‑1 (MQ=255)
                                                                      aTAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGAt                                                                                                                                                                   >  1:339343/1‑251 (MQ=255)
                                                                        aaGTTACACACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACGGCGCATCCGGCAGACATTTACTGTTCCCAGCGCGGGACTTTATGATtc                                                                                                                                                                 >  2:333359/1‑251 (MQ=255)
                                                                             aCAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTc                                                                                                                                                            <  2:362989/251‑1 (MQ=255)
                                                                                aaCGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGc                                                                                                                                                          >  1:107014/1‑250 (MQ=255)
                                                                                aaCGGCGTCGGGTCCGGGACGGTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAAATTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGc                                                                                                                                                          >  2:407249/1‑250 (MQ=255)
                                                                                        cGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCGCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGcc                                                                                                                                                  <  2:343556/250‑1 (MQ=255)
                                                                                         gggTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGccac                                                                                                                                                >  1:303887/1‑251 (MQ=255)
                                                                                              cGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCg                                                                                                                                           >  1:42035/1‑251 (MQ=255)
                                                                                                 gACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCTCGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGtga                                                                                                                                        <  1:466927/251‑1 (MQ=255)
                                                                                                           gACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACTGCAACCGAATGCGTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACGCGTGCCATGTGCATCGCCATGACCGCTCATCCTGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCgg                                                                                                                              <  1:458826/251‑1 (MQ=255)
                                                                                                                       ttCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGACTGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCg                                                                                                                                           >  1:69894/1‑226 (MQ=255)
                                                                                                                       ttCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCg                                                                                                                                           <  2:69894/226‑1 (MQ=255)
                                                                                                                           gCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCAt                                                                                                             <  2:208707/251‑1 (MQ=255)
                                                                                                                                   gTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGTCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGc                                                                                                      >  2:81145/1‑251 (MQ=255)
                                                                                                                                         aGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTc                                                                                                <  1:383957/251‑1 (MQ=255)
                                                                                                                                          ggTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCa                                                                                                              <  2:19280/236‑1 (MQ=255)
                                                                                                                                          ggTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCa                                                                                                              >  1:19280/1‑236 (MQ=255)
                                                                                                                                              gagaAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAg                                                                                           <  1:265630/251‑1 (MQ=255)
                                                                                                                                                  aaaCCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCttt                                                                                       <  1:272232/251‑1 (MQ=255)
                                                                                                                                                     ccGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCTCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGtc                                                                                    <  1:97707/251‑1 (MQ=255)
                                                                                                                                                      cGACAATTGTATACTGATCAGGGCAACCCAATGCTTCACGCACTAGCGGAATTTCAGTTACTGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCAGTTGCATCGCCATGACCGCGCATCCGCCAGACATATCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTTATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCATTGTCGATAGCTTTGtc                                                                                    >  2:470726/1‑250 (MQ=255)
                                                                                                                                                       gaccaTTTTATCCTGCTCACGGCAACCGAATGCTTCCCGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGtcat                                                                                  <  2:330260/247‑1 (MQ=255)
                                                                                                                                                              ttATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTAcc                                                                           >  2:285620/1‑251 (MQ=255)
                                                                                                                                                                        cACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATc                                                                 >  1:414447/1‑251 (MQ=255)
                                                                                                                                                                        cACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATc                                                                 >  2:461061/1‑251 (MQ=255)
                                                                                                                                                                             cAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGcccc                                                            >  2:151045/1‑251 (MQ=255)
                                                                                                                                                                                  gAATGCTTCACGCAATACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTAGATCGCTTTGTCATCACTCCCGGCCGCAATCGCCCCCTGtt                                                       >  2:69937/1‑251 (MQ=255)
                                                                                                                                                                                                  aCCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAGCGGCGCATTACGGGCGTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCTCCCTGTTCCAGTACGTCGCTGaa                                       <  2:59406/251‑1 (MQ=255)
                                                                                                                                                                                                       aCTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCt                                  <  2:42035/251‑1 (MQ=255)
                                                                                                                                                                                                          ttCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTc                                 <  1:251304/249‑1 (MQ=255)
                                                                                                                                                                                                           tCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCAGCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTcgc                               <  2:59704/250‑1 (MQ=255)
                                                                                                                                                                                                              gTTAATACGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTCCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTACGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGTTGAAACGATCGCGccc                           >  2:308858/1‑251 (MQ=255)
                                                                                                                                                                                                                 aaTGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGcccttc                        <  2:66178/251‑1 (MQ=255)
                                                                                                                                                                                                                              gcCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTTTACCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCATCAATCGCGCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATcaca           <  1:313349/251‑1 (MQ=255)
                                                                                                                                                                                                                              gcCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATcac            >  1:274058/1‑250 (MQ=255)
                                                                                                                                                                                                                                       cgccAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAAAAATGc  >  2:214213/1‑251 (MQ=255)
                                                                                                                                                                                                                                       cgccAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAAAAATGc  >  1:285376/1‑251 (MQ=255)
                                                                                                                                                                                                                                       cgccAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAAAAATGc  >  1:378312/1‑251 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   
GATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAAAAATGC  >  NC_000913/1848061‑1848542

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