Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,442,632	0	T	G	63.6%	14.6 / 21.4	22	V290G (GTA→GGA)	mnmC	fused 5‑methylaminomethyl‑2‑thiouridine‑forming enzyme methyltransferase and FAD‑dependent demodification enzyme
Reads supporting (aligned to +/- strand):  ref base T (3/5);  new base G (14/0);  total (17/5)
Fisher's exact test for biased strand distribution p-value = 2.13e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.12e-04

GAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAGGCCCCCGCACTGGGTGCTTCCGGCAATCGCCAGGGGGCGCTGTATCCGTTATTAAGCAAACACGATGAGGCGCTAAACCGCTTTTTCTCTAAT  >  NC_000913/2442496‑2442747                                                                                                                                         |                                                                                                                    gAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGGAACGCTTTATTtcg                                                                                                       >  2:158347/1‑151 (MQ=255) gAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTGCCAGCGCGTTGTTGTCGCCGGCGCTATTTCGGCGCGGCTGGCAGGGAACGCCTttttttg                                                                                                       >  1:156359/1‑148 (MQ=255)  aaCAGACATTACCGCCCCCATGCTCCGCGCCGTGGTTTAACCGCACGGACAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCCGCGCGTTGTTGTCGCTGGGGCTATTACGGCGCTGCTGGCTGCTAACGCTTTATTgcag                                                                                                      <  1:601139/151‑3 (MQ=255)     aGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTGCCAGCGCGGTGTTGTCGCTGGCGCTATTAAGGCGCGGCTGGCAGGGAACGCTTTATTGCGCggg                                                                                                   >  2:216814/1‑150 (MQ=255)        cATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGGAACGCTTTATTGCGCGGGTGa                                                                                                >  2:503697/1‑151 (MQ=255)        cATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTGCCCGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGGTGGCAGGGAAAGCTTTATTGGGCGGATgg                                                                                                >  2:101241/1‑150 (MQ=255)        cATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTGCCAGCGCGTTGTTGTTGCTGGGGCCATTAAGGGGCGGGTGGGCGGGGACGCTTTTTTgttcgggggg                                                                                                >  1:443344/1‑142 (MQ=255)            aCCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAGGcc                                                                                            <  2:96398/151‑1 (MQ=255)                cTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGGAACGCTTTATTTCGCGGATGAGGccccc                                                                                         >  2:353092/1‑150 (MQ=255)                  ccccTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGTGGTGGGATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGGAACGCCTTATTTCGCGGATGAGGccccccc                                                                                       >  2:224148/1‑148 (MQ=255)                  ccccTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGGAACGCTTTTTTTCGCGGGTGAGGgcccccc                                                                                       >  2:558513/1‑148 (MQ=255)                  ccccTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGAGTTCCAGCGCGGTTGTGTCGCAGGCGGTGTTAAGACGCGGGCGGGGGGGAAAGGCTTAATGCGCGGGGGAGGccccccc                                                                                       >  1:652094/1‑148 (MQ=255)                   cccTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAGGCCCCCGCAc                                                                                     <  2:140657/151‑1 (MQ=255)                        cTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAGGCCCCCGCACTGGGt                                                                                <  2:463675/151‑1 (MQ=255)                          ccgcgccgcgGTAGAACCGCCCGGGCAGCAGCAAGCGGGAAGCGGCGATTATCGGCGCTGGTATTGCCAGCGCGTTGTTGTCGCCGGCGCTATTACGGCGCGGCTGGCAGGTACCGCTTTATTGCGCGGCTGAGGCCCCCGCACTGGGTg                                                                               <  1:334280/150‑1 (MQ=255)                                    gTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGGAAAGCTTTTTTTCGCGGGTGAGGCCCCCGCACTGGGGGCTTCCGGCaaa                                                                    >  2:677129/1‑150 (MQ=255)                                          ccGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTGCCAGCGCGGTGTTGTCGCTGGGGCTATTACGGCGCGGGTGGGAGGGAAAGCTTTTTTGCGCGGGTGAGGCCCCCGCCCTGGGGGGTTTCGGCAATCGcccg                                                              >  2:324862/1‑149 (MQ=255)                                             cACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGGAACGCTTTTTTTCGCGGGTGAGGCCCCCCCACTGGGGGCTTCCGGCAATCGCCCgggg                                                           >  2:491166/1‑151 (MQ=255)                                                                            aTCGGCGGTGGTATTGCCAGCGCGTTGTTTTCGCTGGCGCTATTACGGCGCGGCTGGCAGGGAAAGCTTTTTTTCGCGGGTGAGGGCCCCGCCCCGGGGGCTTTCCGGAATCGCCAGGGGGCGCCGTATCCCTTATTTAGCAAACACGATg                            >  1:591720/1‑151 (MQ=255)                                                                                   gtggtATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAGGCCCCCGCACTGGGTGCTTCCGGCAATCGCCAGGGGGCGCTGTATCCGTTATTAAGCAAACACGATGAGgcg                       >  1:230624/1‑149 (MQ=255)                                                                                               gcgcGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAGGCCCCCGCACTGGGTGCTTCCGGCAATCGCCAGGGGGCGCTGTATCCGTTATTAAGCAAACACGATGAGGCGCTAAACCGCttttt         >  2:372435/1‑151 (MQ=255)                                                                                                      gttgtCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAGGCCCCCGCACTGGGTGCTTCCGGCAATCGCCAGGGGGCGCTGTATCCGTTATTAAGCAAACACGATGAGGCGCTAAACCGCTTTTTCTCTAAt  >  2:593607/1‑151 (MQ=255)                                                                                                                                         |                                                                                                                    GAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAGGCCCCCGCACTGGGTGCTTCCGGCAATCGCCAGGGGGCGCTGTATCCGTTATTAAGCAAACACGATGAGGCGCTAAACCGCTTTTTCTCTAAT  >  NC_000913/2442496‑2442747

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