Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	213,341	0	A	G	47.3%	64.4 / 21.3	55	G337G (GGA→GGG)	tilS	tRNA(Ile)‑lysidine synthetase
Reads supporting (aligned to +/- strand):  ref base A (2/27);  new base G (26/0);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 1.14e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGA  >  NC_000913/213195‑213478                                                                                                                                                   |                                                                                                                                           gCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAgtg                                                                                                                                        <  1:269540/151‑1 (MQ=255)     cACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACa                                                                                                                                    <  2:445679/151‑1 (MQ=255)        ccTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCt                                                                                                                                 <  2:399032/151‑1 (MQ=255)         cTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGAGGACAGCtt                                                                                                                                >  1:525644/1‑151 (MQ=255)           gTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTa                                                                                                                               >  2:141289/1‑150 (MQ=255)               aCGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCCTAATGCg                                                                                                                          >  2:429196/1‑151 (MQ=255)               aCGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAGCTTaa                                                                                                                              >  1:591534/1‑147 (MQ=255)                   ttGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGCggg                                                                                                                        >  2:349209/1‑149 (MQ=255)                    tGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACAGCTTAATGCgggggg                                                                                                                     >  1:139140/1‑151 (MQ=255)                    tGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAg                                                                                                                      <  1:197432/150‑1 (MQ=255)                     gggCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAGCTTAATGCgggggg                                                                                                                     >  2:434059/1‑150 (MQ=255)                      ggCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCg                                                                                                                   <  2:483313/151‑1 (MQ=255)                            tttGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACAGGTTAAAGCGGGAGGGGAGAtt                                                                                                              >  2:456084/1‑150 (MQ=255)                             ttGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTc                                                                                                             <  1:230263/150‑1 (MQ=255)                             ttGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTc                                                                                                             <  2:201692/150‑1 (MQ=255)                               gAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAGCTTAATGGGGGGGGGGGAATTCGcc                                                                                                          >  1:636078/1‑151 (MQ=255)                                aaaTCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGTAAAGCTTTATGCGGGGGGCGGTATTTGccc                                                                                                         >  2:16707/1‑151 (MQ=255)                                    ccGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTcc                                                                                                      <  2:337337/150‑1 (MQ=255)                                    ccGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCcg                                                                                                     <  1:714758/151‑1 (MQ=255)                                            aTCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGGGGACCGCCTAATACGGGAGGCGATAttttcccccccc                                                                                                    >  2:703666/1‑134 (MQ=255)                                                gTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGaa                                                                                         <  1:154375/151‑1 (MQ=255)                                                 tCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGTGGAGGCTTATGCGGGGGGCGGTATTCGCCCTCCCCGTGCCGACg                                                                                           >  1:573508/1‑148 (MQ=255)                                                    cAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCTTAATACGGGGGGCGGTATTTGCCCCCCCCGTGCAGACGAAGCgg                                                                                     >  1:303548/1‑151 (MQ=255)                                                    cAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCgg                                                                                     <  1:643559/151‑1 (MQ=255)                                                         gtggtgGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAgc                                                                                <  1:719141/151‑1 (MQ=255)                                                         gtggtgGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACAGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCa                                                                                  <  2:247218/149‑1 (MQ=255)                                                          tggtggATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAGCTTAAATCGGGAGGGGATATTTGCCCCCCCCGTGCAGACGGAgggg                                                                                     >  1:431082/1‑145 (MQ=255)                                                                   aaaTCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACAACTTAATTCGGGGGGGGGTATTTTCCCCCCGCGTGGAGACGAAAAGGGCAGCGTGG‑Gttt                                                                       >  2:489615/1‑149 (MQ=255)                                                                   aaaTCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTc                                                                      <  1:421037/150‑1 (MQ=255)                                                                     aTCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCaaa                                                                   <  2:362113/151‑1 (MQ=255)                                                                       ccGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACAGCTTAATGCGGGGGGGGGGATTTGCCCTCCGCGTGCAGAAGAAGCGGGCAGCGGGC‑GTTTCaaa                                                                   >  1:695872/1‑149 (MQ=255)                                                                       ccGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCTTAATGCGGGGGGGGGTATTCGCCCTCCGCGGGCAGACGAAGCGGGCAGCGGGC‑GGTTCaaa                                                                   >  2:15957/1‑149 (MQ=255)                                                                       ccGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAg                                                                  <  2:458913/150‑1 (MQ=255)                                                                        cGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAgcg                                                                <  1:70493/151‑1 (MQ=255)                                                                              cGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGCGGGGGGGGGTATTTTCCCCCCGCGGGCAGAAGAAGCGGGCAGCGGGC‑GGTTCAAAACGCCAgg                                                           >  1:486507/1‑150 (MQ=255)                                                                                 gCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATtg                                                       <  2:339322/151‑1 (MQ=255)                                                                                      gCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGCGGGGGGGGGTATTTGCCCCCCGCGTGCAGACGAAGCGGGCAGCGGGC‑GTTTTAAAACGCCAGGATtg                                                       >  2:500717/1‑146 (MQ=255)                                                                                         aaaaCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCatat                                               <  1:187935/151‑1 (MQ=255)                                                                                          aaaCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACAACTTAATGCGGGGGGGGGTATTTGCCCCCCCCGGGCCGACGAAGCGGGCAGCGGGC‑GGTTTAAAACGCCCGGGTTGCTGCATAtt                                              >  2:65313/1‑151 (MQ=255)                                                                                          aaaCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATAtt                                              <  1:62608/151‑1 (MQ=255)                                                                                            aCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGt                                            <  2:367459/151‑1 (MQ=255)                                                                                                  tgCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAATGTACAACGTTATTCGGGGGGCGATATTTGCCCCCCCCCGGCAGAAGAAGCGGGCCGCGGGCGGTTTCAAA‑CGGCCGGATTGCTGCATAttttacg                                         >  2:547584/1‑143 (MQ=255)                                                                                                  tgCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGCGGGGGGGGGTATTCGCCCCCCGCGTGCAGACGAAGCGGTCCGCGTGC‑GTTTTAAAGCGCCAGGGTTGCTGCATATTTTCGGGCg                                      >  1:278273/1‑151 (MQ=255)                                                                                                   gCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACAGCTTTAAGCGGGGGGGGGTATTCGCCCCCCCCGTGGAGACGAAGCGGGCAGCGGGC‑GGTTTAAAACGCCCGGATTGCTGCCTATTTTCggggg                                      >  2:330148/1‑148 (MQ=255)                                                                                                     cGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGATAGCTTAATGCGGGGGGGGATATTCGCCCCCCGCGGGCAGACCCAGCGGGGCGGGGGC‑GTTTTCAAACGCCCGGGTTGGTGCATATTGTCcgg                                        >  1:166195/1‑146 (MQ=255)                                                                                                      gTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACCGCGTTATGCGGGGGGGGGTATTCGCCCCCCGCGTGGAGCCAGAGGGGGCCGCGGGG‑GGTTTTAAACGCCCGGGATGGTGCATATTGTCggggggca                                   >  1:159816/1‑145 (MQ=255)                                                                                                             aCGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGAc                           <  2:14217/151‑1 (MQ=255)                                                                                                                tGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAGCTTAATGCGGGGGGGGATATTCGCCCTCCGCGGGGAGACGAAGCGGGCAGCGGGC‑GGTTCAAAACGCCAGGACTGCTGCATATTGTCGGGCGGAAaggggg                             >  1:60171/1‑146 (MQ=255)                                                                                                                 ggCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTaa                       <  2:243166/151‑1 (MQ=255)                                                                                                                    ttCAACCGCTGGAATTACCGGCGGGGCTGGGGAGAGGACAGCTTAATGCGGGAGGGGATATTCGCCCTCCGCGTGCAGACGAAGCGGGCAGCGGGC‑GGTTCAAAACGCCAGGATTGCTGCATATTGTTGGGCGGAACGGCGGGCGGAAGCt                    >  2:673194/1‑151 (MQ=255)                                                                                                                     tCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGc                     <  1:727029/149‑1 (MQ=255)                                                                                                                      cAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTaa                  <  2:575412/151‑1 (MQ=255)                                                                                                                                 aTTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTgg       <  2:340451/151‑1 (MQ=255)                                                                                                                                     ccGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCa     >  2:346750/1‑149 (MQ=255)                                                                                                                                      cggcggGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAaga  <  2:83516/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                           GCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGA  >  NC_000913/213195‑213478

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