Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,943,488	0	T	G	51.7%	12.5 / 26.4	29	noncoding (54/76 nt)	gltU	tRNA‑Glu
Reads supporting (aligned to +/- strand):  ref base T (8/6);  new base G (15/0);  total (23/6)
Fisher's exact test for biased strand distribution p-value = 6.32e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.63e-02

CCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTT  >  NC_000913/3943342‑3943543                                                                                                                                                   |                                                        cctcctTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCg                                                      >  2:912578/1‑150 (MQ=14)    ccttaccttaAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGACCTTTCACGGCGGTAAAAGGGGGTCGAAt                                                   >  2:662055/1‑150 (MQ=9)      ttaccttaAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGGAAAAGGGGGTCGAATccc                                                >  2:56081/1‑151 (MQ=9)       taccttaAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATcccc                                               >  2:955010/1‑151 (MQ=14)        accttaAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCt                                              <  2:738719/151‑1 (MQ=31)          cttaAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTTGAATCCCCt                                              >  2:257604/1‑149 (MQ=14)          cttaAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAAAAAGGGGTCGAAAcccc                                               >  1:895795/1‑148 (MQ=9)           ttaAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCCAgg                                           >  1:373320/1‑151 (MQ=14)            taAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTTGAATCCCCTTggg                                          >  1:1019357/1‑151 (MQ=14)               agaagCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCTAGGGGa                                        >  2:555184/1‑150 (MQ=14)                     gTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACt                                 <  1:860017/151‑1 (MQ=31)                     gTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCTAGGGGACGCCAc                                  >  2:529001/1‑150 (MQ=14)                     gTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGCAACAGGGGGTTGCATCCCCCAggggg                                        >  1:614866/1‑143 (MQ=11)                           tGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAAGGGGTCGAATCCCCCAGGGGGAGCCACTTGCTgg                           >  1:802967/1‑151 (MQ=255)                             cAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCTAGGGGACGCCACTTGCTGGt                          >  1:907114/1‑150 (MQ=14)                              aGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGttt                        <  1:582559/151‑1 (MQ=31)                              aGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTTGAATCCCCCAGGGGACGCCACTTGCTGGttt                        >  1:13950/1‑151 (MQ=9)                               gTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGttt                        >  1:33575/1‑150 (MQ=31)                                tGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAAGGGTTCGAATCCCCTAGGGGGCGCCACCTGCCGGtt                         >  2:278718/1‑148 (MQ=9)                                  cTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCCGGTTtgtg                     >  2:79872/1‑150 (MQ=14)                                   tCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCCAAGGGGCGCCACTTTCTGGCTtttgcg                   >  1:204696/1‑149 (MQ=11)                                      cacaGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTgtgagtga                <  2:448460/151‑1 (MQ=31)                                         aGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCTAGGGGGCGCCCCTTGCTGGTTTTTGAGGGaa               >  1:472232/1‑149 (MQ=255)                                         aGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGGAAACGGGGGTCGAAACCCCTAAGGGATGCCACCTGCTGGGTTgtgagggca               >  2:415352/1‑145 (MQ=11)                                                 tggtAGGAAGTGTAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCACTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGAAACAGGGGTGTGGAGCACGCAGGGGACGCCTCTTCACGGATTATGAGCGAAAgagagct       >  1:224747/4‑143 (MQ=2)                                                 tGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACCGGGGGTCGAATCCCCCAGGGGACGCCACTTTCTGGTTTGTGAGGGAAAGTCAcc        >  1:1058694/1‑148 (MQ=255)                                                  gATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTg      <  1:527474/149‑1 (MQ=31)                                                  gATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGcc    <  1:964600/151‑1 (MQ=31)                                                    tAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAAAAAGGGTTCGAGTCCCCTAGGGGGCGCCACCTTCCGGTTTTTTGATGAAAGACAACTGCCtt  >  1:925531/1‑151 (MQ=11)                                                                                                                                                   |                                                        CCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTT  >  NC_000913/3943342‑3943543

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: ‑