Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	541,284	0	A	G	44.7%	50.3 / 10.3	38	G240G (GGA→GGG)	ybbY	putative uracil/xanthine transporter
Reads supporting (aligned to +/- strand):  ref base A (4/17);  new base G (17/0);  total (21/17)
Fisher's exact test for biased strand distribution p-value = 2.20e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGCAACGTTTTGCCCGTTATGGCCTGCTGGTCGGCACCATAACCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGT‑GCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGC‑GATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGA  >  NC_000913/541134‑541432                                                                                                                                                        |                                                                                                                                                      gCCGCAACGTTTTGCCCGTTATGGCCTGCTGGTCGGCACCATAACCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGT‑GCTTTGTCGCCGGGa                                                                                                                                                       <  1:738447/151‑1 (MQ=255)   cGCAACGTTTTGCCCGTTATGGCCTGCTGGTCGGCACCATAACCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGT‑GCTTTGTCGCCGGGAat                                                                                                                                                     <  2:380153/151‑1 (MQ=255)         gTTTTGCCCGTTATGGCCTGCTGGTCGGCACCATAACCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGT‑GCTTTTTCGCCGGGGATTATTCt                                                                                                                                               >  2:109187/1‑151 (MQ=255)         gTTTTGCCCGTTATGGCCTGCTGGTCGGCACCATAACCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGT‑GCTTTGTCGCCGGGGattatt                                                                                                                                                 >  2:613966/1‑149 (MQ=255)                  gTTATGGCCTGCTGGTCGGCACCATAACCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGTGGCTTTTTCCCGGGGATTTTTTTTgcggggg                                                                                                                                        >  2:921992/1‑144 (MQ=255)                     aTGGCCTGCTGGTCGGCACCATAACCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGT‑GCTTTGTCGCCGGGAATTATTCTGACGGCGGTGAt                                                                                                                                   <  2:461503/151‑1 (MQ=255)                       ggCCTGCTGGTCGGCACCATAACCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGT‑GGTTTTTCGCCGGGGATTATTTTTAcggcgggggttt                                                                                                                                 >  1:492786/1‑145 (MQ=255)                           tgctgGTCGGCACCATAACCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGT‑GCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACa                                                                                                                               <  1:372835/149‑1 (MQ=255)                               ggTCGGCACCATAACCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGT‑GCTTTTTCGCCGGGGATTATTCTGAcggcggggggtatagggct                                                                                                                          >  1:102562/1‑137 (MQ=255)                                 tCGGCACCATAACCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGT‑GCTTTTTTGCCGGGGATTATTTTGACGGCGGGGGGTACAggggt                                                                                                                          >  1:726412/1‑145 (MQ=255)                                  cGGCACCATAACCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGT‑GCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGt                                                                                                                       <  2:920939/150‑1 (MQ=255)                                    gCACCATAACCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGT‑GCTTTTTCGCCCGGGATTTTTTTGACGGCGGGGATTACAGGTGTgggta                                                                                                                     >  2:578851/1‑147 (MQ=255)                                     cACCATAACCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGT‑GCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAat                                                                                                                   <  1:477985/151‑1 (MQ=255)                                            aCCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGT‑GCTTTTTTGCCGGGGATTATTTTGACGGCGGGGGTTTCAGGGCTGGGaaa                                                                                                                    >  1:748474/1‑143 (MQ=255)                                            aCCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGT‑GCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCa                                                                                                            <  2:652369/151‑1 (MQ=255)                                                           tgGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGT‑GCTTTTTCGCCGGGGATTATTTTTACGGCGGGGGTTATAGGGCTGGGAAAAATTAACAAAACCCAccggggg                                                                                              >  2:11468/1‑144 (MQ=255)                                                             gTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGT‑GCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGC‑GAtt                                                                                          <  2:348562/151‑1 (MQ=255)                                                               aCTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGT‑GCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGC‑GATTCg                                                                                        <  1:462902/151‑1 (MQ=255)                                                                                 cGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGT‑GCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGC‑GATTCGGGGCACGGATGtttttt                                                                       <  1:790623/150‑1 (MQ=255)                                                                                 cGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGT‑GCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGC‑GATTCGGGGCACGGATGtttttt                                                                       <  2:222285/150‑1 (MQ=255)                                                                                            tcCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGT‑GCTTTGTCGCCGGGGATTATTCTGACGGCGGGGGTTTCAGGGCTGGGAAAAATTAGCAAAACCCACGGGGG‑GGTTTGGGGGACGGGGGgttttttttc                                                                   >  2:747403/1‑139 (MQ=255)                                                                                            tcCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGT‑GCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGC‑GATTCGGGGCACGGATGTTTTTTATCCGcagcag                                                            <  1:747421/150‑1 (MQ=255)                                                                                             ccGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGT‑GCTTTGTCGCCGGGAATTATTCTGGCGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGC‑GATTCGGGGCACGGATGTTTTTTATCCGcagcag                                                            <  1:617106/149‑1 (MQ=255)                                                                                             ccGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGT‑GCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGC‑GATTCGGGGCACGGATGTTTTTTATCCgcagca                                                             <  1:732720/148‑1 (MQ=255)                                                                                             ccGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGT‑GCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGC‑GATTCGGGGCACGGATGTTTTTTATCCGCAGCAgg                                                           <  2:121078/150‑1 (MQ=255)                                                                                                  gAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGT‑GCTTTTTCGCCGGGGATTATTCTGACGGCGGGGGTTACAGGGCTGGGAAAAATCAGCAATACCTACGGGGG‑GATTTGGGGCACGGGTGTTTTTTTTTCCCAGCAggggg                                                        >  1:597370/1‑146 (MQ=255)                                                                                                     ttGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGT‑GCTTTTTCGCCGGGGAATATTTTTACAGCGGGGGTTTCAGGGCTGGGAAAAATCAGCAATAACCACGGTGC‑GGTTCGGGGGACGGGTGTTTTTTTTTCGCCGCCGGGCGGAggg                                                   >  1:338392/1‑150 (MQ=255)                                                                                                      tGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGT‑GCTTTTTTGCCGGGGAATATTCTTACGGGGGGGGTTTCACGGCTGGGTAAAATCAACAATACCCACCGGGCGGATTGGGGGACCGGGTGTTTTTT‑TCCgcagcccggggccggg                                                   >  2:496998/1‑138 (MQ=255)                                                                                                       gCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGT‑GCTTTTTCGCCGGGGATTATTTTTACGGCGGGGGTTTCAGGGCTGGGAAATATCAGCAAAAACCACGGGGC‑GGTTTGGGGGACGGGTGGTTTTTTTTCCCAGCCGGGGGCCgggga                                                 >  1:28176/1‑148 (MQ=255)                                                                                                          cTGGCAGTGGTTCCCGCTCGGCAGTGGCGGT‑GCTTTTTCGCCGGGGATTATTCTTACGGCGGGGATTACAGGGCTGGGAAAAATCAACAAAACCCACGGGGC‑GGTTTGGGGGACGGGTGGTTTTTTTTCGCAGCCGGGGGCAggggaaaa                                              >  1:689593/1‑145 (MQ=255)                                                                                                            ggCAGTGGTTCCCGCTCGGCAGTGGCGGT‑GCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGC‑GATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATAcgc                                            <  2:432102/150‑1 (MQ=255)                                                                                                                      cccGCTCGGCAGTGGCGGT‑GCTTTGTCGCCGGGGATTATTCTGACGGCGGGGGTTATAGGGCTGGGAAATATCAGCAAAACCTACGGGGC‑GGTTCGGGGGACGGGTGTTTTTTTTCCGCCGCAGGGCGCAGGGAAAACGCGtttt                                       >  1:740140/1‑143 (MQ=255)                                                                                                                         gCTCGGCAGTGGCGGT‑GCTTTGTCGCCGGGGATTATTCTGACGGCGGGGATTACAGGTCTTGGAAAAATCAGCAAAAACTACCGGGG‑GGTTCGGGGGACGGGGGTTTTTTTTTCGCCGCAGGGGGCCGGGGAAACGCGtttttgt                                    >  1:679415/1‑140 (MQ=255)                                                                                                                             ggCAGTGGCGGT‑GCTTTTTCGCCGGGGATTATTCTGACGGCGGGGATTACAGGGGTGGGAAATATCAGCAATACCCACGGGGG‑GGTTTCGGGGGCGGGTGTTTTTTTTCCGCAGCAGGGGGCAGGGGATACGCGTTATCGTCGTAGCttt                           >  2:380831/1‑150 (MQ=255)                                                                                                                                      ggT‑GCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGC‑GATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCg                 >  1:563208/1‑151 (MQ=255)                                                                                                                                          gCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGC‑GATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGAt              <  1:813255/151‑1 (MQ=255)                                                                                                                                                  gCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGC‑GATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGAc       >  2:520017/1‑150 (MQ=255)                                                                                                                                                      ggAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGC‑GATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTg   >  1:329799/1‑150 (MQ=255)                                                                                                                                                      ggAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGC‑GATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGa  <  2:123298/151‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                      GCCGCAACGTTTTGCCCGTTATGGCCTGCTGGTCGGCACCATAACCGGCTGGTTGTTGTGGTACTTTTGCTTTCCTTCTTCGCACTCGCTCTCCGGTGAGTTGCACTGGCAGTGGTTCCCGCTCGGCAGTGGCGGT‑GCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGC‑GATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGA  >  NC_000913/541134‑541432

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