Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	396,645	A→T	K3* (AAG→TAG)	sbmA →	peptide antibiotic/peptide nucleic acid transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	396,645	0	A	T	100.0%	126.5 / NA	33	K3* (AAG→TAG)	sbmA	peptide antibiotic/peptide nucleic acid transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (18/15);  total (18/15)

TGACGCTTTGTAGCGGTCATGCGGTTAATACACAGCAGTTGTTAAAACGATAAGAAGTTAGCAGGAGTGCATATGTTTAAGTCTTTTTTCCCAAAGCCGGGAACGTTTTTTCTCTCGGCCTTCGTTTGGGCATTGATTGCCGTTATCTTCTGGCA  >  NC_000913/396567‑396721                                                                               |                                                                             tGACGCTTTGTAGCGGTCATGCGGTTAATACACAGCAGTTGTTAAAACGATAAGAAGTTAGCAGGAGTGCATATGTTTTAg                                                                            >  1:81558/1‑81 (MQ=255)   aCGCTTTGTAGCGGTCATGCGGTTAATACACAGCAGTTGTTAAAACGATAAGAAGTTAGCAGGAGTGCATATGTTTTAg                                                                            <  1:81526/79‑1 (MQ=255)         tGTAGCGGTCATGCGGTTAATACACAGCAGTTGTTAAAACGATAAGAAGTTAGCAGGAGTGCATATGTTTTAGTCtttt                                                                      <  1:81527/79‑1 (MQ=255)         tGTAGCGGTCATGCGGTTAATACACAGCAGTTGTTAAAACGATAAGAAGTTAGCAGGAGTGCATATGTTTTAGTCtttt                                                                      <  2:81528/79‑1 (MQ=255)           tAGCGGTCATGCGGTTAATACACAGCAGTTGTTAAAACGATAAGAAGTTAGCAGGAGTGCATATGTTTTAGTCTTTTTTc                                                                   <  2:81529/80‑1 (MQ=255)             gCGGTCATGCGGTTAATACACAGCAGTTGTTAAAACGATAAGAAGTTAGCAGGAGTGCATATGTTTTAGTCTTTTTTccc                                                                 >  2:81552/1‑80 (MQ=255)                  cATGCGGTTAATACACAGCAGTTGTTAAAACGATAAGAAGTTAGCAGGAGTGCATATGTTTTAGTCTTTTTTCCCAAAGcc                                                           >  2:81560/1‑81 (MQ=255)                      cGGTTAATACACAGCAGTTGTTAAAACGATAAGAAGTTAGCAGGAGTGCATATGTTTTAGTCTTTTTTCCCAAAGCCggg                                                        >  1:81536/1‑80 (MQ=255)                      cGGTTAATACACAGCAGTTGTTAAAACGATAAGAAGTTAGCAGGAGTGCATATGTTTTAGTCTTTTTTCCCAAAGCCggg                                                        <  2:81530/80‑1 (MQ=255)                         ttAATACACAGCAGTTGTTAAAACGATAAGAAGTTAGCAGGAGTGCATATGTTTTAGTCTTTTTTCCCAAAGCCGGGAACg                                                    >  2:81559/1‑81 (MQ=255)                            aTACACAGCAGTTGTTAAAACGATAAGAAGTTAGCAGGAGTGCATATGTTTTAGTCTTTTTTCCCAAAGCCGGGAACGttt                                                 >  1:81565/1‑81 (MQ=255)                               cacaGCAGTTGTTAAAACGATAAGAAGTTAGCAGGAGTGCATATGTTTTAGTCTTTTTTCCCAAAGCCGGGAACGtttttt                                              >  2:81561/1‑81 (MQ=255)                                 cagcagTTGTTAAAACGATAAGAAGTTAGCAGGAGTGCATATGTTTTAGTCTTTTTTCCCAAAGCCGGGAACGTTTTTtct                                            <  2:81531/81‑1 (MQ=255)                                   gcagTTGTTAAAACGATAAGAAGTTAGCAGGAGTGCATATGTTTTAGTCTTTTTTCCCAAAGCCGGGAACGTTTTTtctct                                          <  2:81532/81‑1 (MQ=255)                                         gttAAAACGATAAGAAGTTAGCAGGAGTGCATATGTTTTAGTCTTTTTTCCCAAAGCCGGGAACGTTTTTTCTCTCGGCCt                                    >  2:81578/1‑81 (MQ=255)                                            aaaaCGATAAGAAGTTAGCAGGAGTGCATATGTTTTAGTCTTTTTTCCCAAAGCCGGGAACGTTTTTTCTCTCGGCCTTc                                  <  2:81533/80‑1 (MQ=255)                                            aaaaCGATAAGAAGTTAGCAGGAGTGCATATGTTTTAGTCTTTTTTCCCAAAGCCGGGAACGTTTTTTCTCTCGGCCTTc                                  <  1:81534/80‑1 (MQ=255)                                              aaCGATAAGAAGTTAGCAGGAGTGCATATGTATTAGTCTTTTTTCCCAAAGCCGGGAACGTTTTTTCTCTCGGCCTTCGtt                               <  2:81535/81‑1 (MQ=255)                                               aCGATAAGAAGTTAGCAGGAGTGCATATGTTTTAGTCTTTTTTCCCAAAGCCGGGAACGTTTTTTCTCTCGGCCTTCGttt                              <  2:81536/81‑1 (MQ=255)                                                  aTAAGAAGTTAGCAGGAGTGCATATGTTTTAGTCTTTTTTCCCAAAGCCGGGAACGTTTTTTCTCTCGGCCTTCGTTTggg                           >  2:81566/1‑81 (MQ=255)                                                   tAAGAAGTTAGCAGGAGTGCATATGTTTTAGTCTTTTTTCCCAAAGCCGGGAACGTTTTTTCTCTCGGCCTTCGTTTGGGc                          <  2:81537/81‑1 (MQ=255)                                                      gaagTTAGCAGGAGTGCATATGTTTTAGTCTTTTTTCCCAAAGCCGGGAACGTTTTTTCTCTCGGCCTTCGTTTGGGCatt                       >  2:81583/1‑81 (MQ=255)                                                        agTTAGCAGGAGTGCATATGTTTTAGTCTTTTTTCCCAAAGCCGGGAACGTTTTTTCTCTCGGCCTTCGTTTGGGCattga                     >  1:81567/1‑81 (MQ=255)                                                        agTTAGCAGGAGTGCATATGTTTTAGTCTTTTTTCCCAAAGCCGGGAACGTTTTTTCTCTCGGCCTTCGTTTGGGCattga                     >  2:81564/1‑81 (MQ=255)                                                         gTTAGCAGGAGTGCATATGTTTTAGTCTTTTTTCCCAAAGCCGGGAACGTTTTTTCTCTCGGCCTTCGTTTGGGCattgat                    <  2:81538/81‑1 (MQ=255)                                                            aGCAGGAGTGCATATGTTTTAGTCTTTTTTCCCAAAGCCGGGAACGTTTTTTCTCTCGGCCTTCGTTTGGGCattgattg                  <  2:81539/80‑1 (MQ=255)                                                              cAGGAGTGCATATGTTTTAGTCTTTTTTCCCAAAGCCGGGAACGTTTTTTCTCTCGGCCTTCGTTTGGGCATTGATTGCCg               >  2:81572/1‑81 (MQ=255)                                                               aGGAGTGCATATGTTTTAGTCTTTTTTCCCAAAGCCGGGAACGTTTTTTCTCTCGGCCTTCGTTTGGGCATTGATTGCCGt              >  1:81579/1‑81 (MQ=255)                                                                ggAGTGCATATGTTTTAGTCTTTTTTCCCAAAGCCGGGAACGTTTTTTCTCTCGGCCTTCGTTTGGGCATTGATTGCCGtt             <  1:81540/81‑1 (MQ=255)                                                                ggAGTGCATATGTTTTAGTCTTTTTTCCCAAAGCCGGGAACGTTTTTTCTCTCGGCCTTCGTTTGGGCATTGATTGCCGt              >  1:81570/1‑80 (MQ=255)                                                                   gTGCATATGTTTTAGTCTTTTTTCCCAAAGCCGGGAACGTTTTTTCTCTCGGCCTTCGTTTGGGCATTGATTGCCGTTAtc          >  1:81571/1‑81 (MQ=255)                                                                     gCATATGTTTTAGTCTTTTTTCCCAAAGCCGGGAACGTTTTTTCTCTCGGCCTTCGTTTGGGCATTGATTGCCGTTAtctt        >  2:81574/1‑81 (MQ=255)                                                                           gTTTTAGTCTTTTTTCCCAAAGCCGGGAACGTTTTTTCTCTCGGCCTTCGTTTGGGCATTGATTGCCGTTATCTTCTGGc   >  2:81586/1‑80 (MQ=255)                                                                            ttttAGTCTTTTTTCCCAAAGCCGGGAACGTTTTTTCTCTCGGCCTTCGTTTGGGCATTGATTGCCGTTATCTTCTGGCa  <  1:81541/80‑1 (MQ=255)                                                                               |                                                                             TGACGCTTTGTAGCGGTCATGCGGTTAATACACAGCAGTTGTTAAAACGATAAGAAGTTAGCAGGAGTGCATATGTTTAAGTCTTTTTTCCCAAAGCCGGGAACGTTTTTTCTCTCGGCCTTCGTTTGGGCATTGATTGCCGTTATCTTCTGGCA  >  NC_000913/396567‑396721

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