Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,389,355	T→A	T6S (ACT→TCT)	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,355	0	T	A	100.0%	114.3 / NA	30	T6S (ACT→TCT)	mscM	miniconductance mechanosensitive channel MscM
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (18/12);  total (18/12)

CTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCG  >  NC_000913/4389281‑4389431                                                                           |                                                                             cTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGAGATAAt                                                                        <  1:829345/81‑1 (MQ=255)    aGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGAGATAATc                                                                       <  1:829346/79‑1 (MQ=255)     gTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGAGATAATCAgg                                                                    >  1:829396/1‑81 (MQ=255)        aTTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGAGATAATCAGgcg                                                                  >  2:829380/1‑80 (MQ=255)        aTTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGAGATAATCAGgcg                                                                  >  1:829382/1‑80 (MQ=255)                gCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGAGATAATCAGGCGCACGTCAGc                                                         >  2:829401/1‑81 (MQ=255)                 cTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGAGATAATCAGGCGCACGTCAGCg                                                        >  1:829397/1‑81 (MQ=255)                      gggggCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGAGATAATCAGGCGCACGTCAGCGTTTcc                                                   >  1:829391/1‑81 (MQ=255)                       ggggCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGAGATAATCAGGCGCACGTCAGCGTTTCCt                                                  >  1:829390/1‑81 (MQ=255)                        gggCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGAGATAATCAGGCGCACGTCAGCGTTTCCtt                                                 >  2:829386/1‑81 (MQ=255)                        gggCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGAGATAATCAGGCGCACGTCAGCGTTTCCtt                                                 <  2:829347/81‑1 (MQ=255)                         ggCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGAGATAATCAGGCGCACGTCAGCGTTTCCttt                                                >  1:829392/1‑81 (MQ=255)                              tCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGAGATAATCAGGCGCACGTCAGCGTTTCCTTTGATgg                                           <  1:829348/81‑1 (MQ=255)                                gcggcgTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGAGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGAt                                         >  1:829389/1‑81 (MQ=255)                                     gTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGAGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGAc                                    >  2:829399/1‑81 (MQ=255)                                      tACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGAGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGAcc                                   >  1:829385/1‑81 (MQ=255)                                           cccccAACTGAGGCACCAGGCCATCAGAAAAGAGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTc                              >  1:829406/1‑81 (MQ=255)                                            ccccAACTGAGGCACCAGGCCATCAGAAAAGAGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCt                             <  2:829349/81‑1 (MQ=255)                                                 aCTGAGGCACCAGGCCATCAGAAAAGAGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCtttttt                        >  2:829398/1‑81 (MQ=255)                                                   tGAGGCACCAGGCCATCAGAAAAGAGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTgt                      >  1:829400/1‑81 (MQ=255)                                                     aGGCACCAGGCCATCAGAAAAGAGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTgtcg                    <  1:829352/81‑1 (MQ=255)                                                     aGGCACCAGGCCATCAGAAAAGAGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTgtc                     <  2:829351/80‑1 (MQ=255)                                                     aGGCACCAGGCCATCAGAAAAGAGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTgt                      <  1:829350/79‑1 (MQ=255)                                                           cAGGCCATCAGAAAAGAGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAAc               >  2:829393/1‑80 (MQ=255)                                                            aGGCCATCAGAAAAGAGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCa             <  2:829353/81‑1 (MQ=255)                                                                   cAGAAAAGAGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCt      >  2:829415/1‑81 (MQ=255)                                                                    aGAAAAGAGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCt      <  1:829354/80‑1 (MQ=255)                                                                       aaaGAGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGcg  <  2:829355/81‑1 (MQ=255)                                                                       aaaGAGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGcg  >  1:829402/1‑81 (MQ=255)                                                                       aaaGAGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGcg  <  1:829356/81‑1 (MQ=255)                                                                           |                                                                             CTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCG  >  NC_000913/4389281‑4389431

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