Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,943,488	0	T	G	29.5%	66.9 / 13.4	44	noncoding (54/76 nt)	gltU	tRNA‑Glu
Reads supporting (aligned to +/- strand):  ref base T (13/18);  new base G (13/0);  total (26/18)
Fisher's exact test for biased strand distribution p-value = 4.09e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCA  >  NC_000913/3943349‑3943535                                                                                                                                            |                                                accttaAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGt                                                 >  1:1250828/1‑140 (MQ=32)   cttaAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTc                                               >  1:1149386/1‑140 (MQ=33)   cttaAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTc                                               >  1:328794/1‑140 (MQ=18)   cttaAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACACGGGTTc                                               >  2:1389198/1‑140 (MQ=17)     taAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGa                                             <  1:1330232/140‑1 (MQ=32)     taAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAAGGGTTCGa                                             >  2:856502/1‑140 (MQ=18)      aaaGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGaa                                            <  2:687689/140‑1 (MQ=32)      aaaGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGaa                                            <  1:796904/140‑1 (MQ=32)      aaaGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGaa                                            <  1:782149/140‑1 (MQ=32)       aagaagCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGaa                                            <  2:1543839/139‑1 (MQ=32)       aagaagCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGaa                                            >  1:1543839/1‑139 (MQ=32)        agaagCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATc                                          >  1:724940/1‑140 (MQ=18)         gaagCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATcc                                         <  1:50492/140‑1 (MQ=32)            gCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCt                                      >  2:193746/1‑140 (MQ=33)            gCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCt                                      >  2:920361/1‑140 (MQ=18)             cGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTa                                     >  1:1136896/1‑140 (MQ=33)               ttCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAgg                                   >  2:1119277/1‑140 (MQ=33)                 cTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAgggg                                 <  2:1501004/140‑1 (MQ=32)                  tttGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCTAggggg                                >  2:1257610/1‑139 (MQ=18)                    tGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACg                              >  2:1524120/1‑140 (MQ=33)                    tGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACg                              <  1:124206/140‑1 (MQ=32)                    tGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCTAGGGGACg                              >  1:191924/1‑140 (MQ=18)                      cAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGcc                            <  1:533897/140‑1 (MQ=32)                        gTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCAc                          <  1:862041/140‑1 (MQ=32)                         tGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACt                         >  1:1416506/1‑140 (MQ=33)                          gCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACtt                        >  1:927900/1‑140 (MQ=33)                          gCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACtt                        <  2:1224224/140‑1 (MQ=32)                           cTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTg                       >  1:1089557/1‑140 (MQ=33)                           cTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCTAGGGGACGCCACTTg                       >  1:608321/1‑140 (MQ=18)                            tCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGc                      >  1:1225337/1‑140 (MQ=33)                              acacaGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCTAGGGGACGCCACTTGCTg                    >  1:189645/1‑140 (MQ=18)                              acacaGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAAAACGGGGTCGAATCCCCTAGGGGGCGCCACTTGCTg                    >  2:643827/1‑140 (MQ=255)                                acaGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCTAGGGGACGCCACTTGCTGGt                  >  2:118384/1‑140 (MQ=18)                                   gATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCAc                          >  1:469714/1‑129 (MQ=33)                                   gATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCAc                          <  2:469714/129‑1 (MQ=32)                                   gATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCCTATCCCCTAGGGGACGCCTCCTGCTGGTTtg               >  2:514650/1‑140 (MQ=11)                                    aTTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTtgt              <  2:1511256/140‑1 (MQ=32)                                    aTTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAAGGGGTCGAATCCCCCAGGGGGCGCCCCCTGCTGGttttt              >  2:777388/1‑138 (MQ=11)                                       gTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTgtgag           <  2:927900/140‑1 (MQ=32)                                        tCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTgtgagg          >  2:1210855/1‑139 (MQ=18)                                          tGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTgtgagtga        <  2:328794/140‑1 (MQ=32)                                          tGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTgtgagtga        <  2:1149386/140‑1 (MQ=32)                                            aTAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGaaa      >  2:506222/1‑140 (MQ=18)                                            aTAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCCAGGGGGCGCCACTTGCTGGGTTGTGAGTGaaa      >  2:891224/1‑140 (MQ=255)                                              aGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGt    >  1:28117/1‑140 (MQ=32)                                              aGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCTAGGGGACGCCCCTTGCTGGTTTGTGAGTGAAAGt    >  1:610337/1‑140 (MQ=17)                                                gAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCa  <  2:28117/140‑1 (MQ=11)                                                gAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCa  <  1:1119277/140‑1 (MQ=11)                                                                                                                                            |                                                ACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCA  >  NC_000913/3943349‑3943535

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