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.

TTTTAAGATTTATATGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGC‑AAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACC  >  NC_000913/1121089‑1121365
                                                                                                                                       |                                                                                                                                              
ttttAAGATTTATATGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGC‑AAAAAAAATCACTg                                                                                                                                  >  3:70184/1‑149 (MQ=255)
  ttAAGATTTATATGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGC‑AAAAAAAATCa                                                                                                                                     <  6:38022/144‑1 (MQ=255)
  ttAAGATTTATATGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGC‑AAAAAAAATCa                                                                                                                                     >  5:38022/1‑144 (MQ=255)
       aTTTATATGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGC‑AAAAAAAATCACTGCAtatta                                                                                                                           <  4:70184/149‑1 (MQ=255)
            tatGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGC‑AAAAAAAATCACTGCATATTATTCGc                                                                                                                      >  2:34714/1‑149 (MQ=255)
                                    ttcatATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGC‑AAAAAAAATCACTGCATATTATTCGCTGACAAACCa                                                                                                            >  8:4224/1‑135 (MQ=255)
                                    ttcatATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGC‑AAAAAAAATCACTGCATATTATTCGCTGACAAACCa                                                                                                            <  7:4224/135‑1 (MQ=255)
                                          tttttttAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGC‑AAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCAc                                                                                         <  7:50227/148‑1 (MQ=255)
                                                   cATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGC‑AAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTg                                                                                 >  5:6995/1‑147 (MQ=255)
                                                             aaTTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGC‑AAAAAAAATCACTGCATATTATTCGCTGACAAAc                                                                                                              >  4:99141/1‑108 (MQ=255)
                                                                               tAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGC‑AAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTAtcttc                                                   <  6:6995/149‑1 (MQ=255)
                                                                                                  aTTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCt                                                                                  <  3:104220/100‑1 (MQ=255)
                                                                                                  aTTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCt                                                                                  >  4:104220/1‑100 (MQ=255)
                                                                                                       taagtaaAATAAAAAAGCCGGGGCGACCCGGC‑AAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATa                            >  1:66564/1‑148 (MQ=255)
                                                                                                              aataaaaAAGCCGGGGCGACCCGGC‑AAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTc                    <  1:34714/149‑1 (MQ=255)
                                                                                                                taaaaAAGCCGGGGCGACCCGGCAAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATcttctt                                                 <  3:98044/119‑1 (MQ=255)
                                                                                                                taaaaAAGCCGGGGCGACCCGGCAAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATcttctt                                                 >  4:98044/1‑119 (MQ=255)
                                                                                                                  aaaaaGCCGGGGCGACCCGGC‑AAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTg                                              >  2:61540/1‑119 (MQ=255)
                                                                                                                  aaaaaGCCGGGGCGACCCGGC‑AAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTg                                              <  1:61540/119‑1 (MQ=255)
                                                                                                                                 cccGGC‑AAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTAcc  >  2:104292/1‑148 (MQ=255)
                                                                                                                                       |                                                                                                                                              
TTTTAAGATTTATATGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGC‑AAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACC  >  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: ‑