Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	111,230	(C)6→7	coding (187/390 nt)	mutT →	8‑oxo‑dGTP diphosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	111,224	1	.	C	88.2%	49.0 / 0.7	17	coding (181/390 nt)	mutT	8‑oxo‑dGTP diphosphatase
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base C (2/13);  total (2/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGA‑TTA‑C‑CCCC‑CAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGC  >  NC_000913/111092‑111341                                                                                                                                       |                                                                                                                        gAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAA‑AAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGGGCGTGAATTTCGGGGAGAAGTAGGGA‑TTA‑c                                                                                                                        >  1:159565/1‑133 (MQ=255)  aaaTCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGA‑TTA‑c                                                                                                                        <  1:596378/133‑1 (MQ=255)                 gTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGGGAAACGCCGGAACAGGCGGTGGCGGGTGAACCTCCGGAAAAAGGAGGGG‑TTAGC‑CCCC‑CAACCatttt                                                                                                        >  2:91999/1‑133 (MQ=255)                               gcgcACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGGGAAACGCCGGAACAGGCGGGGGTGGGGGGACTTCAGGGAGAAAGCGGGG‑GTA‑CACTCT‑CCACATTTTTTGCCATTTTaaaaa                                                                                          >  2:292867/1‑133 (MQ=255)                               gcgcACATAGCGAATCAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGA‑TTACC‑CCCC‑CAACATTTTTCGCTATTTGaaaaa                                                                                          <  1:292867/133‑1 (MQ=255)                                cgcACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGA‑TTT‑C‑CCCCGCAACATTTTTCGCTATTTGAAAAAc                                                                                         <  2:98729/133‑1 (MQ=255)                                         cGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGA‑TTACC‑CCCC‑CAAGATTTTTCGCTAGTTGAAAAACCGGAATATg                                                                                >  1:169973/1‑133 (MQ=255)                                         cGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGA‑TTACC‑CCCC‑CAACATTTTTCGCTATTTGAAAAACTGGAATATg                                                                                <  1:493132/133‑1 (MQ=255)                                            aTAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGA‑TTACC‑CCCC‑CAACATTTTTCGCTATTTGAAAAACTGGAATATGAAt                                                                             <  1:583152/133‑1 (MQ=255)                                                 cTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGA‑TTACC‑CCCC‑CAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCg                                                                        <  2:333330/133‑1 (MQ=255)                                                    gAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGA‑TTACC‑CCCC‑CAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGAc                                                                     <  2:511595/133‑1 (MQ=255)                                                      gTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGA‑TTACC‑CCCC‑CAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAg                                                                   <  1:605413/133‑1 (MQ=255)                                                                  tAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGA‑TTACC‑CCCC‑CAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACt                                                       <  2:371804/133‑1 (MQ=255)                                                                   aaaaTTGAAATGGGTTAAACGCCGTAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGA‑TTACC‑CCCC‑CAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACtg                                                      <  2:258242/133‑1 (MQ=255)                                                                                  gAAACGCCGGAACAGCCGGTGGTGCGTGTACTCCAGGAAGAAGTCGGGA‑TTACC‑CCCC‑CAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTc                                       <  2:223432/133‑1 (MQ=255)                                                                                     aCGCCGGAACAGGCGGTGGGGCGTGAACTTCAGGAAGAAGTCGGGA‑TTACC‑CCCC‑CAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCGaa                                    >  1:109136/1‑133 (MQ=255)                                                                                                  cGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGA‑TTACC‑CCCC‑CAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCGAACGCTGGGAAgggg                       <  1:423426/133‑1 (MQ=255)                                                                                                       gtgCGTGAACTTCAGGAAGAAGTCGGGA‑TTACC‑CCCC‑CAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCGAACGCTGGGAAGGGGAGCCg                  <  1:56586/133‑1 (MQ=255)                                                                                                         gCGTGAACTTAAGGACGAAGTCGGTATTTC‑C‑CCCC‑CAACATCTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCGAACGCTGGGAAGGGGAGCCGTg                <  2:159565/133‑1 (MQ=255)                                                                                                               aCTTCAGGAAGAAGTCGGGA‑TTACC‑CCCC‑CAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTaa          <  1:91999/133‑1 (MQ=255)                                                                                                                       aagaagTCGGGA‑TTACC‑CCCC‑CAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGc  <  2:109136/133‑1 (MQ=255)                                                                                                                                       |                                                                                                                        GAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGA‑TTA‑C‑CCCC‑CAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGC  >  NC_000913/111092‑111341

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