Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	489,541	T→G	R1002R (CGT→CGG)	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	489,541	0	T	G	100.0%	156.4 / NA	43	R1002R (CGT→CGG)	mscK	potassium dependent, small conductance mechanosensitive channel
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (23/20);  total (23/20)

CGGTAAGTAAGATCCGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCA  >  NC_000913/489414‑489667                                                                                                                                |                                                                                                                               cGGTAAGTAAGATCCGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGGCTGGt                                                                                                                           >  1:596576/1‑133 (MQ=255)         aagATCCGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGGCTGGTGATCCGtc                                                                                                                   <  2:532927/133‑1 (MQ=255)         aagATCCGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGGCTGGTGATCCGtc                                                                                                                   <  1:411628/133‑1 (MQ=255)            aTCCGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGGCTGGTGATCCGTCTCg                                                                                                                >  2:382093/1‑133 (MQ=255)            aTCCGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGGCTGGTGATCCGTCTCg                                                                                                                >  1:80280/1‑133 (MQ=255)            aTCCGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCAT‑CAGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGGCTGGTGATCCGTCTCgg                                                                                                               >  2:161772/1‑133 (MQ=255)              ccGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGGCTGGTGATCCGTCTCGGc                                                                                                              >  2:533104/1‑133 (MQ=255)               cGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTTGTTGACTGACACAAGTACGCGGCTGGTGATAGGTCTCGGCg                                                                                                             >  2:283127/1‑133 (MQ=255)               cGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGGCTGGTGATCCGTCTCGGCg                                                                                                             <  2:80280/133‑1 (MQ=255)                 tATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGCTGACTGACACTAGTACCCGGCTGGTGATCCATCTCGGCGTg                                                                                                           >  2:280081/1‑133 (MQ=255)                      gTGCGACAACGATTATCGATTTAGATCGCATAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGGCTGGTGATCCGTCTCGGCGTGGCCTa                                                                                                      <  1:283127/133‑1 (MQ=255)                      gTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGGCTGGTGATCCGTCTCGGCGTGGCCTa                                                                                                      <  1:352000/133‑1 (MQ=255)                      gTGCGACAACAATTACCGATTTCGTTCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGGCTGGTGATCCGTCTCGGCGTGGCCTa                                                                                                      <  2:243068/133‑1 (MQ=255)                               cGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGGCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGa                                                                                             <  2:596576/133‑1 (MQ=255)                                    aCCGATTTCGATCGCAAAGAAGTGATCATCCCGCACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGGCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTgg                                                                                        >  1:626831/1‑133 (MQ=255)                                           tCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGGCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGt                                                                                 >  2:462956/1‑133 (MQ=255)                                                cGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGGCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTa                                                                            <  1:508691/133‑1 (MQ=255)                                                  cAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGGCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTaaa                                                                          >  2:411443/1‑133 (MQ=255)                                                  cAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGGCTGGTGATCCGTCTCGGCGTGGCCTATGGCTACGATCTGGAAAAAGTGCGTa                                                                            >  1:499622/1‑131 (MQ=255)                                                  cAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGGCTGGTGATCCGTCTCGGCGTGGCCTATGGCTACGATCTGGAAAAAGTGCGTa                                                                            <  2:499622/131‑1 (MQ=255)                                                  cAAAGAAGGGATCATGCCCAACCAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGGCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTaaa                                                                          <  1:280081/133‑1 (MQ=255)                                                             tcatcCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGGCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGaa                                                               >  1:39062/1‑133 (MQ=255)                                                             tcatcCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGGCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGaa                                                               >  2:455779/1‑133 (MQ=255)                                                                tcCCGAACAAAGCGTCTGTTACCGCGCGTCTGCTCAACTGGCCGTTGACTGACACTACTACGCGGCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAggc                                                            <  2:179303/133‑1 (MQ=255)                                                                 cccGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGGCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGGTCTGGAAAAAGTGCGTAAAGTGTTACTGAAggcg                                                           >  1:427286/1‑133 (MQ=255)                                                                 cccGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGGCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAggcg                                                           >  2:528792/1‑133 (MQ=255)                                                                       cAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGGCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACt                                                     <  1:161772/133‑1 (MQ=255)                                                                         aaGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGGCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGa                                                   <  1:462956/133‑1 (MQ=255)                                                                         aaGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGGCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGa                                                   <  1:382093/133‑1 (MQ=255)                                                                            cGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGGCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCa                                                <  2:427286/133‑1 (MQ=255)                                                                                tGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGGCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAg                                                  <  2:349736/127‑1 (MQ=255)                                                                                tGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGGCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAg                                                  >  1:349736/1‑127 (MQ=255)                                                                                tGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGGCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCa                                            <  2:39062/133‑1 (MQ=255)                                                                                      cGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGGCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTg                                      >  1:146729/1‑133 (MQ=255)                                                                                          cGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGGCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGc                                  >  2:29086/1‑133 (MQ=255)                                                                                                       ggTCGTTGACTGACACTACTACGCGGCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGcc                     >  2:498332/1‑133 (MQ=255)                                                                                                         tCGTTGACTGACACTACTACGCGGCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCgg                   <  1:528792/133‑1 (MQ=255)                                                                                                          cGTTGACTGACACTACTACGCGGCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCGCCCAAGGTTGATGCCCGAACCAATGCCGGa                  >  1:211139/1‑133 (MQ=255)                                                                                                             tgactgacACTACTACGCGGCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGt               <  1:533104/133‑1 (MQ=255)                                                                                                                   acacTACTACGCGGCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTcttctt         >  1:340902/1‑133 (MQ=255)                                                                                                                   acacTACTACGCGGCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTcttctt         >  1:567657/1‑133 (MQ=255)                                                                                                                    cacTACTACGCGGCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCttt        <  1:411443/133‑1 (MQ=255)                                                                                                                          tacGCGGCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCa  >  2:27656/1‑133 (MQ=255)                                                                                                                                |                                                                                                                               CGGTAAGTAAGATCCGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCA  >  NC_000913/489414‑489667

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