Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,122,017	T→C	E321E (GAA→GAG)	cpsG ←	phosphomannomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,122,017	0	T	C	100.0%	147.0 / NA	43	E321E (GAA→GAG)	cpsG	phosphomannomutase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/23);  total (20/23)

GTACCAGTTCGCCCAGCGTTTTATCTTTCAGGCACACCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGCCACCTGCGGCAGTCACCACATCAACGGTGTTCCAGGAGAGACGTGGATCGTGGATGATCTTCGCGCCGGGATT  >  NC_000913/2121881‑2122148                                                                                                                                         |                                                                                                                                    gTACCAGTTCGCCCAGCGTTTTATCTTTCAGGCACACCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGtcctcc                                                                                                                                  <  1:576568/140‑1 (MQ=255)  tACCAGTTCGCCCAGCGTTTTATCTTTCAGGCACACCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTcctcct                                                                                                                                 >  2:790640/1‑140 (MQ=255)   aCCAGTTCGTCCAGCGTTTTATCTTTCAGGCACACCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCtt                                                                                                                                >  2:616661/1‑140 (MQ=255)       gTTCGCCCAGCGTTTTATCTTTCAGGCACACCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTgcgc                                                                                                                            >  2:434198/1‑140 (MQ=255)            cccAGCGTTTTATCTTTCAGGCACACCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACg                                                                                                                       <  2:488574/140‑1 (MQ=255)            cccAGCGTTTTATCTTTCAGGCACACCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACg                                                                                                                       <  2:301424/140‑1 (MQ=255)            cccAGCGTTTTATCTTTCAGGCACACCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACg                                                                                                                       >  1:362997/1‑140 (MQ=255)             ccAGCGTTTTATCTTTCAGGCACACCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGt                                                                                                                      >  1:89973/1‑140 (MQ=255)                 cGTTTTATCTTTCAGGCACACCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCtt                                                                                                                  >  2:702228/1‑140 (MQ=255)                 cGTTTTATCTTTCAGGCACACCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCtt                                                                                                                  <  2:611019/140‑1 (MQ=255)                 cGTTTTATCTTTCAGGCACACCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCtt                                                                                                                  <  2:362997/140‑1 (MQ=255)                         cTTTCAGGCACACCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCTTTAATAAAg                                                                                                          >  2:806551/1‑140 (MQ=255)                             cAGGCACACCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCTTTAATAAAGGCgt                                                                                                      >  2:45430/1‑140 (MQ=255)                                    acCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGt                                                                                               >  2:395377/1‑140 (MQ=255)                                        gTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTc                                                                                           <  1:64570/140‑1 (MQ=255)                                           cGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGAc                                                                                        >  2:707986/1‑140 (MQ=255)                                              cGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACAtt                                                                                     >  1:252153/1‑140 (MQ=255)                                                 ccAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTAcc                                                                                  <  1:45430/140‑1 (MQ=255)                                                       gcCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTg                                                                            >  1:320701/1‑140 (MQ=255)                                                              gATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGCCACCTg                                                                     <  1:800398/140‑1 (MQ=255)                                                                   tGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGCCACCTGCGGCa                                                                <  2:406433/140‑1 (MQ=255)                                                                   tGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGCCACCTGCGGCa                                                                <  2:320701/140‑1 (MQ=255)                                                                   tGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGCCACCTGCGGCa                                                                >  2:794029/1‑140 (MQ=255)                                                                        cTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGCCACCTGCgg                                                                  <  1:345838/133‑1 (MQ=255)                                                                        cTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGCCACCTGCgg                                                                  >  2:345838/1‑133 (MQ=255)                                                                        cTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGCCACCTGCGGCAGTcac                                                           <  2:89973/140‑1 (MQ=255)                                                                          gTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGCCACCTGCGGCAGTcacca                                                         >  1:609301/1‑140 (MQ=255)                                                                           tCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGCCACCTGCGGCAGTcaccac                                                        >  2:594998/1‑140 (MQ=255)                                                                                  aaGCGAAATCACGGAAGTAATGTTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGCCACCTGCGGCAGTCACCACATCAACg                                                 <  1:694314/140‑1 (MQ=255)                                                                                   aGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGCCACCTGCGGCAGTCACCACATCAACgg                                                >  2:804396/1‑140 (MQ=255)                                                                                   aGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGCCACCTGCGGCAGTCACCACATCAACgg                                                <  1:616661/140‑1 (MQ=255)                                                                                   aGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGCCACCTGCGGCAGTCACCACATCAACgg                                                <  1:395377/140‑1 (MQ=255)                                                                                           cACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGCCACCTGCGGCAGTCACCACATCAACGGTGTTCCAg                                        <  1:804396/140‑1 (MQ=255)                                                                                             cGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGCCACCTGCGGCAGTCACCACATCAACGGTGTTCCAGga                                      <  2:609301/140‑1 (MQ=255)                                                                                             cGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGCCACCTGCGGCAGTCACCACATCAACGGTGTTCCAGga                                      <  1:434198/140‑1 (MQ=255)                                                                                                aaGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGCCACCTGCGGCAGTCACCACATCAACGGTGTTCCAGgagag                                   <  1:707986/140‑1 (MQ=255)                                                                                                  gTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGCCACCTGCGGCAGTCACCACATCAACGGTGTTCCAGGAGAGAc                                 <  1:702228/140‑1 (MQ=255)                                                                                                   tAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTTCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGCCACCTGCGGCAGTCACCACATCAACGGTGTTCCAGGAGAGAcg                                <  1:380556/140‑1 (MQ=255)                                                                                                   tAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGCCACCTGCGGCAGTCACCACATCAACGGTGTTCCAGGAGAGAcg                                <  1:794029/140‑1 (MQ=255)                                                                                                   tAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGCCACCTGCGGCAGTCACCACATCAACGGTGTTCCAGGAGAGAcg                                <  1:790640/140‑1 (MQ=255)                                                                                                           ggCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGCCACCTGCGGCAGTCACCACATCAACGGTGTTCCAGga                                      <  2:29881/126‑1 (MQ=255)                                                                                                           ggCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGCCACCTGCGGCAGTCACCACATCAACGGTGTTCCAGga                                      >  1:29881/1‑126 (MQ=255)                                                                                                           ggCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGCCACCTGCGGCAGTCACCACATCAACGGTGTTCCAGGAGAGAcgtggatcg                         <  2:90562/139‑1 (MQ=255)                                                                                                           ggCGCTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGCCACCTGCGGCAGTCACCACATCAACGGTGTTCCAGGAGAGAcgtggatcg                         >  1:90562/1‑139 (MQ=255)                                                                                                             cgcTCATTTCGCCACCATAGATGGCGTCCTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGCCACCTGCGGCAGTCACCACATCAACGGTGTTCCAGGAGAGAcgtggatcgtgg                      >  1:696313/1‑140 (MQ=255)                                                                                                                                 aTGGCGTCCTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGCCACCTGCGGCAGTCACCACATCAACGGTGTTCCAGGAGAGACGTGGATCGTGGATGATCTTCGCGCCGGGAtt  >  1:114312/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                    GTACCAGTTCGCCCAGCGTTTTATCTTTCAGGCACACCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGCCACCTGCGGCAGTCACCACATCAACGGTGTTCCAGGAGAGACGTGGATCGTGGATGATCTTCGCGCCGGGATT  >  NC_000913/2121881‑2122148

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