Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,389,327	(C)6→7	coding (44/3324 nt)	mscM ←	miniconductance mechanosensitive channel MscM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,389,321	1	.	C	100.0%	148.4 / NA	44	A17G (GCG→GGG)	mscM	miniconductance mechanosensitive channel MscM
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (19/25);  total (19/25)

GAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACG‑CCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCGTGTTCTGCTTCGAT  >  NC_000913/4389186‑4389448                                                                                                                                         |                                                                                                                                gAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACG‑cccc                                                                                                                             <  2:1245365/140‑1 (MQ=255)   aCCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACG‑cccccc                                                                                                                           >  1:588538/1‑140 (MQ=255)    ccTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACG‑ccccccc                                                                                                                          >  2:389094/1‑139 (MQ=255)    ccTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACG‑ccccccc                                                                                                                          <  1:547182/140‑2 (MQ=255)    ccTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACG‑ccccccc                                                                                                                          <  1:893840/140‑2 (MQ=255)    ccTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACG‑ccccccc                                                                                                                          >  1:303939/1‑139 (MQ=255)    ccTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACG‑ccccccc                                                                                                                          >  2:893840/1‑139 (MQ=255)    ccTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACG‑c                                                                                                                                <  2:660722/134‑1 (MQ=255)    ccTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACG‑c                                                                                                                                >  1:660722/1‑134 (MQ=255)     cTTTTCGTTCCTCAAGCGCATTCAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACG‑ccccccca                                                                                                                         <  1:1041481/140‑3 (MQ=255)          cGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGa                                                                                                                     <  2:503886/140‑1 (MQ=255)           gTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACG‑c                                                                                                                                <  1:305175/127‑1 (MQ=255)           gTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACG‑c                                                                                                                                >  2:305175/1‑127 (MQ=255)               cTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACG‑cccc                                                                                                                             >  1:81228/1‑126 (MQ=255)               cTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACG‑cccc                                                                                                                             <  2:81228/126‑1 (MQ=255)               cTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCAGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCAc                                                                                                                >  2:954102/1‑140 (MQ=255)                tCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCAcc                                                                                                               >  2:987332/1‑140 (MQ=255)                          ttAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAg                                                                                                             >  1:877502/1‑132 (MQ=255)                          ttAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAg                                                                                                             <  2:877502/132‑1 (MQ=255)                                  gaatGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGAGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGa                                                                                             >  2:392679/4‑140 (MQ=255)                                   aCTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGAt                                                                                            >  1:1421393/1‑140 (MQ=255)                                   aCTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGAt                                                                                            >  2:719459/1‑140 (MQ=255)                                        aGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCa                                                                                       >  2:309608/1‑140 (MQ=255)                                           gcCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGgc                                                                                    <  1:389094/140‑1 (MQ=255)                                            ccTCTACGACTTCCGGCTGTGCGGGTTTCGCAGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGgcg                                                                                   <  1:954102/140‑1 (MQ=255)                                             ctctACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGgcgc                                                                                  >  1:1167071/1‑140 (MQ=255)                                             ctctACGACTTCCGGCTGTGCGGGTTTCCCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGgcgc                                                                                  >  1:782636/1‑140 (MQ=255)                                               ctACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCAc                                                                                >  2:820774/1‑140 (MQ=255)                                                  cGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTc                                                                             <  1:719459/140‑1 (MQ=255)                                                      ttCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCg                                                                         <  1:820774/140‑1 (MQ=255)                                                       tctggCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGt                                                                        >  1:708467/4‑140 (MQ=255)                                                       tCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGt                                                                        >  2:583212/1‑140 (MQ=255)                                                       tCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGt                                                                        >  2:1232400/1‑140 (MQ=255)                                                       tCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGt                                                                        <  2:303939/140‑1 (MQ=255)                                                       tCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGt                                                                        <  2:487547/140‑1 (MQ=255)                                                          ggCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTc                                                                     >  2:947216/1‑140 (MQ=255)                                                               tgCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTg                                                                <  2:238720/140‑1 (MQ=255)                                                                gCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGa                                                               <  2:798371/140‑1 (MQ=255)                                                                gCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGa                                                               <  2:56200/140‑1 (MQ=255)                                                                gCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGa                                                               <  2:912760/140‑1 (MQ=255)                                                                gCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGa                                                               <  2:1167071/140‑1 (MQ=255)                                                                gCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGa                                                               <  1:987332/140‑1 (MQ=255)                                                                gCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGa                                                               <  1:947216/140‑1 (MQ=255)                                                                     tttCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGAt                                                          <  1:1164632/140‑1 (MQ=255)                                                                                     tGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCttttt                                          <  1:543200/140‑1 (MQ=255)                                                                                        tCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTgt                                       >  1:711192/1‑140 (MQ=255)                                                                                               ccTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAAc                                >  1:224745/1‑140 (MQ=255)                                                                                                 tGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCa                              <  2:1421393/140‑1 (MQ=255)                                                                                                 tGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCa                              >  1:1082937/1‑140 (MQ=255)                                                                                                 tGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCa                              >  2:897837/1‑140 (MQ=255)                                                                                                 tGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCa                              <  2:1082937/140‑1 (MQ=255)                                                                                                  gAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCaa                             >  1:787610/1‑140 (MQ=255)                                                                                                         ttgtttgCTATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTg                      <  2:787610/140‑1 (MQ=255)                                                                                                                 tATCGGGGGCCGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATAGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCgtg              <  1:1232400/140‑1 (MQ=255)                                                                                                                           cGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCGTGTTCTGCTTCg    <  2:711192/140‑1 (MQ=255)                                                                                                                           cGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCGTGTTCTGCTTCg    <  2:480403/140‑1 (MQ=255)                                                                                                                           cGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTCTTTGTCGTCAACCAATGGGCTg                      <  2:112782/122‑1 (MQ=255)                                                                                                                           cGTCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTCTTTGTCGTCAACCAATGGGCTg                      >  1:112782/1‑122 (MQ=255)                                                                                                                             tCGCGGCGTACGCCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCGTGTTCTGCTTCGAt  <  2:1177986/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                GAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACG‑CCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCGTGTTCTGCTTCGAT  >  NC_000913/4389186‑4389448

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