Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	485,143	A→G	D1051G (GAT→GGT)	mscK →	mechanosensitive channel MscK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	485,143	0	A	G	100.0%	188.4 / NA	53	D1051G (GAT→GGT)	mscK	mechanosensitive channel MscK
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/29);  total (24/29)

CGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATCAACATTGCCTTTAACCAGCTTGAAGTGCATCTGC  >  NZ_CP009273/485007‑485277                                                                                                                                         |                                                                                                                                       cggctcGGCGTTGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTc                                                                                                                                      <  1:700171/136‑1 (MQ=255)       ggCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAGTGCCGCCATCCCTCTTCACGGCATTTGGTGCCAGCACGTTGGGTCATGAg                                                                                                                                 <  1:339713/138‑1 (MQ=255)        gCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGGCTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCt                                                                                                                               <  1:646248/139‑1 (MQ=255)             gCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCt                                                                                                                               >  1:227520/1‑134 (MQ=255)             gCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCt                                                                                                                               <  2:227520/134‑1 (MQ=255)                  tGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTa                                                                                                                      <  2:775978/138‑1 (MQ=255)                  tGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGC‑CGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTAtg                                                                                                                    >  1:141905/1‑139 (MQ=255)                        cGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGt                                                                                                               >  1:597912/1‑139 (MQ=255)                        cGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGt                                                                                                               >  1:511717/1‑139 (MQ=255)                          aTCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTg                                                                                                              >  2:324771/1‑138 (MQ=255)                               gAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGc                                                                                                        >  1:820013/1‑139 (MQ=255)                                aaaaaGTGCGTAAAGGGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCg                                                                                                       <  1:229077/139‑1 (MQ=255)                                                     ctgaGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCg                                                                                   <  2:686544/135‑1 (MQ=255)                                                       aaGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGAt                                                                                 >  1:24254/1‑138 (MQ=255)                                                         ggcggcGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAg                                                                              >  2:311442/1‑139 (MQ=255)                                                         ggcggcGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAg                                                                              >  2:553759/1‑139 (MQ=255)                                                         ggcggcGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAg                                                                              >  1:272656/1‑139 (MQ=255)                                                           cggcgACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCt                                                                            >  1:422793/1‑139 (MQ=255)                                                           cggcgACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCt                                                                            <  1:578595/139‑1 (MQ=255)                                                            ggcgACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGTCAGCCCGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTg                                                                           <  1:148785/139‑1 (MQ=255)                                                            ggcgACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCGCGTTGGGTCATGAGCTGCGTCTGTAGGAGCGTGACCTGCGTGACCGTAGTAGTGCTGTCGATGAACTg                                                                           <  2:543848/139‑1 (MQ=255)                                                            ggcgACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTg                                                                           <  1:18566/139‑1 (MQ=255)                                                              cgACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGaa                                                                         <  2:597912/139‑1 (MQ=255)                                                              cgACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGaa                                                                         >  1:293613/1‑139 (MQ=255)                                                                         ccAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCg                                                              <  1:553759/139‑1 (MQ=255)                                                                         ccAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCg                                                              <  1:324771/139‑1 (MQ=255)                                                                         ccAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCg                                                              <  2:62270/139‑1 (MQ=255)                                                                              ggTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAg                                                         >  2:545250/1‑139 (MQ=255)                                                                                   tGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCt                                                                                                                               <  1:554412/64‑1 (MQ=255)                                                                                   tGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCt                                                                                                                               >  2:554412/1‑64 (MQ=255)                                                                                   tGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCtgtg                                                    >  2:542378/1‑139 (MQ=255)                                                                                   tGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCtgtg                                                    >  1:414693/1‑139 (MQ=255)                                                                                   tGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCtgtg                                                    >  1:630169/1‑139 (MQ=255)                                                                                     cACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGcc                                                  >  2:757565/1‑139 (MQ=255)                                                                                           ccAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGaaa                                            <  2:195409/139‑1 (MQ=255)                                                                                               tGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTa                                                                                                                      >  2:568789/1‑61 (MQ=255)                                                                                               tGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTa                                                                                                                      <  1:568789/61‑1 (MQ=255)                                                                                               tGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGAt                                                            >  2:234311/1‑119 (MQ=255)                                                                                               tGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGAt                                                            <  1:234311/119‑1 (MQ=255)                                                                                                 cctgAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGAc                                       <  2:414693/135‑1 (MQ=255)                                                                                                      aGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTg                                                                                                              >  2:382416/1‑62 (MQ=255)                                                                                                      aGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTg                                                                                                              <  1:382416/62‑1 (MQ=255)                                                                                                        tcttctTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATCAACAt                               <  1:435897/139‑1 (MQ=255)                                                                                                        tcttctTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATCAACAt                               <  1:542378/139‑1 (MQ=255)                                                                                                                 cGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATCAACATTGCCTTTaa                      <  1:311442/139‑1 (MQ=255)                                                                                                                 cGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATCAACATTGCCTTTaa                      <  1:861303/139‑1 (MQ=255)                                                                                                                  ggCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATCAACATTGCCTTTAAc                     >  1:75218/1‑139 (MQ=255)                                                                                                                    cATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATCAACATTGCCTTTAAc                     >  2:567779/1‑137 (MQ=255)                                                                                                                          gTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATCAACATTGCCTTTAACCAGCTTGa             <  2:293613/139‑1 (MQ=255)                                                                                                                                gCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGaaaa                                           >  1:699248/1‑103 (MQ=255)                                                                                                                                gCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGaaaa                                           <  2:699248/103‑1 (MQ=255)                                                                                                                                 cACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATCAACATTGCCTTTAACCAGCTTGAAGTGCAt      <  2:734185/139‑1 (MQ=255)                                                                                                                                     ttGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATCAACATTGCCTTTAACCAGCTTGAAGTGCATCTGc  <  1:545250/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                       CGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATCAACATTGCCTTTAACCAGCTTGAAGTGCATCTGC  >  NZ_CP009273/485007‑485277

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