Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,209,827	0	T	G	56.7%	15.0 / 36.5	31	noncoding (54/76 nt)	gltV	tRNA‑Glu
Reads supporting (aligned to +/- strand):  ref base T (2/11);  new base G (17/0);  total (20/11)
Fisher's exact test for biased strand distribution p-value = 1.43e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.15e-06

ATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTT  >  NC_000913/4209677‑4209882                                                                                                                                                       |                                                        aTCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGt                                                         <  2:214060/151‑1 (MQ=31)          ttaccttaAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACCGGGGGTCGAATccc                                                >  2:683936/1‑151 (MQ=14)           taccttaAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATcc                                                 <  2:1062148/149‑1 (MQ=31)            accttaAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGACTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGCAACAGGGGTTCCAATccc                                                >  2:412917/1‑149 (MQ=9)               ttaAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCTacg                                           >  2:796905/1‑149 (MQ=14)                taAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAg                                            <  2:926067/149‑1 (MQ=31)                taAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCTAggg                                          >  2:1118114/1‑151 (MQ=14)                  aagaagCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCTAgggg                                         >  2:788924/1‑150 (MQ=14)                     aagCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGc                                     <  2:1170797/151‑1 (MQ=31)                     aagCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCTAGGGGa                                        >  2:114965/1‑148 (MQ=14)                           tCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACtt                                <  2:1126813/150‑1 (MQ=31)                           tCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCTAGGGGACGCCACtt                                >  2:664279/1‑150 (MQ=14)                            cTTTGCAGTGCTCACACAGATTGTCTGATCGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTACCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGAAAAAGGGGTTCGAATCCCCCAGGGGAccgcacctg                               >  2:452576/1‑142 (MQ=11)                             tttGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAAAAGGGGGTCGGATCCCCTAGGGGACGCCACtt                                >  1:627459/1‑148 (MQ=9)                               tGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTg                            <  1:1022407/150‑1 (MQ=31)                                    tGCTCACGCAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTtg                       <  1:686948/150‑1 (MQ=14)                                       tCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCCAGGGGACGCCACTTGCTGGTTtgt                      >  1:1109825/1‑148 (MQ=14)                                        cacacaGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTgtgag                   <  1:453483/150‑1 (MQ=31)                                         acacaGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTCACAGGGGCTCGAAACCCCTAGGGGACCCCACTTGCCGGTTTgtgggt                  >  2:830852/1‑150 (MQ=11)                                         acacaGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTgtga                    <  2:269962/148‑1 (MQ=31)                                             aGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCTAGGGGACGCCACTTGCTGGGTTGTGAGTGaa               >  2:627645/1‑149 (MQ=14)                                             aGATTGTCTGATAGAAAGCGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAAAAGGGGGTCGAATCCCCCAGGGGACGCCACTTGCTGGTTTGTGAGTGaa               >  2:1001874/1‑149 (MQ=9)                                              gATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCTAGGGGGCGCCACTTGCTGGTTTGTGAGTGaaa              >  1:604729/1‑149 (MQ=14)                                                ttGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGaaa              >  2:672123/1‑147 (MQ=14)                                                 tGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCTTGGGGGCGCCACTTGCTGGTTTGTGAGTGAAAGTCa          >  2:470269/1‑150 (MQ=9)                                                   tCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAAAAGGGGGTCGGATCCCCCTGGGGACGCCACCTGCTGGGTTTTGGGGGGAAGTCAcc        >  1:776545/1‑150 (MQ=11)                                                      gATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGc     <  1:335273/150‑1 (MQ=31)                                                       aTAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCt   <  1:95811/151‑1 (MQ=31)                                                       aTAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCCAGGGGACGACACTTGCTGGGTTTTGGGGGAAAGTCACCTGcc    >  1:723192/1‑150 (MQ=11)                                                         aGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCtt  <  1:440713/150‑1 (MQ=31)                                                          gAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCTAGGGGAAGCCACTTGCTGGTTTGTGGGGGAAAATCACCACCCtt  >  1:177660/1‑149 (MQ=11)                                                          gAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGGAACAGGGGGTCGAATCCCCCAGGGGGCGCCACTTGCTGGTTTGTGAGGGAAAATCACCTGCCtt  >  1:358920/1‑149 (MQ=11)                                                                                                                                                       |                                                        ATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTT  >  NC_000913/4209677‑4209882

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