Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,212,597	T→C	G497G (GGT→GGC)	dnaG →	DNA primase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,212,597	0	T	C	100.0%	135.5 / NA	40	G497G (GGT→GGC)	dnaG	DNA primase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/20);  total (20/20)

TTATAGGGTTGCTGGTGCAAAATCCAGAATTAGCGACGTTGGTCCCGCCGCTTGAGAATCTGGATGAAAATAAGCTCCCTGGACTTGGCTTATTCAGAGAACTGGTCAACACTTGTCTCTCCCAGCCAGGTCTGACCACCGGGCAACTTTTAGAGCACTATCGTGGTACAAATAATGCTGCCACCCTTGAAAAACTGTCGATGTGGGACGATATAGCAGATAAGAATATTGCTGAGCAAACCTTCACCGACTCACTC  >  NC_000913/3212467‑3212723                                                                                                                                   |                                                                                                                               ttATAGGGTTGCTGGTGCGAAACCCAGAATAAGCGACGTTGGTCCCGCCGCTTGAGAATCTGGATGAAAATAAGCTCCCTGGACTTGGCTTATTCAGAGAACTGGTCAACACTTGTCTCTCCCAGCCAGGCCt                                                                                                                              <  1:223553/133‑1 (MQ=255)   ataGGGTTGCTGGTGCAAAATCCAGAATTAGCGACGTTGGTCCCGCCGCTTGAGAATCTGGATGAAAATAAGCTCCCTGGACTTGGCTTATTCAGAGAACTGGTCACCGCTTGTCTCGCCCATCCAGGCCTGa                                                                                                                            >  2:223553/1‑133 (MQ=255)   ataGGGTTGCTGGTGCAAAATCCAGAATTAGCGACGTTGGTCCCGCCGCTTGAGAATCTGGATGAAAATAAGCTCCCTGGACTTGGCTTATTCAGAGAACTGGTCAACACTTGTCTCTCCCAGCCAGGCCTGa                                                                                                                            >  1:443337/1‑133 (MQ=255)   ataGGGTTGCTGGTGCAAAATCCAGAATTAGCGACGTTGGTCCCGCCGCTTGAGAATCTGGATGAAAATAAGCTCCCTGGACTTGGCTTATGCAGAGAACTGGTGAACACTTGTCTCGCCCATCCAGGCCTGa                                                                                                                            >  2:144657/1‑133 (MQ=255)   ataGGGTTGCTGGGGCAAAATCCAGAATTAGCGACGTTGGTCCCGCCGCTTGAGAATCTGGATGAAAATAAGCTCCCTGGACTTGGCTTATTCAGAGAACTGGTCAACACTTGTCTCTCCCAGCCAGGCCTGa                                                                                                                            >  2:395550/1‑133 (MQ=255)     aGGGTTGCTGGTGCAAAATCCAGAATTAGCGACGTTGGTCCCGCCGCTTGAGAATCTGGATGAAAATAAGCTCCCTGGACTTGGCTTATTCAGAGAACTGGTCAACACTTGTCTCTCCCAGCCAGGCCTGacc                                                                                                                          >  1:18661/1‑133 (MQ=255)      gggTTGCTGGTGCAAAATCCAGAATTAGCGACGTTGGTCCCGCCGCTTGAGAATCTGGATGAAAATAAGCTCCCTGGACTTGGCTTATTCAGAGAACTGGTCAACACTTGTCTCTCCCAGCCAGGCCTGacca                                                                                                                         >  1:265567/1‑133 (MQ=255)       ggTTGCTGGTGCAAAATCCAGAATTAGCGACGTTGGTCCCGCCGCTTGAGAATCTGGATGAAAATAAGCTCCCTGGACTTGGCTTATTCAGAGAACTGGTCAACACTTGTCTCTCCCAGCCAGGCCTGaccac                                                                                                                        >  1:6234/1‑133 (MQ=255)         ttGCTGGTGCAAAATCCAGAATTAGCGACGTTGGTCCCGCCGCTTGAGAATCTGGATGAAAATAAGCTCCCTGGACTTGGCTTATTCAGAGAACTGGTCAACTCTTGTCTCTCCAAGACAGGCCGGACCACCg                                                                                                                      >  1:191616/1‑133 (MQ=255)         ttGCTGGTGCAAAATCCAGAATTAGCGACGTTGGTCCCGCCGCTTGAGAATCTGGATGAAAATAAGCTCCCTGGACTTGGCTTATTCAGAGAACTGGTCAACACTTGTCTCTCCCAGCCAGGCCTGACCACCg                                                                                                                      >  2:443651/1‑133 (MQ=255)         ttGCTGGTGCAAAATCCAGAATTAGCGACGTTGGTCCCGCCGCTTGAGAATCTGGATGAAAATAAGCTCCCTGGACTTGGCTTATTCAGAGAACTGGTCAACACTTGTCTCTCCCAGCCAGGCCTGACCACCg                                                                                                                      >  2:47454/1‑133 (MQ=255)           gCTGGTGCAAAATCCAGAATTAGCGACGTTGGTCCCGCCGCTTGAGAATCTGGATGAAAATAAGCTCTCTGGACTTGGCTTATTCAGAGAACTGGTCAACACTTGTCTCTCCCAGCCAGGCCTGACCACCggg                                                                                                                    >  2:46949/1‑133 (MQ=255)           gCTGGTGCAAAATCCAGAATTAGCGACGTTGGTCCCGCCGCTTGAGAATCTGGATGAAAATAAGCTCCCTGGACTTGGCTTATTCAGAGAACTGGTCAACACTTGTCTCTCCCAGCCAGGCCTGACCACCggg                                                                                                                    >  1:42648/1‑133 (MQ=255)           gCTGGTGCAAAATCCAGAATTAGCGACGTTGGTCCCGCCGCTTGAGAATCTGGATGAAAATAAGCTCCCTGGACTTGGCTTATTCAGAGAACTGGTCAACACTTGTCTCTCCCAGCCAGGCCTGACCACCggg                                                                                                                    >  1:367483/1‑133 (MQ=255)           gCTGGTGCAAAATCCAGAATTAGCGACGTTGGTCCCGCCGCTTGAGAATCTGGATGAAAATAAGCTCCCTGGACTTGGCTTATTCAGAGAAATGGTCAACGCTTGTCTCTCCCAGCCTGGCCAGACCACCggg                                                                                                                    >  2:162354/1‑133 (MQ=255)                   aaaaTCCAGAATTAGCGACGTTGGTCCCGCCGCTTGAGAATCTGGATGAAAATAAGCTCCCTGGACTTGGCTTATTCAGAGAACTGGTCAACACTTGTCTCTCCCAGCCAGGCCTGACCACCGGGCAACtttt                                                                                                            <  1:401020/133‑1 (MQ=255)                    aaaTCCAGCATTAGCTACATTGCTCCCGCCACTTGAGAATCTGGATGAAAATAAGCTCCCTGGACTTGGCTTATTCAGAGAACTGGTCAACACTTGTCTCTCCCAGCCAGGCCTGACCACCGGGCAACTTTTa                                                                                                           <  1:242219/133‑1 (MQ=255)                        ccAGAATTAGCGACGTTGGTCCCGCCGCTTGAGAATCTGGATGAAAATAAGCTCCCTGGACTTGGCTTATTCAGAGAACTGGTCAACACTTGTCTCTCCCAGCCAGGCCTGACCACCGGGCAACTTTTAGAGc                                                                                                       <  2:477811/133‑1 (MQ=255)                              ttAGCGACGTTGGTCCCGCCGCTTGAGAATCTGGATGAAAATAAGCTCCCTGGACTTGGCTTATTCAGAGAACTGGTCAACACTTGTCTCTCCCAGCCAGGCCTGACCACCGGGCAACTTTTAGAGCACTATc                                                                                                 <  1:47454/133‑1 (MQ=255)                               tAGCGACGTTGGTCCCGCCGCTTGAGAATCTGGATGAAAATAAGCTCCCTGGACTTGGCTTATTCAGAGAACTGGTCAACACTTGTCTCTCCCAGCCAGGCCTGACCACCGGGCAACTTTTAGAGCACTATCg                                                                                                <  2:367483/133‑1 (MQ=255)                               tAGCGACGATGGTCCTGCTGCTTGAGAATCTGGATGAAAATAAGCTCCCTGGACTTGGCTTATTCAGAGAACTGGTCAACACTTGTCTCTCCCAGCCAGGCCTGACCACCGGGCAACTTTTAGAGCACTATCg                                                                                                <  1:220046/133‑1 (MQ=255)                                                    ttGAGAATCTGGATGAAAATAAGCTCCCTGGACTTGGCTTATTCAGAGAACTGGTCAACACTTGTCTCTCCCAGCCAGGCCTGACCACCGGGCAACTTTTAGAGCACTATCGTGGTACAAATAATGCTGCCAc                                                                           <  2:443337/133‑1 (MQ=255)                                                    ttGAGAATATGGATGCAACTAACCTCCCTGCACTTGGCTTATTCAGAGAACTGGTCAACACTTGTCTCTCCCAGCCAGGCCTGACCACCGGGCAACTTTTAGAGCACTATCGTGGTACAAATAATGCTGCCAc                                                                           <  1:144657/133‑1 (MQ=255)                                                            cTGGATGAAAATAAGCTCCCTGGACTTGGCTTATTCAGAGAACTGGTCAACACTTGTCTCTCCCAGCCAGGCCTGACCACCGGGCAACTTTTAGAGCACTATCGTGGTACAAATAATGCTGCCACCCTTGaaa                                                                   <  2:42648/133‑1 (MQ=255)                                                                       tAAGCTCCCTGGACTTGGCTTATTCAGAAAACCGGTCAACACTTGTCTCTCCCAGCCAGGCCTGACCACCGGGCAACTTTTAGAGCACTATCGTGGTACAAATAATGCTGCCACCCTTGAAAAACTGTCGAtg                                                        <  2:191616/133‑1 (MQ=255)                                                                            tCCCTGGACTTGGCTTATTCAGAGAACTGGTCAACACTTGTCTCTCCCAGCCAGGCCTGACCACCGGGCAACTTTTAGAGCACTATCGTGGTACAAATAATGCTGCCACCCTTGAAAAACTGTCGATGTGGGa                                                   <  1:443651/133‑1 (MQ=255)                                                                            tCCCTGGACTTGGCTTATTCAGAGAACTGGTCAACACTTGTCTCTCCCAGCCAGGCCTGACCACCGGGCAACTTTTAGAGCACTATCGTGGTACAAATAATGCTGCCACCCTTGAAAAACTGTCGATGTGGGa                                                   <  1:395550/133‑1 (MQ=255)                                                                                    cTTGGCTTATTCAGAGAACTGGTCAACACTTGTCTCTCCCAGCCAGGCCTGACCACCGGGCAACTTTTAGAGCACTATCGTGGTACAAATAATGCTGCCACCCTTGAAAAACTGTCGATGTGGGACGATATAg                                           <  2:265567/133‑1 (MQ=255)                                                                                      tGGCTTATTCAGAGAACTGGTCAACACTTGTCTCTCCCAGCCAGGCCTGACCACCGGGCAACTTTTAGAGCACTATCGTGGTACAAATAATGCTGCCACCCTTGAAAAACTGTCGATGTGGGACGATATAGCa                                         <  2:18661/133‑1 (MQ=255)                                                                                       ggCTTATTCAGAGAACTGGTCAACACTTGTCTCTCCCAGCCAGGCCTGACCACCGGGCAACTTTTAGAGCACTATCGTGGTACAAATAATGCTGCCACCCTTGAAAAACTGTCGATGTGGGACGATATAGCAg                                        >  2:458282/1‑133 (MQ=255)                                                                                         cTTATTCAGAGAACTGGTCAACACTTGTCTCTCCCAGCCAGGCCTGACCACCGGGCAACTTTTAGAGCACTATCGTGGTACAAATAATGCTGCCACCCTTGAAAAACTGTCGATGTGGGACGATATAGCAGAt                                      <  2:6234/133‑1 (MQ=255)                                                                                         cTTATTCAGAGAACTGGTCAACACTTGTCTCTCCCAGCCAGGCCTGACCACCGGGCAACTTTTAGAGCACTATCGTGGTACAAATAATGCTGCCACCCTTGAAAAACTGTCGATGTGGGACGATATAGCAGAt                                      <  1:46949/133‑1 (MQ=255)                                                                                                   gaACTGGTCAACACTTGTCTCTCCCAGCCAGGCCTGACCACCGGGCAACTTTTAGAGCACTATCGTGGTACAAATAATGCTGCCACCCTTGAAAAACTGTCGATGTGGGACGATATAGCAGATAAGAATATTg                            <  1:430014/133‑1 (MQ=255)                                                                                                        ggTCAACACTTGTCTCTCCCAGCCAGGCCTGACCACCGGGCAACTTTTAGAGCACTATCGTGGTACAAATAATGCTGCCACCCTTGAAAAACTGTCGATGTGGGACGATATAGCAGATAAGAATATTGCTGAg                       >  2:126919/1‑133 (MQ=255)                                                                                                                   gTCTCTCCCAGCCAGGCCTGACCACCGGGCAACTTTTAGAGCACTATCGTGGTACAAATAATGCTGCCACCCTTGAAAAACTGTCGATGTGGGACGATATAGTAGATAAGGATAGTGGTGAGCAAACCTTCAc            >  1:157021/1‑133 (MQ=255)                                                                                                                   gTCTCTCCCAGCCAGGCCTGACCACCGGGCAACTTTTAGAGCACTATCGTGGTACAAATAATGCTGCCACCCTTGAAAAACTGTCGATGTGGGACGATATAGCAGATAAGAATATTGCTGAGCAAACCTTCAc            >  2:49193/1‑133 (MQ=255)                                                                                                                   gTCTCTCCCAGCCAGGCCTGACCACCGGGCAACTTTTAGAGCACTATCGTGGTACAAATAATGCTGCCACCCTTGAAAAACTGTCGATGTGCGACGATATTGCAGAAACCAATATTGCTGAGCTAACCTTctc            >  2:239048/1‑131 (MQ=255)                                                                                                                        tcCCAGCCAGGCCTGACCACCGGGCAACTTTTAGAGCACTATCGTGGTACAAATAATGCTGCCACCCTTGAAAAACTGTCGATGTGGGACGATATAGCAGATGAGAATATGGCTCAGTAAACCTTCACCGact       >  1:250932/1‑133 (MQ=255)                                                                                                                          ccagccagGCCTGACCACCGGGCAACTTTTAGAGCACTATCGTGGTACAAATAATGCTGCCACCCTTGAAAAACTGTCGATGTGGGACGATATAGCAGATAAGAATATTGCTGAGCAAACCTTCACCGactca     >  1:371506/1‑133 (MQ=255)                                                                                                                             gccagGCCTGACCACTGGTCAACTTTTAGACCACTATCGTGGTACAAATAATGCTGCCACCCTTGAAAAACTGTCGATGTGGGACGATATAGCAGATAAGAATATTGCTGAGCAAACCTTCACCGactcactc  <  1:239048/133‑1 (MQ=255)                                                                                                                             gccagGCCTGACCACCGGGCAACTTTTAGAGCACTATCGTGGTACAAATAATGCTGCCACCCTTGAAAAACTGTCGATGTGGGACGATATAGCAGATAAGAATATTGCTGAGCAAACCTTCACCGactcactc  <  2:371506/133‑1 (MQ=255)                                                                                                                             gccagGCCTGACCACCGGGCAACTTTTAGAGCACTATCGTGGTACAAATAATGCTGCCACCCTTGAAAAACTGTCGATGTGGGACGATATAGCAGATAAGAATATTGCTGAGCAAACCTTCACCGactcactc  <  1:49193/133‑1 (MQ=255)                                                                                                                                   |                                                                                                                               TTATAGGGTTGCTGGTGCAAAATCCAGAATTAGCGACGTTGGTCCCGCCGCTTGAGAATCTGGATGAAAATAAGCTCCCTGGACTTGGCTTATTCAGAGAACTGGTCAACACTTGTCTCTCCCAGCCAGGTCTGACCACCGGGCAACTTTTAGAGCACTATCGTGGTACAAATAATGCTGCCACCCTTGAAAAACTGTCGATGTGGGACGATATAGCAGATAAGAATATTGCTGAGCAAACCTTCACCGACTCACTC  >  NC_000913/3212467‑3212723

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