Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	488,313	C→A	A593E (GCG→GAG)	mscK →	potassium dependent, small conductance mechanosensitive channel

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	488,313	0	C	A	100.0%	199.0 / NA	56	A593E (GCG→GAG)	mscK	potassium dependent, small conductance mechanosensitive channel
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (32/24);  total (32/24)

AACACACCAAAAGCGATCCTTATCGACCTGATCCGTGCGCTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGCGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGACCAGCCACTGGCGTCGGCAAATTGTCCGCATCAGTCTCGCATTG  >  NC_000913/488177‑488446                                                                                                                                         |                                                                                                                                      aaCACACCAAAAGCGATCCTTATCGACCTGATCCGTGCGCTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGat                                                                                                                                    <  1:603561/140‑1 (MQ=255)      acCAAAAGCGATCCTTATCGACCTGATCCGTGCGCTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTTGTTCAGCAAAAAACTGGAGATATTct                                                                                                                               >  1:732969/1‑140 (MQ=255)        cAAAAGCGATCCTTATCGACCTGATCCGTGCGCTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTctgg                                                                                                                             >  1:1028543/1‑140 (MQ=255)          aaaGCGATCCTTATCGACCTGATCCGTGCGCTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTctggct                                                                                                                           <  2:982223/140‑1 (MQ=255)                aTCCTTATCGACCTGATCCGTGCGCTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGtt                                                                                                                     >  1:37903/1‑140 (MQ=255)                  ccTTATCGACCTGATCCGTGCGCTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTg                                                                                                                   >  2:194206/1‑140 (MQ=255)                    ttATCGACCTGATCCGTGCGCTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGc                                                                                                                 <  1:118817/140‑1 (MQ=255)                    ttATCGACCTGATCCGTGCGCTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGc                                                                                                                 <  2:711806/140‑1 (MQ=255)                      aTCGACCTGATCCGTGCGCTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCt                                                                                                               <  1:711359/140‑1 (MQ=255)                       tCGACCTGATCCGTGCGCTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCtg                                                                                                              <  2:37903/140‑1 (MQ=255)                        cGACCTGATCCGTGCGCTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCtgt                                                                                                             >  2:258290/1‑140 (MQ=255)                          aCCTGATCCGTGCGCTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGc                                                                                                           <  1:953730/140‑1 (MQ=255)                           ccTGATCCGTGCGCTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCt                                                                                                               >  1:847975/1‑135 (MQ=255)                           ccTGATCCGTGCGCTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCt                                                                                                               <  2:847975/135‑1 (MQ=255)                              gATCCGTGCGCTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGa                                                                                                       <  1:258290/140‑1 (MQ=255)                               aTCCGTGCGCTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGc                                                                                                                 <  1:1000820/129‑1 (MQ=255)                               aTCCGTGCGCTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGc                                                                                                                 >  2:1000820/1‑129 (MQ=255)                                tCCGTGCGCTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAg                                                                                                     >  1:822410/1‑140 (MQ=255)                                   gTGCGCTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTa                                                                                                  <  1:577416/140‑1 (MQ=255)                                    tGCGCTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTAc                                                                                                 >  2:854520/1‑140 (MQ=255)                                    tGCGCTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTAc                                                                                                 >  2:149996/1‑140 (MQ=255)                                      cgcTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTg                                                                                               >  2:937689/1‑140 (MQ=255)                                                       attattCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGc                                                                              >  2:368835/1‑140 (MQ=255)                                                           ttCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTAGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTa                                                                          >  2:742218/1‑140 (MQ=255)                                                             ctcGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACg                                                                        >  2:380132/1‑140 (MQ=255)                                                                 cgGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCAc                                                                    >  1:354273/1‑140 (MQ=255)                                                                     tGGCCTGATTCTGTTGACCATGCCGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGCTATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCg                                                                >  2:918387/1‑140 (MQ=255)                                                                      ggCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCgg                                                               <  2:608106/140‑1 (MQ=255)                                                                        ccTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCa                                                             <  1:380132/140‑1 (MQ=255)                                                                             ttCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCg                                                        <  2:142366/140‑1 (MQ=255)                                                                              tCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCgg                                                       <  1:475684/140‑1 (MQ=255)                                                                              tCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATAT‑‑TGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGaa                                                     >  1:642017/1‑140 (MQ=255)                                                                                  ttGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAAca                                                   >  1:428709/1‑140 (MQ=255)                                                                                   tGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAAcag                                                  <  1:937689/140‑1 (MQ=255)                                                                                      ccATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAAcagcag                                               >  1:282984/1‑140 (MQ=255)                                                                                      ccATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAAcagcag                                               >  2:685586/1‑140 (MQ=255)                                                                                      ccATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAAcagcag                                               >  2:240818/1‑140 (MQ=255)                                                                                         tGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGAcc                                            >  1:960101/1‑140 (MQ=255)                                                                                          gCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGACCa                                           >  2:1123901/1‑140 (MQ=255)                                                                                          gCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGACCa                                           >  1:397090/1‑140 (MQ=255)                                                                                          gCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGACCa                                           >  1:874955/1‑140 (MQ=255)                                                                                               tCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGACCAGCCAc                                      >  1:1072797/1‑140 (MQ=255)                                                                                                   cATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGACCAGCCACTGGc                                  >  2:205556/1‑140 (MQ=255)                                                                                                   cATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGAGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGACCAGCCACTGGc                                  >  1:59660/1‑140 (MQ=255)                                                                                                            aCTGCTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGACCAGCCACTGGCGTCGGCaaa                         <  1:194206/140‑1 (MQ=255)                                                                                                                cTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTCCGTCACTTCGGCATGCGGGAACAGCAGACCAGCCACTGGCGTCGGCAAATTGt                     <  2:428709/140‑1 (MQ=255)                                                                                                                cTATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGACCAGCCACTGGCGTCGGCAAATTGt                     >  1:610061/1‑140 (MQ=255)                                                                                                                 tATGGTCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGACCAGCCACTGGCGTCGGCAAATTGTc                    <  1:685586/140‑1 (MQ=255)                                                                                                                      tCGTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGACCAGCCACTGGCGTCGGCAAATTGTCCGCAt               <  1:149996/140‑1 (MQ=255)                                                                                                                        gTTCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGACCAGCCACTGGCGTCGGCAAATTGTCCGCATCa             >  2:529968/1‑140 (MQ=255)                                                                                                                         ttCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGACCAGCCACTGGCGTCGGCAAATTGTCCGCATCAg            <  2:1068153/140‑1 (MQ=255)                                                                                                                          tCAGCAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGACCAGCCACTGGCGTCGGCAAATTGTCCGCATCAGt           >  1:9377/1‑140 (MQ=255)                                                                                                                              cAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGACCAGCCACTGGCGTCGGCAAATTGTCCGCATCAGTCTCg       >  1:830546/1‑140 (MQ=255)                                                                                                                              cAAAAAACTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGACCAGCCACTGGCGTCGGCAAATTGTCCGCATCAGTCTCg       <  1:1123901/140‑1 (MQ=255)                                                                                                                                  aaaCTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGACCAGCCACTGGCGTCGGCAAATTGTCCGCATCAGTCTCGCAtt   <  2:610061/140‑1 (MQ=255)                                                                                                                                   aaCTGGAGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGACCAGCCACTGGCGTCGGCAAATTGTCCGCATCAGTCTCGCATtg  <  2:732969/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                      AACACACCAAAAGCGATCCTTATCGACCTGATCCGTGCGCTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGCGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGACCAGCCACTGGCGTCGGCAAATTGTCCGCATCAGTCTCGCATTG  >  NC_000913/488177‑488446

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