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%	93.2 / NA	27	D622G (GAC→GGC)	mscM	miniconductance mechanosensitive channel MscM
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/12);  total (15/12)

GTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGGTCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGATGCTCATCAGGTAGTAGCGCATCCCACGGGAAACACGTTCGCGCGGCCAGCCAAAATGAGCGATAAAC  >  NZ_CP009273/4377188‑4377445                                                                                                                                        |                                                                                                                           gtAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGgccgt                                                                                                                         >  1:100885/1‑139 (MQ=255) gtAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGgccgt                                                                                                                         >  2:688977/1‑139 (MQ=255) gtAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGgccgt                                                                                                                         >  2:130805/1‑139 (MQ=255)        tgtCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGt                                                                                                                   >  1:483929/1‑138 (MQ=255)          tCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTa                                                                                                                 <  2:783520/138‑1 (MQ=255)            gccgcTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATc                                                                                                               >  2:29434/1‑138 (MQ=255)               gcTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAAc                                                                                                           <  2:351483/139‑1 (MQ=255)               gcTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAAc                                                                                                           <  1:51486/139‑1 (MQ=255)                                 aTACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAAt                                                                                          >  2:132381/1‑138 (MQ=255)                                      gCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCgg                                                                                    <  1:130805/139‑1 (MQ=255)                                      gCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCgg                                                                                    <  1:393775/139‑1 (MQ=255)                                         gAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGca                                                                                  >  1:258924/1‑138 (MQ=255)                                                 gCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCgg                                                                                    <  1:108691/128‑1 (MQ=255)                                                 gCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCgg                                                                                    >  2:108691/1‑128 (MQ=255)                                                         cAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGATGCt                                                                  <  2:61528/138‑1 (MQ=255)                                                             cTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGATGCtcatca                                                             <  1:29434/139‑1 (MQ=255)                                                                aGGGTGACCACCGCCAGCGCACCGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGATGCTCATCAgg                                                           >  1:723625/1‑138 (MQ=255)                                                                aGGGTGACCACCGCCAGCGCACCGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGATGCTCATCAgg                                                           >  2:348946/1‑138 (MQ=255)                                                                      accaccGCCAGCGCACCGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGATGCtca                                                                >  2:157914/1‑127 (MQ=255)                                                                      accaccGCCAGCGCACCGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGATGCtca                                                                <  1:157914/127‑1 (MQ=255)                                                                            gCCAGCGCACCGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGATGCTCATCAGGTAGTAGCGCATcc                                              >  1:694468/1‑139 (MQ=255)                                                                                     ccGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGATGCTCATCAGgta                                                         >  1:503396/1‑119 (MQ=255)                                                                                     ccGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGATGCTCATCAGgta                                                         <  2:503396/119‑1 (MQ=255)                                                                                     ccGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGATGCTCATCAGGTAGTAGCGCATccc                                             <  1:64610/131‑1 (MQ=255)                                                                                     ccGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGATGCTCATCAGGTAGTAGCGCATccc                                             >  2:64610/1‑131 (MQ=255)                                                                                                             cGTCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGATGCTCATCAGGTAGTAGCGCATCCCACGGGAAACACGTTCGCGCGGCCAGCCAAAAt             <  1:786642/139‑1 (MQ=255)                                                                                                               tCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGATGCTCATCAGGTAGTAGCGCATCCCACGGGAAACACGTTCGCGCGGCCAGCCAAAAt             >  1:857647/1‑137 (MQ=255)                                                                                                               tCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGATGCTCATCAGGTAGTAGCGCATCCCACGGGAAACACGTTCGCGCGGCCAGCCAAAAt             >  1:686837/1‑137 (MQ=255)                                                                                                               tCCCAGCGAACCGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGATGCTCATCAGGTAGTAGCGCATCCCACGGGAAACACGTTCGCGCGGCCAGCCAAAAt             >  2:636022/1‑137 (MQ=255)                                                                                                                         ccGGAGAATTCACGGCCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGATGCTCATCAGGTAGTAGCGCATCCCACGGGAAACACGTTCGCGCGGCCAGCCAAAATGAGCGATAAAc  >  2:112358/1‑138 (MQ=255)                                                                                                                                        |                                                                                                                           GTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATAAAGCAAAGCCGTCCCAGCGAACCGGAGAATTCACGGTCGTCGAGGTTATCGAACATCATCAGCGCCATAATCAGCGGCACAATAAGCCCGATGCTCATCAGGTAGTAGCGCATCCCACGGGAAACACGTTCGCGCGGCCAGCCAAAATGAGCGATAAAC  >  NZ_CP009273/4377188‑4377445

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