Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,552,776	C→A	R381S (CGC→AGC)	uxuA →	D‑mannonate dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,552,776	0	C	A	100.0%	157.4 / NA	45	R381S (CGC→AGC)	uxuA	D‑mannonate dehydratase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (19/26);  total (19/26)

TCGTAAAGCGGAAGGCAAAGAAGACCTGATCCCGATGCGTCCGGACCACGGTCATCAGATGCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGCGCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTA  >  NC_000913/4552643‑4552907                                                                                                                                      |                                                                                                                                    tcgtAAAGCGGAAGGCAAAGAAGACCTGATCCCGATGCATCCGGACCACGGTCATCAGATGCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGgtg                                                                                                                               >  1:263461/1‑140 (MQ=255)           gAAGGCAAAGAAGACCTGATCCCGATGCGTCCGGACCACGGTCATCAGATGCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAACTGGc                                                                                                                     >  1:575672/1‑140 (MQ=255)              ggCAAAGAAGACCTGATCCCGATGCGTCCGGACCACGGTCATCAGATGCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAAc                                                                                                                         >  2:39177/1‑133 (MQ=255)              ggCAAAGAAGACCTGATCCCGATGCGTCCGGACCACGGTCATCAGATGCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAAc                                                                                                                         <  1:39177/133‑1 (MQ=255)                cAAAGAAGACCTGATCCCGATGCGTCCGGACCACGGTCATCAGATGCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGGTAGCGGTGTCGAACTGGCGATcc                                                                                                                >  2:39580/1‑140 (MQ=255)                   agaagaCCTGATCCCGATGCGTCCGGACCACGGTCATCAGATGCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAgc                                                                                                             <  2:367491/140‑1 (MQ=255)                      agaCCTGATCCCGATGCGTCCGGACCACGGTCATCAGATGCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAgcgcg                                                                                                          <  1:86671/140‑1 (MQ=255)                               cccGATGCGTCCGGACCACGGTCATCAGATGCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTa                                                                                                 <  1:18595/140‑1 (MQ=255)                                ccGATGCGTCCGGACCACGGTCATCAGATGCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAg                                                                                                <  2:575672/140‑1 (MQ=255)                                ccGATGCGTCCGGACCACGGTCATCAGATGCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAg                                                                                                <  2:341464/140‑1 (MQ=255)                                       gTCCGGACCACGGTCATCAGATGCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTaa                                                                                         >  2:398668/1‑140 (MQ=255)                                       gTCCGGACCACGGTCATCAGATGCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTaa                                                                                         <  1:361477/140‑1 (MQ=255)                                       gTCCGGACCACGGTCATCAGATGCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTaa                                                                                         >  2:203358/1‑140 (MQ=255)                                        tCCGGACCACGGTCATCAGATGCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAg                                                                                                <  1:40615/132‑1 (MQ=255)                                        tCCGGACCACGGTCATCAGATGCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAg                                                                                                >  2:40615/1‑132 (MQ=255)                                                  ggTCATCAGATGCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGc                                                                              >  2:500221/1‑140 (MQ=255)                                                  ggTCATCAGATGCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGc                                                                              >  2:533094/1‑140 (MQ=255)                                                  ggTCATCAGATGCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGc                                                                              >  1:714951/1‑140 (MQ=255)                                                   gTCATCAGATGCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCa                                                                             <  2:263461/140‑1 (MQ=255)                                                   gTCATCAGATGCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCa                                                                             <  1:498781/140‑1 (MQ=255)                                                              cTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTgcgcgc                                                                  >  1:70700/1‑140 (MQ=255)                                                              cTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTgcgcgc                                                                  >  2:669031/1‑140 (MQ=255)                                                                 gacgacCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGCGCCGt                                                               <  1:533094/140‑1 (MQ=255)                                                                  acgacCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGCGCCGTg                                                              >  1:332690/1‑140 (MQ=255)                                                                      ccTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGCGCCGTGCCgg                                                          >  2:609958/1‑140 (MQ=255)                                                                      ccTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGCGCCGTGCCgg                                                          <  2:714951/140‑1 (MQ=255)                                                                                aaaaCCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGCGCCGTGCCGGTTCCttcttc                                                >  1:610890/1‑140 (MQ=255)                                                                                 aaaCCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGCGCCGTGCCGGTTCCTTCTTcc                                               <  2:332690/140‑1 (MQ=255)                                                                                    ccaaccCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGCGCCGTGCCGGTTCCTTCTTCCTTg                                            >  2:492996/1‑140 (MQ=255)                                                                                        cccAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGCGCCGTGCCGGTTCCTTCTTCCTTGCCGt                                        <  1:492996/140‑1 (MQ=255)                                                                                         ccAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGCGCCGTGCCGGTTCCTTCTTCCTTGCCGTc                                       <  2:92321/140‑1 (MQ=255)                                                                                         ccAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGCGCCGTGCCGGTTCCTTCTTCCTTGCCGTc                                       <  1:609958/140‑1 (MQ=255)                                                                                          cAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGCGCCGTGCCGGTTCCTTCTTCCTTGCCGTCa                                      <  2:271156/140‑1 (MQ=255)                                                                                                  tCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGCGCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGa                              >  2:319186/1‑140 (MQ=255)                                                                                                  tCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGCGCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGa                              >  1:35212/1‑140 (MQ=255)                                                                                                  tCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGCGCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGa                              >  1:171236/1‑140 (MQ=255)                                                                                                      cAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGCGCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGAc                          <  1:319186/140‑1 (MQ=255)                                                                                                                       ggCCTGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGCGCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACg         <  2:319902/140‑1 (MQ=255)                                                                                                                       ggCCTGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGCGCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACg         <  1:500221/140‑1 (MQ=255)                                                                                                                           tGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGCGCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGttt            >  1:62809/1‑133 (MQ=255)                                                                                                                           tGGCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGCGCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGttt            <  2:62809/133‑1 (MQ=255)                                                                                                                            ggCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGCGCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGac    <  1:669031/140‑1 (MQ=255)                                                                                                                            ggCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGCGCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGac    <  1:429048/140‑1 (MQ=255)                                                                                                                            ggCCGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGCGCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGac    <  1:139287/140‑1 (MQ=255)                                                                                                                              ccGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGCGCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGacta  <  1:398668/140‑1 (MQ=255)                                                                                                                              ccGAAGTTAGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGCGCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGacta  <  1:203358/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                    TCGTAAAGCGGAAGGCAAAGAAGACCTGATCCCGATGCGTCCGGACCACGGTCATCAGATGCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGCGCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTA  >  NC_000913/4552643‑4552907

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