Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	213,341	0	A	G	50.0%	39.5 / 11.7	38	G337G (GGA→GGG)	tilS	tRNA(Ile)‑lysidine synthetase
Reads supporting (aligned to +/- strand):  ref base A (1/18);  new base G (19/0);  total (20/18)
Fisher's exact test for biased strand distribution p-value = 1.13e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCT  >  NC_000913/213192‑213481                                                                                                                                                      |                                                                                                                                              gATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGaa                                                                                                                                              <  1:34675/151‑1 (MQ=255)       tCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGTCGGGGGAg                                                                                                                                             >  1:779216/1‑146 (MQ=255)             tGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGCGGGGAGGGTACAGCtt                                                                                                                                   >  2:66712/1‑150 (MQ=255)                ttACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAA‑t                                                                                                                               <  2:259355/150‑1 (MQ=255)                    gTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACCGCTTAA‑TGCggg                                                                                                                          >  2:750307/1‑151 (MQ=255)                     tttGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAACTACCGGCGGGGGTGGGGAGGGTACAGCTTAA‑TGCgggg                                                                                                                         >  1:125018/1‑150 (MQ=255)                       tGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAA‑TGCGGGAg                                                                                                                        <  1:63030/150‑1 (MQ=255)                            gcgTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAA‑TGCGGGAGGCGata                                                                                                                  <  1:334994/151‑1 (MQ=255)                             cgTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGGGGACAGCTTAA‑AGCGGGGGGCg                                                                                                                     >  1:354551/1‑147 (MQ=255)                                    aaTCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGTAAAGCTTAA‑TACGGGGGGGGAGATTCGccc                                                                                                           >  2:497695/1‑150 (MQ=255)                                      tCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGCGGGGAGTGGTAAGCGTAATTCCGGGGGGGGAAATTCGccccc                                                                                                         >  1:659901/1‑149 (MQ=255)                                              tATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGGGGGGAGGGGAAAACTTTA‑TTCGGGGGGGGGTATTTGCCCCCCCCGTg                                                                                                   >  1:131713/1‑148 (MQ=255)                                              tATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGGAAAGCTTAA‑TGCGGGGGGCGAAATTTGCCCCCcgc                                                                                                      >  2:299970/1‑145 (MQ=255)                                                 cAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGAc                                                                                              <  1:169860/150‑1 (MQ=255)                                                  aGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGAAAGCTTAA‑TGCGGGGGGGGATATTTGCCCCCCGCGCGCAGa                                                                                               >  2:624176/1‑148 (MQ=255)                                                  aGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGa                                                                                            <  1:295261/151‑1 (MQ=255)                                                        aaCTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCg                                                                                        <  1:727343/149‑1 (MQ=255)                                                            gtggtgGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCa                                                                                    <  1:788695/149‑1 (MQ=255)                                                            gtggtgGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCa                                                                                    <  2:590868/149‑1 (MQ=255)                                                                 ggATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAgcgtgcg                                                                             <  1:146128/151‑1 (MQ=255)                                                                   aTTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGAGTACAGCTTAA‑TGCGGGGGGGGGTATTTTCCCCCCGCGTGCAGACGAAGCGGGCAgcgggcggt                                                                           >  2:62759/1‑149 (MQ=255)                                                                          ccGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACAGCTTAA‑TGCGGGAGGCGATATTTGCCCTCCCCGGGCAGACGAAGAGGGCAgcgtggggtttaa                                                                       >  2:205338/1‑140 (MQ=255)                                                                          ccGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAgc                                                                    <  1:341625/151‑1 (MQ=255)                                                                          ccGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAgc                                                                    <  1:248664/151‑1 (MQ=255)                                                                            gTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGTACAGCTTAA‑TGCGGGAGGCGATATTTGCCCTCCCCGCGCAGACGAAGCGGGCAGCGTGCGTTTCAAAgc                                                                    >  2:376502/1‑149 (MQ=255)                                                                            gTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGAAAGCTTTA‑TGCGGGGGGCGGTATTTCCCCTCCCCCCGCCGACGAAGCGGGCAGCGGGGGGTTCaaa                                                                      >  1:144736/1‑147 (MQ=255)                                                                              cACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACAGGCTAA‑TGCGGGAGGCGATATTCGCCCCCCGCGGGCGGACGAAGCGGGCAGCGGGCGGTTCAAAGCGcc                                                                 >  1:443163/1‑150 (MQ=255)                                                                                      aaaGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTAAAGCCTTA‑TGCGGGGGGCGGGATTTTCCCCCCCCGTGCAGACGAAAGGGGCAGCGGGCGGTTTAAAGCGCCCggggtggt                                                        >  1:511605/1‑145 (MQ=255)                                                                                                cATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTc                                              <  2:809563/151‑1 (MQ=255)                                                                                                    gtgCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAA‑TGCGGGGGGGGATATTCGCCCTCCGCGGGCAGACGAAGCGGGCAGCGGGCGGTTTAAAACGCCAGGATTGCTGCATATTTTTgggg                                          >  2:863830/1‑150 (MQ=255)                                                                                                       ccGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTACAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTa                                       <  2:171233/151‑1 (MQ=255)                                                                                                           ggCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGTTTGCTGCATCGTCGCGGGCGTAAcg                                    <  1:677010/150‑1 (MQ=255)                                                                                                               gACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTAAAGCTTAA‑TTCGGGGGGCGGTATTTGCCCCCCCCGTGCAGACGAAGCGGGCAGCGGGGGGTTCAAAGCGCCCGGAATGCTGCCTATTTTTGGGCgggaaggggg                                >  2:871856/1‑141 (MQ=255)                                                                                                                    ggCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTaa                          <  2:419849/151‑1 (MQ=255)                                                                                                                     gCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACAGCTTAA‑TGCGGGGGGGGGTATTTGCCCTCCGCGGGGCGACGAAGCGGTGCGCGGGCGGTTTTAAACGCCAGGATTGCTGGATATTTTTGGGGGGAAAGGCGGGCGGaa                          >  1:262619/1‑150 (MQ=255)                                                                                                                          aaCCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAA‑TGGGGGAGGCGGTATTCGCCCTCCCCCTGCCGGCGAAGAGGGCCGCGGGGGGTTCAAAGGGCCAAGAGTGCTGCATATTTTTGGGGGTAAACGGGGACGGaaaccca                     >  2:83680/1‑145 (MQ=255)                                                                                                                                       aCCGGCGGGGCTGGGAAGTGTACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCaa       <  1:708686/151‑1 (MQ=255)                                                                                                                                            cggGGCTGGGAAGTGTACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCt  <  2:796046/151‑1 (MQ=255)                                                                                                                                             ggggCTGGGAAGTGTACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCt  <  1:66702/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                              GATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCT  >  NC_000913/213192‑213481

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