Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,973,5350TC54.5% 4.7 / 12.7 11D597G (GAT→GGT) cheAfused chemotactic sensory histidine kinase in two‑component regulatory system with CheB and CheY: sensory histidine kinase/signal sensing protein
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base C (0/6);  total (2/9)
Fisher's exact test for biased strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-03

TTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCC  >  NC_000913/1973387‑1973669
                                                                                                                                                      |                                                                                                                                      
ttCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATc                                                                                                                                       <  2:315805/150‑1 (MQ=255)
                  cAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAgcgg                                                                                                                    <  2:245470/151‑1 (MQ=255)
                  cAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAgcg                                                                                                                     <  2:383884/150‑1 (MQ=255)
                                         gTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATca                                                                                              >  1:369860/1‑150 (MQ=255)
                                                        cgcgAAGATGGGTAGCAGCAAAAATGCGGGGTCTTTTGGGATATTTCATTCAAAGTTTTTAACCCCAACCCGGGTGTGAAAGAATAATTTGCCCCCCCAGAAGGGGGTGGGGGGGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGt                                                                              <  1:350488/151‑1 (MQ=255)
                                                         gcaaaGAATGGCGCAAGCAAAAACGCGGGGGAGTTGGCAATAGTTTCTTTCAAGGTTTTAACCCACACCCGGTTGTTGCCAAATTAATGGACCCCCCAGCAGGGGGTAGGGGGGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGt                                                                              <  1:218927/147‑1 (MQ=255)
                                                                                                 atatttttttaaaaGTTTTTAAACCACACCCGGTTGTTAACCAATTAATACACCCCCCAACAAGGCGTGGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCcac                                      <  2:215543/139‑1 (MQ=255)
                                                                                                 aTAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCca                                       >  2:131455/1‑150 (MQ=255)
                                                                                                                        ccACACCCGTGTGTTGACAAATTAATTAACCCCCCAGCAGGGGGTGGGGGCGCCCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTCGACGATGGGCAGATATTc               <  2:76164/151‑1 (MQ=255)
                                                                                                                            aaCTTGGTGTTAACAAATAAATTCACCCACCAGCAGGGCGGAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTCGACGATGGGCAGATATTCAccc           <  1:14534/151‑1 (MQ=255)
                                                                                                                                      ttaacaaataaatgGACCCCCCGGAAAGGGGGAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTcc  <  2:283816/145‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                      
TTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCC  >  NC_000913/1973387‑1973669

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