Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,196,428	A→G	V109A (GTA→GCA)	rfaE ←	fused heptose 7‑phosphate kinase/heptose 1‑phosphate adenyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,196,428	0	A	G	100.0%	73.7 / NA	26	V109A (GTA→GCA)	rfaE	fused heptose 7‑phosphate kinase/heptose 1‑phosphate adenyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/11);  total (15/11)

GTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGATCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTCGGCCAGAGATTTACTCAGCGCGCGCGCTGCATCGTCAATGCCCGTCAACCCGACCAGGCG  >  NC_000913/3196301‑3196552                                                                                                                                |                                                                                                                             gTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGATCAACACCTTCGAAACCTTCTTCAAAATCCAGAGGGATCTGTTGTTGGTTTCGGGAAAGTGCCCGt                                                                                                                         >  2:227084/1‑133 (MQ=255) gTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGATCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTGCCCGt                                                                                                                         >  1:417231/1‑133 (MQ=255) gTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGATCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTATTGGTTCCGGGAAAGTGCCCGt                                                                                                                         >  1:235004/1‑133 (MQ=255)     gAAAGCACCAGCGCGCCAATCGAACTCACCGCCTGATTAATCCGCTCATGCAGCGGCTGCGGATCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTGCTCGTAAtt                                                                                                                     <  2:438976/133‑1 (MQ=255)            ccAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGATCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTGCCCGTAATTTGGTAAt                                                                                                              >  2:428941/1‑133 (MQ=255)                   gcCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGATCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTGCCCGTAATTTGGTAATGGTCGGa                                                                                                       <  2:110461/133‑1 (MQ=255)                    ccAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGATCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTCTTGGTTGCGGGAAAGTGCCCGTAATTTGGTGATGGTCGGAt                                                                                                      >  2:205295/1‑133 (MQ=255)                        tCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGATCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTGCCCGTAATTTGGTAATGGTCGGATGCGt                                                                                                  <  2:287256/133‑1 (MQ=255)                         cGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGATCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTGCCCGTAATTTGGTAATGGTCGGATGCGTc                                                                                                 <  2:186484/133‑1 (MQ=255)                             cTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGATCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTGCCCGTAATTTGGTAATGGTCGGATGCGTCGGTa                                                                                             <  1:366727/133‑1 (MQ=255)                                  cgcCTGATTAATCCGCTCGTGCAGCGGCTGCGGATCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTGCCCGTCATTTGGTAATGGTCGGTTGCGTCGGTACAGaa                                                                                        >  2:215938/1‑133 (MQ=255)                                               cGCTCGTGTAGGGGCCGCTGATTAACACCTGCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTGCCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCAtt                                                                           <  1:220437/133‑1 (MQ=255)                                                  tCGTGCAGCGGCTGCGGATCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTGCCCGTAATTTGGTAATGGTCGGAGGCGTCGGTACAGAAGCGAAGTAGCATttga                                                                        >  1:132183/1‑133 (MQ=255)                                                            gCTGCGGATCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTGCCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTc                                                              >  2:12314/1‑133 (MQ=255)                                                                             tCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTGCCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTCGGCCAGAGATTTACTCa                                             >  2:290958/1‑133 (MQ=255)                                                                                               tCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTGCCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTCGGCCAGAGATTTACTCAGCGCGCGCGCTGCATCGt                           <  1:380570/133‑1 (MQ=255)                                                                                                                ttggttgCGGGAAAGTGCCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTCGGCCAGAGATTTACTCAGCGCGCGCGCTGCATCGTCAATGCCCGTCAACCCg          >  2:281911/1‑133 (MQ=255)                                                                                                                  ggttgCGGGACAGTGCCCGTAATTTGGTAATGTTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTCGGCCAGAGATTTACTCAGCGCGCGCGCTGCATCGTCAATGCCCGTCAACCCGAc        <  2:132183/133‑1 (MQ=255)                                                                                                                  ggttgCGGGAAAGTGCCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTCGGCCAGAGATTTACTCAGCGCGCGCGCTGCATCGTCAATGCCCGTCAACCCGAc        >  1:88315/1‑133 (MQ=255)                                                                                                                  ggttgCGGGAAAGTGCCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTCGGCCAGAGATTTACTCAGCGCGCGCGCTGCATCGTCAATGCCCGTCAACCCGAc        >  1:475278/1‑133 (MQ=255)                                                                                                                  ggttgCGGGAAAGTGCCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATCTGACGTTGACGTCGGCCAGAGATTTACTCAGCGCGCGCGCTGCATCGTCAATGCCCGTCAACCCGAc        >  1:289398/1‑133 (MQ=255)                                                                                                                  gattgCGGGACAGTGCCCGTAATTCGATAATGTTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTCGGCCAGAGATTTACTCAGCGCGCGCGCTGCATCGTCAATGCCCGTCAACCCGAc        <  2:235004/131‑1 (MQ=255)                                                                                                                    ttgCGGGAAAGTGCCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTCGGCCAGAGATTTACTCAGCGCGCGCGCTGCATCGTCAATGCCCGTCAACCCGACCa      >  1:30351/1‑133 (MQ=255)                                                                                                                       cGGGAAAGTGCCCGTTATTTGGAAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTCGGCCAGAGATTTACTCAGCGCGCGCGCTGCATCGTCAATGCCCGTCAACCCGACCAGGc   <  1:227084/133‑1 (MQ=255)                                                                                                                       cGGGAAAGTGCCCGTAATTTGGTAATGGTCGGGTGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTCGGCCAGAGATTTACTCAGCGCGCGCGCTGCATCGTCAATGCCCGTCAACCCGACCAGGc   >  2:356578/1‑133 (MQ=255)                                                                                                                        gggAAAGTGCCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTCGGCCAGAGATTTACTCAGCGCGCGCGCTGCATCGTCAATGCCCGTCAACCCGACCAGGCg  <  1:12314/133‑1 (MQ=255)                                                                                                                                |                                                                                                                             GTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGATCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTCGGCCAGAGATTTACTCAGCGCGCGCGCTGCATCGTCAATGCCCGTCAACCCGACCAGGCG  >  NC_000913/3196301‑3196552

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