Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,217,2430GA25.0% 37.2 / 10.8 24S55S (AGC→AGTyehXputative ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (7/11);  new base A (3/3);  total (10/14)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCAGCGAGCGAATTTCTTCTCCGGCAAAGCGGATCTCTCCGCTGTCATGCTCCACCAGGCGGTTAATCATTTTCAGGGTGGTGGATTTGCCGGAGCCAGATGTGCCAATCAGCACCGAAAAACTCCCTTCCTGAAAATTGAGATTGAGATCGTTAACGGCTTTTTGTGC  >  NC_000913/2217104‑2217371
                                                                                                                                           |                                                                                                                                
aTTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCAGCGAGCGAATTTCTTCTCCGGCAAAGCGGATCTCTCCGCTGTCATg                                                                                                                          <  5:3628/148‑1 (MQ=255)
              gCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCAGCGAGCGAATTTCTTCTCCGGCAAAGCGGATCTCTCCGCTGTCATGCTCCACCAGGCGGt                                                                                                            <  3:61254/148‑1 (MQ=255)
                    ttttGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCAGCGAGCGAATTTCTTCTCCGGCAAAGCGGATCTCTCCGCTGTCATGCTCCACCAGGCGGTTAATCa                                                                                                      <  5:138422/148‑1 (MQ=255)
                        gcgcCACGCTCCAGTGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCAGCGAGCGAATTTCTTCTCCGGCAAAGCGGATCTCTCCGCTGTCATGCTCCACCAGGCGGTTAATCAtttt                                                                                                  <  2:41708/148‑1 (MQ=255)
                                                   cAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCAGCGAGCGAATTTCTTCTCCGGCAAAGCGGATCTCTCCGCTGTCATGCTCCACCAGGCGGTTAATCATTTTCAGGGTGGTGGATTTGCCGGAGCCAGAt                                                                      >  5:14857/1‑149 (MQ=255)
                                                     aTAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCAGCGAGCGAATTTCTTCTCCGGCAAAGCGGATCTCTCCGCTGTCATGCTCCACCAGGCGGTTAATCATTTTCAGggt                                                                                            >  4:65842/1‑125 (MQ=255)
                                                      tAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCAGCGAGCGAATTTCTTCTCCGGCAAAGCGGATCTCTCCGCTGTCATGCTCCACCAGGCGGTTAATCATTTTCAGGGTGGTGGATTTGCCGGAGCCAGAtgtg                                                                   >  4:58377/1‑149 (MQ=255)
                                                        gATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCAGCGAGCGAATTTCTTCTCCGGCAAAGCGGATCTCTCCGCTGTCATGCTCCACCAGGCGGTTAATCATTTTCAGGGTGGTGGATTTGCCGGAGCCAGATGTGc                                                                  <  4:127051/148‑1 (MQ=255)
                                                                        ccATCCGGCGGCGCAACTCCAGTACTGGCAGCGAGCGAATTTCTTCTCCGGCAAAGCGGATCTCTCCGCTGTCATGCTCCACCAGGCGGTTAATCATTTTCAGGGTGGTGGATTTGCCGGAGCCAGATGTGCCAATCAGCACCGaaaaa                                                 >  5:5172/1‑149 (MQ=255)
                                                                        ccATCCGGCGGCGCAACTCCAGTACTGGCAGCGAGCGAATTTCTTCTCCGGCAAAGCGGATCTCTCCGCTGTCATGCTCCACCAGGCGGTTAATCATTTTCAGGGTGGTGGATTTGCCGGAGCCAGATGTGCCAATCAGCACCGaaaaa                                                 >  6:100740/1‑149 (MQ=255)
                                                                           tCCGGCGGCGCAACTCCAGTACTGGCAGCGAGCGAATTTCTTCTCCGGCAAAGCGGATCTCTCCGCTGTCATGCTCCACCAGGCGGTTAATCATTTTCAGGGTGGTGGATTTGCCGGAGCCAGATGTGCCAATCAGCACCGAAAAACTc                                              <  3:58377/149‑1 (MQ=255)
                                                                             cggcggCGCAACTCCAGTACTGGCAGCGAGCGAATTTCTTCTCCGGCAAAGCTGATCTCTCCGCTGTCATGCTCCACCAGGCGGTTAATCATTTTCAGGGTGGTGGATTTGCCGGAGCCAGATGTGCCAATCAGCACCGAAAAACTcc                                             <  6:5172/148‑1 (MQ=255)
                                                                                         tCCAGTACTGGCAGCGAGCGAATTTCTTCTCCGGCAAAGCGGATCTCTCAGCTGTCATGCTCCACCAGGCGGTTAATCATTTTCAGGGTGGTGGATTTGCCGGAGCCAGATGTGCCAATCAGCACCGAAAAACTCCCTTCCTGAAAatt                                >  8:26578/1‑149 (MQ=255)
                                                                                                   gCAGCGAGCGAATTTCTTCTCCGGCAAAGCGGATCTCTCCGCTGTCATGCTCCACCAGGCGGTTAATCATTTTCAGGGTGGTGGATTTGCCGGAGCCAGATGTGCCAATCAGCACCGAAAAACTCCCTTCCTGAAAattgagattga                        <  7:26578/147‑1 (MQ=255)
                                                                                                    cAGCGAGCGAATTTCTTCTCCGGCAAAGCGGATCTCTCCGCTGTCATGCTCCACCAGGCGGTTAATCATTTTCAGGGTGGTGGATTTGCCGGAGCCAGATGTGCCAATCAGCACCGAAAAACTCCCTTCCTGAAAATtgagattgagat                     <  5:100740/149‑1 (MQ=255)
                                                                                                                        cGGCAAAGCGGATCTCTCCGCTGTCATGCTCCACCAGGCGGTTAATCATTTTCAGGGTGGTGGATTTGCCGGAGCCAGATGTGCCAATCAGCACCGAAAAACTCCCTTCCTGAAAATTGAGATTGAGATCGTTAACGGCTTTTTGTgc  <  6:14857/148‑1 (MQ=255)
                                                                                                                          gCAAAGCGGATCTCTCCGCTGTCATGCTCCACCAGGCGGTTAATCATTTTCAGGGTGGTGGATTTGCCGGAGCCAGATGTGCCAATCAGCACCGAAAAACTCCCTTCCTGAAAATTGAGATTGAGATCGTTAACg             >  8:5166/1‑135 (MQ=255)
                                                                                                                          gCAAAGCGGATCTCTCCGCTGTCATGCTCCACCAGGCGGTTAATCATTTTCAGGGTGGTGGATTTGCCGGAGCCAGATGTGCCAATCAGCACCGAAAAACTCCCTTCCTGAAAATTGAGATTGAGATCGTTAACg             <  7:5166/135‑1 (MQ=255)
                                                                                                                                 ggATCTCTCCACTGTCATGCTCCACCAGGCGGTTAATCATTTTCAGGGTGGTGGATTTGCCGGAGCCAGATGTGCCAATCAGCACCGAAAAACTCCCTTCCTGAAAatt                                <  6:89080/109‑1 (MQ=255)
                                                                                                                                 ggATCTCTCCACTGTCATGCTCCACCAGGCGGTTAATCATTTTCAGGGTGGTGGATTTGCCGGAGCCAGATGTGCCAATCAGCACCGAAAAACTCCCTTCCTGAAAatt                                <  6:91358/109‑1 (MQ=255)
                                                                                                                                 ggATCTCTCCACTGTCATGCTCCACCAGGCGGTTAATCATTTTCAGGGTGGTGGATTTGCCGGAGCCAGATGTGCCAATCAGCACCGAAAAACTCCCTTCCTGAAAatt                                >  5:91358/1‑109 (MQ=255)
                                                                                                                                 ggATCTCTCCACTGTCATGCTCCACCAGGCGGTTAATCATTTTCAGGGTGGTGGATTTGCCGGAGCCAGATGTGCCAATCAGCACCGAAAAACTCCCTTCCTGAAAatt                                >  5:89080/1‑109 (MQ=255)
                                                                                                                                 ggATCTCTCCACTGTCATGCTCCACCAGGCGGTTAATCATTTTCAGGGTGGTGGATTTGCCGGAGCCAGATGTGCCAATCAGCACCGAAAAACTCCCTTCCTGAAAatt                                <  4:31445/109‑1 (MQ=255)
                                                                                                                                 ggATCTCTCCACTGTCATGCTCCACCAGGCGGTTAATCATTTTCAGGGTGGTGGATTTGCCGGAGCCAGATGTGCCAATCAGCACCGAAAAACTCCCTTCCTGAAAatt                                >  3:31445/1‑109 (MQ=255)
                                                                                                                                           |                                                                                                                                
ATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCAGCGAGCGAATTTCTTCTCCGGCAAAGCGGATCTCTCCGCTGTCATGCTCCACCAGGCGGTTAATCATTTTCAGGGTGGTGGATTTGCCGGAGCCAGATGTGCCAATCAGCACCGAAAAACTCCCTTCCTGAAAATTGAGATTGAGATCGTTAACGGCTTTTTGTGC  >  NC_000913/2217104‑2217371

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