Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,245,862	(C)6→5	coding (1515/1698 nt)	treA ←	periplasmic trehalase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,245,857	0	C	.	100.0%	101.9 / NA	25	coding (1520/1698 nt)	treA	periplasmic trehalase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (15/10);  total (15/10)

ACTTAACGGTCGGACGCGTCGCCGGAACATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGT  >  NC_000913/1245726‑1245985                                                                                                                                    |                                                                                                                                 aCTTAACGGTCGGACGCGTCGCCGGAACATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCgccgcc                                                                                                                                 <  2:397474/133‑1 (MQ=255)   ttAACGGTCGGACGCGTCGCCGGAACATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGcccc                                                                                                                               <  1:480777/133‑1 (MQ=255)                 cgTCGCCGGAACATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑CCCCCGTTCCGGTGGTGc                                                                                                                <  1:469844/133‑1 (MQ=255)                   tCGCCGGAACATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑CCCCCGTTCCGGTGGTGCTg                                                                                                              >  2:56098/1‑133 (MQ=255)                              ttGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑CCCCCGTTCCGGTGGTGCTGACATCATAttt                                                                                                   >  2:104281/1‑133 (MQ=255)                                      cGGTTGCTCTTTCGGGCATATCTAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTccacca                                                                                           <  1:237959/133‑1 (MQ=255)                                      cGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑CCCCCGATCCTGTGTTGCTGACATCATACTTTTccacca                                                                                           >  2:169009/1‑133 (MQ=255)                                         ttGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑CCCCCGTTCCGCTGGTGCTGACATCATATTTTTGCACCAGCt                                                                                        >  1:251006/1‑133 (MQ=255)                                           gCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCttt                                                                                      >  2:415302/1‑133 (MQ=255)                                           gCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑CCACCGTTCCGGTGGAGTTTCCATCATAATTTTCGACCAGCttt                                                                                      >  2:243341/1‑133 (MQ=255)                                           gCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑CCACCGTTCCGCTGGCGTTTACATCATAGTTTTGAACCAGCtat                                                                                      >  2:237959/1‑131 (MQ=255)                                                ttCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTtctc                                                                                 <  2:373761/133‑1 (MQ=255)                                                ttCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTtctc                                                                                 <  1:6338/133‑1 (MQ=255)                                                ttCGGGCAGAGCAAACCCAGCATTTTCAGCGTGACGCCATTTGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTtctc                                                                                 <  2:164396/133‑1 (MQ=255)                                                 tCGGGCAGATCAAATTCAACATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTcc                                                                                <  1:169009/133‑1 (MQ=255)                                                 tCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTcc                                                                                >  1:483903/1‑133 (MQ=255)                                                  ctggCAGATTCAATCAAGCATTTTCACCATCACGCCATTGGCCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTccc                                                                               <  2:252457/131‑1 (MQ=255)                                                         aTCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑CCCCCGTTCCGGTGGTTCTGACATCTTATATTTACACCAGCTGTTTCTCCCGGTCATa                                                                        >  2:135428/1‑133 (MQ=255)                                                                   gCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGa                                                              >  2:303468/1‑133 (MQ=255)                                                                               taaCGCCATTAGTCCAGCCAAAGCCATCTTGTAATGGATATTCGCCACCGCCG‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAgg                                                  <  2:251006/131‑1 (MQ=255)                                                                                            ccAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTaa                                     >  1:67060/1‑133 (MQ=255)                                                                                               gCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGt                                  >  1:49363/1‑133 (MQ=255)                                                                                                    aGCCATCCTGTAATGGATATTCGCCACCGCCG‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCACAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATc                             >  1:492548/1‑133 (MQ=255)                                                                                                                      atTCGCCATCGCCG‑CGCCCGTTCCGGTGGTTCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCg           <  1:135428/133‑1 (MQ=255)                                                                                                                      atTCGCCACCGCCG‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCg           >  1:257994/1‑133 (MQ=255)                                                                                                                        tCGCCACCGCCG‑CTCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCg           >  2:468713/1‑131 (MQ=255)                                                                                                                        tCGCCACCGCCG‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTa         <  2:49363/133‑1 (MQ=255)                                                                                                                        tCGCCACCGCCG‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTACCAGCTAATGTCCATCGCCACCTCTTTTTGCCCg           <  1:468713/131‑1 (MQ=255)                                                                                                                               cgccg‑cCCGCGTTCCTGTAGTGATGACATCTTATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGt  <  1:217347/133‑1 (MQ=255)                                                                                                                                    |                                                                                                                                 ACTTAACGGTCGGACGCGTCGCCGGAACATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGT  >  NC_000913/1245726‑1245985

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