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%	154.2 / NA	41	D1051G (GAT→GGT)	mscK	mechanosensitive channel MscK
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/28);  total (13/28)

AAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCC  >  NZ_CP009273/485061‑485229                                                                                   |                                                                                       aaGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGc                                                                                <  1:751083/91‑1 (MQ=255) aaGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGc                                                                                <  2:488590/91‑1 (MQ=255)  aGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCt                                                                               >  2:540908/1‑91 (MQ=255)  aGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCt                                                                               >  1:535566/1‑91 (MQ=255)  aGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCt                                                                               >  2:779975/1‑91 (MQ=255)   ggcggcGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTg                                                                              <  2:209712/91‑1 (MQ=255)     cggcgACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCg                                                                            <  1:328628/91‑1 (MQ=255)     cggcgACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCg                                                                            <  1:391259/91‑1 (MQ=255)            tGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTAt                                                                     <  1:94416/91‑1 (MQ=255)            tGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTAt                                                                     <  1:674879/91‑1 (MQ=255)                    cAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGa                                                             <  2:383060/91‑1 (MQ=255)                    cAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGa                                                             <  1:779975/91‑1 (MQ=255)                           gATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGt                                                      <  2:349198/91‑1 (MQ=255)                           gATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGt                                                      >  2:420251/1‑91 (MQ=255)                           gATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGt                                                      <  2:127131/91‑1 (MQ=255)                            aTGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCAt                                                                                    <  2:789356/60‑1 (MQ=255)                            aTGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCAt                                                                                    >  1:789356/1‑60 (MQ=255)                            aTGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTg                                                     >  1:188946/1‑91 (MQ=255)                            aTGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTg                                                     >  2:870125/1‑91 (MQ=255)                                aCGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCg                                                 >  1:699046/1‑91 (MQ=255)                                  gAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTa                                               <  2:645010/91‑1 (MQ=255)                                        aTGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTa                                         >  1:372319/1‑91 (MQ=255)                                         tgtcGGAAGTATTCTCTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTTCGTCTGTATGTGGGTGAACTGCGTGACCGTAGTCGTAc                                        <  2:403265/88‑1 (MQ=255)                                          gCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACt                                       <  2:659608/91‑1 (MQ=255)                                               aaGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGa                                  <  1:485939/91‑1 (MQ=255)                                                    ttcttTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGc                             <  2:593949/91‑1 (MQ=255)                                                    ttcttTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGc                             <  1:540908/91‑1 (MQ=255)                                                     tcttTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCt                            >  1:149568/1‑91 (MQ=255)                                                      cttTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTg                           <  1:383558/91‑1 (MQ=255)                                                         tACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAAc                        <  1:870125/91‑1 (MQ=255)                                                          aCGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAAc                        >  2:244680/1‑90 (MQ=255)                                                          aCGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAAc                        <  1:244680/90‑1 (MQ=255)                                                            ggCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGt                     <  2:382640/91‑1 (MQ=255)                                                               aTTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACt                  <  2:302471/91‑1 (MQ=255)                                                               aTTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACt                  >  1:518441/1‑91 (MQ=255)                                                               aTTTGGCGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACt                  <  1:846871/91‑1 (MQ=255)                                                                    gTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGa             <  2:699046/91‑1 (MQ=255)                                                                    gTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGa             <  2:822444/91‑1 (MQ=255)                                                                      gCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATc           <  1:632589/91‑1 (MQ=255)                                                                              gTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGc   <  1:495981/91‑1 (MQ=255)                                                                               ttGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGGACTATCGAACAGGtgtggc  >  1:153524/1‑89 (MQ=255)                                                                                   |                                                                                       AAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCC  >  NZ_CP009273/485061‑485229

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