Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,069,351	T→C	T153A (ACC→GCC)	mscS ←	mechanosensitive channel protein, small conductance

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,069,351	0	T	C	100.0%	137.0 / NA	39	T153A (ACC→GCC)	mscS	mechanosensitive channel protein, small conductance
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (23/16);  total (23/16)

CCTGATCGATATCGGAATCATACGCCACGCCAATAATAAATTCGTTACGGCGAACTGGCTCGCGGGAGAAGTTAATAATATTTCCGGCAATAATTTTACCGTTCGGAATAACGATAATTTTACCGTCTGCAGTACGCATGGTGGTGGAGAAAATCTGCACACTCAGCACAGTACCGGCTACGCCGCCCAGGTCAACATATTCTCCGGCACGGAACGGGCGGAACATGACAAGTAACACGCCAGCGGCCAGGTTAGACAGTGACCCCTGCAAAGCCAGACCAACAGCTAAGC  >  NC_000913/3069207‑3069497                                                                                                                                                 |                                                                                                                                                   ccTGATCGATATCGGAATCATACGCCACGCCAATAATAAATTCGTTACGGCGAACTGGCTCGCGGGAGAAGTTAATAATATTTCCGGCAATAATTTTACCGTTCGGAATAACGATAATTTTACCGTCTGCAGTACGCATGGTGGCGgaga                                                                                                                                               >  2:699606/1‑150 (MQ=255) ccTGATCGATATCGGAATCATACGCCACGCCAATAATAAATTCGTTACGGCGAACTGGCTCGCGGGAGAAGTTAATAATATTTCCGGCAATAATTTTACCGTTCGGAATAACGATAATTTTACCGTCTGCAGTACGCATGGTGGCGgaga                                                                                                                                               >  1:54150/1‑150 (MQ=255)                 aTCATACGCCACGCCAATAATAAATTCGTTACGGCGAACTGGCTCGCGGGAGAAGTTAATAATATTTCCGGCAATAATTTTACCGTTCGGAATAACGATAATTTTACCGTCTGCAGTACGCATGGTGGCGGAGAAAATCTGCACACTCAGc                                                                                                                              <  2:146871/151‑1 (MQ=255)                 aTCATACGCCACGCCAATAATAAATTCGTTACGGCGAACTGGCTCGCGGGAGAAGTTAATAATATTTCCGGCAATAATTTTACCGTTCGGAATAACGATAATTTTACCGTCTGCAGTACGCATGGTGGCGGAGAAAATCTGCACACTCAGc                                                                                                                              <  1:415055/151‑1 (MQ=255)                   cATACGCCACGCCAATAATAAATTCGTTACGGCGAACTGGCTCGCGGGAGAAGTTAATAATATTTCCGGCAATAATTTTACCGTTCGGAATAACGATAATTTTACCGTCTGCAGTACGCATGGTGGCGGAGAAAATCTGCACACTCAGca                                                                                                                             <  1:518458/150‑1 (MQ=255)                    aTACGCCACGCCAATAATAAATTCGTTACGGCGAACTGGCTCGCGGGAGAAGTTAATAATATTTCCGGCAATAATTTTACCGTTCGGAATAACGATAATTTTACCGTCTGCAGTACGCATGGTGGCGGAGAAAATCTGCACACTCAGcac                                                                                                                            >  2:658091/1‑150 (MQ=255)                        gccacgccaATAATAAATTCGTTACGGCGAACTGGCTCGCGGGAGAAGTTAATAATATTTCCGGCAATAATTTTACCGTTCGGAATAACGATAATTTTACCGTCTGCAGTACGCATGGTGGCGGAGAAAATCTGCACACTCAGCACAGTAc                                                                                                                       <  1:511642/151‑1 (MQ=255)                          cacgccaATAATAAATTCGTTACGGCGAACTGGCTCGCGGGAGAAGTTAATAATATTTCCGGCAATAATTTTACCGTTCGGAATAACGATAATTTTACCGTCTGCAGTACGCATGGTGGCGGAGAAAATCTGCACCCTCAGCACAGTACCg                                                                                                                     <  1:635387/151‑1 (MQ=255)                           acgccaATAATAAATTCGTTACGGCGAACTGGCTCGCGGGAGAAGTTAATAATATTTCCGGCAATAATTTTACCGTTCGGAATAACGATAATTTTACCGTCTGCAGTACGCATGGTGGCGGAGAAAATCTGCACACTCAGCACAGTACCgg                                                                                                                    >  1:766863/1‑151 (MQ=255)                           acgccaATAATAAATTCGTTACGGCGAACTGGCTCGCGGGAGAAGTTAATAATATTTCCGGCAATAATTTTACCGTTCGGAATAACGATAATTTTACCGTCTGCAGTACGCATGGTGGCGGAGAAAATCTGCACACTCAGCACAGTACCgg                                                                                                                    >  2:367070/1‑151 (MQ=255)                                     taaATTCGTTACGGCGAACTGGCTCGCGGGAGAAGTTAATAATATTTCCGGCAATAATTTTACCGTTCGGAATAACGATAATTTTACCGTCTGCAGTACGCATGGTGGCGGAGAAAATCTGCACACTCAGCACAGTACCGGCTACgccgcc                                                                                                          >  1:321922/1‑151 (MQ=255)                                       aaTTCGTTACGGCGAACTGGCTCGCGGGAGAAGTTAATAATATTTCCGGCAATAATTTTACCGTTCGGAATAACGATAATTTTACCGTCTGCAGTACGCATGGTGGCGGAGAAAATCTGCACACTCAGCACAGTACCGGCTACGCCGccc                                                                                                         >  2:63965/1‑150 (MQ=255)                                          tCGTTACGGCGAAATGGCTCGCGGGAGAAGTTAATAATTTTTCCGGCAATAATTTTACCGTGCGGAATAACGATAATTTTACCGTCGTCAGGACTCATGGTGGCGGAGAAAATGTGCACACTCAGCACAGTACAGGCAAGGCCGCCCAGGt                                                                                                     >  1:426175/1‑151 (MQ=255)                                            gTTACGGCGAACTGGCTCGCGGGAGAAGTTAATAATATTTCCGGCAATAATTTTACCGTTCGGAATAACGATAATTTTACCGTCTGCAGTACGCATGGTGGCGGAGAAAATCTGCACACTCAGCACAGTACCGGCTACGCCGCCCAGGTCa                                                                                                   <  1:699532/151‑1 (MQ=255)                                             ttACGGCGAACTGGCTCGCGGGAGAAGTTAATAATATTTCCGGCAATAATTTTACCGTTCGGAATAACGATAATTTTACCGTCTGCAGTACGCATGGTGGCGGAGAAAATCTGCACACTCAGCACAGTACCGGCTACGCCGCCCAGGTCa                                                                                                   >  1:632019/1‑150 (MQ=255)                                              tACGGCGAACTGGCTCGCGGGAGAAGTTAATAATATTTCCGGCAATAATTTTACCGTTCGGAATAACGATAATTTTACCGTCTGCAGTACGCATGGTGGCGGAGAAAATCTGCACACTCAGCACAGTACCGGCTACGCCGCCCAGGTCaa                                                                                                  <  1:378572/150‑1 (MQ=255)                                              tACGGCGAACTGGCTCGCGGGAGAAGTTAATAATATTTCCGGCAATAATTTTACCGTTCGGAATAACGATAATTTTACCGTCTGCAGTACGCATGGTGGCGGAGAAAATCTGCACACTCAGCACAGTACCGGCTACGCCGCCCAGGTCAAc                                                                                                 >  1:52635/1‑151 (MQ=255)                                                      aCTGGCTCGCGGGAGAAGTTAATAATATTTCCGGCAATAATTTTACCGTTCGGAATAACGATAATTTTACCGTCTGCAGTACGCATGGTGGCGGAGAAAATCTGCACACTCAGCACAGTACCGGCTACGCCGCCCAGGTCAACATATtctc                                                                                         <  1:620065/151‑1 (MQ=255)                                                         ggCTCGCGGGAGAAGTTAATAATATTTCCGGCAATAATTTTACCGTTCGGAATAACGATAATTTTACCGTCTGCAGTACGCATGGTGGCGGAGAAAATCTGCACACTCAGCACAGTACCGGCTACGCCGCCCAGGTCAACATATTCTCCgg                                                                                      <  2:310551/151‑1 (MQ=255)                                                                       gTTAATAATATTTCCGGCAATAATTTTACCGTTCGGAATAACGATAATTTTACCGTCTGCAGTACGCATGGTGGCGGAGAAAATCTGCACACTCAGCACAGTACCGGCTACGCCGCCCAGGTCAACATATTCTCCGGCACGGAACGGGCgg                                                                        <  2:52666/151‑1 (MQ=255)                                                                           ataataTTTCCGGCAATAATTTTACCGTTCGGAATAACGATAATTTTACCGTCTGCAGTACGCATGGTGGCGGAGAAAATCTGCACACTCAGCACAGTACCGGCTACGCCGCCCAGGTCAACATATTCTCCGGCACGGAACGGGCGGAAc                                                                     >  2:679217/1‑150 (MQ=255)                                                                             aataTTTCCGTCAATAATTTTACCGTTCGGAATAACGATAATTTTACCGTCTGCAGTACGCATGGTGGCGGAGAAAATCTGCACACTCAGCACAGTACCGGCTACGCCGCCCAGGTCAACATATTCTCCGGCACGGAACGGGCGGAACATg                                                                  >  1:81342/1‑151 (MQ=255)                                                                              atatTTCCGGCAATAATTTTCCCGCTCGGAATAACGATAATTTTACCGCCTGCAGTACGCATGGTCGCGGAGAAAATCTGCCCGCTCAGCACAGTTCCGCCTACGCCGCCCAGGTCAATATATGCTCCGGCACGGAACGGGCGTAACATg                                                                  <  1:577292/150‑1 (MQ=255)                                                                                 tttCCGGCAATAATTTTACCGTTCGGAATAACGATAATTTTACCGTCTGCAGTACGCATGGTGGCGGAGAAAATCTGCACACTCAGCACAGTACCGGCTACGCCGCCCAGGTCAACATATTCTCCGGCACGGAACGGGCGGAACATGACa                                                               >  1:72575/1‑150 (MQ=255)                                                                                 tttCCGGCAATAATTTTACCGTTCGGAATAACGATAATTTTACCGTCTGCAGTACGCATGGTGGCGGAGAAAATCTGCACACTCAGCACAGTACCGGCTACGCCGCCCAGGTCAACATATTCTCCGGCACGGAACGGGCGGAACATGACa                                                               >  1:292504/1‑150 (MQ=255)                                                                                        cAATAATTTTACCGTTCGGAATAACGATAATTTTACCGTCTGCAGTACGCATGGTGGCGGAGAAAATCTGCACACTCAGCACAGTACCGGCTACGCCGCCCAGGTCAACATATTCTCCGGCACGGAACGGGCGGAACATGACAAGTAacac                                                       <  2:60013/151‑1 (MQ=255)                                                                                          ataatTTTACCGTTCGGAATAACGATAATTTTACCGTCTGCAGTACGCATGGTGGCGGAGAAAATCTGCACACTCAGCACAGTACCGGCTACGCCGCCCAGGTCAACATATTCTCCGGCACGGAACGGGCGGAACATGACAAGTAACACGc                                                     <  1:480382/151‑1 (MQ=255)                                                                                              ttttACCGTTCGGAATAACGATAATTTTACCGTCTGCAGTACGCATGGTGGCGGAGAAAATCTGCACACTCAGCACAGTACCGGCTACGCCGCCCAGGTCAACATATTCTCCGGCACGGAACGGGCGGAACATGACAAGTAACACGCCAGc                                                 >  1:709396/1‑151 (MQ=255)                                                                                               tttACCGTTCGGAATAACGATAATTTTACCGTCTGCAGTACGCATGGTGGCGGAGAAAATCTGCACACTCAGCACAGTACCGGCTACGCCGCCCAGGTCAACATATTCTCCGGCACGGAACGGGCGGAACATGACAAGTAACACGCCAGc                                                 <  1:90530/150‑1 (MQ=255)                                                                                                    cGTTCGGAATAACGATAATTTTACCGTCTGCAGTACGCATGGTGGCGGAGAAAATCTGCACACTCAGCACAGTACCGGCTACGCCGCCCAGGTCAACATATTCTCCGGCACGGAACGGGCGGAACATGACAAGTAACACGCCAGCGGcc                                             >  1:284333/1‑149 (MQ=255)                                                                                                          gAATAACGATAATTTTACCGTCTGCAGTACGCATGGTGGCGGAGAAAATCTGCACACTCAGCACAGTACCGGCTACGCCGCCCAGGTCAACATATTCTCCGGCACGGAACGGGCGGAACATGACAAGTAACACGCCAGCGGCCAGGTTAGa                                     <  2:24659/151‑1 (MQ=255)                                                                                                             tAACGATAATTTTACCGTCTGCAGTACGCATGGTGGCGGAGAAAATCTGCACACTCAGCACAGTACCGGCTACGCCGCCCAGGTCAACATATTCTCCGGCACGGAACGGGCGGAACATGACAAGTAACACGCCAGCGGCCAGGTTAGAc                                    >  1:332613/1‑149 (MQ=255)                                                                                                               aCGATAATTTTACCGTCTGCAGTACGCATGGTGGCGGAGAAAATCTGCACACTCAGCACAGTACCGGCTACGCCGCCCAGGTCAACATATTCTCCGGCACGGAACGGGCGGAACATGACAAGTAACACGCCAGCGGCCAGGTTAGACAGTg                                >  2:769090/1‑151 (MQ=255)                                                                                                                    aaTTTTACCGTCTGCAGTACGCATGGTGGCGGAGAAAATCTGCACACTCAGCACAGTACCGGCTACGCCGCCCAGGTCAACATATTCTCCGGCACGGAACGGGCGGAACATGACAAGTAACACGCCAGCGGCCAGGTTAGACAGTGAcccc                           >  1:450663/1‑151 (MQ=255)                                                                                                                           ccGTCTGCAGTACGCATGGTGGCGGAGAAAATCTGCACACTCAGCACAGTACCGGCTACGCCGCCCAGGTCAACATATTCTCCGGCACGGAACGGGCGGAACATGACAAGTAACACGCCAGCGGCCAGGTTAGACAGTGACCCCTGCAAAg                    <  2:55629/151‑1 (MQ=255)                                                                                                                                  cAGTACGCATGGTGGCGGAGAAAATCTGCACACTCAGCACAGTACCGGCTACGCCGCCCAGGTCAACATATTCTCCGGCACGGAACGGGCGGAACATGACAAGTAACACGCCAGCGGCCAGGTTAGACAGTGACCCCTGCAAAGCCAGAcc             >  1:53496/1‑151 (MQ=255)                                                                                                                                   aGTACGCATGGTGGCGGAGAAAATCTGCACACTCAGCACAGTACCGGCTACGCCGCCCAGGTCAACATATTCTCCGGCACGGAACGGGCGGAACATGACAAGTAACACGCCAGCGGCCAGGTTAGACAGTGACCCCTGCAAAGCCAGAcc             >  2:24735/1‑150 (MQ=255)                                                                                                                                        gCATGGTGGCGGAGAAAATCTGCACACTCAGCACAGTACCGGCTACGCCGCCCAGGTCAACATATTCTCCGGCACGGAACGGGCGGAACATGACAAGTAACACGCCAGCGGCCAGGTTAGACAGTGACCCCTGCAAAGCCAGACCAACAGc       >  1:107788/1‑151 (MQ=255)                                                                                                                                         cATGGTGGCGGAGAAAATCTGCACACTCAGCACAGTACCGGCTACGCCGCCCAGGTCAACATATTCTCCGGCACGGAACGGGCGGAACATGACAAGTAACACGCCAGCGGCCAGGTTAGACAGTGACCCCTGCAAAGCCAGACCAACAGCt      >  2:699401/1‑151 (MQ=255)                                                                                                                                             gtggcggAGAAAATCTGCACACTCAGCACAGTACCGGCTACGCCGCCCAGGTCAACATATTCTCCGGCACGGAACGGGCGGAACATGACAAGTAACACGCCAGCGGCCAGGTTAGACAGTGACCCCTGCAAAGCCAGACCAACAGCTAAGc  <  2:214559/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                   CCTGATCGATATCGGAATCATACGCCACGCCAATAATAAATTCGTTACGGCGAACTGGCTCGCGGGAGAAGTTAATAATATTTCCGGCAATAATTTTACCGTTCGGAATAACGATAATTTTACCGTCTGCAGTACGCATGGTGGTGGAGAAAATCTGCACACTCAGCACAGTACCGGCTACGCCGCCCAGGTCAACATATTCTCCGGCACGGAACGGGCGGAACATGACAAGTAACACGCCAGCGGCCAGGTTAGACAGTGACCCCTGCAAAGCCAGACCAACAGCTAAGC  >  NC_000913/3069207‑3069497

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