Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,958,309	T→C	A192A (GCA→GCG)	ptrA ←	protease III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,958,309	0	T	C	100.0%	172.2 / NA	48	A192A (GCA→GCG)	ptrA	protease III
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/24);  total (24/24)

AATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTGCTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGTGCCGGGTTAATGGTTTCTGCGCTGACCTGTGCCATGCGCATCCCGTCACGCGTACGCGCCATGGTTAATTCAGCGTTCACCGCATTACGCTCACGTTCGGCATATTTCTTGTCGAGCAAAGGT  >  NC_000913/2958180‑2958432                                                                                                                                  |                                                                                                                            aaTAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTGCTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGCGCCGGGTTa                                                                                                                    >  1:746501/1‑139 (MQ=255)    aaCCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTGCTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATg                                                                                                                 >  1:343626/1‑139 (MQ=255)      ccGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTGCTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGt                                                                                                               >  2:437906/1‑139 (MQ=255)      ccGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTGCTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGt                                                                                                               >  2:422698/1‑139 (MQ=255)       cGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTGCTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGtt                                                                                                              >  2:441346/1‑139 (MQ=255)       cGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTGCTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGtt                                                                                                              >  1:318301/1‑139 (MQ=255)         ccTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTGCTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTgcgc                                                                                                                           >  1:56487/1‑124 (MQ=255)             catcaAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTGCTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTgcg                                                                                                        <  2:1350806/139‑1 (MQ=255)             catcaAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTGCTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTgcg                                                                                                        <  1:778687/139‑1 (MQ=255)             catcaAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTGCTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTgcg                                                                                                        <  2:82959/139‑1 (MQ=255)              atcaAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTGCTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTgcgc                                                                                                       >  1:193659/1‑139 (MQ=255)                       gCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTGCTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTGCGCTGACCtgtg                                                                                              <  2:664372/139‑1 (MQ=255)                                tACTTCTCGTGGAAATCTTTCAGCGCCTGCTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTGCGCTGACCTGTGCCATGCGCa                                                                                     <  2:1371172/139‑1 (MQ=255)                                  cTTCTCGTGGAAATCTTTCAGCGCCTGCTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTgcg                                                                                                        <  1:892662/118‑1 (MQ=255)                                  cTTCTCGTGGAAATCTTTCAGCGCCTGCTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTgcg                                                                                                        >  2:892662/1‑118 (MQ=255)                                     ctcGTGGAAATCTTTCAGCGCCTGCTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTGCGCt                                                                                                      >  1:234695/1‑117 (MQ=255)                                     ctcGTGGAAATCTTTCAGCGCCTGCTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTGCGCt                                                                                                      <  2:234695/117‑1 (MQ=255)                                     ctcGTGGAAATCTTTCAGCGCCTGCTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTGCGCTGACCTGTGCCATGCGCATCCCg                                                                                >  1:426035/1‑139 (MQ=255)                                         tGGAAATCTTTCAGCGCCTGCTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTGCGCTGACCTGTGCCATGCGCATCCCGTCAc                                                                            <  2:343626/139‑1 (MQ=255)                                         tGGAAATCTTTCAGCGCCTGCTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTGCGCTGACCTGTGCCATGCGCATCCCGTCAc                                                                            >  1:275141/1‑139 (MQ=255)                                         tGGAAATCTTTCAGCGCCTGCTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTGCGCTGACCTGTGCCATGCGCATCCCGTCAc                                                                            <  2:426035/139‑1 (MQ=255)                                              aTCTTTCAGCGCCTGCTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTGCGCTGACCTGTGCCATGCGCATCCCGTCACGCGTa                                                                       <  2:275141/139‑1 (MQ=255)                                              aTCTTTCAGCGCCTGCTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTGCGCTGACCTGTGCCATGCGCATCCCGTCACGCGTa                                                                       <  2:193659/139‑1 (MQ=255)                                                  ttCAGCGCCTGCTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTGCGCTGACCTGTGCCATGCGCATCCCGTCACGCGTAcgcg                                                                   <  2:746501/139‑1 (MQ=255)                                                    cAGCGCCTGCTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTGCGCTGACCTGTGCCATGCGCATCCCGTCACGCGTACGCGcc                                                                 <  1:422698/139‑1 (MQ=255)                                                      gcgcCTGCTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTGCGCTGACCTGTGCCATGCGCATCCCGTCACGCGTACGCGCCAt                                                               <  1:441346/139‑1 (MQ=255)                                                         ccTGCTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACTAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTGCGCTGACCTGTGCCATGCGCATCCCGTCACGCGTACGCGcc                                                                 >  2:930581/1‑134 (MQ=255)                                                         ccTGCTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAGACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTGCGCTGACCTGTGCCATGCGCATCCCGTCACGCGTACGCGcc                                                                 <  1:930581/134‑1 (MQ=255)                                                                               tttGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTGCGCTGACCTGTGCCATGCGCATCCCGTCACGCGTACGCGCCATGGTTAATTCAGCGTTCACCGc                                          >  2:998864/1‑135 (MQ=255)                                                                               tttGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTGCGCTGACCTGTGCCATGCGCATCCCGTCACGCGTACGCGCCATGGTTAATTCAGCGTTCACCGCATTa                                      >  2:1487242/1‑139 (MQ=255)                                                                               tttGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGCTGCGCCGGGTTAATGGTTTCTGCGCTGACCTGTGCCATGCGCATCCCGTCACGCGTACGCGCCATGGTTAATTCAGCGTTCACCGc                                          <  1:998864/135‑1 (MQ=255)                                                                                 tGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTGCGCTGACCTGTGCCATGCGCATCCCGTCACGCGTACGCGCCATGGTTAATTCAGCGTTCACCGCATTACg                                    <  2:318301/139‑1 (MQ=255)                                                                                  gTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTGCGCTGACCTGTGCCATGCGCATCCCGTCACGCGTa                                                                       <  2:789817/103‑1 (MQ=255)                                                                                  gTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTGCGCTGACCTGTGCCATGCGCATCCCGTCACGCGTa                                                                       >  1:789817/1‑103 (MQ=255)                                                                                  gTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTGCGCTGACCTGTGCCATGCGCATCCCGTCACGCGTACGCGCCATGGTTAATTCAGCGTTCACCGCATTACGc                                   >  2:808872/1‑139 (MQ=255)                                                                                   tCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTGCGCTGACCTGTGCCATGCGCATCCCGTCACGCGTACGCGCCATGGTTAATTCAGCGTTCACCGCATTACGCt                                  >  1:43198/1‑139 (MQ=255)                                                                                             tttCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTGCGCTGACCTGTGCCATGCGCATCCCGTCACGCGTACGCGCCATGGTTAATTCAGCGTTCACCGCATTACGCTCACGTTCGGc                        <  1:437906/139‑1 (MQ=255)                                                                                                 gAGGTTACCACCAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTGCGCTGACCTGTGCCATGCGCATCCCGTCACGCGTACGCGCCATGGTTAATTCAGCGTTCACCGCATTACGCTCACGTTCGGCatat                    >  2:87629/1‑139 (MQ=255)                                                                                                       accaccAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTGCGCTGACCTGTGCCATGCGCATCCCGTCACGCGTAc                                                                      <  2:1433063/83‑1 (MQ=255)                                                                                                       accaccAGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTGCGCTGACCTGTGCCATGCGCATCCCGTCACGCGTAc                                                                      >  1:1433063/1‑83 (MQ=255)                                                                                                             aGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTgcg                                                                                                        <  2:1367076/43‑1 (MQ=255)                                                                                                             aGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTgcg                                                                                                        >  1:1367076/1‑43 (MQ=255)                                                                                                             aGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTGCGCTGACCTGTGCCATGCGCATCCCGTCACGCGTACGCGCCATGGTTAATTCAGCGTTCACCGCAtt                                       >  2:1042650/1‑108 (MQ=255)                                                                                                             aGAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTGCGCTGACCTGTGCCATGCGCATCCCGTCACGCGTACGCGCCATGGTTAATTCAGCGTTCACCGCAtt                                       <  1:1042650/108‑1 (MQ=255)                                                                                                              gAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTGCGCTGACCTGTGCCATGCGCATCCCGTCACGCGTACGCGCCATGGTTAATTCAGCGTTCACCGCATTACGCTCACGTTCGGCATATTTCTTGTCGAGCa       <  1:371414/139‑1 (MQ=255)                                                                                                              gAAAACTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTGCGCTGACCTGTGCCATGCGCATCCCGTCACGCGTACGCGCCATGGTTAATTCAGCGTTCACCGCATTACGCTCACGTTCGGCATATTTCTTGTCGAGCa       >  2:371414/1‑139 (MQ=255)                                                                                                                   cTTTGAACCGGGGTGCGCCGGGTTAATGGTTTCTGCGCTGACCTGTGCCATGCGCATCCCGTCACGCGTACGCGCCATGGTTAATTCAGCGTTCACCGCATTACACTCACGTTCGGCATATTTCTTGTCGAGCAAAGGt  >  2:126586/1‑139 (MQ=255)                                                                                                                      tGAACCGGGGTGCGCCGGGTTAATGGTTTCTGCGCTGACCTGTGCCATGCGCATcc                                                                                  >  1:266875/1‑56 (MQ=255)                                                                                                                      tGAACCGGGGTGCGCCGGGTTAATGGTTTCTGCGCTGACCTGTGCCATGCGCATcc                                                                                  <  2:266875/56‑1 (MQ=255)                                                                                                                                  |                                                                                                                            AATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTGCTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGTGCCGGGTTAATGGTTTCTGCGCTGACCTGTGCCATGCGCATCCCGTCACGCGTACGCGCCATGGTTAATTCAGCGTTCACCGCATTACGCTCACGTTCGGCATATTTCTTGTCGAGCAAAGGT  >  NC_000913/2958180‑2958432

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