Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,387,000	G→A	L791L (CTG→TTG)	mscM ←	miniconductance mechanosensitive channel MscM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,387,000	0	G	A	100.0%	123.2 / NA	35	L791L (CTG→TTG)	mscM	miniconductance mechanosensitive channel MscM
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/18);  total (17/18)

GGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTCCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCGCCCGATCAAACGCCAGCCTGCGGCGCTGGATCAGCATCCAGCGGCGGATAACGTGATAAACCACCAGTAGC  >  NC_000913/4386755‑4387241                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  ggCGTAACCCGTACCCGGCGTTAAATCAAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAAGGCg                                                                                                                                                                                                                                               >  1:71552/1‑250 (MQ=255)  gCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAAGGCGAt                                                                                                                                                                                                                                             <  1:238163/251‑1 (MQ=255)   cGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAAGGCGATc                                                                                                                                                                                                                                            >  1:3224/1‑251 (MQ=255)         ccGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTGCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAATCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAAGGCGATCAACATc                                                                                                                                                                                                                                      <  1:164749/251‑1 (MQ=255)                 ggCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTGTCGAGGAAGCCGAAAGCGGAATGGATGTCTGACCACAGCACAATGACAGAAAGCAAGGCGATCAACATCAAAATTGa                                                                                                                                                                                                                              <  2:165063/251‑1 (MQ=255)                   cGTTAAATCCAGGGGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCGCCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATGTCTGACCACAGCACAATGACAGAAAGCAAGGCGATCAACATCAAAAGTGAgc                                                                                                                                                                                                                            <  2:119226/251‑1 (MQ=255)                           ccAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAAGGCGATCAACATCAAAATTGAGCGCACCAGc                                                                                                                                                                                                                    <  1:208140/251‑1 (MQ=255)                                    gcAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAAGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAgg                                                                                                                                                                                                             <  1:187390/249‑1 (MQ=255)                                                gTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAAGGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAAGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATg                                                                                                                                                                                               <  2:134679/251‑1 (MQ=255)                                                   ccAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAAGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCa                                                                                                                                                                                            <  1:161715/251‑1 (MQ=255)                                                     agcagcGCGGGCAAGTTGCGCCCCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAAGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATc                                                                                                                                                                                          >  2:111701/1‑251 (MQ=255)                                                           gcgGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAAGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATc                                                                                                                                                                                    <  2:191719/251‑1 (MQ=255)                                                                                              aaaCACCAGACTGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAACGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAAGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCCATTGCTCCTTCCg                                                                                                                                                 >  2:39293/1‑251 (MQ=255)                                                                                                acacCAGAATGGCCATCAGCACCGCACAGAGGGAAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAAGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCggg                                                                                                                                               >  2:105342/1‑251 (MQ=255)                                                                                                  acCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACGGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAAGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCt                                                                                                                                             >  1:157804/1‑251 (MQ=255)                                                                                                              aTCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAAGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGatgatg                                                                                                                                 <  1:146648/251‑1 (MQ=255)                                                                                                                cAGCACCGCGCCGAGGTTAAGGTGCTCCAGACTTTCAACGCCCTGAACCGTGGAGGTTACATCCCACATCGAAATATTTTCGAGGAAGCGGAAAGCGGAATGGATTTCAGACCACAGCACAAAGACAGAAAGCAAGGCGATCAACAGCAAAATTGAGCGCCCCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACGTCAATTGCGCCTTCCGGGCTACTGTGAGGATgc                                                                                                                                <  2:50917/250‑1 (MQ=255)                                                                                                                             aGGGAAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGCACTATTTTCGAGGAAGCCGACAGCGGAAAGGAATTCTGACCACAGCACAATGACAGAAAGCAAGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGATTCCTCTTCGcc                                                                                                                  >  1:161489/1‑251 (MQ=255)                                                                                                                                     tGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAAGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTCCTCTTCGCCACGCGCa                                                                                                           <  1:116008/250‑1 (MQ=255)                                                                                                                                                                  tGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAAGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCa                                                                                                                                                                       <  1:4831/161‑1 (MQ=255)                                                                                                                                                                  tGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAAGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCa                                                                                                                                                                       >  2:4831/1‑161 (MQ=255)                                                                                                                                                                     aGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAAGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTCCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTc                                                                          >  2:98413/1‑251 (MQ=255)                                                                                                                                                                         gACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAAGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTCCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCGccc                                                                      <  1:91074/251‑1 (MQ=255)                                                                                                                                                                              cccACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAAGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTCCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCGCCCGATc                                                                  <  1:231838/250‑1 (MQ=255)                                                                                                                                                                                        aaTATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAAGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTCCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCGCCCGATCAAACGCCAGcc                                                       >  1:57633/1‑251 (MQ=255)                                                                                                                                                                                          tatTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAAGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTCCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCGCCCGATCAAACGCCAGCCTg                                                     >  2:243288/1‑251 (MQ=255)                                                                                                                                                                                             tttCGAGCAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAAGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTCCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGCTGCTTCGCCCGATCAAACGCCACACTgcgg                                                  >  2:249284/1‑251 (MQ=255)                                                                                                                                                                                                  aGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAAGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTCCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCGCCCGATCAAACGCCAGCCTGCGGCGCTg                                             >  1:20834/1‑251 (MQ=255)                                                                                                                                                                                                                   aTGGATTTCTGACCACAGCACAATGACAGAAAGCAAGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTCCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCGCCCGATCAAACGCCAGCCTGCGGCGCTGGATCAGCATCCAgcggc                            >  2:172650/1‑251 (MQ=255)                                                                                                                                                                                                                   aTGGATTTCTGACCACAGCACAATGACAGAAAGCAAGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTCCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCGCCCGATCAAACGCCAGCCTGCGGCGCTGGATCAGCATCCAgcggc                            <  1:86184/251‑1 (MQ=255)                                                                                                                                                                                                                     ggATTTCTGACCACAGCACAATGACAGAAAGCAAGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGTGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTCCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCGCCCGATCAAACGCCAGCCTGCGGCGCTGGATCAGCATCTAgcggcg                           <  1:84933/250‑1 (MQ=255)                                                                                                                                                                                                                          tcttaCCACAGCACAATGACAGAAAGCAAGGCGTTCAACATCAAAATTGAGCCCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTCCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCGCCCGATCAAACGCCATCCTGCGGCGCTGGATCAGCATCCAGCGGCGGATaa                      <  2:246772/246‑1 (MQ=255)                                                                                                                                                                                                                                        aaTGACAGAAAGCAAGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTCCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCGCCCGATCAAACGCCAGCCTGCGGCGCTGGATCAGCATCCAGCGGCGGATAACGTGATAAaccacc        >  2:24478/1‑250 (MQ=255)                                                                                                                                                                                                                                         aTGACAGAAAGCAAGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTCCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCGCCCGATCAAACGCCAGCCTGCGGCGCTGGATCAGCATCCAGCGGCGGATAACGTGATAAAccacca       >  2:216040/1‑250 (MQ=255)                                                                                                                                                                                                                                             cAGAAAGCAAGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTCCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCGCCCGATCAAACGCCAGCCTGCGGCGCTGGATCAGCATCCAGCGGCGGATAACGTGATAAACCACCAGTAGc  >  2:151037/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  GGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTCCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCGCCCGATCAAACGCCAGCCTGCGGCGCTGGATCAGCATCCAGCGGCGGATAACGTGATAAACCACCAGTAGC  >  NC_000913/4386755‑4387241

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