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

CTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCA  >  NC_000913/3943351‑3943535                                                                                                                                          |                                                cttaAAGGAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTc                                               <  2:1183311/140‑1 (MQ=17)   taAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGa                                             <  2:269350/140‑1 (MQ=32)   taAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGa                                             >  1:1208023/1‑140 (MQ=33)   taAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGa                                             <  2:1668059/140‑1 (MQ=32)   taAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGa                                             >  2:1530118/1‑140 (MQ=33)   taAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGa                                             <  2:553218/140‑1 (MQ=32)   taAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGa                                             <  1:1530118/140‑1 (MQ=32)    aaaGAAGCGTTCTTTGCAGTTCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGa                                             <  1:1180137/139‑1 (MQ=17)    aaaGAAGCGTTCTTTGCAGTTCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGa                                             >  2:1180137/1‑139 (MQ=18)    aaaGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGaa                                            <  2:1656773/140‑1 (MQ=32)     aagaagCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAAt                                           <  2:1659619/140‑1 (MQ=32)     aagaagCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGaaa                                           >  2:46625/1‑139 (MQ=18)      agaagCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATc                                          >  2:628205/1‑140 (MQ=33)      agaagCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATc                                          >  2:1007697/1‑140 (MQ=18)      agaagCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATc                                          >  1:307036/1‑140 (MQ=18)      agaagCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGGAACAGGGGATCGaaac                                          >  1:820765/1‑138 (MQ=12)      agaagCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGAAACAGGGGGTCggaac                                          >  2:386510/1‑136 (MQ=255)        aagCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATccc                                        <  1:585575/140‑1 (MQ=32)        aagCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATccc                                        <  2:202186/140‑1 (MQ=32)          gCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATccccc                                      >  2:923844/1‑139 (MQ=18)          gCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGGAACAGGGGGTCGGAGccccc                                      >  2:206301/1‑139 (MQ=255)           cGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCTa                                     >  2:1540903/1‑140 (MQ=18)           cGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCTa                                     >  2:206419/1‑140 (MQ=18)             ttCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCTAgg                                   >  2:1605601/1‑140 (MQ=18)             ttCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAAAAGGGGGTCGAATCCCCCAgg                                   >  1:1228099/1‑140 (MQ=16)                tttGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGa                                >  1:909912/1‑140 (MQ=33)                   gCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGc                             <  1:1540903/140‑1 (MQ=32)                      gTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCAc                          <  2:1604814/140‑1 (MQ=32)                      gTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCAc                          <  1:1605601/140‑1 (MQ=32)                        gCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACtt                        >  1:940621/1‑140 (MQ=33)                        gCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAAACCCCTAGGGGACGCCACtt                        >  1:871058/1‑140 (MQ=18)                         cTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTCACAGGGGGTCGAATCCCCCACGGGACGccccctg                       >  1:576103/1‑135 (MQ=11)                         cTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTg                       >  1:28616/1‑140 (MQ=32)                         cTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTg                       >  1:650641/1‑140 (MQ=33)                         cTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAAGGGGTCGAATCCCCCTGGGGACGCCACTTg                       >  2:1525580/1‑140 (MQ=255)                          tCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGc                      <  2:909912/140‑1 (MQ=32)                          tCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGc                      <  2:28616/140‑1 (MQ=32)                          tCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGc                      <  2:307036/140‑1 (MQ=32)                            acacaGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCTAGGGGACGCCACTTGCTg                    >  2:493243/1‑140 (MQ=18)                                aGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGttt                <  1:367000/140‑1 (MQ=32)                                aGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGttt                <  1:1095954/140‑1 (MQ=32)                                  aTTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTtgt              <  1:1120082/140‑1 (MQ=32)                                  aTTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTtgt              >  1:1435329/1‑140 (MQ=33)                                    tGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTgtga            <  1:628205/140‑1 (MQ=32)                                     gTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTgtgag           <  1:179768/140‑1 (MQ=32)                                      tCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGGCGCCACTTGCTGGGTTgtgagg          >  2:1091351/1‑139 (MQ=16)                                      tCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTgtgag           >  2:1319648/1‑139 (MQ=33)                                      tCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTgtgag           <  1:1319648/139‑1 (MQ=32)                                      tCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCTAAGGGGCGCCACCTGCTGGCTTttgagg          >  2:179768/1‑139 (MQ=11)                                          aTAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCTAGGGGGCGCCACTTGCTGGGTTGTGAGTGaaa      >  2:1095954/1‑140 (MQ=16)                                          aTAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCTAGGGGACGCCACTTGCTGGGTTGTGAGTGaaa      >  1:910712/1‑140 (MQ=18)                                          aTAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAAACCCCTAGGGGGCGCCACCTGCTGGTTTGTGAGTGaaa      >  2:1607832/1‑140 (MQ=255)                                            aGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGGCGCCACTTGCTGGTTTgtgag           >  1:1612517/1‑133 (MQ=17)                                            aGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTgtgag           <  2:1612517/133‑1 (MQ=32)                                            aGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTgtgag           <  2:1482621/133‑1 (MQ=32)                                            aGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTgtgag           >  1:1482621/1‑133 (MQ=32)                                              gAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCa  <  1:85475/140‑1 (MQ=11)                                                                                                                                          |                                                CTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCA  >  NC_000913/3943351‑3943535

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