Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	541,275	0	G	T	46.4%	33.0 / 10.7	28	L237F (TTG→TTT)	ybbY	putative uracil/xanthine transporter
Reads supporting (aligned to +/- strand):  ref base G (4/11);  new base T (13/0);  total (17/11)
Fisher's exact test for biased strand distribution p-value = 6.72e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

CCGCAACGTTTTGCCCGTTATGGCCTGCTGGTCGGCACCATAACCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGTGCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGC  >  NC_000913/541135‑541412                                                                                                                                             |                                                                                                                                          ccGCAACGTTTTGCCCGTTATGGCCTGCTGGTCGGCACCATAACCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGTGCTTTGTCGCCGGGaa                                                                                                                                 <  1:533416/151‑1 (MQ=255)     aaCGTTTTGCCCGTTATGGCCTGCTGGTCGGCACCATAACCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGTGCTTTTTCGCCggggagt                                                                                                                               >  2:37658/1‑145 (MQ=255)       cGTTTTGCCCGTTATGGCCTGCTGGTCGGCACCATAACCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGTGCTTTGTTGCCGGGGatttttt                                                                                                                           >  1:615850/1‑150 (MQ=255)          tttGCCCGTTATGGCCTGCTGGTCGGCACCATAACCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGTGGTTTGTCGCCGGGGATTTTTCTg                                                                                                                         >  2:131786/1‑150 (MQ=255)             gCCCGTTATGGCCTGCTGGTCGGCACCATAACCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGTGCTTTTTCGCCGGGGATTATTCTGAcgg                                                                                                                     >  2:118091/1‑151 (MQ=255)                 gTTATGGCCTGCTGGTCGGCACCATAACCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGTGCTTTGTCGCCGGGAATTATTCTGACGGCGGt                                                                                                                 <  2:662657/151‑1 (MQ=255)                     tGGCCTGCTGGTCGGCACCATAACCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGGGCCTTTTCGCCGGGGGTTATTTTGTcggcgg                                                                                                                  >  1:262712/1‑146 (MQ=255)                       gCCTGCTGGTCGGCACCATAACCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGTGCTTTTTCGCCGGGGATTTTTTTGAcggcgggggtt                                                                                                             >  2:341213/1‑144 (MQ=255)                          tgctgGTCGACACCATAACCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGTGCTTTTGCGCCCGGGCTTATTCTCGCAGCGGTGATTACAg                                                                                                         >  2:131191/1‑150 (MQ=255)                           gctgGTCGGCACCATAACCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGTGCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAg                                                                                                         <  2:7345/149‑1 (MQ=255)                                    cACCATAACCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGTGCTTTTTCGCCGGGGATTATTTTTACGGCGGGGGTTTCAAGGCTGGTaaa                                                                                               >  2:149144/1‑150 (MQ=255)                                       cATAACCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGTGCTTTTTTGCCGGGGATTATTCTGACGGCGGGGGTTTAAGGGCTGGGAAata                                                                                             >  2:660042/1‑149 (MQ=255)                                        aTAACCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGTGCTTTGTCGCCGGGGATTATTTTGACGGCGGGGGTTACAGGGCTGGTAAata                                                                                             >  2:342056/1‑148 (MQ=255)                                           aCCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGTGCTTTTTCGCCGGGGATTTTTTTTACGGCGGGGGTTACAGGGCTgggtaa                                                                                               >  2:475947/1‑139 (MQ=255)                                           aCCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGTGCTTTTTCGCCGGGAATTATTTTGACGGCGGGGATTACAGGGCTGGTAAata                                                                                             >  2:378814/1‑145 (MQ=255)                                                cTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGTGCTTTGTCGCCGGGGATTTTTTTGACGGCGGGGGTTACAGGGCTGGGAAATATCAGCaaaaac                                                                                  >  1:90272/1‑147 (MQ=255)                                                     tgttgtGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGTGCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACg                                                                              <  1:447735/150‑1 (MQ=255)                                                           ggTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGTGCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGCGAt                                                                       <  1:208997/151‑1 (MQ=255)                                                                ttttGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGTGCTTTTTTGCCGGGGATTTTTTTTTCGGCGGGGGTTTCAGGGGTGGGGAATATTAGCAAAAACTTCgggggggtttggg                                                                  >  1:196437/1‑140 (MQ=255)                                                                       ttCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGTGCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGCGATTCGGGGCACgg                                                            <  2:457056/150‑1 (MQ=255)                                                                               tcGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGTGCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGCGATTCGGGGCACGGATGtttt                                                     <  1:514766/149‑1 (MQ=255)                                                                               tcGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGTGCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGCGATTCGGGGCACGGATGtttt                                                     <  1:157765/149‑1 (MQ=255)                                                                                cGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGTGCTTTTTCGCCGGGGATTATTCTTACGGCGGGGGTTACAGGGCTGGGAAAAAATAGAAAAAACTACGGGGGGGTTTGGGGGCCGGGGGtttttt                                                   >  2:578336/1‑150 (MQ=255)                                                                                            ccGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGTGCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAgg                                       <  1:286724/150‑1 (MQ=255)                                                                                                                  gTTCCCGCTCGGCAGTGGCGGTGCTTTTTCGCCGGGGAATATTCTGACGGCGGGGATTACAGGGCTGGGAAAAATCAACAAAACCTACGGGGCGGTTTGGGGGGCGGGTGGTTTTTTTTCGCAGCAGGGGGGCGGGGATACGCGGTAt                   >  2:629348/1‑148 (MQ=255)                                                                                                                    tCCCGCTCGGCAGTGGCGGTGCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTAtcgtc               <  1:629109/150‑1 (MQ=255)                                                                                                                        gCTCGGCAGTGGCGGTGCTTTTTCGCCGGGGATTATTTTTACGGCGGGGGTTACAGGGGTTGGAAATATTAGCAAAAACTACCGTGCGGTTCGGGGGGCGGGTGTTTTTTTTTTGCAGCAGGGCGCCGGGGAAAAGCGGTTtcgt                >  1:428273/1‑145 (MQ=255)                                                                                                                                gTGGCGGTGCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGc  <  2:67823/151‑1 (MQ=255)                                                                                                                                             |                                                                                                                                          CCGCAACGTTTTGCCCGTTATGGCCTGCTGGTCGGCACCATAACCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGTGCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGC  >  NC_000913/541135‑541412

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