Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	428,188	G→T	R181L (CGC→CTC)	secD →	SecYEG protein translocase auxillary subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	428,188	0	G	T	100.0%	147.8 / NA	43	R181L (CGC→CTC)	secD	SecYEG protein translocase auxillary subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (20/23);  total (20/23)

CAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCGCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCTGCGTG  >  NC_000913/428091‑428280                                                                                                  |                                                                                             cAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCctcg                                                                                            <  1:1206760/100‑1 (MQ=255) cAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCctcg                                                                                            <  1:723455/100‑1 (MQ=255)     atatCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGc                                                                                        <  1:1133804/100‑1 (MQ=255)     atatCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGc                                                                                        <  2:85754/100‑1 (MQ=255)     atatCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGc                                                                                        <  1:943729/100‑1 (MQ=255)       atCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTa                                                                                      >  2:855966/1‑100 (MQ=255)        tCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTaa                                                                                     <  2:177791/100‑1 (MQ=255)          gATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAg                                                                                   >  1:14825/1‑100 (MQ=255)              gCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCg                                                                               >  1:129708/1‑100 (MQ=255)                  aCGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGACGCTAAAGCTCGtgat                                                                           <  2:862413/100‑1 (MQ=255)                   cGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGtgatg                                                                          >  1:1160164/1‑100 (MQ=255)                     cAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGaa                                                                        >  2:505829/1‑100 (MQ=255)                        tGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGcc                                                                     <  1:505833/100‑1 (MQ=255)                         gACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCa                                                                    >  1:607806/1‑100 (MQ=255)                         gACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCa                                                                    >  2:1241846/1‑100 (MQ=255)                          aCCTGCGCGAAAAGGGCATCCCGTATCCCACTGTTCGTAAAGAAAGCAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGGGATGAAGCCAt                                                                   <  1:437307/100‑1 (MQ=255)                            cTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTg                                                                 <  2:225062/100‑1 (MQ=255)                              gcgcGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCg                                                               <  2:183977/100‑1 (MQ=255)                                  gAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATc                                                           <  1:531676/100‑1 (MQ=255)                                   aaaaGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCt                                                          >  1:965821/1‑100 (MQ=255)                                    aaaGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTg                                                         >  1:80783/1‑100 (MQ=255)                                     aaGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGa                                                        <  1:702954/100‑1 (MQ=255)                                     aaGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGa                                                        >  2:59016/1‑100 (MQ=255)                                      aGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAg                                                       >  1:698491/1‑100 (MQ=255)                                       gggCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGc                                                      <  1:1212984/100‑1 (MQ=255)                                          cATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAg                                                   >  2:35047/1‑100 (MQ=255)                                                   taCCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCg                                          >  1:566076/1‑100 (MQ=255)                                                     cgaCTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGGGATGAAGCCATTGCGTATCTGAGCGAGCGCCa                                              >  2:1142224/3‑94 (MQ=255)                                                        cTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGAcc                                      <  1:1241857/99‑1 (MQ=255)                                                         tGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTg                                    <  1:419349/100‑1 (MQ=255)                                                           ttCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGt                                  <  1:814479/100‑1 (MQ=255)                                                            tCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTg                                 <  1:855972/100‑1 (MQ=255)                                                             cGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGa                                >  1:1000483/1‑100 (MQ=255)                                                                  aGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTagc                           <  1:784346/100‑1 (MQ=255)                                                                   gAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTagca                          <  2:607802/100‑1 (MQ=255)                                                                     aaaCAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTagcagc                        >  1:224523/1‑100 (MQ=255)                                                                             aCGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGgcag                <  2:129706/100‑1 (MQ=255)                                                                             aCGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGgcag                <  1:1154277/100‑1 (MQ=255)                                                                              cGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGgcagc               <  2:190148/100‑1 (MQ=255)                                                                                gCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCaa             >  1:374293/1‑100 (MQ=255)                                                                                atcTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCaa             <  1:319470/98‑1 (MQ=255)                                                                                    gAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCaaagat         >  2:319467/1‑96 (MQ=255)                                                                                      gCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCt       >  1:214225/1‑100 (MQ=255)                                                                                         taaCTTTCCTCGATGCGAAAGCGCGGGAGGAAGCCAGTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCtgcg    <  2:616902/98‑1 (MQ=255)                                                                                           aCTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCtgcgtg  >  1:972504/1‑100 (MQ=255)                                                                                                  |                                                                                             CAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCGCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCTGCGTG  >  NC_000913/428091‑428280

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