Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,296,268	0	T	C	28.1%	33.0 / 20.2	34	intergenic (+474/+168)	gltP/yjcO	glutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base T (6/17);  new base C (3/6);  total (9/25)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.27e-02
Rejected as consensus: Frequency below/above cutoff threshold.

GTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGA‑CATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCA  >  NC_000913/4296021‑4296507                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   gTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGcc                                                                                                                                                                                                                                                 <  1:363122/251‑1 (MQ=21) gTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGcc                                                                                                                                                                                                                                                 >  1:485791/1‑251 (MQ=21)       cggcaacggCGTTGAAGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTAGCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATc                                                                                                                                                                                                                                            <  2:447613/243‑1 (MQ=255)            aaGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGA‑CAt                                                                                                                                                                                                                                     <  1:264151/251‑1 (MQ=38)              ggCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGA‑CATCa                                                                                                                                                                                                                                   <  1:525054/251‑1 (MQ=38)                  tttACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAACGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGA‑CATCaatgc                                                                                                                                                                                                                               >  1:224820/1‑248 (MQ=12)                     aCGCCGCAGCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGA‑CATCAATGCCTg                                                                                                                                                                                                                            <  1:377702/251‑1 (MQ=21)                            aTCCGACATCAACGCCTGATGAGACGCTTGCCGCATCTTATCAGGCCTACACCAGACAGCCCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGACCTACGCCAGACAGCGCCATAGCCTGATTTAGCGTGATTTTGTAGGTCGTATAAGGCGTTTACGCCGCATCCGA‑CATCAATGCCTGATGCGAc                                                                                                                                                                                                                     >  1:542690/1‑251 (MQ=11)                                gACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGA‑CATCAATGCCTGATGCGACGCtt                                                                                                                                                                                                                 <  1:511161/251‑1 (MQ=14)                                                                           tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCAGA‑CATCAACGCCTGGTGCGACGCTTGACGCGTCTTATCAGGCCTACGCCAGACGGGGCAATAGACTg                                                                                                                                                                       >  2:331919/1‑250 (MQ=255)                                                                                                                                                                             gcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGA‑CATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                                                                    <  1:556993/251‑1 (MQ=11)                                                                                                                                                                                   tgccAGGCCTACGCCAGACAGCGCAAGAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGAATCCGA‑CATCAACGCCTGATGCGACGCTTAACGCGTCGTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGtt                                                                 <  2:311789/245‑1 (MQ=11)                                                                                                                                                                                    atcaggtatcccgcAGACAGCGCAATAGACTGATTTAGCGTGATTGTGTAGGTCGGATAAGGCGTTTATGGCGCATCGGG‑CATCAACGCCTGATGCGACGCTTGTCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCAGGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                                                              <  2:286846/237‑1 (MQ=11)                                                                                                                                                                                    aTCAGTCCTATGCCATACATGGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTTGGCCTCACCCGC‑CATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTGTGTAGGTCGGATAAGGCGTTTACGCCGCACCCGACATCAATGCCTGATGCGGCGCTTGCCGCGTCGTATCAGGCCTATCTTAACCGTTGGt                                                              <  2:177942/250‑1 (MQ=11)                                                                                                                                                                                    aTCAGGCAGACGCCAGACGTCTCAATAGCCTGATTTAGCGTGATTTTGTATGTCGGATAAGGCATTTATGCCGCATCCGA‑CATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGtt                                                             <  2:312156/251‑1 (MQ=11)                                                                                                                                                                                        ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGA‑CATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtt                                                         <  1:366877/251‑1 (MQ=25)                                                                                                                                                                                           ggatgtcaggcggcgcaattgtaggATTGAGCGTGATTTTGTAGGTAGGAGAAGGCGTTTAGGCCGCATC‑GAGCATCAACGCCGTATGCGACGCTAACCTCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCTATATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                                                      <  2:620404/227‑1 (MQ=11)                                                                                                                                                                                             aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGA‑CATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTATCCGTTGGTTAATTTTTCa                                                    <  1:108959/251‑1 (MQ=255)                                                                                                                                                                                                 ccgccAGCGCAATAACCTGATTTAGCGTGATTGTGTAGGTCGGATAAGGAGTTTATGCCGCATCCGA‑CATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGGGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTTCCGCGTCTTATCAGGCATATTTTAACCGTTGGTTAATTTTTCAAAct                                                <  1:371137/247‑1 (MQ=11)                                                                                                                                                                                                     cagtgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAGGGCGTTTATGCCGGATCCGA‑CATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGGCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc                                            <  1:273886/247‑1 (MQ=14)                                                                                                                                                                                                           aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGA‑CATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAccc                                      >  1:473429/1‑251 (MQ=255)                                                                                                                                                                                                                 catctgtAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGA‑CATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATc                                <  2:94256/245‑1 (MQ=255)                                                                                                                                                                                                                   gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGA‑CATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCa                               >  2:605762/1‑250 (MQ=255)                                                                                                                                                                                                                      ttAGCGTGATTTTGTAGTGCGGATAAGGGGTTTATGCCGCATCCGC‑CATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCTACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACTCCCGGTATCAAAc                            <  2:522753/250‑1 (MQ=14)                                                                                                                                                                                                                         gCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCAACCGA‑CATCATCGCCTGATGCGACGCTTAACGCGTCTTAACAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCtt                        <  2:351868/251‑1 (MQ=255)                                                                                                                                                                                                                          cGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGAATCCGA‑CATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                       >  1:469252/1‑251 (MQ=255)                                                                                                                                                                                                                                 ttGTAGGTCGGATAAGGCGTTTATGCCGCATCCGA‑CATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACa                 <  1:163184/250‑1 (MQ=255)                                                                                                                                                                                                                                  tggtcgTCGGATAAGGCGTTCGTGCCGCTCCCGA‑CATCAACGCCTGATGCGACGATTAACGCGTCTTATCAGGCCTACGCCAGAAAGCGCAATAGCCTGATGTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACTCCCGGTATCAATCCCTTCCAGACAg                <  2:117302/245‑1 (MQ=17)                                                                                                                                                                                                                                     aGGTCGGATAAGGCGTTTATGCCGCATCCGA‑CATCAACGCCTGATGCGTCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTATGTCGGATAAGGCGTTTACGCCGCATCCGATATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAGCTTATCCGGTTGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTc             <  1:169468/250‑1 (MQ=17)                                                                                                                                                                                                                                         cGGATAAGGCGTTTATGCCGCATCCGA‑CATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGt         <  1:318305/250‑1 (MQ=255)                                                                                                                                                                                                                                          ggATAAGCCGTTTATACCGCATCCGC‑CATCAACGCCTGATGCGACGCTTAACGCGTCTTATCCGGCCTACGCCAGCCCGCGCAATAGCCTGATTTAGAGTGATTTTGTAGATCGGATACGTCGTTTAAGCCGCATCGGACATCAATTACTGATGCGACTCTTGCCGCGTCTTCTCGGGCCTATCTAAACCGTCAGTTATTTTTTAAAATTCTTCACACCCGTTCTCAAACCCTTGCATAccattccggttc       >  2:270514/1‑240 (MQ=9)                                                                                                                                                                                                                                           gATAAGGCGTTTATGCCGCAGCCGA‑CATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCa      <  1:357866/251‑1 (MQ=255)                                                                                                                                                                                                                                            aTAAGGCGTTTATGCCGCATCCGA‑CATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCa      >  1:481834/1‑250 (MQ=255)                                                                                                                                                                                                                                               aGGCGTTTATGCCGCATCCGA‑CATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCa  <  1:605762/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   GTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGA‑CATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCA  >  NC_000913/4296021‑4296507

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