Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,474,725	A→G	7.7%	A159A (GCT→GCC)	tal ←	transaldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,474,725	0	A	G	7.7%	146.0 / 2.5	53	A159A (GCT→GCC)	tal	transaldolase
Reads supporting (aligned to +/- strand):  ref base A (22/26);  new base G (2/2);  total (25/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

ATTGGCCTTGTAGTAGTTGTAGATGCGGGTGACCGACTGCACGCCTGGGTC‑CTCGGCGCCCACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATACCTTCCTTCTCCAGCTGCTCGGCGGCGCGAATGCCTTCCCAGGTGGAGGCCAGCTTGATCAGTACACGCTCTTTTGCCACTCCGGCCGCCTCGTACAGG  >  NC_002947/2474578‑2474872                                                                                                                                                     |                                                                                                                                                     aTTGGCCTTGTAGTAGTTGTAGATGCGGGTGACCGACTGCACGCCTGGGTC‑CTCGGCGCCCACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAg                                                                                                                                                     >  6:309732/1‑149 (MQ=255) aTTGGCCTTGTAGTAGTTGTAGATGCGGGTGACCGACTGCACGCCTGGGTC‑CTCGGCGCCCACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAg                                                                                                                                                     >  3:7813/1‑149 (MQ=255) aTTGGCCTTGTAGTAGTTGTAGATGCGGGTGACCGACTGCACGCCTGGGTC‑CTCGGCGCCCACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAg                                                                                                                                                     >  2:226750/1‑149 (MQ=255)                 ttGTAGATGCGGGTGACCGACTGCACGCCTGGGTC‑CTCGGCGCCCACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCCCAGgcttgtgccgatgc                                                                                                                                        <  1:89708/146‑6 (MQ=255)                 ttGTAGATGCGGGTGACCGACTGCACGCCTGGGTC‑CTCGGCGCCCACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCCCAGgcttgtgccgatgc                                                                                                                                        <  1:93478/146‑6 (MQ=255)                 ttGTAGATGCGGGTGACCGACTGCACGCCTGGGTC‑CTCGGCGCCCACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCCCAGgcttgtgccgatgc                                                                                                                                        >  2:93478/1‑141 (MQ=255)                 ttGTAGATGCGGGTGACCGACTGCACGCCTGGGTC‑CTCGGCGCCCACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCCCAGgcttgtgccgatgc                                                                                                                                        >  2:89708/1‑141 (MQ=255)                    tAGATGCGGGTGACCGACTGCACGCCTGGGTC‑CTCGGCGCCCACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAgg                                                                                                                                   <  7:80891/148‑1 (MQ=255)                      gATGCGGGTGACCGACTGCACGCCTGGGTC‑CTCGGCGCCCACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGgag                                                                                                                                 <  8:100590/148‑1 (MQ=255)                       aTGCGGGTGACCGACTGCACGCCTGGGTC‑CTCGGCGCCCACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGaa                                                                                                                               >  5:180502/1‑149 (MQ=255)                       aTGCGGGTGACCGACTGCACGCCTGGGTC‑CTCGGCGCCCACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGaa                                                                                                                               <  1:112756/149‑1 (MQ=255)                         gCGGGTGACCGACTGCACGCCTGGGTCGCTCGGCGCCCACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAgcctgggcctg                                                                                                                                           >  1:13285/1‑136 (MQ=255)                         gCGGGTGACCGACTGCACGCCTGGGTC‑CTCGGCGCCCACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAgcctgggcctg                                                                                                                                           <  2:13285/135‑1 (MQ=255)                                  cGACTGCACGCCTGGGTC‑CTCGGCGCCCACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑t                                                                                                                     >  7:25818/1‑147 (MQ=255)                                  cGACTGCACGCCTGGGTC‑CTCGGCGCCCACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑t                                                                                                                     >  1:258767/1‑147 (MQ=255)                                           gCCTGGGTC‑CTCGGCGCCCACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACACGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCt                                                                                                          >  4:221153/1‑149 (MQ=255)                                             cTGGGTC‑CTCGGCGCCCACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTgg                                                                                                        <  2:62773/149‑1 (MQ=255)                                             cTGGGTC‑CTCGGCGCCCACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTg                                                                                                         >  7:105148/1‑148 (MQ=255)                                                ggTC‑CTCGGCGCCCACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGa                                                                                                       <  5:309732/147‑1 (MQ=255)                                                ggTC‑CTCGGCGCCCACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGAt                                                                                                      <  4:7813/148‑1 (MQ=255)                                                 gTC‑CTCGGCGCCCACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATAc                                                                                                    >  6:133256/1‑149 (MQ=255)                                                   c‑cTCGGCGCCCACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAAAAAAAAACACCCCGGCATCGCCACAAGCCTGGGCATGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGAATAact                                                                                                  >  6:7951/1‑149 (MQ=255)                                                        ggCGCCCACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATAccttcc                                                                                               <  3:34260/148‑1 (MQ=255)                                                        ggCGCCCACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATAccttcc                                                                                               <  5:131266/148‑1 (MQ=255)                                                             ccACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATACCTTCCTtctc                                                                                          >  4:106323/1‑148 (MQ=255)                                                                                 ctctTCTTGTACCCGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATACCTTCCTTCTCCAGCTGCTCGGCGGCGCGAAt                                                                     <  3:20687/149‑1 (MQ=255)                                                                                 ctctTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGGCTCAGGTTGGTCTGGATACCTTCCTTCTCCAGCTGCTCGGCGGCGCGaa                                                                      <  1:226750/149‑1 (MQ=255)                                                                                 ctctTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑GCAGGTTGGTATGGATACCTTCCTTCTCCAGCTGCGCGGCGGCGCGAAt                                                                     >  8:194120/1‑149 (MQ=255)                                                                                    ttcttGTAGCAGTCGTAGATACGGCCCACTAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGTG‑TCAGGTTGGTCTGGATACCTTCCTTCTCCAGCTGCTCGGCGGCGCGAAt                                                                     <  5:7951/146‑1 (MQ=255)                                                                                      cttGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATACCTTCCTTCTCCAGCTGCTCGGCGGCGCGAATGCCtt                                                                >  3:167118/1‑149 (MQ=255)                                                                                      cttGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATACCTTCCTTCTCCAGCTGCTCGGCGGCGCGAATGCCtt                                                                >  1:161812/1‑149 (MQ=255)                                                                                         gTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATACCTTCCTTCTCCAGCTGCTCGGCGGCGCGAATGCCTTcc                                                              <  1:174036/148‑1 (MQ=255)                                                                                             cAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATACCTTCCTTCTCCAGCTGCTCGGCGGCGCGAATGCCTTCCCAGGt                                                         <  5:133256/149‑1 (MQ=255)                                                                                              aGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATACCTTCCTTCTCCAGCTGCTCGGCGGCGCGAATGCCTTCCCAGGTg                                                        <  8:105148/149‑1 (MQ=255)                                                                                                  gTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATACCTTCCTTCTCCAGCTGCTCGGCGGCGCGAATGCCTTCCCAGGTGGAgg                                                    <  4:289120/149‑1 (MQ=255)                                                                                                  gTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATACCTTCCTTCTCCAGCTGCTCGGCGGCGCGAATGCCTTCCCAGGTGGAg                                                     <  8:25818/148‑1 (MQ=255)                                                                                                        aCGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATACCTTCCTTCTCCAGCTGCTCggcggc                                                                           >  8:144682/1‑120 (MQ=255)                                                                                                        aCGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATACCTTCCTTCTCCAGCTGCTCggcggc                                                                           <  7:144682/120‑1 (MQ=255)                                                                                                            cccACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATACCTTCCTTCTCCAGCTGCTCGGCGGCGCGAATGCCTTCCCAGGTGGAGGCCAGCTTGAt                                          <  2:258767/149‑1 (MQ=255)                                                                                                                 gAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATACCTTCCTTCTCCAGCTGCTCGGCGGCGCGAATGCCTTCCCAGGTGGAGGCCAGCTTGATCAGTa                                     >  4:1394/1‑149 (MQ=255)                                                                                                                  aaCGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATACCTTCCTTCTCCAGCTGCTCGGCGGCGCGAATGCCTTCCCAGGTGGAGGCCAGCTTGATCAGTac                                    <  3:1394/149‑1 (MQ=255)                                                                                                                   aCGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATACCTTCCTTCTCCAGCTGCTCGGCGGCGCGAATGCCTTCCCAGGTGGAGGCCAGCTTGATCAGTaca                                   <  6:180502/149‑1 (MQ=255)                                                                                                                   aCGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATACCTTCCTTCTCCAGCTGCTCGGCGGCGCGAATGCCTTCCCAGGTGGAGGCCAGCTTGATCAGTaca                                   <  3:106323/149‑1 (MQ=255)                                                                                                                     ggCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATACCTTCCTTCTCCAGCTGCTCGGCGGCGCGAATGCCTTCCCAGGTGGAGGCCAGCTTGATCAGTACACg                                 >  1:213349/1‑149 (MQ=255)                                                                                                                      gCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATACCTTCCTTCT‑CACCAGCTAGGCCGCCCCAATGCCTTCCCACGGGGGGGGCAGCGTGAAGCGTACACGca                               >  4:256525/1‑148 (MQ=255)                                                                                                                            tCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATACCTTCCTTCTCCAGCTGCTCGGCGGCGCGAATGCCTTCCCAGGTGGAGGCCAGCTTGATCAGTACACGCTCtttt                          <  7:194120/149‑1 (MQ=255)                                                                                                                            tCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATACCTTCCTTCTCCAGCTGCTCGGCGGCGCGAATGCCTTCCCAGGTGGAGGCCAGCTTGATCAGTACACGCTCtttt                          >  6:100186/1‑149 (MQ=255)                                                                                                                              aaaaaCACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATACCTTCCTTCTCCAGCTGCTCGGCGGCGCGAATGCCTTCCCAGGTGGAGGCCAGCTTGATCAGTACACGCTCTTTTGc                        >  7:148045/1‑149 (MQ=255)                                                                                                                              aaaaaCACCCCGGCATCGGCACAAGCATGGGCCTGGGCGAAGTAGAACAGCAGGG‑TCAGGTTGGTCTGGATACCTTCCTTCTCCAGCTGCTCGGCGGCGCGAATGCCTTCCCAGGTGGAGGCCAGCTTGATCAGTACACGCTCTTTTGc                        <  7:243764/149‑1 (MQ=255)                                                                                                                                aaaCACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATACCTTCCTTCTCCAGCTGCTCGGCGGCGCGAATGCCTTCCCAGGTGGAGGCCAGCTTGATCAGTACACGCTCTTTTGCCa                      >  3:97976/1‑149 (MQ=255)                                                                                                                                            aTCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATACCTTCCTTCTCCAGCTGCTCGGCGGCGCGAATGCCTTCCCAGGTGGAGGCCAGCTTGATCAGTACACGCTCTTTTGCCACTCCGGCCGCCt          <  2:213349/149‑1 (MQ=255)                                                                                                                                                 cacaAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATACCTTCCTTCTCCAGCTGCTCGGCGGCGCGAATGCCTTCCCAGGTGGAGGCCAGCTTGATCAGTACACGCTCTTTTGCCACTCCGGCCGCCTCGTAc     <  8:148045/149‑1 (MQ=255)                                                                                                                                                 cacaAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATACCTTCCTTCTCCAGCTGCTCGGAGGCGCGAATGCCTTCCCAGGTGGAGGCCAGCTTGATCAGTACACGCTCTTTTGCCACTCCGGCCGCCTCGTa      >  8:142970/1‑148 (MQ=255)                                                                                                                                                    aaGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATACCTTCCTTCTCCAGCTGCTCGGCGGCGCGAATGCCTTCCCAGGTGGAGGCCAGCTTGATCAGTACACGCTCTTTTGCCACTCCGGCCGCCTCGTACAgg  <  7:142970/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     ATTGGCCTTGTAGTAGTTGTAGATGCGGGTGACCGACTGCACGCCTGGGTC‑CTCGGCGCCCACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATACCTTCCTTCTCCAGCTGCTCGGCGGCGCGAATGCCTTCCCAGGTGGAGGCCAGCTTGATCAGTACACGCTCTTTTGCCACTCCGGCCGCCTCGTACAGG  >  NC_002947/2474578‑2474872

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