Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	487,837	G→C	10.5%	T155T (ACC→ACG)	pdxA ←	4‑hydroxythreonine‑4‑phosphate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	487,837	0	G	C	10.5%	157.7 / 9.6	57	T155T (ACC→ACG)	pdxA	4‑hydroxythreonine‑4‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (24/27);  new base C (2/4);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01

TACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAG‑GCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGTGATCATCCCGGCGAAGTGCCCGTCCAGGC  >  NC_002947/487695‑487979                                                                                                                                                |                                                                                                                                                tACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACCCGCGGCCTG‑CGCGTGGcc                                                                                                                                           <  5:249320/149‑1 (MQ=255) tACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACCCGCGGCCTG‑CGCGTGGcc                                                                                                                                           <  2:321598/149‑1 (MQ=255)         ttGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGCGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGcatcat                                                                                                                                   >  3:107947/1‑149 (MQ=255)                    aaaCTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACCCGCGGCCTG‑CGCGTGGccctggt                                                                                                                                      <  4:100445/135‑6 (MQ=255)                    aaaCTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACCCGCGGCCTG‑CGCGTGGccctggt                                                                                                                                      >  3:100445/1‑130 (MQ=255)                     aaCTTGTCGCGCATGTCGGCATGCATGATACGGGTTACACGTTCGCCGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCCGGGTGGCCAGCATCATTACCACTtgt                                                                                                                         <  1:198917/148‑1 (MQ=255)                     aaCTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTtgt                                                                                                                         >  2:198917/1‑147 (MQ=255)                             gcgcATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGtgcggtgcgg                                                                                                                <  4:261105/148‑1 (MQ=255)                             gcgcATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGtgcggtgcgg                                                                                                                <  3:139414/148‑1 (MQ=255)                              tgcATGCCGACATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTgtgcg                                                                                                                      <  6:101036/140‑1 (MQ=255)                              tgcATGCCGACATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTgtgcg                                                                                                                      >  5:101036/2‑141 (MQ=255)                              cgcATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGcc                                                                                                                                           >  5:104329/1‑120 (MQ=255)                              cgcATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGcc                                                                                                                                           <  6:104329/120‑1 (MQ=255)                               gcATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGt                                                                                                              <  4:63994/148‑1 (MQ=255)                               gcATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGt                                                                                                              >  3:63994/1‑148 (MQ=255)                                       gCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCgg                                                                                                      >  7:166980/1‑148 (MQ=255)                                       gCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCgg                                                                                                      >  2:164229/1‑148 (MQ=255)                                         aTGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAAAATCACGCAGTGGCAGGTGCGTGGTAAACAGGGCCACGCGCAGGCCG‑AGGGTGGCCAGCATCATTACCACTTGTCCGGTGCGGGTCAGTTCGGc                                                                                                     >  3:280966/1‑147 (MQ=255)                                          tGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAg                                                                                                  <  4:261334/149‑1 (MQ=255)                                                  aCGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGCGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCg                                                                                           <  4:107947/148‑1 (MQ=255)                                                   cGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCgg                                                                                          <  4:238741/148‑1 (MQ=255)                                                         acacGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGcc                                                                                    <  5:276490/148‑1 (MQ=255)                                                             gTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGaa                                                                               <  8:166980/149‑1 (MQ=255)                                                                  aCGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAG‑Gc                                                                         >  7:126718/1‑149 (MQ=255)                                                                      gcTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGCCAAGGAATTCGGTATTACCGGAAAAG‑GAGATa                                                                     >  6:252662/1‑149 (MQ=255)                                                                          tGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAG‑GCGATACCGc                                                                 >  6:31734/1‑149 (MQ=255)                                                                           gCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAG‑GCGATACCGCt                                                                >  4:290925/1‑149 (MQ=255)                                                                            cAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCgg                                                                                          <  1:64939/123‑1 (MQ=255)                                                                            cAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCgg                                                                                          >  2:64939/1‑123 (MQ=255)                                                                            cAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAG‑GCGATACCGCt                                                                <  5:31734/148‑1 (MQ=255)                                                                             aGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAG‑GCGATACCGCttt                                                              >  2:345293/1‑149 (MQ=255)                                                                               tGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAG‑GCGATACCGCTTTCg                                                            >  5:354501/1‑149 (MQ=255)                                                                               tGATAGCGTCGGCAACATAACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAG‑GCGATACCGCTTTCg                                                            >  4:308956/1‑149 (MQ=255)                                                                                 atatCGTCGGCAACATCAAGCAGTGGCAGGTGCGTGGTCACCAGGGCCTCGCGCAGGGCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAATAG‑GCGATACCGCTTTCGtt                                                          <  4:280966/149‑1 (MQ=255)                                                                                    gCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAG‑GCGATACCGCTTTCGTTGAt                                                       <  8:126718/149‑1 (MQ=255)                                                                                    gCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGAGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAG‑GCGATACCGCTTTCGTTGa                                                        >  2:261462/1‑148 (MQ=255)                                                                                      gTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAG‑Gc                                                                         >  5:131300/1‑129 (MQ=255)                                                                                      gTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAG‑Gc                                                                         <  6:131300/129‑1 (MQ=255)                                                                                         ggCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGTGCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAG‑GCGATACCGCTTTCGTTGATCacac                                                  <  5:252662/149‑1 (MQ=255)                                                                                           cAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAG‑GCGATACCGCTTTCGTTGATCACAccc                                                <  1:164229/149‑1 (MQ=255)                                                                                            aaCATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTAGGTATGCCCGGAAAAG‑GCGATACCGCTTTCGTTGATCACACCCt                                               >  2:256190/1‑149 (MQ=255)                                                                                             aCATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAG‑GCGATACCGCTTTCGTTGATCACACCCtt                                              <  3:225718/149‑1 (MQ=255)                                                                                               aTCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAG‑GCGATACCGCTTTCGTTGATCACACCCTtgt                                            <  4:257021/149‑1 (MQ=255)                                                                                                tCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCCCAGGCCG‑CGGGTGGCCAGCATCAATACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAAATCCGTATGCCCGGAAAAG‑GCGATACCGCTTTCGTTGATCACACCCTtgtc                                           >  6:67165/1‑148 (MQ=255)                                                                                                               tGCGTGGTCACCAGGGCCACGCGCAGGCCT‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGTGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAG‑GCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGTGat                            <  3:308956/149‑1 (MQ=255)                                                                                                               tGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGCGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGTGa                             <  1:261462/149‑1 (MQ=255)                                                                                                                     gTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAG‑GCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGCGATCATccc                      >  8:44839/1‑149 (MQ=255)                                                                                                                          cAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAG‑Gc                                                                         >  3:213654/1‑93 (MQ=255)                                                                                                                          cAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAG‑Gc                                                                         <  4:213654/93‑1 (MQ=255)                                                                                                                             ggCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCAGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAG‑GCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGTGATCATCCCGGCGAAGt              >  5:113153/1‑149 (MQ=255)                                                                                                                             ggCCACGCGCAGGCCG‑CGGGGGGCCAGCATCATTACCACTTGTGCGGGGCGGGTGAGTTGGGCAAGGATTTCGGTATGGCCGGAAAAG‑GCGTTACCGCTTTCGTTGATCAGACCCTTGTGCACCGGGGCCGTGATCATCCCGGCGAAGt              <  5:67165/149‑1 (MQ=255)                                                                                                                             ggCCACCCGCGGCCTG‑CGCGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAG‑GCGATACCGCTTTCGTTGATCACACCCTTGTGCACCggg                                    >  2:241250/1‑127 (MQ=255)                                                                                                                             ggCCACCCGCGGCCTG‑CGCGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAG‑GCGATACCGCTTTCGTTGATCACACCCTTGTGCACCggg                                    <  1:241250/127‑1 (MQ=255)                                                                                                                               ccACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAG‑GCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGCGATCATCCCGGCGAAGTGc            <  7:44839/149‑1 (MQ=255)                                                                                                                                   gcgcAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCg                                                                                           <  7:273622/67‑1 (MQ=255)                                                                                                                                   gcgcAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCg                                                                                           >  8:273622/1‑67 (MQ=255)                                                                                                                                         gCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCTGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAG‑GCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGTGATCATCCCGGCGAAGTGCCCGTCCAGGc  <  6:113153/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                TACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCG‑CGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAG‑GCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGTGATCATCCCGGCGAAGTGCCCGTCCAGGC  >  NC_002947/487695‑487979

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