Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,122,707	G→A	8.9%	G91G (GGC→GGT)	cpsG ←	phosphomannomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,122,707	0	G	A	8.9%	166.3 / 5.8	56	G91G (GGC→GGT)	cpsG	phosphomannomutase
Reads supporting (aligned to +/- strand):  ref base G (27/24);  new base A (2/3);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCATCGACGGGAGGAAAGTCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTAAACCTTTCGCCAGCGCCAGTTTTAAGGTTTCGCTGGTGAGGCGGACATCACCG  >  NC_000913/2122559‑2122848                                                                                                                                                     |                                                                                                                                              tttCATCGACGGGAGGAAAGTCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATgccg                                                                                                                                               >  6:263547/1‑149 (MQ=255)      tCGACGGGAGGAAAGTCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTACAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATc                                                                                                                                          <  7:321888/149‑1 (MQ=255)           gggAGGAAAGTCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGCGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGc                                                                                                                                     >  6:248272/1‑149 (MQ=255)              aggaAAGTCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCga                                                                                                                                  <  7:126207/149‑1 (MQ=255)               ggaAAGTCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGcc                                                                                                                                    >  8:240247/1‑146 (MQ=255)                gaAAGTCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCgaga                                                                                                                                <  7:240247/149‑1 (MQ=255)                gaAAGTCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCgaga                                                                                                                                <  1:252988/149‑1 (MQ=255)                     tCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGaa                                                                                                                           <  4:314799/149‑1 (MQ=255)                       gTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACg                                                                                                                         >  7:52492/1‑149 (MQ=255)                          tgattCAGCCAGACGCTGGACGTCGCGCAGTCCGGGATCTCCGCTGAGCGGGCGAGCCCCCTCGCGAACCAGCGTCCTGCCGGTATCATACATCGGAGTATGGCTGGCGGTCACTTCAATGCCGCCATCCACGCCGAGATGGAACGTg                                                                                                                       <  2:147249/145‑1 (MQ=255)                           gCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTgg                                                                                                                      <  1:282613/148‑1 (MQ=255)                             ttCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATTGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCg                                                                                                                    >  8:337972/1‑148 (MQ=255)                             ttCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCg                                                                                                                    >  3:318087/1‑148 (MQ=255)                              tCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGa                                                                                                                   >  8:241566/1‑148 (MQ=255)                                    aGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGc                                                                                                                     >  8:177198/1‑140 (MQ=255)                                    aGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGc                                                                                                                     <  7:177198/140‑1 (MQ=255)                                    aGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGaaa                                                                                                                 >  2:120418/1‑144 (MQ=255)                                    aGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGAt                                                                                                            <  4:318087/149‑1 (MQ=255)                                            gACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCgg                                                                                                    >  2:369773/1‑149 (MQ=255)                                                tCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGATGGCTGTAACTTCAATGCCACCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTg                                                                                                  >  2:294817/1‑147 (MQ=255)                                                     cAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACa                                                                                           <  5:263547/149‑1 (MQ=255)                                                        tCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCACCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATAc                                                                                        >  2:22168/1‑149 (MQ=255)                                                          cGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATAcc                                                                                       <  1:114071/148‑1 (MQ=255)                                                                   cGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATcc                                                                                <  3:171670/146‑1 (MQ=255)                                                                   cGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCACCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGc                                                                             <  1:22168/149‑1 (MQ=255)                                                                        aTCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAg                                                                              >  7:281257/1‑143 (MQ=255)                                                                        aTCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAg                                                                              <  8:281257/143‑1 (MQ=255)                                                                            ggCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATc                                                                                 <  3:151883/136‑1 (MQ=255)                                                                            ggCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATc                                                                                 >  4:151883/1‑136 (MQ=255)                                                                             gCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGc                                                                   <  1:120418/149‑1 (MQ=255)                                                                                      cTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCACCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCt                                                          <  2:362407/149‑1 (MQ=255)                                                                                         gcgAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTa                                                       <  7:29012/149‑1 (MQ=255)                                                                                              ccAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTAAACCt                                                  >  2:238486/1‑149 (MQ=255)                                                                                               cAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGaaa                                                                                                                 >  3:196183/1‑85 (MQ=255)                                                                                               cAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGaaa                                                                                                                 <  4:196183/85‑1 (MQ=255)                                                                                                aGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTAAACCttt                                                >  7:217373/1‑149 (MQ=255)                                                                                                       tGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTAAACCTTTcgccagc                                         >  6:372714/1‑149 (MQ=255)                                                                                                         ccGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTAAACCTTTcgccagcgc                                       >  2:202506/1‑149 (MQ=255)                                                                                                           gTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTAAACCTTTc                                               >  7:22374/1‑139 (MQ=255)                                                                                                           gTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTAAACCTTTc                                               <  8:22374/139‑1 (MQ=255)                                                                                                           gTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAAGATCCAGCACGTCAACGCCCGCATCCTGTAAACCTTTc                                               >  4:242853/1‑139 (MQ=255)                                                                                                                  tCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTAAACCTTTCGCCAGCGCCAGTTTTaa                              <  8:52492/149‑1 (MQ=255)                                                                                                                  tCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTAAACCTTTCGCCAGCGCCAGTTTTaa                              <  7:337972/149‑1 (MQ=255)                                                                                                                  tCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTAAACCTTTCGCCAGCGCCAGTTTTaa                              <  5:248272/149‑1 (MQ=255)                                                                                                                     aTCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTAAACCTTTCGCCAGCGCCAGTTTTAAGGt                           >  6:256075/1‑149 (MQ=255)                                                                                                                     aTCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTAAACCTTTCGCCAGCGCCAGTTTTAAGGt                           >  1:242345/1‑149 (MQ=255)                                                                                                                        ggATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTAAACCTTTCGCCAGCGCCAGTTTTaa                              <  2:124130/143‑1 (MQ=255)                                                                                                                        ggATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTAAACCTTTCGCCAGCGCCAGTTTTaa                              >  1:124130/1‑143 (MQ=255)                                                                                                                             aTGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTAAACCTTTCGCCAGCGCCAGTTTTAAGGTTTCGCTgg                   >  1:364934/1‑149 (MQ=255)                                                                                                                                       gTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTAAACCTTTCGCCAGCGCCAGTTTTAAGGTTTCGCTGGTGAGGCGGAc         <  1:202506/149‑1 (MQ=255)                                                                                                                                       gTAACTTCAATGCCACCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTAAACCTTTCGCCAGCGCCAGTTTTAAGGTTTCGCTGGTGAGGCGGAc         <  1:294817/149‑1 (MQ=255)                                                                                                                                        tAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTAAACCTTTCGCCAGCGCCAGTTTTAAGGTTTCGCTGGTGAGGCGGACa        >  7:342964/1‑149 (MQ=255)                                                                                                                                         aaCTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTAAACCTTTCGCCAGCGCCAGTTTTAAGGTTTCGCTGGTGAGGCGGACAt       >  1:119466/1‑149 (MQ=255)                                                                                                                                           cTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTAAACCTTTCGCCAGCGCCAGTTTTAAGGTTTCGCTGGTGAGGCGGACATCa     >  4:256475/1‑149 (MQ=255)                                                                                                                                            ttCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTAAACCTTTCGCCAGCGCCAGTTTTAAGGTTTCGCTGGTGAGGCGGACATCAc    >  6:86831/1‑149 (MQ=255)                                                                                                                                            ttCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTAAACCTTTCGCCAGCGCCAGTTTTAAGGTTTCGCTGGTGAGGCGGACATCAc    >  1:259360/1‑149 (MQ=255)                                                                                                                                               aaTGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTAAACCTTTCGCCAGCGCCAGTTTTAAGGTTTCGCTGGTGAGGCGGACATCACCg  <  1:369773/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                              TTTCATCGACGGGAGGAAAGTCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTAAACCTTTCGCCAGCGCCAGTTTTAAGGTTTCGCTGGTGAGGCGGACATCACCG  >  NC_000913/2122559‑2122848

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