Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,474,885	G→C	intergenic (+29/‑99)	yfdQ → / → yfdR	DUF2303 domain‑containing protein YfdQ/5'‑deoxynucleotidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,474,885	0	G	C	89.5%	45.6 / 0.3	20	intergenic (+29/‑99)	yfdQ/yfdR	DUF2303 domain‑containing protein YfdQ/5'‑deoxynucleotidase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base C (10/7);  total (12/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01

ATTTAAATGCGTTCCGTTTGAAGGTCTGAAAGAACGTCCGTTTAAATTACGCCTCAGCATTATCACTGGCGATCGTCCTGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCAC‑CGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATCGCGGTT  >  NC_000913/2474649‑2475070                                                                                                                                                                                                                                             |                                                                                                                                                                                           aTTTAAATGCGTTCCGTTTGAAGGTCTGAAAGAACGTCCGTTTAAATTACGCCTCAGCATTATCACTGGCGATCGTCCTGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGGGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCTCCGCCTGATTTCATTACTGCAAATGCCCCTGCTGGGCCATTTATGGAAACg                                                                                                                                                                              <  1:441209/251‑1 (MQ=255)    tAAATGCGTTCCGTTTGAAGGTCTGAAAGAACGTCCGTTTAAATTACGCCTCAGCATTATCACTGGCGATCGTCCTGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCCGGAAGATATGGCTAACGAATTTCGTGCTCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTa                                                                                                                                                                            >  2:431907/1‑250 (MQ=255)             tCCGTTTGAAGGTCTGAAAGAACGTCCGTTTAAATTACGCCTCAGCATTATCACTGGCGATCGTCCTGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGCCATTTATGGAAACGTAATTAACTCa                                                                                                                                                                   >  1:232178/1‑250 (MQ=255)                  ttGAAGGTCTGAAAGAACGTCCGTTTAAATTACGCCTCAGCTTTATCACTGGCGATCGTCCTGTACTGGTTATGCGCATTATTCAGATGGAAGCGGTGCAGGAAGATATGGCTAACGCATTTCGTGATCTGCTTGTTGAGAAATTCACAGACAGCACAGTAGACACCTTTATTGGTACTTTCACCTCATGATTTCATTACTGCAAATGCCCCTGCGGGGCCATTTACTGAAACGTACTTAACTCaataat                                                                                                                                                              >  2:11894/1‑250 (MQ=255)                        gTCTGAAAGAACGTCCGTTTAAATTACGCCTCAGCATTATCACTGGCGATCGTACTGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAACGACAGCAAAGTAGAAACCTTTATTGGTACTTTCAACGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGCCATTTATGGCAACGTAATTAACTCAATAATCAC‑Cgg                                                                                                                                                       >  2:425907/1‑250 (MQ=255)                                           tAAATTACGCCTCAGCATTATCACTGGCGATCGTCCTGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGCCATTTATGGAAACGTAATTAACTCAATAATCAC‑CGGATGGTGAGGGCTTCCtttt                                                                                                                                    >  1:153747/1‑250 (MQ=255)                                                          cATTATCACTGGCGATCGTCCTGTACTGGTTCGGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGTTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGCCATTTATGGAAACGTAATTAACTCAATAATCAC‑CGGATGGTGAGGGCTGCCTTGTACCCAAACTCAGcgcg                                                                                                                    <  2:232178/251‑1 (MQ=255)                                                           aTTATCACTGGCGATCGTCCTGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAACTGCCCCTGCGGGGCCATTTATGGAAACGTAATTAACTCACTAATCAC‑CGGATGGTGAGGGCTTCCTTTTACCCAACCTAAGCGCgg                                                                                                                   >  2:23043/1‑251 (MQ=255)                                                                      cGATCGTCCTGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGCCATTTATGGAAACGTAATTAACTCAATAATCAC‑CGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCata                                                                                                        <  1:5055/251‑1 (MQ=255)                                                                                gTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGCCATTTATGGAAACGTAATTAACTCAATAATCAC‑CGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGgaga                                                                                              <  1:23043/251‑1 (MQ=255)                                                                                  aCTGTTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTACGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTTCGGGGCCATTTATGGAAACGTAATTAACTCAATAATCAC‑CGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCATCGCATATACGTGGAGaa                                                                                             <  2:165125/250‑1 (MQ=255)                                                                                       ttCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTTCGGGGCCATTTATGGAAACGTAATTAACTCAATAATCAC‑CGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAAt                                                                                       >  1:11184/1‑251 (MQ=255)                                                                                               attattCAGCTGGAAGCGATGCAGGAAGATATGGCTAACGAATTTCGTGATCTGATTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGATCATTTCACCGCCTGATTTCATTACTGCAAATGCACCTGAGGGGACATTTATGGAACCGTAATTAACTCAATACT‑ACACGGATGGATAGGGCTTCCTTTTACCCAAACTTAGCGCTGTGCAGTGCATATACGCGGAGAACAAAATGTCAtttt                                                                               >  2:407534/1‑250 (MQ=255)                                                                                                          ggAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGCCATTTATGGAAACGTAATTAACTCAATAATCAC‑CGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACtttt                                                                    <  1:396517/251‑1 (MQ=255)                                                                                                                                  tAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGCCATTTATGGAAACGTAATTAACTCaata                                                                                                                                                                >  2:190943/1‑136 (MQ=255)                                                                                                                                   aaCGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTAACCGCCTGATTTCATTACTGCAAATGCCCCTGAGGGGCCATTTATGGAAACGAAATTAACTAACTAATAAC‑CGGTTGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGAAAGCATTTTTATTATGAc                                            >  1:440921/1‑250 (MQ=255)                                                                                                                                           tCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCAC‑CGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTTTATTATGACAAGataaa                                    <  1:431907/250‑1 (MQ=255)                                                                                                                                                       tgttgAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTAATTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGCCATTTATGGAAACGTAATTAACTCAATAATCAC‑CGGATGGTGAGGGCTTCATTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTTTATTATGACAAGATAAATAAAGACGACAt                        <  1:271083/250‑1 (MQ=255)                                                                                                                                                                         cAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGCCATTTATGGAAACGTAATTAACTCAATAATCAC‑CGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATcgcg     >  1:435231/1‑251 (MQ=255)                                                                                                                                                                            cAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGCCATTTATGGAAACGTAATTAACTCAATAATCAC‑CGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATCGCGGtt  >  1:84457/1‑251 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                           ATTTAAATGCGTTCCGTTTGAAGGTCTGAAAGAACGTCCGTTTAAATTACGCCTCAGCATTATCACTGGCGATCGTCCTGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCAC‑CGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATCGCGGTT  >  NC_000913/2474649‑2475070

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