Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,937,676	A→T	Q128L (CAA→CTA)	rbsK →	ribokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,937,676	0	A	T	100.0%	135.1 / NA	35	Q128L (CAA→CTA)	rbsK	ribokinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (17/18);  total (17/18)

TTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGC  >  NC_000913/3937598‑3937751                                                                               |                                                                            ttAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCTAc                                                                           <  2:719824/81‑1 (MQ=255)     tGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCTACGTGa                                                                       <  1:719825/81‑1 (MQ=255)        cGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCTACGTGAGCg                                                                    <  1:719826/81‑1 (MQ=255)        cGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCTACGTGAGCg                                                                    >  1:719871/1‑81 (MQ=255)               gagaATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCTACGTGAGCGTATTGc                                                              >  1:719861/1‑80 (MQ=255)                agaATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCTACGTGAGCGTATTGCCa                                                            <  1:719828/81‑1 (MQ=255)                agaATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCTACGTGAGCGTATTGCCa                                                            <  1:719827/81‑1 (MQ=255)                 gaATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCTACGTGAGCGTATTGCCaa                                                           <  2:719829/81‑1 (MQ=255)                  aaTGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCTACGTGAGCGTATTGCCAAc                                                          <  1:719831/81‑1 (MQ=255)                  aaTGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCTACGTGAGCGTATTGCCAAc                                                          >  2:719863/1‑81 (MQ=255)                  aaTGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCTACGTGAGCGTATTGCCAAc                                                          <  1:719830/81‑1 (MQ=255)                   aTGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCTACGTGAGCGTATTGCCAAcg                                                         >  2:719869/1‑81 (MQ=255)                      tCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCTACGTGAGCGTATTGCCAACGCGt                                                      >  1:719870/1‑81 (MQ=255)                        aTCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCTACGTGAGCGTATTGCCAACGCGTc                                                     <  2:719832/80‑1 (MQ=255)                           ggTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCTACGTGAGCGTATTGCCAACGCGTCAGc                                                  <  1:719833/80‑1 (MQ=255)                                 cATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCTACGTGAGCGTATTGCCAACGCGTCAGCattatt                                            >  1:719876/1‑80 (MQ=255)                                 cATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCTACGTGAGCGTATTGCCAACGCGTCAGCattatt                                            >  1:719872/1‑80 (MQ=255)                                 cATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCTACGTGAGCGTATTGCCAACGCGTCAGCattatt                                            <  2:719834/80‑1 (MQ=255)                                  aTGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCTACGTGAGCGTATTGCCAACGCGTCAGCAttatta                                           <  2:719835/80‑1 (MQ=255)                                  aTGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCTACGTGAGCGTATTGCCAACGCGTCAGCATTATTaa                                          >  2:719866/1‑81 (MQ=255)                                       ggCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCTACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAg                                     <  2:719836/81‑1 (MQ=255)                                       ggCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCTACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAg                                     >  2:719881/1‑81 (MQ=255)                                         cgcTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCTACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCt                                   <  1:719837/81‑1 (MQ=255)                                          gcTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCTACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTg                                  <  1:719838/81‑1 (MQ=255)                                           cTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCTACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTgg                                 <  2:719839/81‑1 (MQ=255)                                             aaTGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCTACGTGAGCGTATTGTCAACGCGTCAGCATTATTAATGCAGCTGGaa                               >  2:719883/1‑81 (MQ=255)                                               tgctgcCCTTTCCCCGGCGCTGGTGGAAGCGCTACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATc                             <  2:719840/81‑1 (MQ=255)                                                    cccTTTCCCCGGCGCTGGTGGAAGCGCTACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATcacc                          <  2:719841/79‑1 (MQ=255)                                                         tCCCCGGCGCTGGTGGAAGCGCTACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGa                    >  1:719887/1‑80 (MQ=255)                                                                cgcTGGTGGAAGCGCTACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGa            >  1:719886/1‑81 (MQ=255)                                                                 gcTGGTGGAAGCGCTACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGAt           >  1:719873/1‑81 (MQ=255)                                                                 gcTGGTGGAAGCGCTACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGAt           >  2:719888/1‑81 (MQ=255)                                                                  cTGGTGGAAGCGCTACGTGAGCGTATTGCCAACGCGTCGGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATg          >  2:719879/1‑81 (MQ=255)                                                                    ggtggAAGCGCTACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATgg         >  2:719882/1‑80 (MQ=255)                                                                          aGCGCTACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAgcggc  >  2:719884/1‑81 (MQ=255)                                                                               |                                                                            TTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGC  >  NC_000913/3937598‑3937751

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