Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,371,109	G→A	V29I (GTT→ATT)	groL →	chaperonin GroEL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,371,109	0	G	A	100.0%	128.4 / NA	35	V29I (GTT→ATT)	groL	chaperonin GroEL
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/15);  total (20/15)

GCTAAAGACGTAAAATTCGGTAACGACGCTCGTGTGAAAATGCTGCGCGGCGTAAACGTACTGGCAGATGCAGTGAAAGTTACCCTCGGTCCAAAAGGCCGTAACGTAGTTCTGGATAAATCTTTCGGTGCACCGACCATCACCAAAGAT  >  NC_000913/4371031‑4371180                                                                               |                                                                        gCTAAAGACGTAAAATTCGGTAACGACGCTCGTGTGAAAATGCTGCGCGGCGTAAACGTACTGGCAGATGCAGTGAAAAtt                                                                       <  1:878667/81‑1 (MQ=255)        aCGTAAAATTCGGTAACGACGCTCGTGTGAAAATGCTGCGCGGCGTAAACGTACTGGCAGATGCAGTGAAAATTACCCTCg                                                                >  2:878727/1‑81 (MQ=255)          gTAAAATTCGGTAACGACGCTCGTGTGAAAATGCTGCGCGGCGTAAACGTACTGGCAGATGCAGTGAAAATTACCCTCGGt                                                              <  1:878668/81‑1 (MQ=255)           tAAAATTCGGTAACGACGCTCGTGTGAAAATGCTGCGCGGCGTAAACGTACTGGCAGATGCAGTGAAAATTACCCTCGGTc                                                             >  2:878735/1‑81 (MQ=255)            aaaaTTCGGTAACGACGCTCGTGTGAAAATGCTGCGCGGCGTAAACGTACTGGCAGATGCAGTGAAAATTACCCTCGGTc                                                             >  2:878716/1‑80 (MQ=255)               aTTCGGTAACGACGCTCGTGTGAAAATGCTGCGCGGCGTAAACGTACTGGCAGATGCAGTGAAAATTACCCTCGGTCCaaa                                                         <  2:878669/81‑1 (MQ=255)                ttCGGTAACGACGCTCGTGTGAAAATGCTGCGCGGCGTAAACGTACTGGCAGATGCAGTGAAAATTACCCTCGGTCCaa                                                          >  1:878722/1‑79 (MQ=255)                  cGGTAACGACGCTCGTGTGAAAATGCTGCGCGGCGTAAACGTACTGGCAGATGCAGTGAAAATTACCCTCGGTCCAAAAgg                                                      >  2:878720/1‑81 (MQ=255)                  cGGTAACGACGCTCGTGTGAAAATGCTGCGCGGCGTAAACGTACTGGCAGATGCAGTGAAAATTACCCTCGGTCCAAAAgg                                                      >  1:878710/1‑81 (MQ=255)                     tAACGACGCTCGTGTGAAAATGCTGCGCGGCGTAAACGTACTGGCAGATGCAGTGAAAATTACCCTCGGTCCAAAAGGCCg                                                   >  1:878719/1‑81 (MQ=255)                     tAACGACGCTCGTGTGAAAATGCTGCGCGGCGTAAACGTACTGGCAGATGCAGTGAAAATTACCCTCGGTCCAAAAGGCCg                                                   <  1:878670/81‑1 (MQ=255)                      aaCGACGCTCGTGTGAAAATGCTGCGCGGCGTAAACGTACTGGCAGATGCAGTGAAAATTACCCTCGGTCCAAAAGGCCg                                                   >  2:878705/1‑80 (MQ=255)                       acgacgCTCGTGTGAAAATGCTGCGCGGCGTAAACGTACTGGCAGATGCAGTGAAAATTACCCTCGGTCCAAAAGGCCGt                                                  >  2:878730/1‑80 (MQ=255)                       acgacgCTCGTGTGAAAATGCTGCGCGGCGTAAACGTACTGGCAGATGCAGTGAAAATTACCCTCGGTCCAAAAGGCCGTa                                                 <  2:878671/81‑1 (MQ=255)                          acgCTCGTGTGAAAATGCTGCGCGGCGTAAACGTACTGGCAGATGCAGTGAAAATTACCCTCGGTCCAAAAGGCCGTAACg                                              >  2:878731/1‑81 (MQ=255)                               cGTGTGAAAATGCTGCGCGGCGTAAACGTACTGGCAGATGCAGTGAAAATTACCCTCGGTCCAAAAGGCCGTAACGTAGtt                                         >  2:878726/1‑81 (MQ=255)                                gtgtGAAAATGCTGCGCGGCGTAAACGTACTGGCAGATGCAGTGAAAATTACCCTCGGTCCAAAAGGCCGTAACGTAGt                                          <  2:878672/79‑1 (MQ=255)                                gtgtGAAAATGCTGCGCGGCGTAAACGTACTGGCAGATGCAGTGAAAATTACCCTCGGTCCAAAAGGCCGTAACGTAGTTc                                        >  2:878741/1‑81 (MQ=255)                                 tgtgAAAATGCTGCGCGGCGTAAACGTACTGGCAGATGCAGTGAAAATTACCCTCGGTCCAAAAGGCCGTAACGTAGTTCt                                       <  1:878673/81‑1 (MQ=255)                                       aaTGCTGCGCGGCGTAAACGTACTGGCAGATGCAGTGAAAATTACCCTCGGTCCAAAAGGCCGTAACGTAGTTCTGGATaa                                 <  1:878674/81‑1 (MQ=255)                                        aTGCTGCGCGGCGTAAACGTACTGGCAGATGCAGTGAAAATTACCCTCGGTCCAAAAGGCCGTAACGTAGTTCTGGATaaa                                <  1:878675/81‑1 (MQ=255)                                        aTGCTGCGCGGCGTAAACGTACTGGCAGATGCAGTGAAAATTACCCTCGGTCCAAAAGGCCGTAACGTAGTTCTGGATaaa                                <  1:878676/81‑1 (MQ=255)                                              cgcgGCGTAAACGTACTGGCAGATGCAGTGAAAATTACCCTCGGTCCAAAAGGCCGTAACGTAGTTCTGGATAAATCTTTc                          <  1:878677/81‑1 (MQ=255)                                                    gTAAACGTACTGGCAGATGCAGTGAAAATTACCCTCGGTCCAAAAGGCCGTAACGTAGTTCTGGATAAATCTTTCGGTGCa                    >  2:878721/1‑81 (MQ=255)                                                     tAAACGTACTGGCAGATGCAGTGAAAATTACCCTCGGTCCAAAAGGCCGTAACGTAGTTCTGGATAAATCTTTCGGTGCAc                   <  2:878678/81‑1 (MQ=255)                                                      aaaCGTACTGGCAGATGCAGTGAAAATTACCCTCGGTCCAAAAGGCCGTAACGTAGTTCTGGATAAATCTTTCGGTGCAcc                  >  1:878737/1‑81 (MQ=255)                                                      aaaCGTACTGGCAGATGCAGTGAAAATTACCCTCGGTCCAAAAGGCCGTAACGTAGTTCTGGATAAATCTTTCGGTGCAcc                  >  2:878740/1‑81 (MQ=255)                                                       aaCGTACTGGCAGATGCAGTGAAAATTACCCTCGGTCCAAAAGGCCGTAACGTAGTTCTGGATAAATCTTTCGGTGCAcc                  >  1:878729/1‑80 (MQ=255)                                                       aaCGTACTGGCAGATGCAGTGAAAATTACCCTCGGTCCAAAAGGCCGTAACGTAGTTCTGGATAAATCTTTCGGTGCACCg                 <  1:878679/81‑1 (MQ=255)                                                          gTACTGGCAGATGCAGTGAAAATTACCCTCGGTCCAAAAGGCCGTAACGTAGTTCTGGATAAATCTTTCGGTGCACCGAc               >  1:878728/1‑80 (MQ=255)                                                            aCTGGCAGATGCAGTGAAAATTACCCTCGGTCCAAAAGGCCGTAACGTAGTTCTGGATAAATCTTTCGGTGCACCGACCAt            >  1:878736/1‑81 (MQ=255)                                                               ggCAGATGCAGTGAAAATTACCCTCGGTCCAAAAGGCCGTAACGTAGTTCTGGATAAATCTTTCGGTGCACCGACCATCAc         >  1:878734/1‑81 (MQ=255)                                                                gCAGATGCAGTGAAAATTACCCTCGGTCCAAAAGGCCGTAACGTAGTTCTGGATAAATCTTTCGGTGCACCGACCATCAcc        <  2:878681/81‑1 (MQ=255)                                                                gCAGATGCAGTGAAAATTACCCTCGGTCCAAAAGGCCGTAACGTAGTTCTGGATAAATCTTTCGGTGCACCGACCATCAc         <  2:878680/80‑1 (MQ=255)                                                                      gCAGTGAAAATTACCCTCGGTCCAAAAGGCCGTAACGTAGTTCTGGATAAATCTTTCGGTGCACCGACCATCACCAAAGAt  >  2:878732/1‑81 (MQ=255)                                                                               |                                                                        GCTAAAGACGTAAAATTCGGTAACGACGCTCGTGTGAAAATGCTGCGCGGCGTAAACGTACTGGCAGATGCAGTGAAAGTTACCCTCGGTCCAAAAGGCCGTAACGTAGTTCTGGATAAATCTTTCGGTGCACCGACCATCACCAAAGAT  >  NC_000913/4371031‑4371180

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