Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,162,308	G→A	S121L (TCG→TTG)	ftsP ←	septal ring component that protects the divisome from stress, multicopy suppressor of ftsI(Ts)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,162,308	0	G	A	100.0%	113.1 / NA	36	S121L (TCG→TTG)	ftsP	septal ring component that protects the divisome from stress, multicopy suppressor of ftsI(Ts)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/20);  total (16/20)

TCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAATAAGCTTAACATCGTCGCCCTTC  >  NC_000913/3162173‑3162424                                                                                                                                        |                                                                                                                     tCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCAAca                                                                                                                   <  2:263044/139‑1 (MQ=255)   aTCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCAAcatc                                                                                                                 >  1:621566/1‑139 (MQ=255)                          cAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCAACATCATCCGTGCCGGACCGCCCATCAg                                                                                          <  1:158802/139‑1 (MQ=255)                          cAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCAACATCATCCGTGCCGGACCGCCCATCAg                                                                                          <  1:317269/139‑1 (MQ=255)                             ggCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCAACATCATCCGTGCCGGACCGCCCATCAGCgg                                                                                       >  2:688740/1‑139 (MQ=255)                               ccGTTATAGACCTGCTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCAACATCATCCGTGCCGGACCGCCCATCAGCGGGc                                                                                     <  1:688740/139‑1 (MQ=255)                                             ctgAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCAACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGc                                                                       >  2:565053/1‑139 (MQ=255)                                                   cGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCAACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGcc                                                                 <  1:134147/139‑1 (MQ=255)                                                      gCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCAACATCATCCGTGTCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCAcc                                                              <  2:87261/139‑1 (MQ=255)                                                            gggAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCAACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCAtt                                                        >  1:790414/1‑139 (MQ=255)                                                                gTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCAACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTgaga                                                    <  1:682686/139‑1 (MQ=255)                                                                    tGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCAACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGcc                                                                      <  1:424872/117‑1 (MQ=255)                                                                    tGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCAACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGcc                                                                      >  2:424872/1‑117 (MQ=255)                                                                           aTACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCAACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTc                                         >  1:611457/1‑139 (MQ=255)                                                                                 caGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCAACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGt                                                            <  2:780148/114‑1 (MQ=255)                                                                                 caGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCAACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGt                                                            >  1:780148/1‑114 (MQ=255)                                                                                 caGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCAACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCgg                                   <  2:738662/139‑1 (MQ=255)                                                                                    agTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCAACATCATCCGTGCCGGAc                                                                                                    >  1:58311/1‑71 (MQ=255)                                                                                    agTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCAACATCATCCGTGCCGGAc                                                                                                    <  2:58311/71‑1 (MQ=255)                                                                                      tAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCAACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCt                              >  1:583955/1‑139 (MQ=255)                                                                                      tAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCAACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCt                              <  2:611457/139‑1 (MQ=255)                                                                                      tAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCAACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCt                              >  1:147035/1‑139 (MQ=255)                                                                                        gCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCAACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGt                            >  2:778798/1‑139 (MQ=255)                                                                                        gCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCAACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGt                            >  1:156600/1‑139 (MQ=255)                                                                                          tGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCAACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTaa                          <  2:621566/139‑1 (MQ=255)                                                                                            cGTTTTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCAACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACTTGTAGCCCGGCTACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAAt                        <  1:737722/139‑1 (MQ=255)                                                                                                 tGGCGAATGGGCAGCACGGGTGCCCAGTCAGCGTTTGGCAACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAATAAGCt                   >  2:592984/1‑139 (MQ=255)                                                                                                  ggCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCAACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAATAAGCtt                  <  2:147035/139‑1 (MQ=255)                                                                                                   gCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCAACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAATAAGCTTa                 <  2:156600/139‑1 (MQ=255)                                                                                                           gCAGTACGGGTGCCCAGTCAGCGTTTGGCAACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAATAAGCTTAACAtcgtc         >  2:99719/1‑139 (MQ=255)                                                                                                            cAGTACGGGTGCCCAGTCAGCGTTTGGCAACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAATAAGCTTAACAtcgtcg        <  1:565053/139‑1 (MQ=255)                                                                                                                  gggTGCCCAGTCAGCGTTTGGCAACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAATAAGCTTAACATCGTCGCCCTTc  <  1:778798/139‑1 (MQ=255)                                                                                                                    gTGCCCAGTCAGCGTTTGGCAACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAAt                        <  1:801244/115‑1 (MQ=255)                                                                                                                    gTGCCCAGTCAGCGTTTGGCAACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAAt                        >  2:801244/1‑115 (MQ=255)                                                                                                                       cccAGTCAGCGTTTGGCAACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGt                                  >  2:565275/1‑102 (MQ=255)                                                                                                                       cccAGTCAGCGTTTGGCAACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGt                                  <  1:565275/102‑1 (MQ=255)                                                                                                                                        |                                                                                                                     TCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAATAAGCTTAACATCGTCGCCCTTC  >  NC_000913/3162173‑3162424

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