Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,983,945	G→A	7.0%	A46T (GCC→ACC)	PP_1776 →	bifunctional mannose‑6‑phosphate isomerase/mannose‑1‑phosphate guanylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,983,945	0	G	A	7.0%	167.4 / 2.2	57	A46T (GCC→ACC)	PP_1776	bifunctional mannose‑6‑phosphate isomerase/mannose‑1‑phosphate guanylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (30/23);  new base A (2/2);  total (32/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.84e-01

AAATATGGAATTGATCCCGGTAATTTTATCCGGTGGCGTTGGTAGCCGTCTGTGGCCAGTATCACGTGAGGCTCATCCCAAGCCGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGAATTGCTGTTCAAGACCGAGGACGAATATCGCACCATCAACAAGCAGAACCTCAGCCAGGGTTACATTCTCGAACCCTTCGGTCGCAATACCGCC  >  NC_002947/1983806‑1984085                                                                                                                                             |                                                                                                                                             aaaTATGGAATTGATCCCGGTAATTTTATCCGGTGGCGTTGGTAGCTGTCTGTGGCCAGTATCACGTGAGGCTCATCCCAAGCCGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTg                                                                                                                                      >  4:248954/1‑148 (MQ=255) aaaTATGGAATTGATCCCGGTAATTTTATCCGGTGGCGTTGGTAGCCGTCTGTGGCCAGTATCACGTGAGGCTCATCCCAAGCCGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTg                                                                                                                                      >  6:253440/1‑148 (MQ=255) aaaTATGGAATTGATCCCGGTAATTTTATCCGGTGGCGTTGGTAGCCGTCTGTGGCCAGTATCACGTGAGGCTCATCCCAAGCCGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGa                                                                                                                                     <  5:253440/149‑1 (MQ=255)   atatGGAATTGATCCCGGTAATTTTATCCGGTGGCGTTGGTAGCCGTCTGTGGCCAGTATCACGTGAGGCTCATCCCAAGCCGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAAc                                                                                                                                   >  2:166862/1‑149 (MQ=255)    tatGGAATTGATCCCGGTAATTTTATCCGGTGGCGTTGGTAGCTGTCTGTGGCCAGTATCACGTGAGGCTCATCCCAAGCCGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACg                                                                                                                                  <  3:248954/149‑1 (MQ=255)    tatGGAATTGATCCCGGTAATTTTATCCGGTGGCGTTGGTAGCCGTCTGTGGCCAGTATCACGTGAGGCTCATCCCAAGCCGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACg                                                                                                                                  >  4:321906/1‑149 (MQ=255)              aTCCCGGTAATTTTATCCGGTGGCGTTGGTAGCCGTCTGTGGCCAGTATCACGTGAGGCTCATCCCAAGCCGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTg                                                                                                                         >  2:50573/1‑148 (MQ=255)                 ccGGTAATTTTATCCGGTGGCGTTGGTAGCCGTCTGTGGCCAGTATCACGTGAGGCTCATCCCAAGCCGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGaaa                                                                                                                      <  1:5221/148‑1 (MQ=255)                    gTAATTTTATCCGGTGGCGTTGGTAGCCGTCTGTGGCCAGTATCACGTGAGGCTCATCCCAAGCCGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGAc                                                                                                                                         <  1:275573/126‑1 (MQ=255)                    gTAATTTTATCCGGTGGCGTTGGTAGCCGTCTGTGGCCAGTATCACGTGAGGCTCATCCCAAGCCGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGAc                                                                                                                                         >  2:275573/1‑126 (MQ=255)                     tAATTTTATCCGGTGGCGTTGGTAGCCGTCTGTGGCCAGTATCACGTGAGGCTCATCCCAAGCCGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCtt                                                                                                                 <  3:232265/149‑1 (MQ=255)                     tAATTTTATCCGGTGGCGTTGGTAGCCGTCTGTGGCCAGTATCACGTGAGGCGCATCCCAAGCCGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCtt                                                                                                                 <  5:383101/149‑1 (MQ=255)                       aTTTTATCCGGTGGCGTTGGTAGCCGTCTGTGGCCAGTATGACGTTCGGCTCATCCCAAGCCGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTa                                                                                                                <  5:220211/148‑1 (MQ=255)                        ttttATCCGGTGGCGTTGGTAGCCGTCTGTGGCCAGTATCACGTGAGGCTCATCCCAAGCCGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACg                                                                                                              <  2:116141/149‑1 (MQ=255)                        ttttATCCGGTGGCGTTGGTAGCCGTCTGTGGCCAGTATCACGTGAGGCTCATCCCAAGCCGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACg                                                                                                              <  3:321906/149‑1 (MQ=255)                          ttATCCGGTGGCGTTGGTAGCCGTCTGTGGCCAGTATCACGTGAGGCTCATCCCAAGCCGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGt                                                                                                            >  1:19608/1‑149 (MQ=255)                                       ttGGTAGCCGTCTGTGGCCAGTATCACGTGAGGCTCATCCCAAGCCGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGaa                                                                                               <  2:280300/149‑1 (MQ=255)                                          gTAGCCGTCTGTGGCCAGTATCACGTGAGGCTCATCCCAAGCCGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGt                                                                                                            <  4:166156/133‑1 (MQ=255)                                          gTAGCCGTCTGTGGCCAGTATCACGTGAGGCTCATCCCAAGCCGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGt                                                                                                            >  3:166156/1‑133 (MQ=255)                                           tAGCCGTCTGTGGCCAGTATCACGTGAGGCTCATCCCAAGCCGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGAATTGc                                                                                           >  3:3944/1‑149 (MQ=255)                                             gCCGTCTGTGGCCAGTATCACGTGAGGCTCATCCCAAGCCGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGAATTGCt                                                                                          <  2:304016/148‑1 (MQ=255)                                                  cTGTGGCCAGTATCACGTGAGGCTCATCCCAAGCCGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGAATTGCTGTTCaa                                                                                    <  3:73841/149‑1 (MQ=255)                                                     tgGCCAGTATCACGTGAGGCTCATCCCAAGCCGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGAATTGCTGTTCAAGa                                                                                  <  7:305216/148‑1 (MQ=255)                                                         aaGTATCACGTGAGGCTCATCCCAAGCCGTTCAT‑TACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGAATTGCTGTTCAAGACCGAg                                                                             <  3:284510/148‑1 (MQ=255)                                                          aGTATCACGTGAGGCTCATCCCAAGCCGTTCATAGCCCCTGCCAGACGGTCCGAACCTCATCCAGAAGACTTACATCCGAGCCGCCGACCTCAACGGTGTGGTTGAAATCATTACGGTCACCAACCGCGAATTGCTg                                                                                         >  1:44257/1‑137 (MQ=255)                                                          aGTATCACGTGAGGCTCATCCCAAGCCGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGAATTGCTg                                                                                         <  2:44257/136‑1 (MQ=255)                                                           gTATCACGTGAGGCTCATCCCAAGCCGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCACCGACCTGAACGGTGTGGTTGaaa                                                                                                                      <  7:227186/106‑1 (MQ=255)                                                           gTATCACGTGAGGCTCATCCCAAGCCGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCACCGACCTGAACGGTGTGGTTGaaa                                                                                                                      >  8:227186/1‑106 (MQ=255)                                                            tATCACGTGAGGCTCATCCCAAGCCGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGAATTGCTGTTCAAGACCGAGGAc                                                                          >  2:185266/1‑149 (MQ=255)                                                              tCACGTGAGGCTCATCCCAAGCCGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGAATTGCTGTTCAAGACCGAGGACg                                                                         <  1:166862/148‑1 (MQ=255)                                                                aCGTGAGGCTCATCCCAAGCCGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGAATTGCTGTTCAAGACCGAGGACGaa                                                                       >  2:263583/1‑148 (MQ=255)                                                                           aTCCCAAGCCGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGAATTGCTGTTCAAGACCGAGGACGAATATCGCACCAt                                                            <  1:50573/148‑1 (MQ=255)                                                                           aTCCCAAGCCGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGAATTGCTGTTCAAGACCGAGGACGAATATCGCACCAt                                                            >  6:317923/1‑148 (MQ=255)                                                                                aaGCCGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCACCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGAATTGCTGTTCAAGACCGAGGACGAATATCGCACCATcaac                                                        >  2:209631/1‑147 (MQ=255)                                                                                   ccGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGAATTGCTGTTCAAGACCGAGGACGAATATCGCACCATCAACAAGCa                                                   >  7:274423/1‑149 (MQ=255)                                                                                         aT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGAATTGCTGTTCAAGACCGAGGACGAATATCGCACCATCAACAAGCAGAACCt                                             >  8:333190/1‑149 (MQ=255)                                                                                            gACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGAATTGCTGTTCAAGACCGAGGACGAATATCGCACCATCAACAAGCAGAACCTCa                                           >  1:352173/1‑149 (MQ=255)                                                                                             aCCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGAATTGCTGTTCAAGACCGAGGACGAATATCGCACCATCAACAAGCAGAACCTCa                                           <  6:377939/148‑1 (MQ=255)                                                                                              cccTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGAATTGCTGTTCAAGACCGAGGACGAATATCGCACCATCAACAAGCAGAACCTCAGc                                         >  6:212074/1‑149 (MQ=255)                                                                                                   ccAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGAATTGCTGTTCAAGACCGAGGACGAATATCGCACCATCAACAAGCAGAACCTCAGCCAggg                                    <  6:89955/149‑1 (MQ=255)                                                                                                    cAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGAATTGCTGTTCAAGACCGAGGACGAATATCGCACCATCAACAAGCAGAACCTCAGCCAGGGt                                   >  4:359461/1‑149 (MQ=255)                                                                                                       aCGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGAATTGCTGTTCAAGACCGAGGACGAATATCGCACCATCAACAAGCAGAACCTCAGCCAGGGTTa                                 >  7:356135/1‑148 (MQ=255)                                                                                                       aCGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGAATTGCTGTTCAAGACCGAGGACGAATATCGCACCATCAACAAGCAGAACCTCAGCCAGGGTTAc                                <  5:212074/149‑1 (MQ=255)                                                                                                          gTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGCACTCCTTACGGTCACCAACCGCGAATTGCTGTTCAAGACCGAGGACGAATATCGCACCATCATCACTCAGAACCTCAGCGAGGGTTACAtt                             >  6:267005/1‑149 (MQ=255)                                                                                                             aGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGAATTGCTGTTCAAGACCGAGGACGAATATCGCACCATCAACAAGCAGAACCTCAGCCAGGGTTACATtctc                          >  1:137923/1‑149 (MQ=255)                                                                                                                 ccTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACCGTCACCAACCGCGAATTGCTGTTCAAGACCGAGGACGAATATCGCACAATCAACAAGCAGAACCTCAGCCAGGGTTACATTCTCGAAc                      >  7:53913/1‑149 (MQ=255)                                                                                                                  cTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGAATTGCTGTTCAAGACCGAGGACGAATATCGCACCATCAACAAGCAGAACCTCAGCCAGGGTTACATTCTCGAAcc                     >  1:202420/1‑149 (MQ=255)                                                                                                                     aTCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCTCGAATTGCTGATCAAGACCGAGGACGAATATCGCACCATCAACAAGCAGAACCTCAGCCAGGGTTACTTTCTCGAACCCtt                  >  5:167341/1‑149 (MQ=255)                                                                                                                        cAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGAATTGCTGTTCAAGACCGAGGACGAATATCGCACCATCAACAAGCAGAACCTCAGCCAGGGTTACATTCTCGAACCCTTCgg               >  6:289778/1‑149 (MQ=255)                                                                                                                         agaagaCTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGAATTGCTGTTCAAGACCGAGGACGAATATCGCACCATCAACAAGCAGAACCTCAGCCAGGGTTACATTCTCGAACCCTTCGGt              >  8:218833/1‑149 (MQ=255)                                                                                                                         agaagaCTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGAATTGCTGTTCAAGACCGAGGACGAATATCGCACCATCAACAAGCAGAACCTCAGCCAGGGTTACATTCTCGAACCCTTCGGt              >  3:128427/1‑149 (MQ=255)                                                                                                                          gaagaCTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGAATTGCTGTTCAAGACCGAGGACGAATATCGCACCATCAACAAGCAGAACCTCAGCCAGGGTTACATTCTCGAACCCTTCGGTc             <  5:317923/149‑1 (MQ=255)                                                                                                                           aagaCTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGAATTGCTGTTCAAGACCGAGGACGAATATCGCACCATCAACAAGCAGAACCTCAGCCAGGGTTACATTCTCGAACCCTTCGGTCg            <  4:128427/149‑1 (MQ=255)                                                                                                                                tttCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGAATTGCTGTTCAAGACCGAGGACGAATATCGCACCATCAACAAGCAGAACCTCAGCCAGGGTTACATTCTCGAACCCTTCGGTCGCAATa       >  2:118201/1‑149 (MQ=255)                                                                                                                                    cTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGAATTGCTGTTCAAGACCGAGGACGAATATCGCACCATCAACAAGCAGAACCTCAGCCAGGGTTACATTCTCGAACCCTTCGGTCGCAATAccgc   >  5:221152/1‑149 (MQ=255)                                                                                                                                    cTGCGAGCCACCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGAATTGCTGTTCAAGACCGAGGACGAATATCGCACCATCAACAAGCAGAACCTCAGCCAGGGTTACATTCTCGAACCCTTCGGTCGCAATAccgc   <  1:209631/149‑1 (MQ=255)                                                                                                                                     tGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGAATTGCTGTTCAAGACCGAGGACGAATATCGCACCATCAACAAGCAGAACCTCAGCCAGGGTTACATTCTCGAACCCTTCGGTCGCAATAccgcc  >  8:103217/1‑149 (MQ=255)                                                                                                                                             |                                                                                                                                             AAATATGGAATTGATCCCGGTAATTTTATCCGGTGGCGTTGGTAGCCGTCTGTGGCCAGTATCACGTGAGGCTCATCCCAAGCCGTTCAT‑GACCCTGCCAGACGGTCAGAACCTCATCCAGAAGACTTTCCTGCGAGCCGCCGACCTGAACGGTGTGGTTGAAATCCTTACGGTCACCAACCGCGAATTGCTGTTCAAGACCGAGGACGAATATCGCACCATCAACAAGCAGAACCTCAGCCAGGGTTACATTCTCGAACCCTTCGGTCGCAATACCGCC  >  NC_002947/1983806‑1984085

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