Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,377,323	T→C	D622G (GAC→GGC)	mscM ←	miniconductance mechanosensitive channel MscM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,377,323	0	T	C	100.0%	81.1 / NA	23	D622G (GAC→GGC)	mscM	miniconductance mechanosensitive channel MscM
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (9/14);  total (9/14)

CCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGGTCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGATGCTCATCAGGTAGTAGCGCATCCCACGGGA  >  NZ_CP009273/4377237‑4377408                                                                                       |                                                                                      ccTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGgccgtc                                                                                   <  2:84571/91‑1 (MQ=255)          ggCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATc                                                                          <  1:515010/91‑1 (MQ=255)            cTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGa                                                                        <  2:161474/91‑1 (MQ=255)                      ccaccGCCAGCGCACCGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAg                                                              <  2:235221/91‑1 (MQ=255)                         ccGCCAGCGCACCGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAgcgc                                                           <  1:268578/91‑1 (MQ=255)                          cGCCAGCCCACCGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGcc                                                          >  1:1046376/1‑91 (MQ=255)                           gCCAGCGCACCGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCa                                                         <  1:64725/91‑1 (MQ=255)                                      gCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGc                                              <  2:545016/91‑1 (MQ=255)                                            gAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATa                                        <  2:1046376/91‑1 (MQ=255)                                                   aagcaaagcCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGa                                 >  2:39926/1‑91 (MQ=255)                                                   aagcaaagcCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGa                                 <  1:165856/91‑1 (MQ=255)                                                    agcaaagcCGCCCCAGCGAACCGGAGAATTCGCGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGAt                                <  2:31318/91‑1 (MQ=255)                                                       aaagcCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGATGCt                             <  1:39926/91‑1 (MQ=255)                                                        aagcCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCTATAATCAGCGGCACAATAAGCCCGATGCtc                            >  2:288428/1‑91 (MQ=255)                                                        aagcCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGATGCtc                            >  1:277264/1‑91 (MQ=255)                                                           cggTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGATGCtcatc                         <  1:425575/89‑1 (MQ=255)                                                            cGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGATGCtcatca                        <  2:323990/91‑1 (MQ=255)                                                             gTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGATGCTCATCAg                       >  1:141896/1‑91 (MQ=255)                                                             gTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGATGCTCATCAg                       >  1:788877/1‑91 (MQ=255)                                                                  aGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGATGCTCATCAGgtagt                  >  2:991251/1‑91 (MQ=255)                                                                    cGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGATGCTCATCAGGtagtag                <  1:739434/91‑1 (MQ=255)                                                                      aaCCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGATGCTCATCAGGTAGTAgcg              >  2:56334/1‑91 (MQ=255)                                                                          ggAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGATGCTCATCAGGTAGTAGCGCATc          >  1:339800/1‑91 (MQ=255)                                                                              aaTTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGATGCTCATCAGGTAGTAGCGCATCCCAc      <  2:339800/91‑1 (MQ=255)                                                                                  cACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGATGCTCATCAGGTAGTAGCGCATCCCACGGGa  <  2:925925/91‑1 (MQ=255)                                                                                       |                                                                                      CCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGGTCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGATGCTCATCAGGTAGTAGCGCATCCCACGGGA  >  NZ_CP009273/4377237‑4377408

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