Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,069,712	C→T	V32V (GTG→GTA)	mscS ←	small conductance mechanosensitive channel MscS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,069,712	0	C	T	100.0%	136.1 / NA	41	V32V (GTG→GTA)	mscS	small conductance mechanosensitive channel MscS
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/24);  total (17/24)

ATACCGTAACGGACTAATGCAGAAAGAAAATCAGCAACAGTGGCATCGATTTTACGGGAGATCATCAGGCGATTCACCGCGTTGGAAATCATCCGCGCGATAATCAAACCAACGATGATGATCGCGAGTGCCGCCACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTCCCGCGCCGTTTATGCTATCGACAACATTCAAATCTTCCATTCAATATTCCTTTTCAACTGACTCAAAATGGAGAAATACA  >  NC_000913/3069578‑3069847                                                                                                                                       |                                                                                                                                        ataCCGTAACGGACTCATGCAGTAACAAACTCAGCAACAGTGGCATCGATTTTACGGGAGATCATCAGGCGATTCACCGCGTTGGAAATCATCCGCGCGATAATCAAACCAACGATGATGATCGCGAGTGCCGCTACGAt                                                                                                                                    <  1:206285/140‑1 (MQ=255)       tAACGGACTAATGCAGAAAGAAAATCAGCAACAGTGGCATCGATTTTACGGGAGATCATCAGGCGATTCACCGCGTTGGAAATCATCCGCGCGATAATCAAACCAACGATGATGATCGCGAGTGCCGCTACGATGTTTAc                                                                                                                              >  2:692040/1‑140 (MQ=255)           ggACTAATGCAGAAAGAAAATCAGCAACAGTGGCATCGATTTTACGGGAGATCATCAGGCGATTCACCGCGTTGGAAATCATCCGCGCGATAATCAAACCAACGATGATGATCGCGAGTGCCGCTACGATGTTTACTGCa                                                                                                                          <  1:680466/140‑1 (MQ=255)             aCTAATGCAGAAAGAAAATCAGCAACAGTGGCATCGATTTTACGGGAGATCATCAGGCGATTCACCGCGTTGGAAATCATCCGCGCGATAATCAAACCAACGATGATGATCGCGAGTGCCGCTACGATGTTTACTGCATa                                                                                                                        >  2:509317/1‑140 (MQ=255)             aCTAATGCAGAAAGAAAATCAGCAACAGTGGCATCGATTTTACGGGAGATCATCAGGCGATTCACCGCGTTGGAAATCATCCGCGCGATAATCAAACCAAAGATGATGATCGCGAGGGCCGCCACGATGTTTACTGCATa                                                                                                                        >  2:293702/1‑140 (MQ=255)              cTAATGCAGAAAGAAAATCAGCAACAGTGGCATCGATTTTACGGGAGATCATCAGGCGATTCACCGCGTTGGAAATCATCCGCGCGATAATCAAACCAACGATGATGATCGCGAGTGCCGCTACGATGTTTACTGCATaa                                                                                                                       <  2:64985/140‑1 (MQ=255)                 aTGCAGAAAGAAAATCACCAACCGTGGCATCCATTATACAGGAGATCATCAGGCGATTCACCGCGTTGGAAATCATCCGCGCGATAATCAAACCAACGATGATGATCGCGAGTGCCGCTACGATGTTTACTGCATAACtt                                                                                                                    <  2:245156/140‑1 (MQ=255)                            aaaTCAGCAACAGTGGCATCGATTTTACGGGAGATCATCAGGCGATTCACCGCGTTGGAAATCATCCGCGCGATAATCAAACCAACGATGATGATCGCGAGTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAgcgc                                                                                                         <  2:130093/140‑1 (MQ=255)                                     aCAGTGGCATCGATTTTACGGGAGATCATCAGGCGATTCACCGCGTTGGAAATCATCCGCGCGATAATCAAACCAACGATGATGATCGCGAGTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGc                                                                                                >  2:138204/1‑140 (MQ=255)                                         tGGCATCGATTTTACGGGAGATCATCAGGCGATTCACCGCGTTGGAAATCATCCGCGCGATAATCAAACCAACGATGATGATCGCGAGTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTAcc                                                                                            <  1:609581/140‑1 (MQ=255)                                           gCATCGATTTTACGGGAGATCATCAGGCGATTCACCGCGTTGGAAATCATCCGCGCGATAATCAAACCAACGATGATGATCGCGAGTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTAccag                                                                                          >  2:373171/1‑140 (MQ=255)                                                 aTTTTACGGGAGATCATCAGGCGATTCACCGCGTTGGAAATCATCCGCGCGATAATCAAACCAACGATGATGATCGCGAGTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCt                                                                                    >  1:394606/1‑140 (MQ=255)                                                   tttACGGGAGATCATCAGGCGATTCACCGCGTTGGAAATCATCCGCGCGATAATCAAACCAACGATGATGATCGCGAGTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTc                                                                                  >  1:571078/1‑140 (MQ=255)                                                       cGGGAGATCATCAGGCGATTCACCGCGTTGGAAATCATCCGCGCGATAATCAAACCAACGATGATGATCGCGAGTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTCCcgc                                                                              <  2:91366/140‑1 (MQ=255)                                                        gggAGATCATCAGGCGATTCACCGCGTTGGAAATCATCCGCGCGATAATCAAACCAACGATGATGATCGCGAGTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTCCcgcg                                                                             <  2:571078/140‑1 (MQ=255)                                                            gaTCATCAGGCGATTCACCGCGTTGGAAATCATCCGCGCGATAATCAAACCAACGATGATGATCGCGAGTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTCCCGCGCCGt                                                                         >  1:53166/1‑140 (MQ=255)                                                                    ggCGATTCACCGCGTTGGAAATCATCCGCGCGATAATCGAACCAACGATGATGATCGCGAGTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTCCCGCGCCGTTTATGCTa                                                                 >  2:90483/1‑140 (MQ=255)                                                                        aTTCACCGCGTTGGAAATCATCCGCGCGATAATCAAACCAACGATGATGATCGCGAGTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTCCCGCGCCGTTTATGCTATCGa                                                             <  2:64941/140‑1 (MQ=255)                                                                        aTTCACCGCGTTGGAAATCATCCGCGCGATAATCAAACCAACGATGATGATCGCGAGTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTCCCGCGCCGTTTATGCTATCGa                                                             >  1:388463/1‑140 (MQ=255)                                                                         ttCACCGCGTTGGAAATCATCCGCGCGATAATCAAACCAACGATGATGATCGCGAGTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTCCCGCGCCGTTTATGCTATCGac                                                            <  1:138204/140‑1 (MQ=255)                                                                         ttCACCGCGTTGGAAATCATCCGCGCGATAATCAAACCAACGATGATGATCGCGAGTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTCCCGCGCCGTTTATGCTATCGac                                                            >  1:577025/1‑140 (MQ=255)                                                                          tCACCGCGTTGGAAATCATCCGCGCGATAATCAAACCAACGATGATGATCGCGAGTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTCCCGCGCCGTTTATGCTATCGaca                                                           <  1:419809/140‑1 (MQ=255)                                                                          tCACCGCGTTGGAAATCATCCGCGCGATAATCAAACCAACGATGATGATCGCGAGTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTCCCGCGCCGTTTATGCTATCGaca                                                           <  2:53166/140‑1 (MQ=255)                                                                             ccGCGTTGGAAATCATCCGCGCGATAATCAAACCAACGATGATGATCGCGAGTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTCCCGCGCCGTTTATGCTATCGacaaca                                                        <  1:692040/140‑1 (MQ=255)                                                                                  ttGGAAATCATCCGCGCGATAATCAAACCAACGATGATGATCGCGAGTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTCCCGCGCCGTTTATGCTATCGACAACATTCaa                                                   <  1:93642/140‑1 (MQ=255)                                                                                  ttGGAAATCATCCGCGCGATAATCAAACCAACGATGATGATCGCGAGTGCCACTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTCCCGCGCCGTTTATGCTATCGACAACATTCaa                                                   <  2:390326/140‑1 (MQ=255)                                                                                         tcatcCGCGCGATAATCGAACCAACGATGATGATCGCGAGTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTCCCGCGCCGTTTATGCTATCGACAACATTCAAATCTTcc                                            <  1:90483/140‑1 (MQ=255)                                                                                            tcCGCGCGATAATCAAACCAACGATGATGATCGCGAGTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTCCCGCGCCGTTTATGCTATCGACAACATTCAAATCtt                                              <  1:111219/135‑1 (MQ=255)                                                                                            tcCGCGCGATAATCAAACCAACGATGATGATCGCGAGTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTCCCGCGCCGTTTATGCTATCGACAACATTCAAATCtt                                              >  2:111219/1‑135 (MQ=255)                                                                                             ccGCGCGATAATCAAACCAACGATGATGATCGCGAGTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTCCCGCGCCGTTTATGCTATCGACAACATTCAAATCTTCCATTc                                        >  1:475299/1‑140 (MQ=255)                                                                                                cgcgATAATCAAACCAACGATGATGATCGCGAGTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTCCCGCGCCTTTTATGCTATCTACTACATTAAAATCTTCCATTCAat                                     >  1:323099/1‑140 (MQ=255)                                                                                                  cgATAATCAAACCAACGATGATGATCGCGAGTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTCCCGCGCCGTTTATGCTATCGACAACATTCAAATCTTCCATTCAatat                                   <  2:577025/140‑1 (MQ=255)                                                                                                  cgATAATCAAACCAACGATGATGATCGCGAGTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTCCCGCGCCGTTTATGCTATCGACAACATTCAAATCTTCCATTCAatat                                   <  1:373171/140‑1 (MQ=255)                                                                                                  cgATAATCAAACCAACGATGATGATCGCGAGTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTCCCGCGCCGTTTATGCTATCGACAACATTCAAATCTTCCATTCAatat                                   <  2:475299/140‑1 (MQ=255)                                                                                                             ccAACGATGATGATCGCGAGTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTCCCGCGCCGTTTATGCTATCGACAACATTCAAATCTTCCATTCAATATTCCTTTTCAAc                        >  1:672342/1‑140 (MQ=255)                                                                                                             ccAACGATAATGATCGCGAGTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTCCCGCGCCGTTTATGCTATCGACAACATTCAAATCTTCCATTCAATATTCCTTTTCAAc                        <  2:560203/140‑1 (MQ=255)                                                                                                                    tgatgatCGCGAGTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTCCCGCGCCGTTTATGCTATCGACAACATTCAAATCTTCCATTCAATATTCCTTTTCAACTGACTCa                 <  2:622148/140‑1 (MQ=255)                                                                                                                         atCGCGAGTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTCCCGCGCCGTTTATGCTATCGACAACATTCAAATCTTCCATTCAATATTCCTTTTCAACTGACTCAAAATg            <  1:509317/140‑1 (MQ=255)                                                                                                                         atCGCGAGTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTCCCGCGCCGTTTATGCTATCGACAACATTCAAATCTTCCATTCAATATTCCTTTTCAACTGACTCAAAATg            <  2:388463/140‑1 (MQ=255)                                                                                                                            gcgAGTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTCCCGCGCCGTTTATGCTATCGACAACATTCAAATCTTCCATTCAATATTCCTTTTCAACTGACTCAAAATGgag         >  1:515423/1‑140 (MQ=255)                                                                                                                            gcgAGTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTCCCGCGCCGTTTATGCTATCGACAACATTCAAATCTTCCATTCAATATTCCTTTTCAACTGACTCAAAATGgag         <  1:545410/140‑1 (MQ=255)                                                                                                                               aGTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTCCCGCGCCGTTTATGCTATCGACAACATTCAAATCTTCCATTCAATATTCCTTTTCAACTGACTCAAAATGGAGaaa      >  1:532890/1‑140 (MQ=255)                                                                                                                                gTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTCCCGCGCCGTTTATGCTATCGACAACATTCAAATCTTCCATTCAATATTCCTTTTCAACTGACTCAAAATGgag         >  1:682706/1‑136 (MQ=255)                                                                                                                                gTGCCGCTACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTCCCGCGCCGTTTATGCTATCGACAACATTCAAATCTTCCATTCAATATTCCTTTTCAACTGACTCAAAATGgag         <  2:682706/136‑1 (MQ=255)                                                                                                                                   ccgctACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTCCCGCGCCGTTTATGCTATCGACAACATTCAAATCTTCCATTCAATATTCCTTTTCAACTGACTCAAAATGGAGAAATaca  >  2:665978/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                        ATACCGTAACGGACTAATGCAGAAAGAAAATCAGCAACAGTGGCATCGATTTTACGGGAGATCATCAGGCGATTCACCGCGTTGGAAATCATCCGCGCGATAATCAAACCAACGATGATGATCGCGAGTGCCGCCACGATGTTTACTGCATAACTTAGCAGCAGCGCCTGGTTAGCTACCAGCCAGCTTCCCGCGCCGTTTATGCTATCGACAACATTCAAATCTTCCATTCAATATTCCTTTTCAACTGACTCAAAATGGAGAAATACA  >  NC_000913/3069578‑3069847

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