Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,211,471	G→A	G122D (GGT→GAT)	dnaG →	DNA primase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,211,471	0	G	A	100.0%	133.9 / NA	34	G122D (GGT→GAT)	dnaG	DNA primase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/18);  total (16/18)

CTGATGAACTACGACAAGCTCGAGTTCGTCGAAACGGTCGAAGAGCTGGCAGCAATGCACAATCTTGAAGTGCCATTTGAAGCAGGCAGCGGCCCCAGCCAGATCGAGCGCCATCAGAGGCAAACGCTTTATCAGTTGATGGACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGTATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGCCCCCCGGCTGGGACAACGTCCTGAAGCG  >  NC_000913/3211326‑3211618                                                                                                                                                  |                                                                                                                                                    cTGATGAACTACGACAAGCTCGAGTTCGTCGAAACGGTCGAAGAGCTGGCAGCAATGCACAATCTTGAAGTGCCATTTGAAGCAGGCAGCGGCCCCAGCCAGATCGAGCGCCATCAGAGGCAAACGCTTTATCAGTTGATGGACGATCTg                                                                                                                                                 >  2:207783/1‑150 (MQ=255)                aaGCTCGAGTTCGTCGAAACGGTCGAAGAGCTGGCAGCAATGCACAATCTTGAAGTGCCATTTGAAGCAGGCAGCGGCCCCAGCCAGATCGAGCGCCATCAGAGGCAAACGCTTTATCAGTTGATGGACGATCTGAATACGTTTTACcaac                                                                                                                                 >  2:561579/1‑151 (MQ=255)                        gTTCGTCGAAACGGTCGAAGAGCTGGCAGCAATGCACAATCTTGAAGTGCCATTTGAAGCAGGCAGCGGCCCCAGCCAGATCGAGCGCCATCAGAGGCAAACGCTTTATCAGTTGATGGACGATCTGAATACGTTTTACCAACAATCTTTa                                                                                                                         <  1:103478/151‑1 (MQ=255)                           cgtcgAAACGGTCGAAGAGCTGGCAGCAATGCACAATCTTGAAGTGCCATTTGAAGCAGGCAGCGGCCCCAGCCAGATCGAGCGCCATCAGAGGCAAACGCTTTATCAGTTGATGGACGATCTGAATACGTTTTACCAACAATCTTTaca                                                                                                                       <  1:455600/150‑1 (MQ=255)                            gtcgAAACGGTCGAAGAGCTGGCAGCAATGCACAATCTTGAAGTGCCATTTGAAGCAGGCAGCGGCCCCAGCCAGATCGAGCGCCATCAGAGGCAAACGCTTTATCAGTTGATGGACGATCTGAATACGTTTTACCAACAATCTTTacaac                                                                                                                     <  1:266285/151‑1 (MQ=255)                                   cGGTCGAAGAGCTGGCAGCAATGCACAATCTTGAAGTGCCATTTGAAGCAGGCAGCGGCCCCAGCCAGATCGAGCGCCATCAGAGGCAAACGCTTTATCAGTTGATGGACGATCTGAATACGTTTTACCAACAATCTTTACAACAACCTGt                                                                                                              <  2:493934/151‑1 (MQ=255)                                        gAAGAGCTGGCAGCAATGCACAATCTTGAAGTGCCATTTGAAGCAGGCAGCGGCCCCAGCCAGATCGAGCGCCATCAGAGGCAAACGCTTTATCAGTTGATGGACGATCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCa                                                                                                         <  1:300825/151‑1 (MQ=255)                                         aaGAGCTGGCAGCAATGCACAATCTTGAAGTGCCATTTGAAGCAGGCAGCGGCCCCAGCCAGATCGAGCGCCATCAGAGGCAAACGCTTTATCAGTTGATGGACGATCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCAc                                                                                                        >  1:441831/1‑151 (MQ=255)                                         aaGAGCTGGCAGCAATGCACAATCTTGAAGTGCCATTTGAAGCAGGCAGCGGCCCCAGCCAGATCGAGCGCCATCAGAGGCAAACGCTTTATCAGTTGATGGACGATCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCAc                                                                                                        <  1:371266/151‑1 (MQ=255)                                          agagCTGGCAGCAATGCACAATCTTGAAGTGCCATTTGAAGCAGGCAGCGGCCCCAGCCAGATCGAGCGCCATCAGAGGCAAACGCTTTATCAGTTGATGGACGATCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACg                                                                                                       >  1:80715/1‑151 (MQ=255)                                            agCTGGCAGCAATGCACAATCTTGAAGTGCCATTTGAAGCAGGCAGCGGCCCCAGCCAGATCGAGCGCCATCAGAGGCAAACGCTTTATCAGTTGATGGACGATCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTc                                                                                                     <  1:647787/151‑1 (MQ=255)                                            agCTGGCAGCAATGCACAATCTTGAAGTGCCATTTGAAGCAGGCAGCGGCCCCAGCCAGATCGAGCGCCATCAGAGGCAAACGCTTTATCAGTTGATGGACGATCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTc                                                                                                     >  1:347457/1‑151 (MQ=255)                                                 gcagcaATGCACAATCTTGAAGTGCCATTTGAAGCAGGCAGCGGCCCCAGCCAGATCGAGCGCCATCAGAGGCAAACGCTTTATCAGTTGATGGACGATCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTgcgc                                                                                                >  2:182757/1‑151 (MQ=255)                                                  cagcaATGCACAATCTTGAAGTGCCATTTGAAGCAGGCAGCGGCCCCAGCCAGATCGAGCGCCATCAGAGGCAAACGCTTTATCAGTTGATGGACGATCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTgcgcg                                                                                               >  2:80718/1‑151 (MQ=255)                                                                    aaGTGCCATTTGAAGCAGGCAGCGGCCCCAGCCAGATCGAGCGCCATCAGAGGCAAACGCTTTATCAGTTGATGGACGATCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGTATCTGGAAAAAcg                                                                             <  2:533603/151‑1 (MQ=255)                                                                         ccATTTGAAGCAGGCAGCGGCCCCAGCCAGATCGAGCGCCATCAGAGGCAAACGCTTTATCAGTTGATGGACGATCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGTATCTGGAAAAACGCGGAt                                                                        <  2:441820/151‑1 (MQ=255)                                                                           aTTTGAAGCAGGCAGCGGCCCCAGCCAGATCGAGCGCCATCAGAGGCAAACGCTTTATCAGTTGATGGACGATCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGTATCTGGAAAAACGCGGAtt                                                                       >  1:542211/1‑150 (MQ=255)                                                                           aTTTGAAGCAGGCAGCGGCCCCAGCCAGATCGAGCGCCATCAGAGGCAAACGCTTTATCAGTTGATGGACGATCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGTATCTGGAAAAACGCGGATTa                                                                      <  1:286880/151‑1 (MQ=255)                                                                                aaGCAGGCAGCGGCCCCAGCCAGATCGAGCGCCATCAGAGGCAAACGCTTTATCAGTTGATGGACGATCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGTATCTGGAAAAACGCGGATTAAGCCa                                                                 <  2:201314/151‑1 (MQ=255)                                                                                    aggcagCGGCCCCAGCCAGATCGAGCGCCATCAGAGGCAAACGCTTTATCAGTTGATGGACGATCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGTATCTGGAAAAACGCGGATTAAGCCACg                                                               >  2:502017/1‑149 (MQ=255)                                                                                    aggcagCGGCCCCAGCCAGATCGAGCGCCATCAGAGGCAAACGCTTTATCAGTTGATGGACGATCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGTATCTGGAAAAACGCGGATTAAGCCACGAg                                                             <  1:400039/151‑1 (MQ=255)                                                                                             ccccAGCCAGATCGAGCGCCATCAGAGGCAAACGCTTTATCAGTTGATGGACGATCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGTATCTGGAAAAACGCGGATTAAGCCACGAGGTTAtcgct                                                    <  2:565183/151‑1 (MQ=255)                                                                                               ccagccagATCGAGCGCCATCAGAGGCAAACGCTTTATCAGTTGATGGACGATCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGTATCTGGAAAAACGCGGATTAAGCCACGAGGTTAtcgctcg                                                  >  2:165937/1‑151 (MQ=255)                                                                                                 agccagATCGAGCGCCATCAGAGGCAAACGCTTTATCAGTTGATGGACGATCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGTATCTGGAAAAACGCGGATTAAGCCACGAGGTTATcgctcgct                                                >  2:458160/1‑151 (MQ=255)                                                                                                          gAGCGCCATCAGAGGCAAACGCTTTATCAGTTGATGGACGATCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGTATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTg                                       >  2:260236/1‑151 (MQ=255)                                                                                                              gcCATCAGAGGCAAACGCTTTATCAGTTGATGGACGATCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGTATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGttt                                   >  2:406960/1‑151 (MQ=255)                                                                                                               ccATCAGAGGCAAACGCTTTATCAGTTGATGGACGATCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGTATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGtttt                                  <  2:344080/151‑1 (MQ=255)                                                                                                               ccATCAGAGGCAAACGCTTTATCAGTTGATGGACGATCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGTATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGtttt                                  <  1:362965/151‑1 (MQ=255)                                                                                                                   cAGAGGCAAACGCTTTATCAGTTGATGGACGATCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGTATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTgcgc                              <  1:347112/151‑1 (MQ=255)                                                                                                                     gagGCAAACGCTTTATCAGTTGATGGACGATCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGTATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGccc                            >  1:371855/1‑151 (MQ=255)                                                                                                                        gCAAACGCTTTATCAGTTGATGGACGATCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGTATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGcccccc                         <  1:584733/151‑1 (MQ=255)                                                                                                                          aaaCGCTTTATCAGTTGATGGACGATCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGTATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGCCCCCCgg                       >  1:330761/1‑151 (MQ=255)                                                                                                                                      aGTTGATGGACGATCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGTATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGCCCCCCGGCTGGGACAACGt           >  1:355916/1‑151 (MQ=255)                                                                                                                                               aCGATCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGTATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGCCCCCCGGCTGGGACAACGTCCTGAAGCg  <  1:368527/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    CTGATGAACTACGACAAGCTCGAGTTCGTCGAAACGGTCGAAGAGCTGGCAGCAATGCACAATCTTGAAGTGCCATTTGAAGCAGGCAGCGGCCCCAGCCAGATCGAGCGCCATCAGAGGCAAACGCTTTATCAGTTGATGGACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGTATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGCCCCCCGGCTGGGACAACGTCCTGAAGCG  >  NC_000913/3211326‑3211618

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