Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,093,668	(T)5→4	coding (605/825 nt)	rhaD ←	rhamnulose‑1‑phosphate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,093,664	0	T	.	100.0%	234.3 / NA	50	coding (609/825 nt)	rhaD	rhamnulose‑1‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (26/24);  total (26/24)

GTGTCGATTAAACCGAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATGTTTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCAT  >  NC_000913/4093584‑4093740                                                                                 |                                                                             gtgtCGATTAAACCGAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATGt                                                                              <  1:808032/81‑1 (MQ=255) gtgtCGATTAAACCGAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATGt                                                                              >  1:808105/1‑81 (MQ=255)    tCGATTAAACCGAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATGtttt                                                                           <  2:808033/81‑1 (MQ=255)      gATTAAACCGAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATGttttgc                                                                         <  1:808034/81‑3 (MQ=255)      gATTAAACCGAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATGttttg                                                                          >  2:808109/1‑79 (MQ=255)       aTTAAACCGAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATGttttgca                                                                        >  2:808035/1‑78 (MQ=255)        ttAAACCGAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATGttttgcat                                                                       <  2:808036/81‑5 (MQ=255)        ttAAACCGAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATGttttgcat                                                                       >  1:808106/1‑77 (MQ=255)         tAAACCGAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCAtc                                                                     >  2:808097/1‑81 (MQ=255)          aaaCCGAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCAtcc                                                                    >  1:808114/1‑80 (MQ=255)           aaCCGAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCAtctc                                                                   <  1:808037/81‑1 (MQ=255)            aCCGAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCATctct                                                                  >  2:808118/1‑81 (MQ=255)             ccGAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCtt                                                                 <  2:808038/81‑1 (MQ=255)               gAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTtgt                                                               >  1:808116/1‑81 (MQ=255)               gAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTtgt                                                               >  1:808107/1‑81 (MQ=255)               gAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTtgt                                                               >  1:808091/1‑81 (MQ=255)                     ttCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTc                                                         >  1:808102/1‑81 (MQ=255)                      tCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTc                                                         <  2:808039/80‑1 (MQ=255)                      tCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCg                                                        >  2:808110/1‑81 (MQ=255)                       cATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCg                                                        <  2:808040/80‑1 (MQ=255)                          ccAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTg                                                    >  1:808099/1‑81 (MQ=255)                           cAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTgg                                                   <  1:808041/81‑1 (MQ=255)                           cAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTgg                                                   >  2:808113/1‑81 (MQ=255)                             gCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGcc                                                 <  2:808042/81‑1 (MQ=255)                               gTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGCCGa                                               >  1:808112/1‑81 (MQ=255)                                  ggTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGCCGAttt                                            <  1:808044/81‑1 (MQ=255)                                  ggTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGCCGAttt                                            <  2:808045/81‑1 (MQ=255)                                  ggTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGCCGAttt                                            >  2:808124/1‑81 (MQ=255)                                  ggTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCAAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGCCGAttt                                            <  1:808043/81‑1 (MQ=255)                                   gTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTc                                           <  1:808047/81‑1 (MQ=255)                                   gTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTc                                           <  1:808046/81‑1 (MQ=255)                                   gTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTc                                           >  1:808115/1‑81 (MQ=255)                                      cGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTc                                        >  2:808096/1‑81 (MQ=255)                                        cTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCg                                      <  2:808048/81‑1 (MQ=255)                                          gCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTg                                    >  1:808125/1‑81 (MQ=255)                                          gCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTg                                    >  2:808119/1‑81 (MQ=255)                                            cGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGcc                                  <  1:808050/81‑1 (MQ=255)                                            cGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGcc                                  <  2:808049/81‑1 (MQ=255)                                              aaGACGCCGTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCgg                                >  2:808052/1‑81 (MQ=255)                                              aaGACGCCGTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCgg                                <  2:808051/81‑1 (MQ=255)                                                  cGCCGTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCAc                            <  1:808053/81‑1 (MQ=255)                                                    ccGTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCa                          <  1:808054/81‑1 (MQ=255)                                                      gTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATc                        <  2:808055/81‑1 (MQ=255)                                                      gTGGAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATc                        >  2:808100/1‑81 (MQ=255)                                                       tGGAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATcc                       <  1:808056/81‑1 (MQ=255)                                                       tGGAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATc                        >  2:808094/1‑80 (MQ=255)                                                        ggAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCa                      >  1:808089/1‑81 (MQ=255)                                                         gAAGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCAc                     >  1:808103/1‑81 (MQ=255)                                                          aaGGGCCACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACg                    >  1:808098/1‑81 (MQ=255)                                                             ggCCACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGc                  >  1:808104/1‑80 (MQ=255)                                                              gCCACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCaa                <  2:808057/81‑1 (MQ=255)                                                               ccACAACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCaa                >  2:808101/1‑80 (MQ=255)                                                                cacaACACCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCaaaa              >  2:808135/1‑81 (MQ=255)                                                                    acacCAGCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGc           >  1:808121/1‑80 (MQ=255)                                                                          gCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGcc    <  1:808058/81‑1 (MQ=255)                                                                          gCGAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGcc    <  1:808059/81‑1 (MQ=255)                                                                            gAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCAt  <  1:808060/81‑1 (MQ=255)                                                                            gAATG‑TTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCAt  <  1:808061/81‑1 (MQ=255)                                                                                 |                                                                             GTGTCGATTAAACCGAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATGTTTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCAT  >  NC_000913/4093584‑4093740

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