Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,036,527	C→T	7.3%	L353F (CTC→TTC)	algC →	phosphomannomutase/phosphoglucomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,036,527	0	C	T	7.3%	158.2 / 2.5	55	L353F (CTC→TTC)	algC	phosphomannomutase/phosphoglucomutase
Reads supporting (aligned to +/- strand):  ref base C (28/23);  new base T (2/2);  total (30/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 = 6.63e-01

TGGAAAACCGGCCATTCGCTGATCAAGAAGAAGATGAAGCAGACCGGTTCGCTGCTGGCAGGCGAGATGAGCGGTCAC‑ATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCAGCATCATTGATGCACTGCAACGCGACGCCGACTGGGGC  >  NC_002947/6036380‑6036667                                                                                                                                                     |                                                                                                                                               tgGAAAACCGGCCATTCGCTGATCAAGAAGAAGATGAAGCAGACCGGTTCGCTGCTGGCAGGCGAGATGAGCGGTCAC‑ATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATcc                                                                                                                                                <  7:277707/148‑1 (MQ=255)   gAAAACCGGCCATTCGCTTATCAAGAAGAAGATGAATCAGACCGGTTCGGTGCTGGCAGTCGAGATGAGCGGTCAC‑ATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCa                                                                                                                                             <  5:56844/149‑1 (MQ=255)        ccGGCCATTCGCTGATCAAGAAGAAGATGAAGCAGACCGGTTCGCTGCTGGCAGGCGAGATGAGCGGTCAC‑ATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAg                                                                                                                                        <  2:41715/149‑1 (MQ=255)              aTTCGCTGATCAAGAAGAAGATGAAGCAGACCGGTTCGCTGCTGGCAGGCGAGATGAGCGGTCAC‑ATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAg                                                                                                                                  >  1:129250/1‑149 (MQ=255)                      aTCAAGAAGAAGATGAAGCAGACCGGTTCGCTGCTGGCAGGCGAGCTGAGCGGTCAC‑ATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAgc‑gc                                                                                                                         >  5:89550/1‑149 (MQ=255)                       tCAAGAAGAAGATGAAGCAGACCGGTTCGCTGCTGGCAGGCGAGATGAGCGGTCAC‑ATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAgc‑gc                                                                                                                         <  4:225592/148‑1 (MQ=255)                          agaagaagaTGAAGCAGACCGGTTCGCTGCTGGCAGGCGAGATGAGCGGTCAC‑ATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGa                                                                                                                     >  6:309580/1‑148 (MQ=255)                              gaagaTGAAGCAGACCGGTTCGCTGCTGGCAGGCGAGATGAGCGGTCAC‑ATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAAcc                                                                                                                >  2:254133/1‑149 (MQ=255)                                            ccGGTTCGCTGCTGGCAGGCGAGATGAGCGGTCAC‑ATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCtt                                                                                                   <  3:110350/148‑1 (MQ=255)                                                   gctgctGGCAGGCGAGATGAGCGGTCACTATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAAc                                                                                            >  1:18680/1‑149 (MQ=255)                                                   gctgctGGCAGGCGAGATGAGCGGTCAC‑ATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAATACCGAGCAGAGC‑GC‑CGAAAACATGTTTGCCGCCTTCCCGAAc                                                                                            >  7:199160/1‑148 (MQ=255)                                                          gCAGGCGAGATGAGCGGTCAC‑ATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTc                                                                                    >  7:198809/1‑149 (MQ=255)                                                              gCGAGATGAGCGGTCAC‑ATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACg                                                                                >  6:244996/1‑149 (MQ=255)                                                                    tGAGCGGTCAC‑ATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTCCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGa                                                                           >  7:246736/1‑148 (MQ=255)                                                                       gCGGTCAC‑ATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGAATGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCATTCCCGAACGACATTTCCACGCCGGAAAt                                                                        >  8:39836/1‑148 (MQ=255)                                                                        aggTCAC‑ATCTTCATCAATGAACGCTGGTACGGTTTTTACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CTAAAACCTGTTTGCCGCATTCCCTAACGACATTTCCACGCCGGAAATc                                                                       <  4:103069/147‑1 (MQ=255)                                                                          gTCAC‑ATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCaa                                                                     >  8:170352/1‑148 (MQ=255)                                                                           tCAC‑ATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAata                                                                   <  7:170352/149‑1 (MQ=255)                                                                             ac‑aTCTTAATCAAGGTACGCTGGTACGGGTTTGACGACGGCATCTATATCGCCGCGCGCCTGCTGGAGATCCTCATCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAatat                                                                  <  8:199160/148‑1 (MQ=255)                                                                                   ttCATCAACGATCGCTGGTTCGGTTTTGTCGACGGGATCTATATCCCCGCGCGGCTGCTGGAGATCCTCACCAAGACCGATCAGAGCGGCTCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGAt                                                              <  1:254133/149‑1 (MQ=255)                                                                                     catcaAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGAtgtg                                                           <  5:244996/148‑1 (MQ=255)                                                                                        caAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAATTCAATATTGATGTGAcc                                                        >  5:201064/1‑148 (MQ=255)                                                                                        caAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGAcc                                                        >  1:314034/1‑148 (MQ=255)                                                                                        caAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGAGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACAGAGCAGAGC‑GC‑CGAAAACCTGTTTGCAGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATTTGa                                                          >  1:318848/1‑146 (MQ=255)                                                                                         aaGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCa                                                                      >  8:239525/1‑133 (MQ=255)                                                                                         aaGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCa                                                                      <  7:239525/133‑1 (MQ=255)                                                                                         aaGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTGTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGGTACcga                                                      >  3:225970/1‑149 (MQ=255)                                                                                          aGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACcgac                                                     >  1:178172/1‑149 (MQ=255)                                                                                          aGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACcgac                                                     >  8:108959/1‑149 (MQ=255)                                                                                          aGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACcgac                                                     >  2:39292/1‑149 (MQ=255)                                                                                           ggAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACcgac                                                     >  3:64051/1‑148 (MQ=255)                                                                                                  tGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTa                                              >  6:285186/1‑148 (MQ=255)                                                                                                  tGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTa                                              >  2:245581/1‑148 (MQ=255)                                                                                                    gTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTaaa                                            <  5:309580/148‑1 (MQ=255)                                                                                                    gTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAAt                                           >  8:225008/1‑149 (MQ=255)                                                                                                    gTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAAt                                           >  2:106305/1‑149 (MQ=255)                                                                                                       cGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTc                                         <  8:198809/148‑1 (MQ=255)                                                                                                       cGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTc                                         <  7:108959/148‑1 (MQ=255)                                                                                                       cGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCa                                        <  2:18680/149‑1 (MQ=255)                                                                                                       cGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCa                                        <  1:245581/149‑1 (MQ=255)                                                                                                        ggTTTTGACGACGGCATCTATAGCTCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCAg                                       <  6:89550/149‑1 (MQ=255)                                                                                                                cgacgGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCAGCATCATTg                               <  6:78979/149‑1 (MQ=255)                                                                                                                 gacgGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCAGCATCATTGa                              >  4:75915/1‑149 (MQ=255)                                                                                                                           ataGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCAGCATCATTGATGCACTGCaa                    >  1:29028/1‑149 (MQ=255)                                                                                                                            taGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGCCATTTCCACGCCGGAAATCAATATTGATGTGCCCGACGAGGGTAAATTCAGCATCATTGATGCACTGCAAc                   >  5:324171/1‑149 (MQ=255)                                                                                                                              gctccgcGCGCCTGCTGGAGATCCTAAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTTCCGCATTCCCTAACGACAATTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCAGCATCATTGATGCACTGCAAcgc                 <  2:318848/146‑1 (MQ=255)                                                                                                                              gcgcCGCGCGCCTGCTGGAGATCCTCAGCAATACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAATTCAATATTGATGTGACCGACGAGGGTAAATTCAGCATCATTGATGCACTGCAAcgc                 <  6:201064/149‑1 (MQ=255)                                                                                                                              gcgcCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCAGCATCATTGATGCACTGCAAcgc                 <  2:314034/149‑1 (MQ=255)                                                                                                                                      cgcCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCAGCATCATTGATGCACTGCAACGCGACGCCGa         <  7:286545/149‑1 (MQ=255)                                                                                                                                       tcttgctgAGGATCTCCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCGCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCAGCATCATTGATGCACTGCAACGCGAc              <  2:98792/140‑1 (MQ=255)                                                                                                                                       tcttgctgAGGATCTCCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCAGCATCATTGATGCACTGCAACGCGAc              >  1:98792/4‑143 (MQ=255)                                                                                                                                       tcttgctgAGGATCTCCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCAGCATCATTGATGCACTGCAACGCGAc              <  2:75680/140‑1 (MQ=255)                                                                                                                                       tcttgctgAGGATCTCCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTGGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCAGCATCATTGATGCACTGCAACGCGAc              >  1:75680/4‑143 (MQ=255)                                                                                                                                            ctgGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCCGCATCATTGATGCACTGCAACGCGACGCCGACTggg    <  3:75915/148‑1 (MQ=255)                                                                                                                                            ctgGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCAGCATCATTGATGCACTGCAACGCGACGCCGACTgggg   <  2:178172/149‑1 (MQ=255)                                                                                                                                             tgGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCAGCATCAGTGATGCACTGCAACGCGACGCCGACTGGGGc  >  2:273023/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                               TGGAAAACCGGCCATTCGCTGATCAAGAAGAAGATGAAGCAGACCGGTTCGCTGCTGGCAGGCGAGATGAGCGGTCAC‑ATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGC‑GC‑CGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCAGCATCATTGATGCACTGCAACGCGACGCCGACTGGGGC  >  NC_002947/6036380‑6036667

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: ‑