Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,170,791	A→G	R1058R (CGT→CGC)	mfd ←	transcription‑repair coupling factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,170,791	0	A	G	94.1%	45.3 / ‑3.8	17	R1058R (CGT→CGC)	mfd	transcription‑repair coupling factor
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (9/7);  total (9/8)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.25e-01

GGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGATGAAATCATCTG  >  NC_000913/1170671‑1170927                                                                                                                         |                                                                                                                                         ggCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTGCGCGCCGGATCTGgcagca                                                                                                                       >  1:629723/1‑140 (MQ=255)                       tcGGCAAATTCGATCGCCCCGCGTTGCTCTTTACCTTACAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTGCGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTc                                                                                                 <  1:220494/140‑1 (MQ=255)                                tCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGAt                                                                                        <  1:125404/140‑1 (MQ=255)                                 cGATCACCCCGCCTTGCTCATTGCCTTCCAGCTTCCTCATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTGCGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGAtc                                                                                       <  1:291569/140‑1 (MQ=255)                                       ccccGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTGCGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTcc                                                                                 <  2:629723/140‑1 (MQ=255)                                          cGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTGCGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGt                                                                              >  2:387697/1‑140 (MQ=255)                                                 ctcATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTGCGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGtttt                                                                       <  1:387697/140‑1 (MQ=255)                                                              aGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTGCGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGCCTTTGATTTCTTCCAGTTCGTTTTCCGTTTTGGCGCt                                                          >  1:222074/1‑140 (MQ=255)                                                                cTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTGCGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCATTTTGGCGCTgg                                                        >  2:28927/1‑140 (MQ=255)                                                                            cAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTGCGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGtt                                               <  1:626294/137‑1 (MQ=255)                                                                            cAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTGCGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGtt                                               >  2:626294/1‑137 (MQ=255)                                                                                   tGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTGCGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGa                                     >  2:551484/1‑140 (MQ=255)                                                                                                    aCGGGCAATATCCAGCAGGGTGTGCGCCGGAGCTGGCAGTAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGt                    <  2:222074/140‑1 (MQ=255)                                                                                                     cGGGCAATATCCAGCAGGGTGCGCGCCGGATCTGGCAGCAGGCCGAGGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTc                   <  1:468995/140‑1 (MQ=255)                                                                                                                  gcagGGTGCGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGATGAAatca      >  1:549788/1‑140 (MQ=255)                                                                                                                    agGGTGCGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGATGAAatcatc    >  1:441130/1‑140 (MQ=255)                                                                                                                      ggTGCGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGATGAAATCATCTg  >  2:130438/1‑140 (MQ=255)                                                                                                                         |                                                                                                                                         GGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGATGAAATCATCTG  >  NC_000913/1170671‑1170927

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