Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,388,120	C→T	A417A (GCG→GCA)	mscM ←	miniconductance mechanosensitive channel MscM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,388,120	0	C	T	100.0%	97.9 / NA	29	A417A (GCG→GCA)	mscM	miniconductance mechanosensitive channel MscM
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/13);  total (16/13)

AGCGGCGCGAGTAAATACTGCAACCGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGCTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGCGCATCCTCCAGTTGCCCGTTGGAGACTTTCAGCTTGGTCAGTTCCAATAGTAGCGTGTCGCCACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCGTCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCGGTATCA  >  NC_000913/4387888‑4388363                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    aGCGGCGCGAGTAAATACTGCAACCGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGCTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGTGCATCCTCCAGTTGCCCg                                                                                                                                                                                                                                   >  1:18364/1‑251 (MQ=255)    ggcgCGAGTAAATACTGCAACCGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGCTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTGGCTTCGTTCACCTCTTTCAGTGCATCCTCCAGTTGCCCGTTg                                                                                                                                                                                                                                <  1:276752/251‑1 (MQ=255)         gAGTAAATACTGCAACCGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGCTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGTGCATCCTCCAGTTGCCCGTTGGAGAc                                                                                                                                                                                                                           <  2:227363/251‑1 (MQ=255)          aGTAAATCCTGCAACCGACCAGAATCAACGCGCCAACCAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGCTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGTGCATCCTCCAGTTGCCCGTTGGAGACt                                                                                                                                                                                                                          >  2:30221/1‑251 (MQ=255)              aaTACTGCAACCGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGCTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGTGCATCCTCCAGTTGCCCGTTGGAGACTTTCa                                                                                                                                                                                                                      <  2:236432/251‑1 (MQ=255)                        ccGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGCTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGTGCATCCTCCAGTTGCCCGTTGGAGACTTTCAGCTTGGTc                                                                                                                                                                                                              >  1:134296/1‑249 (MQ=255)                           aCCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGCTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGTGCATCCTCCAGTTGCCCGTTGGAGACTTTCAGCTTGGTCAGTTc                                                                                                                                                                                                         >  1:79082/1‑251 (MQ=255)                             cAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGCTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGTGCATCCTCCAGTTGCCCGTTGGAGACTTTCAGCTTGGTCAGTTcc                                                                                                                                                                                                        >  2:51760/1‑250 (MQ=255)                              aGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGCTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGTGCATCCTCCAGTTGCCCGTTGGAGACTTTCAGCTTGGTCAGTTCCa                                                                                                                                                                                                       >  2:118112/1‑250 (MQ=255)                              aGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGCTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCCGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGTGCATCCTCCAGTTGCCCGTTGGAGACTTTCAGCTTGGTCAGTTCCa                                                                                                                                                                                                       >  1:7781/1‑250 (MQ=255)                                   cAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACGGGCTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGTGCATCCTCCAGTGGCCCGTTGGAGACTTTCATCTTGGTCAGTTCCAAtagt                                                                                                                                                                                                  <  2:193555/250‑1 (MQ=255)                                                            gTCTCTTTGCTGGTCAGCATCATCACACTGGCTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGTGCATCCTCCAGTTGCCCGTTGGAGACTTTCAGCTTGGTCAGTTCCAATAGTAGCGTGTCGCCACCCTGCAACAATg                                                                                                                                                                         >  1:287982/1‑250 (MQ=255)                                                              ctctTTGCTGGTCAGCATCATCACACTGGCTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGTGCATCCTCCAGTTGCCCGTTGGAGACTTTCAGCTTGGTCAGTTCCAATAGTAGCGTGTCTCCACCCTGCAACAATGAGt                                                                                                                                                                      >  2:252937/1‑251 (MQ=255)                                                                  ttGCTGGTCAGCATCATCACACTGGCTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGTGCATCCTCCAGTTGCCCGTTGGAGACTTTCAGCTTGGTCAGTTCCACTAGTAGCGTGTCGCCACCCTGCAACAATGAGTTc                                                                                                                                                                    >  2:244881/1‑249 (MQ=255)                                                                   tGCTGGTCAGCATCATCACACTGGCTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGTGCATCCTCCAGTTGCCCGTTGGAGACTTTCAGCTTGGTCAGTTCCAATAGTAGCGTGTCGCCACCCTGCAACAATGAGTTCAg                                                                                                                                                                  <  1:26259/250‑1 (MQ=255)                                                                       ggTCAGCATCATCACACTGGCTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTGCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGTGCATCCTCCAGTTGCCCGTTGGAGACTTTCAGCTTGGTCAGTTCCAATAGTAGCGTGTCGCCACCCTGCAACAATGAGTTCAGCAAc                                                                                                                                                              <  1:81036/250‑1 (MQ=255)                                                                                    ccacTGGCTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGTGCATCCTCCAGTTGCCCGTTGGAGACTTTCAGCTTGGTCAGTTCCAATAGTAGCGTGTCGCCACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTgcg                                                                                                                                                <  2:17601/250‑1 (MQ=255)                                                                                      acTGGCTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGTGCATCCTCCAGTTGCCCGTTGGAGACTTTCAGCGTGGTCAGTTCCAATAGTAGCGTGTCGCCACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCa                                                                                                                                              <  2:100453/251‑1 (MQ=255)                                                                                      acTGGCTTTGCACAAATGACTGAAGGTGTCCAGCGCAATGAGACGACGCAGCTCCTAGGCGCTTTCCAGCGGCCAGGCGATGGTCCTCGGGCGCACGTCAGAGGCCCAGAACAGATAGCGGTGCGTTGCTTCGTTAACCTCTTTCAGTGCATCCTCCAGTTGCCCGTTGGAGACTTTCAGCTTGGTCAGTTCCAATAGTAGCGTGTCGCCACCCTGCAACACTGACTTCAGCAACTCACGCTGAGTgcgc                                                                                                                                               >  2:305029/1‑250 (MQ=255)                                                                                         ggCTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGTGCATCCTCCAGTTGCCCGTTGGAGACTTTCAGCTTGGTCAGTTCCAATAGTAGCGTGTCGCCACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAg                                                                                                                                             <  2:164877/249‑1 (MQ=255)                                                                                                      tgactgaAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGTGCATCCTCCAGTTGCCCGTTGGAGACTTTCAGCTTGGTCAGTTCCAATAGTAGCGTGTCGCCACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAAt                                                                                                                              >  1:195284/1‑251 (MQ=255)                                                                                                                                                 cGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGTGCATCCTCCAGTTGCCCGTTGGAGACTTTCAGCTTGGTCAGTTCCAATAGTAGCGTGTCGCCACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCGTCGGCCtgatg                                                                                    >  1:101299/1‑250 (MQ=255)                                                                                                                                                                 ggCGATGGTCATCGGGCGCACGTCAGATGTCCTGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGTGCATCCTCCATTTGTCCGTTGGAGACTTTCAGCTTGGTCAGTTCCAATAGTAGCGTGTCGCCACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAAGACGTTTTGGCTCTGCAGTCAGCGGCTGTCCGTCGTCCTGATGAATTTTCCGTAGCAGCg                                                                   <  1:270912/251‑1 (MQ=255)                                                                                                                                                                    gATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGTGCATCCTCCAGTTGCCCGTTGGAGACTTTCAGCTTGGTCAGTTCCAATAGTAGCGTGTCGCCACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCGTCGGCCTGATGAATTTGCCGTAGCAGCGGCt                                                                <  1:156666/251‑1 (MQ=255)                                                                                                                                                                                           agGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGTGCATCCTCCAGTTGCCCGTTGGAGACTTTCAGCTTGGTCAGTTCCAATAGTAGCGTGTCGCCACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGTTTTTGCTCGGCAGTCAGCGTCTGACCGTCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATaa                                         <  2:104311/251‑1 (MQ=255)                                                                                                                                                                                             gTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGTGCATCCTCCAGTTGCCCGTTGGAGACTTTCAGCTTGGTCAGTTCCAATAGTAGCGTGTCGCCACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCGTCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACg                                       >  1:127327/1‑251 (MQ=255)                                                                                                                                                                                                       gATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGTGCATCCTCCAGTTGCCCGTTGGAGACTTTCAGCTTGGTCAGTTCCAATAGTAGCGTGTCGCCACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCGTCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGc                             >  1:174471/1‑251 (MQ=255)                                                                                                                                                                                                                gCGTTGCTTCGTTCACCTCTTTCAGTGCATCCTCCAGTTGCCCGTTGGAGACTTTCAGCTTGGTCAGTTCCAATAGTAGCGTGTCGCCACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCGTCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAAc                    <  2:18364/251‑1 (MQ=255)                                                                                                                                                                                                                                    ttCAGTGCATCCTCCAGTTGCCCGTTGGAGACTTTCAGCTTGGTCAGTTCCAATAGTAGCGTGTCGCCACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCGTCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCGGTATCa  >  1:147705/1‑249 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    AGCGGCGCGAGTAAATACTGCAACCGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGCTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGCGCATCCTCCAGTTGCCCGTTGGAGACTTTCAGCTTGGTCAGTTCCAATAGTAGCGTGTCGCCACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCGTCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCGGTATCA  >  NC_000913/4387888‑4388363

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