Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,607,869	G→A	G23S (GGC→AGC)	holD →	DNA polymerase III subunit psi

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,607,869	0	G	A	100.0%	188.9 / NA	53	G23S (GGC→AGC)	holD	DNA polymerase III subunit psi
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/27);  total (26/27)

CGATTCGCTATGACATCCCGACGAGACTGGCAGTTACAGCAACTGGGCATTACCCAGTGGTCGCTGCGTCGCCCTGGCGCGTTGCAGGGGGAGATTGCCATTGCGATCCCGGCACACGTCCGTCTGGTGATGGTGGCAAACGATCTTCCCGC  >  NC_000913/4607794‑4607945                                                                            |                                                                             cGATTCGCTATGACATCCCGACGAGACTGGCAGTTACAGCAACTGGGCATTACCCAGTGGTCGCTGCGTCGCCCTAGcgcg                                                                         >  2:932857/1‑81 (MQ=255) cGATTCGCTATGACATCCCGACGAGACTGGCAGTTACAGCAACTGGGCATTACCCAGTGGTCGCTGCGTCGCCCTAGcgcg                                                                         <  2:932780/81‑1 (MQ=255)  gATTCGCTATGACATCCCGACGAGACTGGCAGTTACAGCAACTGGGCATTACCCAGTGGTCGCTGCGTCGCCCTAGCGCGt                                                                        >  2:932863/1‑81 (MQ=255)  gATTCGCTATGACATCCCGACGAGACTGGCAGTTACAGCAACTGGGCATTACCCAGTGGTCGCTGCGTCGCCCTAGCGCGt                                                                        <  2:932781/81‑1 (MQ=255)  gATTCGCTATGACATCCCGACGAGACTGGCAGTTACAGCAACTGGGCATTACCCAGTGGTCGCTGCGTCGCCCTAGCGCGt                                                                        >  1:932846/1‑81 (MQ=255)   aTTCGCTATGACATCCCGACGAGACTGGCAGTTACAGCAACTGGGCATTACCCAGTGGTCGCTGCGTCGCCCTAGCGCGtt                                                                       <  2:932782/81‑1 (MQ=255)      cGCTATGACATCCCGACGAGACTGGCAGTTACAGCAACTGGGCATTACCCAGTGGTCGCTGCGTCGCCCTAGCGCGTTGc                                                                     <  2:932784/80‑1 (MQ=255)      cGCTATGACATCCCGACGAGACTGGCAGTTACAGCAACTGGGCATTACCCAGTGGTCGCTGCGTCGCCCTAGCGCGTTGc                                                                     <  2:932783/80‑1 (MQ=255)          aTGACATCCCGACGAGACTGGCAGTTACAGCAACTGGGCATTACCCAGTGGTCGCTGCGTCGCCCTAGCGCGTTGCAgggg                                                                >  2:932886/1‑81 (MQ=255)           tGACATCCCGACGAGACTGGCAGTTACAGCAACTGGGCATTACCCAGTGGTCGCTGCGTCGCCCTAGCGCGTTGCAggggg                                                               <  1:932785/81‑1 (MQ=255)            gACATCCCGACGAGACTGGCAGTTACAGCAACTGGGCATTACCCAGTGGTCGCTGCGTCGCCCTAGCGCGTTGCAgggg                                                                >  1:932859/1‑79 (MQ=255)             aCATCCCGACGAGACTGGCAGTTACAGCAACTGGGCATTACCCAGTGGTCGCTGCGTCGCCCTAGCGCGTTGCAgggg                                                                >  1:932870/1‑78 (MQ=255)              cATCCCGACGAGACTGGCAGTTACAGCAACTGGGCATTACCCAGTGGTCGCTGCGTCGCCCTAGCGCGTTGCAgggg                                                                >  1:932856/1‑77 (MQ=255)               aTCCCGACGAGACTGGCAGTTACAGCAACTGGGCATTACCCAGTGGTCGCTGCGTCGCCCTAGCGCGTTGCAGGGGGAGAt                                                           >  2:932842/1‑81 (MQ=255)               aTCCCGACGAGACTGGCAGTTACAGCAACTGGGCATTACCCAGTGGTCGCTGCGTCGCCCTAGCGCGTTGCAGGGGGAGAt                                                           <  2:932786/81‑1 (MQ=255)               aTCCCGACGAGACTGGCAGTTACAGCAACTGGGCATTACCCAGTGGTCGCTGCGTCGCCCTAGCGCGTTGCAGGGGGAGAt                                                           >  1:932866/1‑81 (MQ=255)                tCCCGACGAGACTGGCAGTTACAGCAACTGGGCATTACCCAGTGGTCGCTGCGTCGCCCTAGCGCGTTGCAGGGGGAGAtt                                                          <  2:932787/81‑1 (MQ=255)                   cgacgaGACTGGCAGTTACAGCAACTGGGCATTACCCAGTGGTCGCTGCGTCGCCCTAGCGCGTTGCAGGGGGAGATTGc                                                        >  1:932867/1‑80 (MQ=255)                   cgacgaGACTGGCAGTTACAGCAACTCGGCATTACCCAGTGGTCGCTGCGTCGCCCTAGCGCGTTGCAGGGGGAGATTGc                                                        >  1:932845/1‑80 (MQ=255)                    gacgaGACTGGCAGTTACAGCAACTGGGCATTACCCAGTGGTCGCTGCGTCGCCCTAGCGCGTTGCAGGGGGAGATTGCCa                                                      >  1:932858/1‑81 (MQ=255)                       gagaCTGGCAGTTACAGCAACTGGGCATTACCCAGTGGTCGCTGCGTCGCCCTAGCGCGTTGCAGGGGGAGATTGCCATTg                                                   >  1:932852/1‑81 (MQ=255)                       gagaCTGGCAGTTACAGCAACTGGGCATTACCCAGTGGTCGCTGCGTCGCCCTAGCGCGTTGCAGGGGGAGATTGCCATTg                                                   >  1:932843/1‑81 (MQ=255)                       gagaCTGGCAGTTACAGCAACTGGGCATTACCCAGTGGTCGCTGCGTCGCCCTAGCGCGTTGCAGGGGGAGATTGCCATTg                                                   >  1:932849/1‑81 (MQ=255)                       gagaCTGGCAGTTACAGCAACTGGGCATTACCCAGTGGTCGCTGCGTCGCCCTAGCGCGTTGCAGGGGGAGATTGCCATTg                                                   >  1:932865/1‑81 (MQ=255)                        agaCTGGCAGTTACAGCAACTGGGCATTACCCAGTGGTCGCTGCGTCGCCCTAGCGCGTTGCAGGGGGAGATTGCCATTGc                                                  <  1:932788/81‑1 (MQ=255)                              gCAGTTACAGCAACTGGGCATTACCCAGTGGTCGCTGCGTCGCCCTAGCGCGTTGCAGGGGGAGATTGCCATTGCGATccc                                            <  1:932790/81‑1 (MQ=255)                              gCAGTTACAGCAACTGGGCATTACCCAGTGGTCGCTGCGTCGCCCTAGCGCGTTGCAGGGGGAGATTGCCATTGCGATccc                                            >  2:932873/1‑81 (MQ=255)                              gCAGTTACAGCAACTGGGCATTACCCAGTGGTCGCTGCGTCGCCCTAGCGCGTTGCAGGGGGAGATTGCCATTGCGATcc                                             <  2:932789/80‑1 (MQ=255)                               cAGTTACAGCAACTGGGCATTACCCAGTGGTCGCTGCGTCGCCCTAGCGCGTTGCAGGGGGAGATTGCCATTGCGATCCCg                                           <  2:932791/81‑1 (MQ=255)                                aGTTACAGCAACTGGGCATTACCCAGTGGTCACTGCGTCGCCCTAGCGCGTTGCAGGGGGAGATTGCCATTGCGATCCCgg                                          <  2:932792/81‑1 (MQ=255)                                 gTTACAGCAACTGGGCATTACCCAGTGGTCGCTGCGTCGCCCTAGCGCGTTGCAGGGGGAGATTGCCATTGCGATCCCGGc                                         <  2:932793/81‑1 (MQ=255)                                    aCAGCAACTGGGCATTACCCAGTGGTCGCTGCGTCGCCCTAGCGCGTTGCAGGGGGAGATTGCCATTGCGATCCCGGcac                                       <  2:932794/80‑1 (MQ=255)                                              ggCATTACCCAGTGGTCGCTGCGTCGCCCTAGCGCGTTGCAGGGGGAGATTGCCATTGCGATCCCGGCACACGTCCGTCTg                            <  1:932795/81‑1 (MQ=255)                                                cATTACCCAGTGGTCGCTGCGTCGCCCTAGCGCGTTGCAGGGGGAGATTGCCATTGCGATCCCGGCACACGTCCGTCTGGt                          >  2:932883/1‑81 (MQ=255)                                                 aTTACCCAGTGGTCGCTGCGTCGCCCTAGCGCGTTGCAGGGGGAGATTGCCATTGCGATCCCGGCACACGTCCGTCTGGt                          >  1:932861/1‑80 (MQ=255)                                                   tACCCAGTGGTCGCTGCGTCGCCCTAGCGCGTTGCAGGGGGAGATTGCCATTGCGATCCCGGCACACGTCCGTCTGGTGAt                       >  2:932869/1‑81 (MQ=255)                                                     cccAGTGGTCGCTGCGTCGCCCTAGCGCGTTGCAGGGGGAGATTGCCATTGCGATCCCGGCACACGTCCGTCTGGTGAtg                      <  2:932796/80‑1 (MQ=255)                                                      ccAGTGGTCGCTGCGTCGCCCTAGCGCGTTGCAGGGGGAGATTGCCATTGCGATCCCGGCACACGTCCGTCTGGTGAtggt                    >  1:932864/1‑81 (MQ=255)                                                           ggTCGCTGCGTCGCCCTAGCGCGTTGCAGGGGGAGATTGCCATTGCGATCCCGGCACACGTCCGTCTGGTGATGGTGGCa                >  2:932871/1‑80 (MQ=255)                                                            gTCGCTGCGTCGCCCTAGCGCGTTGCAGGGGGAGATTGCCATTGCGATCCCGGCACACGTCCGTCTGGTGATGGTGGCaaa              <  2:932798/81‑1 (MQ=255)                                                            gTCGCTGCGTCGCCCTAGCGCGTTGCAGGGGGAGATTGCCATTGCGATCCCGGCACACGTCCGTCTGGTGATGGTGGCaa               <  2:932797/80‑1 (MQ=255)                                                             tCGCTGCGTCGCCCTAGCGCGTTGCAGGGGGAGATTGCCATTGCGATCCCGGCACACGTCCGTCTGGTGATGGTGGCAAAc             >  2:932874/1‑81 (MQ=255)                                                               gCTGCGTCGCCCTAGCGCGTTGCAGGGGGAGATTGCCATTGCGATCCCGGCACACGTCCGTCTGGTGATGGTGGCAAACg            <  1:932800/80‑1 (MQ=255)                                                               gCTGCGTCGCCCTAGCGCGTTGCAGGGGGAGATTGCCATTGCGATCCCGGCACACGTCCGTCTGGTGATGGTGGCAAACGa           <  1:932799/81‑1 (MQ=255)                                                               gCTGCGTCGCCCTAGCGCGTTGCAGGGGGAGATTGCCATTGCGATCCCGGCACACGTCCGTCTGGTGATGGTGGCAAACGa           >  2:932862/1‑81 (MQ=255)                                                                  gCGTCGCCCTAGCGCGTTGCAGGGGGAGATTGCCATTGCGATCCCGGCACACGTCCGTCTGGTGATGGTGGCAAACGATCt        >  1:932885/1‑81 (MQ=255)                                                                   cGTCGCCCTAGCGCGTTGCAGGGGGAGATTGCCATTGTGATCCCGGCACACGTCCGTCTGGTGATGGTGGCAAACGATCt        <  2:932801/80‑1 (MQ=255)                                                                    gTCGCCCTAGCGCGTTGCAGGGGGAGATTGCCATTGCGATCCCGGCACACGTCCGTCTGGTGATGGTGGCAAACGATCTTc      <  2:932802/81‑1 (MQ=255)                                                                    gTCGCCCTAGCGCGTTGCAGGGGGAGATTGCCATTGCGATCCCGGCACACGTCCGTCTGGTGATGGTGGCAAACGATCTTc      >  1:932876/1‑81 (MQ=255)                                                                       gCCCTAGCGCGTTGCAGGGGGAGATTGCCATTGCGATCCCGGCACACGTCCGTCTGGTGATGGTGGCAAACGATCTTCCCg   <  2:932805/81‑1 (MQ=255)                                                                       gCCCTAGCGCGTTGCAGGGGGAGATTGCCATTGCGATCCCGGCACACGTCCGTCTGGTGATGGTGGCAAACGATCTTCCCg   <  2:932804/81‑1 (MQ=255)                                                                       gCCCTAGCGCGTTGCAGGGGGAGATTGCCATTGCGATCCCGGCACACGTCCGTCTGGTGATGGTGGCAAACGATCTTCCCg   <  2:932803/81‑1 (MQ=255)                                                                        cccTAGCGCGTTGCAGGGGGAGATTGCCATTGCGATCCCGGCACACGTCCGTCTGGTGATGGTGGCAAACGATCTTCCCGc  <  2:932806/81‑1 (MQ=255)                                                                            |                                                                             CGATTCGCTATGACATCCCGACGAGACTGGCAGTTACAGCAACTGGGCATTACCCAGTGGTCGCTGCGTCGCCCTGGCGCGTTGCAGGGGGAGATTGCCATTGCGATCCCGGCACACGTCCGTCTGGTGATGGTGGCAAACGATCTTCCCGC  >  NC_000913/4607794‑4607945

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