Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	109,960	G→A	R561H (CGT→CAT)	secA →	protein translocation ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	109,960	0	G	A	100.0%	92.9 / NA	31	R561H (CGT→CAT)	secA	protein translocation ATPase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/16);  total (15/16)

CCAACGAAGCGGCGATTGTTGCTCAGGCAGGTTATCCGGCTGCGGTGACTATCGCGACCAATATGGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGT‑ACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCGTCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACTAAAGCGATTGCCAACGCCCAGCGTAAAGTTGAAAGCCGTAACTTCGACATTCG  >  NC_000913/109732‑110203                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    ccAACGATGCGGCGATTGTTGCTCAGGTAGGTTATCCTGCTGCGGTGACTATCGCGACCAATATGGCTGGGCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGT‑ACTGGAAGCAGGTGGCCTGCATTTCATCGGTCCCGAGCATCACGAATCCCGTCGTATCGAt                                                                                                                                                                                                                               <  2:151609/251‑1 (MQ=255)                    tgCTCAGGCAGGTTATCCGGCTGCGGTGACTATCGCGACCAATATGGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCTCGATGCGGT‑ACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGATCATCAGGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGt                                                                                                                                                                                                             <  1:36625/250‑1 (MQ=255)                            caggTTATCCGGCGTCGTTGACTATCGCGACCAAGAGGGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGTTACGTCACGATGCGGT‑ACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTTGTATCGATATCCAGTTGCGCGGTCGTTCTGgtcgt                                                                                                                                                                                                    <  2:291895/251‑1 (MQ=255)                             aggTTATCCGGCTGCGGTGACTATCGCGTCCAATATGGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGT‑ACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGgtcgtc                                                                                                                                                                                                   <  2:311729/251‑1 (MQ=255)                                ttATCCGGCTGCGGTGACTATCGCGACCAATATGGCGGGTCGAGGGACAGAAAATGTGGTCGGGGGGAGCTGGGAGGGAGAAGTTGCCGCGCTGGAAAAACCGACCGCAGAGCAAATTGAAAAAAATAAAGCCGCATGGCAGGGACGGCCCCAATCGGTCAC‑GGGAGAAGGGGGGCTGCATATCATCGGGAACGAGCATCACGAATCACGTCGTGTCGGTAAACAGATGAGCggggggtttgggcgggcgg                                                                                                                                                                                                >  1:281566/1‑234 (MQ=255)                                            ggTGACTATCGCGACCAAGATGGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGT‑ACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGtt                                                                                                                                                                                    <  2:2604/251‑1 (MQ=255)                                                  tATCGCGACCAATATGGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGT‑ACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGt                                                                                                                                                                                     <  1:436664/244‑1 (MQ=255)                                                            aaTATGGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGT‑ACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCt                                                                                                                                                                    <  1:11245/251‑1 (MQ=255)                                                                  gCGGGTCGTGGGACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGT‑ACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGGCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGAt                                                                                                                                                              <  2:335612/251‑1 (MQ=255)                                                                     ggTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGT‑ACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGa                                                                                                                                                           <  1:113985/251‑1 (MQ=255)                                                                           ggTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGT‑ACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGGATCCCGTCGTATCGATAACCAGTTGCGCGGGCGTTCTGGTCGTCCGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGACGATgg                                                                                                                                                     >  1:227039/1‑250 (MQ=255)                                                                                aGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGT‑ACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTgg                                                                                                                                                >  1:1557/1‑250 (MQ=255)                                                                                            cGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGT‑ACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTg                                                                                                                                    >  1:254101/1‑251 (MQ=255)                                                                                                   aGCTGGCAGGCAGCAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCCAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGT‑ACTGGGAGCCGGTGGCCTGCATATCATCGGTACCGGGCATCACGAATCCCGGCGGATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGGTGGTTCTTCCCCTTTCTACCTGTCGCTGGAAGATGCGCTGAGGCGTCTTTTTCCTTccg                                                                                                                              >  2:406352/1‑250 (MQ=255)                                                                                                                                    ccGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGT‑ACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACTAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTa                                                                                             <  1:263612/250‑1 (MQ=255)                                                                                                                                     cGACCGCAGAGAAAATTGAAAACATTACAGCCGACTGGCACGACCGTCACGCTGGCGT‑ACTGGAAGAAGGTGGCCCGCAAATCCCTGGCACCGAGGATCAACAATCCAGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGCCGCCCGGGGGCTGCTGGTTATTCCCCTTTATACTTCTCGGAGGAAGATGCGCTGAAGCGTATTTTTGCTTCCGACCGCGTAGACGCCCTGATGCGTACACTGGGttg                                                                                            >  2:386635/1‑248 (MQ=255)                                                                                                                                                cAAATTGAAAAAATTAAAGCCGACTGGCAGCTACGTCACGATGCGGA‑ACAGGCACCAGGTGGCCTCCACCTCATCGGTACCGAGCACCACGAATCCCGCCGTATCGATCCCCAGTTACGCCATCACCCTGATCGTCCGGGGGATGCTGGTTCCTCCCGTTTCTACCTGTCGATAGAAGATGCGCTGATGCGTATTTTTGCTTCCGCCCGAGTATACGGCATGATCCGTACAATGGGTATGCAGCCAGGCg                                                                                 >  1:414008/1‑250 (MQ=255)                                                                                                                                                           aaTTAAAGCCGACTGGCAGGTACGTCTCGATGCGGT‑ACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAAGCCCGTCGTATCGATAACCAGTTGCGCGGTCGTGCTGGTCTTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGTAAGATTCGCTGATGCGTATTTTTGCTTCCGACCGAGTAGCCTGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAac                                                                     <  1:151850/251‑1 (MQ=255)                                                                                                                                                              tAAAGCTGACTGGCAGGTACGTCACGATGCGGT‑ACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACAcc                                                                  >  1:118509/1‑251 (MQ=255)                                                                                                                                                                   ccGACTGGCAGGTACGTCACGATGCGGT‑ACTGGAAGCAGGTGGCCTGCATATCCTCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCgggg                                                             >  1:26535/1‑248 (MQ=255)                                                                                                                                                                    cGACTGGCAGGGACGTCACGATGCGGT‑ACTGGAAGCAGGGGGCCTGCATATCATCGGTACCGAGGATCACGAATCCCGTCGTAACGATAACCCGTTGCGCGGTCGTTCTGGTCGTCAGGGGGGTGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGGGCTGATGCGTATTTTTGCTTCCGCCCGAGTATCCGGCATGGTGCGTAAACTGGGTATGAAGCCAGGCGGAGCCATTGAACCCCCGgggg                                                            >  2:357023/1‑251 (MQ=255)                                                                                                                                                                          gCAGGGACGTCACGAGGCGGT‑ACGGGAAGCAGGTGGCCTGAATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGGCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGCCCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACTa                                                      <  2:12996/251‑1 (MQ=255)                                                                                                                                                                              gTACGTCCCGATGCGGT‑GCTGGAAGCAGGTGGCCTGCATATCATCTGTACCGCGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGGCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTCTCCGGCATGATGCGTAAACTGGGTATGTAGCCAGGCGAAGCCATTGAACCCCCGTGGGTGACTAAAGc                                                  <  1:222091/251‑1 (MQ=255)                                                                                                                                                                                aCGTCACGATGCGGT‑ACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGGCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCTATTGATCACCCGTGGGTGATTAAAGCg                                                 <  2:282801/250‑1 (MQ=255)                                                                                                                                                                                            ggT‑ACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACTAAAGCGATTGCCAACGccc                                    >  1:356469/1‑251 (MQ=255)                                                                                                                                                                                                   ggAAGCAGGAGGCCTGCATATCATCGGTACAGAACATAACGAATCACGTCGCATCGATCACCCGTTCCGCGGTCGATCTGGTCGTAAGGGGAATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGCGTATCCGGCATGATGCGTAAACTGGATATGAAGCAAGGCGAACCCATTGAACACCTGGGGGTGACTAAAGCAATTGCCAAAGACCAGCGt                               >  2:278123/1‑250 (MQ=255)                                                                                                                                                                                                    gAAGCCGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACTAAAGCGATTGCCAACGCCCAGCGTa                              >  1:96173/1‑250 (MQ=255)                                                                                                                                                                                                       gCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACTAAAGCGATTGCCAACGCCCAGCGTAAAGt                          <  1:37391/251‑1 (MQ=255)                                                                                                                                                                                                                 tGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACTAAAGCGATTGCCAACGCCCAGCGTAAAGTTGAAAGCCGt                >  1:350974/1‑251 (MQ=255)                                                                                                                                                                                                                           tggTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGACGTCCGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACTAAAGCGATTGCAAACGCCCAGCGTAAAGTTGAAAGCCGTAACTTCGAc       >  2:89125/2‑250 (MQ=255)                                                                                                                                                                                                                               aCCGAGCAGCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAATCTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACTAAAGCGATTGCCAACGTCCAGCGTAAAGTTGAAAGCCGTAACTTCGTCATTCg  <  2:355680/251‑1 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    CCAACGAAGCGGCGATTGTTGCTCAGGCAGGTTATCCGGCTGCGGTGACTATCGCGACCAATATGGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGT‑ACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCGTCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACTAAAGCGATTGCCAACGCCCAGCGTAAAGTTGAAAGCCGTAACTTCGACATTCG  >  NC_000913/109732‑110203

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