Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,121,2231.A20.0% 45.7 / 12.6 20intergenic (‑268/+19)bssS/dinIbiofilm regulator/DNA damage‑inducible protein I
Reads supporting (aligned to +/- strand):  ref base . (9/7);  new base A (2/2);  total (11/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTTTAAGATTTATATGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACC  >  NC_000913/1121089‑1121365
                                                                                                                                       |                                                                                                                                              
ttttAAGATTTATATGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTg                                                                                                                                  >  3:70184/1‑149 (MQ=255)
  ttAAGATTTATATGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCa                                                                                                                                     <  6:38022/144‑1 (MQ=255)
  ttAAGATTTATATGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCa                                                                                                                                     >  5:38022/1‑144 (MQ=255)
       aTTTATATGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCAtatta                                                                                                                           <  4:70184/149‑1 (MQ=255)
            tatGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGc                                                                                                                      >  2:34714/1‑149 (MQ=255)
                                    ttcatATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCa                                                                                                            >  8:4224/1‑135 (MQ=255)
                                    ttcatATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCa                                                                                                            <  7:4224/135‑1 (MQ=255)
                                          tttttttAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCAc                                                                                         <  7:50227/148‑1 (MQ=255)
                                                   cATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTg                                                                                 >  5:6995/1‑147 (MQ=255)
                                                             aaTTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAAc                                                                                                              >  4:99141/1‑108 (MQ=255)
                                                                               tAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTAtcttc                                                   <  6:6995/149‑1 (MQ=255)
                                                                                                  aTTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCt                                                                                  <  3:104220/100‑1 (MQ=255)
                                                                                                  aTTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCt                                                                                  >  4:104220/1‑100 (MQ=255)
                                                                                                       taagtaaAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATa                            >  1:66564/1‑148 (MQ=255)
                                                                                                              aataaaaAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTc                    <  1:34714/149‑1 (MQ=255)
                                                                                                                taaaaAAGCCGGGGCGACCCGGCAAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATcttctt                                                 <  3:98044/119‑1 (MQ=255)
                                                                                                                taaaaAAGCCGGGGCGACCCGGCAAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATcttctt                                                 >  4:98044/1‑119 (MQ=255)
                                                                                                                  aaaaaGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTg                                              >  2:61540/1‑119 (MQ=255)
                                                                                                                  aaaaaGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTg                                              <  1:61540/119‑1 (MQ=255)
                                                                                                                                 cccGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTAcc  >  2:104292/1‑148 (MQ=255)
                                                                                                                                       |                                                                                                                                              
TTTTAAGATTTATATGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACC  >  NC_000913/1121089‑1121365

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: