Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	485,143	A→G	D1051G (GAT→GGT)	mscK →	mechanosensitive channel MscK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	485,143	0	A	G	100.0%	97.3 / NA	28	D1051G (GAT→GGT)	mscK	mechanosensitive channel MscK
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/16);  total (12/16)

TCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATCAACATTGCCTTTAA  >  NZ_CP009273/485009‑485257                                                                                                                                       |                                                                                                                   tctcGGCGCGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCAt                                                                                                                >  2:96087/1‑139 (MQ=255)    cGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGTGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAg                                                                                                             <  1:205458/139‑1 (MQ=255)      gCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCt                                                                                                           <  2:431755/139‑1 (MQ=255)                tGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTa                                                                                                  <  1:417469/138‑1 (MQ=255)                       gATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGCACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCCTGAGCTGCGTCTGTATGTGCGt                                                                                           >  1:344463/1‑138 (MQ=255)                        aTCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTg                                                                                          >  1:21391/1‑138 (MQ=255)                              aaaaaGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCg                                                                                   <  1:247336/139‑1 (MQ=255)                                     gCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCg                                                                             >  2:232256/1‑138 (MQ=255)                                            gtgtTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTa                                                                     <  2:388245/139‑1 (MQ=255)                                             tgtTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTAc                                                                    <  2:177208/139‑1 (MQ=255)                                                tACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGt                                                                 >  2:139383/1‑139 (MQ=255)                                                       ggcggcGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAg                                                          >  1:191887/1‑139 (MQ=255)                                                       ggcggcGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCGGCGTAGGGAGGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAg                                                          <  1:1651/139‑1 (MQ=255)                                                            cgACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGaa                                                     <  1:96087/139‑1 (MQ=255)                                                            cgACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGaa                                                     >  2:175378/1‑139 (MQ=255)                                                            cgACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGaa                                                     <  1:139383/139‑1 (MQ=255)                                                                  aGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTa                                                <  2:413429/138‑1 (MQ=255)                                                                  aGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTAc                                               <  2:191887/139‑1 (MQ=255)                                                                  aGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTAc                                               <  1:175378/139‑1 (MQ=255)                                                                      cccAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATc                                           >  2:146542/1‑139 (MQ=255)                                                                         aaGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGGATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGag                                        <  1:338/139‑2 (MQ=255)                                                                            ggTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAg                                     <  1:74556/139‑1 (MQ=255)                                                                                 tGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCtgtg                                >  2:192227/1‑139 (MQ=255)                                                                                     cGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGt                            >  2:120150/1‑139 (MQ=255)                                                                                                  gAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATc                >  2:20057/1‑138 (MQ=255)                                                                                                  gAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGCAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGATGTGCCGTGAAAACGACAta                >  2:49158/1‑137 (MQ=255)                                                                                                      tcttctTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATCAACAt           <  1:257350/139‑1 (MQ=255)                                                                                                               cGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATCAACATTGCCTTTaa  <  2:316378/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                   TCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATCAACATTGCCTTTAA  >  NZ_CP009273/485009‑485257

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