Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,594,1320C.54.5% 12.4 / 20.0 11intergenic (‑1690/+4485)hipB/lsrKantitoxin/DNA‑binding transcriptional repressor HipB/autoinducer‑2 kinase
Reads supporting (aligned to +/- strand):  ref base C (0/5);  new base . (4/2);  total (4/7)
Fisher's exact test for biased strand distribution p-value = 6.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGTTATCAACCCGAACGCGAGTGTTCCCTTGGGTATCGCCTTTCACGTTCATCTTATCGGTGGGCGAATTATCGCCGCCCAGCGTCGCGTTGAACGTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGCTACGCTCAGTGG  >  NC_000913/1593894‑1594380
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                        
tGTTATCAACCCGAACGCGAGTGTTCCCTTGGGTATCGCCTTTCACGTTCATCTTATCGGTGGGCGACTTATCGCCGCCCAGCGTCGCGTTGAACGTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGCTAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTc                                                                                                                                                                                                                                             >  1:177304/1‑251 (MQ=255)
 gTTATCAACCCGAACGCGAGTGTTCCCTTGGGTATCGCCTTTCACGTTCATCTTATCGGTGGGCGAATTATCGCCGCCCAGCGTCGCGTTGAACGTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCTGTTTCCCCCCTTCGCTATTCATTTGAACCGTTCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCTTCACGCAGTGCCCCCCCCAGGCATAACGTTc                                                                                                                                                                                                                                             <  1:131933/251‑1 (MQ=255)
    aTCAACCCGAACGCGAGTGTTCCCTTGGGTATCGCCTTTCACGTTCATCTTATCGGTGGGCGAATTATCGCCGCCCAGCGTCGCGTTGAACGTCATCAGGCCATGGTTTCCGGTAGATTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCata                                                                                                                                                                                                                                          <  1:200682/251‑1 (MQ=255)
            gAACGCGAGTGTTCCCTTGGGTATCGCCTTTCACGTTCATCTTATCGGTGGGCGAATTATCGCCGCCCAGCGTCGCGTTGACCGTCATCAGGCCATTGTTTTCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGc                                                                                                                                                                                                                                  <  1:71024/251‑1 (MQ=255)
                            ttGGGTATCGCCTTTCACGTTCATCTTATCGGTGGGCGAATTATCGCCGCCCAGCGTCGCGTTGAACGTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCg                                                                                                                                                                                                                 >  1:217536/1‑251 (MQ=255)
                              gggTATCGCCTGTCACTTTCATCTTAACTGTGGGCGAATTATCGCCGCCGAGCGTCGCGTTGAACGTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGTTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCCCCAGGTGGCCGCCGATAGTGGTTTTAGCGACGCGTAGTGGCCCCCCCAGGCATGACGTTCATATCTCCTGCAGTGGTGACATTTCCCGa                                                                                                                                                                                                                <  1:202609/251‑1 (MQ=255)
                               gttATCGCCTTTCACTTTCATCTTATCGGTGGTCGGATTTTCGCCGCCAATCGGCTCGTTGAACGTCTTCAGGCCATTGTTTCTGTTTTATTTGCTGTTACCGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCGCTTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTGTATATCACCTGCAGGTGTGACATTTCCCGACa                                                                                                                                                                                                              <  2:281623/249‑1 (MQ=255)
                                                                                                                                  aaCGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTAccc                                                                                                            <  1:40303/250‑1 (MQ=255)
                                                                                                                                          ggatctcccccGGTTTCCCCCCTTCGCTATTCATTTGACCCGGGCCGCCATCCTCCAGGTGGCCGCCGATAGTGGTTTTAGCGACTCGCAGGGCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTTACATTTCCCGACATCACACCGCTGGTATTGATATTAATGTGTTGACTTGCCATATTAATTGCGGAGGCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTTCCCTGATAcgcc                                                                                                   <  2:177304/244‑1 (MQ=255)
                                                                                                                                                                             tGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGCGTCAGGATCCCGGTAcc                                                                >  1:137357/1‑251 (MQ=255)
                                                                                                                                                                                  cGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGCTCCCGGTACCCAGtt                                                           >  2:36199/1‑251 (MQ=255)
                                                                                                                                                                                                                                          cccccAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGCTACGCTCAGt    <  1:36199/251‑1 (MQ=255)
                                                                                                                                                                                                                                            cccAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGCTACGCTCAGTgg  >  1:154835/1‑251 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                        
TGTTATCAACCCGAACGCGAGTGTTCCCTTGGGTATCGCCTTTCACGTTCATCTTATCGGTGGGCGAATTATCGCCGCCCAGCGTCGCGTTGAACGTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGCTACGCTCAGTGG  >  NC_000913/1593894‑1594380

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