Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	541,284	0	A	G	55.6%	14.0 / 34.4	27	G240G (GGA→GGG)	ybbY	putative uracil/xanthine transporter
Reads supporting (aligned to +/- strand):  ref base A (2/10);  new base G (15/0);  total (17/10)
Fisher's exact test for biased strand distribution p-value = 7.82e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.55e-05

CCGCAACGTTTTGCCCGTTATGGCCTGCTGGTCGGCACCATAACCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGTGCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGAT  >  NC_000913/541135‑541433                                                                                                                                                      |                                                                                                                                                      ccGCAACGTTTTGCCCGTTATGGCCTGCTGGTCGGCACCATAACCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGTGCTTTGTCGCCGGGaa                                                                                                                                                      <  1:533416/151‑1 (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:157765/149‑1 (MQ=255)                                                                               tcGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGTGCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGCGATTCGGGGCACGGATGtttt                                                                          <  1:514766/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)                                                                                                                                                     gAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGAt  >  1:640738/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                      CCGCAACGTTTTGCCCGTTATGGCCTGCTGGTCGGCACCATAACCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGTGCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGAT  >  NC_000913/541135‑541433

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