Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	181,536	C→G	A218G (GCC→GGC)	degP →	periplasmic serine endoprotease DegP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	181,536	0	C	G	100.0%	192.3 / NA	53	A218G (GCC→GGC)	degP	periplasmic serine endoprotease DegP
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (23/30);  total (23/30)

GATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGCCGAAAACTACGAAAACTTCATCCAGACCGATGCAGCGATCAACCGTGGTAACTCCGGTGGTGCGCTGGTTAACCTGA  >  NC_000913/181459‑181613                                                                              |                                                                              gATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCg                                                                             >  1:529657/1‑80 (MQ=255)  aTTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCGaa                                                                           <  1:38975/81‑1 (MQ=255)   ttGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCGaaa                                                                          <  1:38977/81‑1 (MQ=255)   ttGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGAGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATAGCGaaa                                                                          <  2:38976/81‑1 (MQ=255)    tGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCGaaaa                                                                         >  2:39037/1‑81 (MQ=255)    tGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCGaaaa                                                                         >  2:39020/1‑81 (MQ=255)      gTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCGAAAAc                                                                        <  2:38980/80‑1 (MQ=255)      gTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCGAAAACt                                                                       <  2:38981/81‑1 (MQ=255)      gTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCGAAAACt                                                                       <  1:38978/81‑1 (MQ=255)      gTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGACTGAATGGCGAAAAc                                                                        <  2:38979/80‑1 (MQ=255)            cGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCGAAAACTACGaaa                                                                 <  2:38983/81‑1 (MQ=255)            cGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCGAAAACTACGaaa                                                                 <  2:38982/81‑1 (MQ=255)               ttGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCGAAAACTACGAAAACt                                                              >  2:39028/1‑81 (MQ=255)                tGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCGAAAACTACGAAAACtt                                                             >  2:39045/1‑81 (MQ=255)                tGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCGAAAACTACGAAAACt                                                              <  2:38984/80‑1 (MQ=255)                 ggTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCGAAAACTACGAAAACtt                                                             >  1:39043/1‑80 (MQ=255)                 ggTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCGAAAACTACGAAAACtt                                                             >  1:39044/1‑80 (MQ=255)                 ggTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCGAAAACTACGAAAACtt                                                             <  2:38985/80‑1 (MQ=255)                     tGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCGAAAACTACGAAAACTTCATc                                                         >  1:39027/1‑80 (MQ=255)                      gggCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCGAAAACTACGAAAACTTCATc                                                         >  1:39047/1‑79 (MQ=255)                       ggCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCGAAAACTACGAAAACTTCATcc                                                        <  2:38986/79‑1 (MQ=255)                        gCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCGAAAACTACGAAAACTTCATCCAGa                                                     <  2:38987/81‑1 (MQ=255)                         cGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCGAAAACTACGAAAACTTCATCCAGAc                                                    <  1:38989/81‑1 (MQ=255)                         cGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCGAAAACTACGAAAACTTCATCCAGAc                                                    <  2:38988/81‑1 (MQ=255)                          gagaCGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCGAAAACTACGAAAACTTCATCCAGAcc                                                   >  2:39046/1‑81 (MQ=255)                          gagaCGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCGAAAACTACGAAAACTTCATCCAGAcc                                                   >  1:39033/1‑81 (MQ=255)                          gagaCGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCGAAAACTACGAAAACTTCATCCAGAcc                                                   >  2:39056/1‑81 (MQ=255)                            gaCGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCGAAAACTACGAAAACTTCATCCAGACCGa                                                 <  2:38990/81‑1 (MQ=255)                             aCGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCGAAAACTACGAAAACTTCATCCAGACCGAt                                                <  1:38991/81‑1 (MQ=255)                               ggTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCGAAAACTACGAAAACTTCATCCAGACCGATGc                                              <  2:38992/81‑1 (MQ=255)                                gTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCGAAAACTACGAAAACTTCATCCAGACCGATg                                               >  1:39038/1‑79 (MQ=255)                                    cTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCGAAAACTACGAAAACTTCATCCAGACCGATGCAGCGa                                         >  1:39039/1‑81 (MQ=255)                                    cTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCGAAAACTACGAAAACTTCATCCAGACCGATGCAGCGa                                         <  1:38993/81‑1 (MQ=255)                                       ccGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCGAAAACTACGAAAACTTCATCCAGACCGATGCAGCGAt                                        <  2:38994/79‑1 (MQ=255)                                       ccGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCGAAAACTACGAAAACTTCATCCAGACCGATGCAGCGAt                                        <  2:38995/79‑1 (MQ=255)                                       ccGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCGAAAACTACGAAAACTTCATCCAGACCGATGCAGCGAt                                        <  2:38996/79‑1 (MQ=255)                                          ggATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGGCGAAAACTACGAAAACTTCATCCAGACCGATGCAGCGATCAAcc                                   <  1:38997/81‑1 (MQ=255)                                                   ctGCGCTGGGGCGTAGCGGCCTGAATGGCGAAAACTACGAAAACTTCATCCAGACCGATGCAGCGATCAACCGTGGTAACt                          >  2:39042/1‑81 (MQ=255)                                                    tGCGCTGGGGCGTAGCGGCCTGAATGGCGAAAACTACGAAAACTTCATCCAGACCGATGCAGCGATCAACCGTGGTaa                            <  1:38998/78‑1 (MQ=255)                                                      cgcTGGGGCGTAGCGGCCTGAATGGCGAAAACTACGAAAACTTCATCCAGACCGATGCAGCGATCAACCGTGGTAACTc                         <  2:38999/79‑1 (MQ=255)                                                      cgcTGGGGCGTAGCGGCCTGAATGGCGAAAACTACGAAAACTTCATCCAGACCGATGCAGCGATCAACCGTGGTAACTCCg                       >  1:39040/1‑81 (MQ=255)                                                      cgcTGGGGCGTAGCGGCCTGAATGGCGAAAACTACGAAAACTTCATCCAGACCGATGCAGCGATCAACCGTGGTAACTCCg                       <  1:39000/81‑1 (MQ=255)                                                       gcTGGGGCGTAGCGGCCTGAATGGCGAAAACTACGAAAACTTCATCCAGACCGATGCAGCGATCAACCGTGGTAACTCCgg                      >  2:39035/1‑81 (MQ=255)                                                          ggggCGTAGCGGCCTGAATGGCGAAAACTACGAAAACTTCATCCAGACCGATGCAGCGATCAACCGTGGTAACTCCggtgg                   >  2:39036/1‑81 (MQ=255)                                                           gggCGTAGCGGCCTGAATGGCGAAAACTACGAAAACTTCATCCAGACCGATGCAGCGATCAACCGTGGTAACTCCggtggt                  <  2:39001/81‑1 (MQ=255)                                                                tAGCGGCCTGAATGGCGAAAACTACGAAAACTTCATCCAGACCGATGCAGCGATCAACCGTGGTAACTCCGGTGGTGCGCt             <  2:39002/81‑1 (MQ=255)                                                                tAGCGGCCTGAATGGCGAAAACTACGAAAACTTCATCCAGACCGATGCAGCGATCAACCGTGGTAACTCCGGTGGTGCGCt             <  2:39003/81‑1 (MQ=255)                                                                    ggCCTGAATGGCGAAAACTACGAAAACTTCATCCAGACCGATGCAGCGATCAACCGTGGTAACTCCGGTGGTGCGCTGGt          >  2:39061/1‑80 (MQ=255)                                                                    ggCCTGAATGGCGAAAACTACGAAAACTTCATCCAGACCGATGCAGCGATCAACCGTGGTAACTCCGGTGGTGCGCTGGt          >  1:39041/1‑80 (MQ=255)                                                                         gAATGGCGAAAACTACGAAAACTTCATCCAGACCGATGCAGCGATCAACCGTGGTAACTCCGGTGGTGCGCTGGTTAACCt    >  1:39062/1‑81 (MQ=255)                                                                         gAATGGCGAAAACTACGAAAACTTCATCCAGACCGATGCAGCGATCAACCGTGGTAACTCCGGTGGTGCGCTGGTTAACCt    >  1:39048/1‑81 (MQ=255)                                                                          aaTGGCGAAAACTACGAAAACTTCATCCAGACCGATGCAGCGATCAACCGTGGTAACTCCGGTGGTGCGCTGGTTAACCt    >  1:39063/1‑80 (MQ=255)                                                                           aTGGCGAAAACTACGAAAACTTCATCCAGACCGATGCAGCGATCATCCGTGGTAACTCCGGTGGTGCGCTGGTTAACCTGa  <  1:39004/81‑1 (MQ=255)                                                                              |                                                                              GATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGCCGAAAACTACGAAAACTTCATCCAGACCGATGCAGCGATCAACCGTGGTAACTCCGGTGGTGCGCTGGTTAACCTGA  >  NC_000913/181459‑181613

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