Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,093,876	C→T	A126V (GCG→GTG)	yqgF →	ribonuclease H‑like domain‑containing nuclease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,093,876	0	C	T	100.0%	45.5 / NA	18	A126V (GCG→GTG)	yqgF	ribonuclease H‑like domain‑containing nuclease
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (10/8);  total (10/8)

GACTGCCAGAGCGCGTAAATTTGCCAACCGTATTCATGGCCGTTTCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGCGGTTATTATTCTCGAAAGCTATTTCGAGCAGGGATATTAAGGCGATTTAAACGCCTGGCGGAGTGTAAATAATATCATCGCGTCTTATTGCCGGATGCGGCGTGAACACCTTATCCAGCACACATCTGGCAGCGGCTATAGGTCTGATAAGACGCGTTAGCGTCGCATCAGACATTTATTGCCTTTGTAGGCCCGATAAGCTTGCGCATCGGGCATGGCAACG  >  NC_000913/3093715‑3094099                                                                                                                                                                  |                                                                                                                                                                                                                                gACTGCCAGAGCGCGTAAATTTGCCAACCGTATTCATGGCCGTTTCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGTGGTTATTATTCTCGAAAGCTATTTCGAGCa                                                                                                                                                                                                   >  1:208823/1‑192 (MQ=255)           gcgcGTAAATTTGCCAACCGTATTCATGGCCGTTTCGGTGGTGAATTAAAGCTCCATGACGAGCGTCTTAGCACTGTTGAAGCCCGTGCCGGTCTGTTTGAACAGGGCTGCTATCGGGCGCTCAACAAAGGCAAAGTGGACTCTTCCTCTGTTTTTATGATTCTCGAAAGCTATTTCGAGCAGAGATATTAAgg                                                                                                                                                                                       <  1:459998/194‑1 (MQ=255)             gcgTAAATTTGCCAACCGTATGCATGGCCGTTTCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCGCAACAAAGGCAAAGTTGACTCTGCCGCTGTGGTTATTATTCTCGAAAGCTATTTCGAGCAGGGATATTAAGGCGATTTAAACGCCTGGCGGAGTGTAAATAATATCATCGCGTCTTATTGCCGGATgcggc                                                                                                                            <  1:271967/251‑1 (MQ=255)                 aaaTTTGCCAACCGTATTCATTGCCGTTTCGGTGTTGAAGGAAAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGTGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGTGGTTATTATTCTCGAAAGCTATTTCGAGCAGGGATATTAAGGCTATTTAAACGCCTGGCGGAGTTTAAATAATATCATCGCGTCTTATTGCCGGATGCGGCGtn                                                                                                                         <  1:53722/250‑2 (MQ=255)                                          gTTTCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGTGGTTATTATTCTCGaaa                                                                                                                                                                                                                >  1:190030/1‑138 (MQ=255)                                               gttgtTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTACGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGTGGTTATTATTCTCGAAAGCTATTTCGAGCAGGGCTATTAAGGCGATTTAAACGCCTGGCGGAGTGTAAATAATATCATCGCGTCTTATTGCCGGATGCGGCGGGAACACCTTATCCAGCACACATCTGGCAGCgg                                                                                          >  1:463939/3‑251 (MQ=255)                                                                    gACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGTGGTTATTATTCTCGAAAGCTATTTCGAGCAGGGATATTAAGGCGATTTAAACGCCTGGCGGAGTGTAAATAATATCATCGCGTCTTATTGCCGGATGCGGCGTGAACACCTTATCCAGCACACATCTGGCAGCGGCTAGAGGTCTGATAAGAcgc                                                                      <  1:175010/250‑1 (MQ=255)                                                                           gTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAACGTTGACTCTGCCTCTGTGGTTATTATTCTCGAAAg                                                                                                                                                                                                               >  1:83061/1‑106 (MQ=255)                                                                                             cccGATCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTATTTGGTAATTATTCTCGAAAGCTATTTCGAGCAGGGATATTAAGGCGAGTTAAACGCCTGGCGGAGNGTANnnnntann                                                                                                                                                     >  2:298725/1‑144 (MQ=255)                                                                                                      gTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCCAAGTTGACTCTGCCTCTGTGGTTATTATTCTCGAAAGCTATTTCGAGCAGGGATATTAAGGCGATTTAAACGCCTGGCGGAGTGTAAATAATATCATCGCGTCTTACTGCCGGATGCGGCGTGAACACCTTATCCAGCACACATCTGGCAGCGGCTATAGGTCTGATAAGACGCGTTAGCg                                                               >  1:506880/1‑223 (MQ=255)                                                                                                                  aNGGNGGCTANCGGGNGCTCAACAAAGGCAANGTTGACTCGGCCTCTGTGGNNANTATTCNCGaaa                                                                                                                                                                                                                <  2:190030/66‑1 (MQ=255)                                                                                                                     gcggcTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGTGGTTATTATTCTCGAAAGCTATATCGAGCAGGGATCTTCAGGCGATTTaac                                                                                                                                                                              >  1:272728/1‑96 (MQ=255)                                                                                                                           aTCGGGCGCTCAACAAAGGCAAAGTCGACTCTGCCTCTGTGGTTATTATCCTCGCAAGCTATTTCGAGCAGGGATCTTAAGGCGACTTAAACGCCTGGCGGAGTGTAAATAATATCATCGCGTCTTATTGCCGGCTGCGGCGTGAACACCTTATCCAGCACACATCTGGCACCGGCTATCGGTCTGCTAAGACGCGTTAGCGTCGCATCAGACACTTATTGCCTTTTCAGGCCCGATAAGCCTGCGAATCg              >  1:442808/1‑251 (MQ=255)                                                                                                                             cGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGTGGTTATTATTCTCGAAAGCTATTTCGAGCAGGGATATTAAGGCGATTTAAACGCCTGGCGGAGTGTAAATAAGATCATCGCGTCGTATTGCCGGATGCGGCGTGAACACCTTATCCAGCACACATCTGGCAGCGGCTATAGGTCTGAGAATACGCGTTTGCGTCGCATCAGACATTTATTGCCTTTGTAGGCCCGATAAGCTTGCGCATCgg             <  1:284330/250‑1 (MQ=255)                                                                                                                                     naacaaAGGCAACGTTGACTCTGCCTCTGTGGTTATTATTCTCGAAAGCTCTTTCGAGCAGGGATATTAAGGCGATTTAAACGCCTGGCGGAGTGTAAATAATATCATCGCGTCTTATTGCCGGATGCGGCGTGAACACCTTATCCAGCACACATCTGGCAGCGGCTATAGGTCTGATa                                                                            >  1:133011/2‑179 (MQ=255)                                                                                                                                     caacaaAGGCAAAGTTGACTCTGCCTCCGTTGTTATTATTCTCGAAAGCTATTTCGAGCAGGGATATTAAGGCGATTTAAACGCCTGGCGGAGTGTAAATAATATCATCGCGTCTTATTGCCGGATGCGGCGTGAACACCTTATCCAGCACACATCTGGCAGCGGCTATAGGTCTGATAAGACGCGTTAGCGTCGCATCAGACATTTATTGCCTTTGTAGGCCCGATAAGCTTGCGCATCTGGCATgg       <  1:135477/248‑1 (MQ=255)                                                                                                                                       acaaAGGCAAAGTTGACTCTGCCTCTGTGGTTATTATTCTCGAAAGCTATTTCGAGCAGGGATATTAAGGCGATTTAAACGCCTGGCGGAGTGTAAATAATATCATCGCGTCTTATTGCCGGATGCGGCGTGAACACCTTATCCAGCACACATCTGGCAGCGGCTATAGGTCTGATAAGACGCGTTAGCGTCGCATCAGACATTTATTGCCTTTGTAGGCCCGATAAGCTTGCGCATCGGGCATGGCAAcn  <  1:298725/251‑2 (MQ=255)                                                                                                                                        naaaGGCAAAGTTGACTCTGCCTCTGTGGTTATTATTCTCGAAAGCTATTTCGAGCAGGGATATTAAGGCGATTTAAACGCCTGGCGGAGTGTAAATAATATCATCGCGTCTTATTGCCGGATGCGGCGTGAACCCCTTATCCAGCACACATCTGGCAGCGGCAATAGGTCTGATAAGACGCGTTAGCGTCGCATCAGAc                                                    >  1:146352/2‑200 (MQ=255)                                                                                                                                                       ctctgcctctgtGGTTATTATTCTCGAAAGCTATTTCGAGCAGGGATATTAAGGCGATTTAAACGCCTGGCGGAGTGTAAATAATATCATCGCGTCTTATTGCCGGATgcg                                                                                                                              >  1:167272/1‑111 (MQ=255)                                                                                                                                                                  |                                                                                                                                                                                                                                GACTGCCAGAGCGCGTAAATTTGCCAACCGTATTCATGGCCGTTTCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGCGGTTATTATTCTCGAAAGCTATTTCGAGCAGGGATATTAAGGCGATTTAAACGCCTGGCGGAGTGTAAATAATATCATCGCGTCTTATTGCCGGATGCGGCGTGAACACCTTATCCAGCACACATCTGGCAGCGGCTATAGGTCTGATAAGACGCGTTAGCGTCGCATCAGACATTTATTGCCTTTGTAGGCCCGATAAGCTTGCGCATCGGGCATGGCAACG  >  NC_000913/3093715‑3094099

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