Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,947,344	0	G	.	21.9%	109.0 / 29.3	41	intergenic (‑166/‑43)	mltA/metZ	membrane‑bound lytic murein transglycosylase A/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base G (14/18);  new base . (5/4);  total (19/22)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AACGAAGATAACAAACCGCCGCGTTTAATGAAATGAGCAGCCTCTCCCTGACGCAAATTTTGCACAAAAAATAGGCTTTAGTGATTTGTTTTTGTTCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGC  >  NC_000913/2947211‑2947486                                                                                                                                      |                                                                                                                                               aaCGAAGATAACAAACCGCCGCGTTTAATGAAATGAGCAGCCTCTCCCTGACGCAAATTTTGCACAAAAAATAGGCTTTAGTGATTTGTTTTTGTTCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAA‑GTTGCATGAAAAcgc                                                                                                                                 >  4:109725/1‑148 (MQ=255)    gAAGATAACAAACCGCCGCGTTTAATGAAATGAGCAGCCTCTCCCTGACGCAAATTTTGCACAAAAAATAGGCTTTAGTGATTTGTTTTTGTTCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAA‑GTTGCATGAAAACGCGAGc                                                                                                                             >  5:102061/1‑149 (MQ=255)            cAAACCGCCGCGTTTAATGAAATGAGCAGCCTCTCCCTGACGCAAATTTTGCACAAAAAATAGGCTTTAGTGATTTGTTTTTGTTCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAA‑GTTGCATGAAAACGCGAGCGGAGtata                                                                                                                     <  3:109725/149‑1 (MQ=255)                        tttAATGAAATGAGCAGCCTCTCCCTGACGCAAATTTTGCACAAAAAATAGGCTTTAGTGATTTGTTTTTGTTCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCa                                                                                                          >  1:174136/1‑149 (MQ=255)                          tAATGAAATGAGCAGCCTCTCCCTGACGCAAATTTTGCACAAAAAATAGGCTTTAGTGATTTGTTTTTGTTCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCAc                                                                                                         <  2:174136/148‑1 (MQ=255)                            atgaaatgaGCAGCCTCTCCCTGACGCAAATTTTGCACAAAAAATAGGCTTTAGTGATTTGTTTTTGTTCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGa                                                                                                      >  1:198376/1‑149 (MQ=255)                            atgaaatgaGCAGCCTCTCCCTGACGCAAATTTTGCACAAAAAATAGGCTTTAGTGATTTGTTTTTGTTCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGa                                                                                                      >  3:17089/1‑149 (MQ=255)                             tgaaatgaGCAGCCTCTCCCTGACGCAAATTTTGCACAAAAAATAGGCTTTAGTGATTTGTTTTTGTTCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGAc                                                                                                     >  4:131583/1‑149 (MQ=255)                                  tgaGCAGCATCTCCCTGACGCAAAATTTGCACAAAAAATAGGCTTTAGTGATTTGTTTTTGTTCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCggg                                                                                                <  7:36812/149‑1 (MQ=255)                                    agcagcCTCTCCCTGACGCAAATTTTGCACAAAAAATAGGCTTTAGTGATTTGTTTTTGTTCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAA‑GTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTg                                                                                             >  1:157161/1‑149 (MQ=255)                                     gcagcCTCTCCCTGACGCAAATTTTGCACAAAAAATAGGCTTTAGTGATTTGTTTTTGTTCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTg                                                                                             <  6:6920/149‑1 (MQ=255)                                      cagcCTCTCCCTGACGCAAATTTTGCACAAAAAATAGGCTTTAGTGATTTGTTTTTGTTCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTgg                                                                                            <  3:131583/149‑1 (MQ=255)                                               ccTGACGCAAATTTTGCACAAAAAATAGGCTTTAGTGATTTGTTTTTGTTCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAAGG‑TGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTg                                                                                   <  2:157161/148‑1 (MQ=255)                                               ccTGACGCAAATTTTGCACAAAAAATAGGCTTTAGTGATTTGTTTTTGTTCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAA‑GTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGcc                                                                                     <  6:102061/146‑1 (MQ=255)                                                   aCGCAAATTTTGCACAAAAAATAGGCTTTAGTGATTTGTTTTTGTTCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAA‑GTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGc                                                                              >  4:59519/1‑149 (MQ=255)                                                     gCAAATTTTGCACAAAAAATAGGCTTTAGTGATTTGTTTTTGTTCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCt                                                                             >  7:106182/1‑149 (MQ=255)                                                         aTTTTGCACAAAAAATAGGCTTTAGTGATTTGTTTTTGTTCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCtcgtc                                                                         <  8:106182/149‑1 (MQ=255)                                                             tGCACAAAAAATAGGCTTTAGTGATTTGTTTTTGTTCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCggg                                                                      <  2:198376/148‑1 (MQ=255)                                                                acaAAAAATAGGCTTTAGTGATTTGTTTTTGTTCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGagcagc                                                                                      <  6:118295/129‑1 (MQ=255)                                                                acaAAAAATAGGCTTTAGTGATTTGTTTTTGTTCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGagcagc                                                                                      >  5:118295/1‑129 (MQ=255)                                                                acaAAAAATAGGCTTTAGTGATTTGTTTTTGTTCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCa                                                                  >  4:51523/1‑149 (MQ=255)                                                                            cTTTAGTGATTTGTTTTTGTTCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGgt                                                      <  4:148439/149‑1 (MQ=255)                                                                                  tGATTTGTTTTTGTTCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAA‑GTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGt                                               >  7:57796/1‑149 (MQ=255)                                                                                    aTTTGTTTTTGTTCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGc                                                                     >  3:106632/1‑126 (MQ=255)                                                                                    aTTTGTTTTTGTTCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGc                                                                     <  4:106632/126‑1 (MQ=255)                                                                                     tttgtttttgttCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCtcgtcg                                                                        >  7:90204/1‑122 (MQ=255)                                                                                     tttgtttttgttCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCtcgtcg                                                                        <  8:90204/122‑1 (MQ=255)                                                                                              gttCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCtcgtc                                                                         <  2:108563/112‑1 (MQ=255)                                                                                              gttCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCtcgtc                                                                         >  1:108563/1‑112 (MQ=255)                                                                                                    aaTCATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCg                               >  2:105372/1‑148 (MQ=255)                                                                                                      tCATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAgg                                                       <  5:33832/122‑1 (MQ=255)                                                                                                      tCATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAgg                                                       >  6:33832/1‑122 (MQ=255)                                                                                                      tCATGCCAAATCCGTGATCGGGGTAAAAAAAA‑GTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAAc                           <  3:59519/149‑1 (MQ=255)                                                                                                        aTGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGa                                                                                           >  2:139765/1‑84 (MQ=255)                                                                                                        aTGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGa                                                                                           <  1:139765/84‑1 (MQ=255)                                                                                                         tGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAAcc                          >  4:39813/1‑148 (MQ=255)                                                                                                                tCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTaaaa                  <  3:51523/149‑1 (MQ=255)                                                                                                                  cGTGATCGGGGTAAAAAAAA‑GTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGTGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAtt                <  8:57796/148‑1 (MQ=255)                                                                                                                     gATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtga             <  4:17089/149‑1 (MQ=255)                                                                                                                               aaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTaaa    <  3:39813/148‑1 (MQ=255)                                                                                                                                aaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGc  <  1:105372/149‑1 (MQ=255)                                                                                                                                      |                                                                                                                                               AACGAAGATAACAAACCGCCGCGTTTAATGAAATGAGCAGCCTCTCCCTGACGCAAATTTTGCACAAAAAATAGGCTTTAGTGATTTGTTTTTGTTCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGC  >  NC_000913/2947211‑2947486

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