Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,037,4960AT71.4% 24.3 / 14.1 29intergenic (+161/‑23)alaT/rrlAtRNA‑Ala/23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (3/5);  new base T (9/11);  total (13/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATT  >  NC_000913/4037248‑4037631
                                                                                                                                                                                                                                                        |                                                                                                                                       
aCTACATGTTATGGAGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAGTTTTCACAACACGATGATGAATCgaa                                                                                                                                        >  2:209851/1‑250 (MQ=255)
  tACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaa                                                                                                                                        >  1:73485/1‑248 (MQ=35)
         tATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCt                                                                                                                              >  1:284453/1‑251 (MQ=21)
          atggggcattcgcccAGCTGGGAGAACTCCTTCTTTTCACGCAGGAGGTCTGCGGTTTGATCCCGCATAGCTCCACCATCTCAGTATTTATTAAATAAAACAGACTTCAGAGCGTACCTGCAAATGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTTGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTTAGTCTCTCAAATTGTCGCAACACAATGATGAATCGTAAGATACATCtt                                                                                                                             <  2:106548/237‑1 (MQ=11)
           tGGGGCTTTAGCTCATCTGGGTGATCGCCTGCATTTCTAGCAGGATGTTTGAGGTTAGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAAAACTTCAGAGTGTATCTGCAAAGGTTCACTTCGAAGTTTTGCGCTTTAAAAATCTGTATCTAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGTGTTTCGTAAATTGTCGCAACACTATGATGAATCTTAAGAACCATCtt                                                                                                                             <  2:216154/250‑1 (MQ=2)
              ggCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGTTCCCGCATAGCTCCACCATCTCTGTAACGATTAAATAACAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAACAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCggg                                                                                                                         >  1:272851/1‑251 (MQ=255)
                         aGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAACTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAACGTTGTTCGTGAGTCACTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATATTCGGGTTGTGAGGtt                                                                                                               >  2:173981/1‑250 (MQ=255)
                                agagCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTACACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTAGCTGAAAAGGTTCACTGCGCAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAAAAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTAGCAACACAATGATGAATCTGAAGAAACATCTGCGGGTTGTTAGGTTAAGcacct                                                                                                       >  2:68423/1‑247 (MQ=11)
                                 gagCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTATCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACt                                                                                                       >  2:180401/1‑250 (MQ=37)
                                     gcCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg                                                                                                  <  2:259024/251‑1 (MQ=21)
                                         gCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTaca                                                                                              <  1:253106/251‑1 (MQ=37)
                                           tttGCACGCAGGAGGTCTGAGGTTCGATCCCGCATCGCTCCACCATATCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTCACTGCCAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGAGCAAGCTGAAATTTGAATCAATGAAAAACGAAAGTTGTTCGTGAGTCTCGCAAATTTACACAACACGAGGATGAATCGACAGATACATATTCGGGTTGTACGGTTAAGCGACTAAGCGTGCACg                                                                                            >  2:70031/1‑251 (MQ=11)
                                                 cGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTTATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGAt                                                                                      <  2:267257/251‑1 (MQ=21)
                                                  gCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAACGTACACGGTGGATg                                                                                     >  1:210477/1‑251 (MQ=37)
                                                             gctggtCGATACCTCATAGCTCCACCATCTCTGTAGTGAATAAATAAAAAATACTTCATAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTc                                                                          <  2:151898/246‑1 (MQ=255)
                                                             gCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGt                                                                           <  2:284453/250‑1 (MQ=25)
                                                             gCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTc                                                                          <  1:152568/251‑1 (MQ=21)
                                                                           cATAGCTACACCATCTCTGTAGTGATTAAATAAAAAATCCTTCAGAGTGTCCCTCCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAGAATCTGGATCAAGCTGAAAATTGAAACACTGACCAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTGCCATCACGATGATGCATCGTAAGAAACATCTTCGGGTTGAGAGGTTACGCGCATACGCGTACACGGTGGAAGCCCTGGCATTCAGGGGCGAAGAAgg                                                            >  2:267592/1‑251 (MQ=11)
                                                                                           ctGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCt                                                <  1:38736/247‑1 (MQ=25)
                                                                                               aGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATa                                          >  1:128099/1‑249 (MQ=25)
                                                                                                   aTTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTc                                    <  1:266561/251‑1 (MQ=35)
                                                                                                    ttAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCg                                   >  1:146951/1‑251 (MQ=18)
                                                                                                     taaataaaAAATACTGCAGCGTGTACCTGCAAAGGTTCATTGCGCAGATTTTCTCTTTAAAAATATGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTAAAATTTTCGAAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTACACGACTACGCGTACACGGTGGATGCCCTGGCAGTGAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGtagg                                  >  2:175817/1‑248 (MQ=11)
                                                                                                      aaataaaAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTTGGt                                 <  1:209008/251‑1 (MQ=21)
                                                                                                          agaaaaTACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTGCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAgg                             <  2:167918/249‑1 (MQ=18)
                                                                                                           aaaaaTTCTTCATAGTTTATCTGCAAACGTTCACTGCGAATTTTTGCACTTTAAAAATCTGGATCAAGCTGAAAATTTAAACACTGAACTACGAAAGTTGTTCGTGCGTCCCTCAAAGTTTCGCATCACGATGATGAATCGTAAGAAACATCTGCGGGTTGTGAGGGTAAGCGACTAAGCGTCCTCGGTGGAGGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTACTCTTCGATAAGCGTCGGTACGGt                            <  2:209455/251‑1 (MQ=11)
                                                                                                             aaaTACTTCAGAGTGTACCTGCCAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGAGATGATGGACGTGCATATCTGCGATAAGCGTCGGTAAGGTg                           >  2:267710/1‑250 (MQ=255)
                                                                                                                         gtttACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAATCATCTTAGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGTAGTCAGAGGCGATTAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTTATATGAACCGTta              <  2:172429/248‑1 (MQ=11)
                                                                                                                           gtACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTtata            <  1:193755/251‑1 (MQ=34)
                                                                                                                                   aaaGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCg     >  1:95714/1‑250 (MQ=18)
                                                                                                                                     aGGTTCATTGCGAAGTTGTGCGCTGTAAAAATCTGGATCAAGCTGAAAATAGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGATTCGAAAGAAACAGCTTCGGGTTGTGAGGTTAAGCGACTAAGCGCACACCGTGGATGCCCTGGCAGTCAGAGGCGATGTAGTACGTGCTAATCTGCGAGAAGCGACGGTAATGTTAGATGAACTGTTATAACCGGCGAtt  <  2:21573/251‑1 (MQ=11)
                                                                                                                                                                                                                                                        |                                                                                                                                       
ACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATT  >  NC_000913/4037248‑4037631

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑