Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	983,623	G→A	*183* (TAG→TAA)	ycbK →	M15A protease‑related family periplasmic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	983,623	0	G	A	100.0%	204.6 / NA	63	*183* (TAG→TAA)	ycbK	M15A protease‑related family periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (31/32);  total (31/32)

TCCATATTGAAGGTATCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAGTAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCTGGCGGCGATGCGGA  >  NC_000913/983496‑983756                                                                                                                                |                                                                                                                                      tCCATATTGAAGGTATCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTaa                                                                                                                            <  1:930049/139‑1 (MQ=255)       ttGAAGGTATCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAAc                                                                                                                      <  1:690985/139‑1 (MQ=255)           aGGTATCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAggg                                                                                                                  <  2:508391/139‑1 (MQ=255)            ggTATCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                                 <  1:267495/139‑1 (MQ=255)                   cgTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATg                                                                                                          <  2:1278555/139‑1 (MQ=255)                     ttAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGaa                                                                                                        <  1:1155340/139‑1 (MQ=255)                        aGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTa                                                                                                     <  2:380681/139‑1 (MQ=255)                                    aaaGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTcc                                                                                         >  1:751216/1‑139 (MQ=255)                                         cgcgTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTaa                                                                                                                            <  2:1347217/99‑1 (MQ=255)                                         cgcgTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTaa                                                                                                                            >  1:1347217/1‑99 (MQ=255)                                                      cgcgcAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGgtaata                                                                                                                                     <  2:875914/77‑1 (MQ=255)                                                        cgcAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                                 >  2:72068/1‑95 (MQ=255)                                                        cgcAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                                 <  1:72068/95‑1 (MQ=255)                                                          cAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCt                                                                                                                               <  1:374365/79‑1 (MQ=255)                                                          cAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCt                                                                                                                               >  2:374365/1‑79 (MQ=255)                                                                          ttatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGtt                                                               >  2:716848/1‑127 (MQ=255)                                                                          ttatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGtt                                                               <  1:716848/127‑1 (MQ=255)                                                                           tatCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGgtg                                                  >  2:13425/1‑139 (MQ=255)                                                                             tCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAAc                                                                                                                           <  2:1361372/64‑1 (MQ=255)                                                                             tCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAAc                                                                                                                           >  1:1361372/1‑64 (MQ=255)                                                                                aCGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTa                                                                                                     >  1:1329069/1‑83 (MQ=255)                                                                                aCGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTa                                                                                                     <  2:1329069/83‑1 (MQ=255)                                                                                aCGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAAc                                                                   >  2:743328/1‑117 (MQ=255)                                                                                aCGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAAc                                                                   <  1:743328/117‑1 (MQ=255)                                                                                    agtaACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAAc                                         >  2:1236164/1‑139 (MQ=255)                                                                                            tgtgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAg                                                                                                              >  2:445629/1‑62 (MQ=255)                                                                                            tgtgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAg                                                                                                              <  1:445629/62‑1 (MQ=255)                                                                                            tgtgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATg                                                                                                          >  1:661997/1‑66 (MQ=255)                                                                                            tgtgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATg                                                                                                          <  2:661997/66‑1 (MQ=255)                                                                                            tgtgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                              <  2:345783/94‑1 (MQ=255)                                                                                            tgtgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                              >  1:345783/1‑94 (MQ=255)                                                                                            tgtgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTgg                                 >  1:550903/1‑139 (MQ=255)                                                                                             gtgCATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGc                                <  1:1236164/139‑1 (MQ=255)                                                                                                cATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                                 >  2:286174/1‑55 (MQ=255)                                                                                                cATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                                 <  1:286174/55‑1 (MQ=255)                                                                                                cATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGt                                     >  1:301207/1‑131 (MQ=255)                                                                                                cATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGt                                     <  2:301207/131‑1 (MQ=255)                                                                                                cATATTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGc                             >  2:237898/1‑139 (MQ=255)                                                                                                 atatTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCa                            <  1:13425/139‑1 (MQ=255)                                                                                                   atTGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACt                          >  1:788623/1‑139 (MQ=255)                                                                                                    ttGATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTg                         >  2:343009/1‑139 (MQ=255)                                                                                                      gATACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGt                       <  2:751216/139‑1 (MQ=255)                                                                                                        tACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCa                                                                                    >  1:162733/1‑76 (MQ=255)                                                                                                        tACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCa                                                                                    <  2:162733/76‑1 (MQ=255)                                                                                                        tACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTcc                                                                         >  2:1255815/1‑87 (MQ=255)                                                                                                        tACCGGGCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTcc                                                                         <  1:1255815/87‑1 (MQ=255)                                                                                                             ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGt                                                                                                             >  1:1041639/1‑46 (MQ=255)                                                                                                             ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGt                                                                                                             <  2:1041639/46‑1 (MQ=255)                                                                                                             ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATg                                                                                                          <  1:566810/49‑1 (MQ=255)                                                                                                             ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATg                                                                                                          >  2:566810/1‑49 (MQ=255)                                                                                                             ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                              <  1:488417/77‑1 (MQ=255)                                                                                                             ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCAt                                                                              >  2:488417/1‑77 (MQ=255)                                                                                                             ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTGTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCt                >  1:836992/1‑139 (MQ=255)                                                                                                             ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGt                       >  1:146967/1‑132 (MQ=255)                                                                                                             ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGt                       <  2:146967/132‑1 (MQ=255)                                                                                                             ggCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCt                >  2:507274/1‑139 (MQ=255)                                                                                                              gCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                                 <  1:798365/41‑1 (MQ=255)                                                                                                              gCCAGCACGGCACTGGTAATAATTGCTTAACGAAACAgggg                                                                                                                 >  2:798365/1‑41 (MQ=255)                                                                                                                 aGCACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCTggcg            >  2:149154/1‑139 (MQ=255)                                                                                                                   cACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCa                                                                                                               <  2:462177/38‑1 (MQ=255)                                                                                                                   cACGGCACTGGTAATAATTGCTTAACGAAACAGGGGCa                                                                                                               >  1:462177/1‑38 (MQ=255)                                                                                                                        cACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGaa                                                                                                        <  1:123224/40‑1 (MQ=255)                                                                                                                        cACTGGTAATAATTGCTTAACGAAACAGGGGCAGTATGaa                                                                                                        >  2:123224/1‑40 (MQ=255)                                                                                                                           tGGTAATAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCTGGCGGCGATGCGGa  <  1:490040/139‑1 (MQ=255)                                                                                                                                |                                                                                                                                      TCCATATTGAAGGTATCGCGTTAAGCAATATTCGCAAAGCCGCGTTATCTATGCGCGCAGGTGGTGTAGGATATTATCCACGTAGTAACTTTGTGCATATTGATACCGGGCCAGCACGGCACTGGTAGTAATTGCTTAACGAAACAGGGGCAGTATGAACTATCGTATTATTCCGGTCACCGCATTCTCCCAGAACTGTTCATTAATCTGGTGTGAACAAACCCGTCTGGCCGCACTGGTCGATCCTGGCGGCGATGCGGA  >  NC_000913/983496‑983756

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