Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,246,015	T→C	A454A (GCA→GCG)	treA ←	periplasmic trehalase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,246,015	0	T	C	100.0%	124.1 / NA	35	A454A (GCA→GCG)	treA	periplasmic trehalase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/11);  total (24/11)

TGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAAC  >  NC_000913/1245870‑1246159                                                                                                                                                  |                                                                                                                                                 tggtgCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGccc                                                                                                                                              >  1:20212/1‑150 (MQ=255)  ggtgCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGccc                                                                                                                                              >  2:52964/1‑149 (MQ=255)  ggtgCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCa                                                                                                                                             <  2:41994/150‑1 (MQ=255)     gCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCa                                                                                                                                         <  1:89077/151‑1 (MQ=255)         aCATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTg                                                                                                                                     >  2:32947/1‑151 (MQ=255)                atTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATc                                                                                                                              >  2:501920/1‑151 (MQ=255)                  ttttCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATccc                                                                                                                            >  2:467971/1‑151 (MQ=255)                   tttCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATccc                                                                                                                            >  2:339292/1‑150 (MQ=255)                       caccaGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCAttgt                                                                                                                       >  1:545990/1‑151 (MQ=255)                        accaGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCAttgtt                                                                                                                      >  1:10986/1‑151 (MQ=255)                              tttttCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCAc                                                                                                                >  1:231618/1‑151 (MQ=255)                               ttttCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACt                                                                                                               >  1:572214/1‑151 (MQ=255)                                   ctcCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTc                                                                                                           >  2:326479/1‑151 (MQ=255)                                    tcCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTc                                                                                                           >  2:327618/1‑150 (MQ=255)                                      ccGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCAcc                                                                                                        <  2:564326/151‑1 (MQ=255)                                                gtgtGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGgtgt                                                                                              >  2:425195/1‑151 (MQ=255)                                                  gtgCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTc                                                                                            >  2:304131/1‑151 (MQ=255)                                                       gAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGgc                                                                                        <  1:624278/150‑1 (MQ=255)                                                        aaCATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGgcc                                                                                       <  2:297658/150‑1 (MQ=255)                                                          cATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGgccgcc                                                                                    <  2:272024/151‑1 (MQ=255)                                                                          ccAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGAtgt                                                                    >  1:69282/1‑151 (MQ=255)                                                                                         cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTcgtcgccgtc                                                     <  2:22658/151‑1 (MQ=255)                                                                                          gCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTcgtcgccgtcg                                                    <  1:395625/151‑1 (MQ=255)                                                                                             aCCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCgtcgccgtcgcc                                                  >  1:613040/1‑150 (MQ=255)                                                                                                   ttttGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTttgt                                           >  2:356742/1‑151 (MQ=255)                                                                                                    tttGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTttgtt                                          >  1:524194/1‑151 (MQ=255)                                                                                                        cccGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGgcg                                      >  1:419970/1‑151 (MQ=255)                                                                                                           gTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGcgcg                                    >  1:207611/1‑150 (MQ=255)                                                                                                                   gTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCg                            <  2:524183/150‑1 (MQ=255)                                                                                                                                  cgACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTa            <  2:159405/151‑1 (MQ=255)                                                                                                                                   gACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTAc           >  1:388337/1‑151 (MQ=255)                                                                                                                                   gACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTAc           >  1:287566/1‑151 (MQ=255)                                                                                                                                         cTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCggg     >  1:60140/1‑151 (MQ=255)                                                                                                                                           gTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGaa   <  2:412275/151‑1 (MQ=255)                                                                                                                                             aaCGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAAc  >  2:131646/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                 TGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAAC  >  NC_000913/1245870‑1246159

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