Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,296,2680TC27.3% 26.0 / 11.2 25intergenic (+474/+168)gltP/yjcOglutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base T (11/5);  new base C (2/4);  total (16/9)
Fisher's exact test for biased strand distribution p-value = 1.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.93e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCA‑TCCGACATCAACG‑CCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCG  >  NC_000913/4296027‑4296515
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                           
cGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCA‑Tcc                                                                                                                                                                                                                                                   <  1:112291/251‑1 (MQ=37)
   aTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTAACCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAAGCCGCA‑TCCGAc                                                                                                                                                                                                                                                >  2:428883/1‑251 (MQ=255)
          gTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCAGACATCAACGCCTGATGCGACGCTTCACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATATCCTGATTTAGCGTGATTTTGTAGGTCGGATCAGGCGATTAAGCCCCA‑TCCGCCATCaatg                                                                                                                                                                                                                                         >  2:409922/1‑249 (MQ=11)
                gccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCA‑TCCGACATCaat                                                                                                                                                                                                                                          >  1:164239/1‑243 (MQ=25)
                gccgcAGCCGACATCAACGCCGGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCA‑TCCGACATCaat                                                                                                                                                                                                                                          <  2:164239/244‑2 (MQ=19)
                    cATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCA‑TCCGACATCAATG‑CCTGATGCGa                                                                                                                                                                                                                              <  1:363177/251‑1 (MQ=18)
                        cGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCA‑TCCGACATCAATG‑CCTGATGCGACGCt                                                                                                                                                                                                                          <  1:36510/251‑1 (MQ=17)
                         gACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTCTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGTGATTTTGTCGCTCGGATCAGGCGTTTACGC‑NCACTCCGACATCAATG‑CCCGATGCGACTCtt                                                                                                                                                                                                                         >  2:237567/1‑251 (MQ=11)
                         gACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTTTGTATATCAGATAATGCGTTTAAGCCGCA‑TCCGA‑ACCAAAGTCCTGATGCTACGCtt                                                                                                                                                                                                                         >  2:281431/1‑251 (MQ=11)
                                                                                                                                                                     cgcgTCTTATCAGCCCTACGCCAGACATCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCA‑TCCGACATCAACG‑CCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc                                                                             <  2:53136/251‑1 (MQ=2)
                                                                                                                                                                             atgaGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCA‑TCCGACATCAACG‑CCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                                                                      <  2:127923/247‑1 (MQ=16)
                                                                                                                                                                                                gcgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCA‑TCCGACATCAACG‑CCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCGTATCCGTTGGTTAATTTTTCAAACTCTTcac                                                  <  1:170313/251‑1 (MQ=255)
                                                                                                                                                                                                   cAATAGCCTGATTTAGCGTGCTTTTGTAGGTCGGATAAGGCGTTTATGCCGCA‑TCCGACATCAACG‑CCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                               <  2:370982/251‑1 (MQ=255)
                                                                                                                                                                                                    aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCA‑TCCGACATCAACG‑CCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTATCGTGATTTTGTAGATCGGATACGGCGATTAAGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGCCCTCTTTTAACCGTTGGTACATTTTTCAAGCTCTTCACAccc                                              >  2:296857/1‑251 (MQ=11)
                                                                                                                                                                                                                              aGGTCGGATAAGGCGTTTATGCCGCA‑TCCGACATCAACG‑CCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCACTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTCACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTc                     >  1:479099/1‑250 (MQ=255)
                                                                                                                                                                                                                              aGGTCGGATAAGGCGTTTATGCCGCA‑TCCGACATCAACG‑CCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTc                     >  1:241222/1‑250 (MQ=255)
                                                                                                                                                                                                                                 gcGGATAAGGCGTTTATGCCGCA‑TCCGACATCAACG‑CCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAACTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGt                 <  1:271615/250‑1 (MQ=255)
                                                                                                                                                                                                                                    gATAAGGCGTTTATGCCGCA‑TCCGCCATCAACG‑CCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCCGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCATATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCCGGtta               >  1:200931/1‑249 (MQ=17)
                                                                                                                                                                                                                                          gCGTTTATGCCGCA‑TCCGACATCAACG‑CCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCNAACCCTTCCATACAGCTCAGGTTCAACCAGt        >  2:246691/1‑251 (MQ=255)
                                                                                                                                                                                                                                            gTTTATGCCGCA‑TCCGACATCAACG‑CCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGCCATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTACCCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAACCCCTTCCATACAGCTCAGGTTCACCCAGTGc      >  2:190639/1‑251 (MQ=18)
                                                                                                                                                                                                                                            gTTTATGCCGCA‑TCCGACATCAACA‑CCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGc      >  1:140302/1‑251 (MQ=255)
                                                                                                                                                                                                                                               tATGCCGCA‑TCCGACATCAACG‑CCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgc   >  1:404245/1‑251 (MQ=255)
                                                                                                                                                                                                                                                aTGCCGCA‑TCCGACATCAACG‑CCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACCTCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCCGTgcggcg  >  2:364284/1‑251 (MQ=255)
                                                                                                                                                                                                                                                aTGCCGCA‑TCCGACATCAACG‑CCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgcg  >  1:393447/1‑251 (MQ=255)
                                                                                                                                                                                                                                                aTGCCGCA‑GCCGCCATCAAAG‑CATGATGAGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGAGTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGCCATCAATGCCTGATGCGACGCTTGCCGCGTCATATCAGGCCGATCTTAAACGTTGGTTAATTTTTCAAACTCCTCACAACCGGGATCACACCCTTCCCTACAGATCAgattaacccatctccgcg  >  2:337607/1‑234 (MQ=14)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                           
CGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCA‑TCCGACATCAACG‑CCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCG  >  NC_000913/4296027‑4296515

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: ‑