Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	13,901	A→G	D580G (GAC→GGC)	dnaK →	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	13,901	0	A	G	100.0%	110.5 / NA	32	D580G (GAC→GGC)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/16);  total (16/16)

CGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAA  >  NC_000913/13824‑13977                                                                              |                                                                             cGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGGCa                                                                            >  2:3306/1‑80 (MQ=255)  gACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGGCa                                                                            >  2:3311/1‑79 (MQ=255)  gACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGGCa                                                                            >  2:3292/1‑79 (MQ=255)     aaaCTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGGCAAAGcc                                                                       <  1:3254/81‑1 (MQ=255)       aCTGCCGGCTGACGACAAAACTGCTATCTAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGGCAAAGCCGc                                                                     >  1:3334/1‑81 (MQ=255)        cTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGGCAAAGCCGCt                                                                    <  1:3256/81‑1 (MQ=255)        cTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGGCAAAGCCGCt                                                                    <  2:3255/81‑1 (MQ=255)           ccGGCTGACGACAAAACTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGGCAAAGCCGCTAt                                                                  <  1:3257/80‑1 (MQ=255)               cTGACGACAAAACTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGGCAAAGCCGCTATCGAAg                                                             >  2:3328/1‑81 (MQ=255)               cTGACGACAAAACTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGGCAAAGCCGCTATCGAAg                                                             >  1:3318/1‑81 (MQ=255)                     acAAAACTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGGCAAAGCCGCTATCGAAGCGaaaa                                                       >  1:3315/1‑81 (MQ=255)                      cAAAACTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGGCAAAGCCGCTATCGAAGCGAAAAt                                                      <  1:3258/81‑1 (MQ=255)                       aaaaCTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGGCAAAGCCGCTATCGAAGCGAAAATg                                                     >  1:3309/1‑81 (MQ=255)                       aaaaCTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGGCAAAGCCGCTATCGAAGCGAAAATg                                                     <  1:3259/81‑1 (MQ=255)                        aaaCTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGGCAAAGCCGCTATCGAAGCGAAAATg                                                     >  1:3326/1‑80 (MQ=255)                        aaaCTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGGCAAAGCCGCTATCGAAGCGAAAATg                                                     >  1:3316/1‑80 (MQ=255)                        aaaCTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGGCAAAGCCGCTATCGAAGCGAAAATGc                                                    <  2:3260/81‑1 (MQ=255)                         aaCTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGGCAAAGCCGCTATCGAAGCGAAAATGCa                                                   <  1:3261/81‑1 (MQ=255)                          aCTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGGCAAAGCCGCTATCGAAGCGAAAATGCAg                                                  <  1:3262/81‑1 (MQ=255)                              cTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGGCAAAGCCGCTATCGAAGCGAAAATGCAGGaa                                               <  2:3263/80‑1 (MQ=255)                                   gAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGGCAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGca                                         >  1:3317/1‑81 (MQ=255)                                   gAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGGCAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGca                                         >  2:3327/1‑81 (MQ=255)                                    aGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGGCAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGcac                                        <  2:3264/81‑1 (MQ=255)                                     gTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGGCAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGcaca                                       <  1:3265/81‑1 (MQ=255)                                               aCTGCACTGGAAACTGCTCTGAAAGGTGAAGGCAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAg                             <  1:3266/81‑1 (MQ=255)                                                cTGCACTGGAAACTGCTCTGAAAGGTGAAGGCAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAg                             >  1:3336/1‑80 (MQ=255)                                                   cACTGGAAACTGCTCTGAAAGGTGAAGGCAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGaa                           <  2:3267/79‑1 (MQ=255)                                                          aaCTGCTCTGAAAGGTGAAGGCAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGa                  <  2:3268/81‑1 (MQ=255)                                                               ctctGAAAGGTGAAGGCAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATc              <  1:3269/80‑1 (MQ=255)                                                                  tGAAAGGTGAAGGCAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGcc           >  1:3305/1‑80 (MQ=255)                                                                          gAAGGCAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCaa  >  2:3325/1‑81 (MQ=255)                                                                           aaGGCAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCaa  >  2:3322/1‑80 (MQ=255)                                                                              |                                                                             CGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAA  >  NC_000913/13824‑13977

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