Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,474,722	A→C	7.8%	C160W (TGT→TGG)	tal ←	transaldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,474,722	0	A	C	7.8%	148.0 / 2.8	51	C160W (TGT→TGG)	tal	transaldolase
Reads supporting (aligned to +/- strand):  ref base A (22/25);  new base C (2/2);  total (24/27)
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.20e-01

AACCATTGGCCTTGTAGTAGTTGTAGATGCGGGTGACCGACTGCACGCCTGGGTC‑CTCGGCGCCCACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATACCTTCCTTCTCCAGCTGCTCGGCGGCGCGAATGCCTTCCCAGGTGGAGGCCAGCTTGATCAGTACACGCTCTTTTGCCACTCCGGCCGCCTCGTAC  >  NC_002947/2474574‑2474869                                                                                                                                                      |                                                                                                                                                     aaCCATTGGCCTTGTAGTAGTTGTAGATGCGGGTGACCGACTGCACGCCTGGGTC‑CTCGGCGCCCACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGca                                                                                                                                                      >  8:80891/1‑149 (MQ=255)     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                                                                                                                            <  1:112756/149‑1 (MQ=255)                           aTGCGGGTGACCGACTGCACGCCTGGGTC‑CTCGGCGCCCACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGaa                                                                                                                            >  5:180502/1‑149 (MQ=255)                             gCGGGTGACCGACTGCACGCCTGGGTCGCTCGGCGCCCACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAgcctgggcctg                                                                                                                                        >  1:13285/1‑136 (MQ=255)                             gCGGGTGACCGACTGCACGCCTGGGTC‑CTCGGCGCCCACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAgcctgggcctg                                                                                                                                        <  2:13285/135‑1 (MQ=255)                                      cGACTGCACGCCTGGGTC‑CTCGGCGCCCACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑t                                                                                                                  >  1:258767/1‑147 (MQ=255)                                      cGACTGCACGCCTGGGTC‑CTCGGCGCCCACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑t                                                                                                                  >  7:25818/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                                                             >  1:161812/1‑149 (MQ=255)                                                                                          cttGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATACCTTCCTTCTCCAGCTGCTCGGCGGCGCGAATGCCtt                                                             >  3:167118/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                                <  3:106323/149‑1 (MQ=255)                                                                                                                       aCGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATACCTTCCTTCTCCAGCTGCTCGGCGGCGCGAATGCCTTCCCAGGTGGAGGCCAGCTTGATCAGTaca                                <  6:180502/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)                                                                                                                                                      |                                                                                                                                                     AACCATTGGCCTTGTAGTAGTTGTAGATGCGGGTGACCGACTGCACGCCTGGGTC‑CTCGGCGCCCACGTATTCCTGGCCGGTGCTCTTCTTGTACCAGTCGTAGATACGGCCCACGAACGGCGAAATCAAAAACACCCCGGCATCGGCACAAGCCTGGGCCTGGGCGAAGGAGAACAGCAGGG‑TCAGGTTGGTCTGGATACCTTCCTTCTCCAGCTGCTCGGCGGCGCGAATGCCTTCCCAGGTGGAGGCCAGCTTGATCAGTACACGCTCTTTTGCCACTCCGGCCGCCTCGTAC  >  NC_002947/2474574‑2474869

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