Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,655,631 T→G 11.0% intergenic (‑86/+66) epd ← / ← tktA D‑erythrose 4‑phosphate dehydrogenase/transketolase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,655,6310TG11.0% 114.7 / 8.8 46intergenic (‑86/+66)epd/tktAD‑erythrose 4‑phosphate dehydrogenase/transketolase A
Reads supporting (aligned to +/- strand):  ref base T (16/23);  new base G (2/3);  total (20/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.25e-01

CAGGACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCG  >  NC_002947/5655483‑5655770
                                                                                                                                                      |                                                                                                                                           
cAGGACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGttt                                                                                                                                             <  5:257994/149‑1 (MQ=255)
 aGGACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTc                                                                                                                                            <  3:142128/149‑1 (MQ=255)
 aGGACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTc                                                                                                                                            <  3:315966/149‑1 (MQ=255)
   gACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCAt                                                                                                                                          <  8:72213/149‑1 (MQ=255)
     cacaGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGc                                                                                                                                        <  1:370462/149‑1 (MQ=255)
           ggCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAAccc                                                                                                                                   <  2:220740/149‑2 (MQ=255)
            tccGATACGGCCATAACCGTTTAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACTAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCTGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACctct                                                                                                                                 <  2:349800/148‑1 (MQ=255)
            gCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACCGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAAccct                                                                                                                                  <  5:356724/149‑3 (MQ=38)
              cGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCtgt                                                                                                                               >  2:357382/1‑149 (MQ=255)
                aTACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTgg                                                                                                                             <  1:357382/149‑1 (MQ=255)
                     gCCATAACCGTTGAGTTCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCg                                                                                                                        <  3:89881/149‑1 (MQ=255)
                          aaCCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGTCATTCGCGGGCAAGCCCCCACCCACAGGGTTTGATGGAACCTCTGTGGGAGCGGGcgt                                                                                                                   >  5:28055/1‑147 (MQ=255)
                                    tGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAAt                                                                                                         <  7:192360/149‑1 (MQ=255)
                                     gCAATTTTGTTGTGACGGTTGGGCATTCGTGTCTAACTAGACAGTTCAATTGTCTTGAGCCTGTGTAATTTTCGGGGTGGCATTCCGGGGCAAGGCCTCTGCATCAGGGTTTCTTTCATTGTCTGTGTGTGCGGGCTTGCCCGCGAAt                                                                                                         <  2:99468/148‑1 (MQ=255)
                                      cAACTTTGTAGGGACGGTGGGGCATTCGCGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGCGCAAGACCGCTACCAAATGGTTACATTCAACCTATGTGTGTGCGGGCTTGCCAGCGAAtt                                                                                                        >  1:99468/1‑147 (MQ=255)
                                       aaCTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATg                                                                                                        >  7:272772/1‑147 (MQ=255)
                                         cTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATg                                                                                                        >  8:212803/1‑145 (MQ=255)
                                         cTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCCGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATg                                                                                                        <  7:212803/145‑1 (MQ=255)
                                               aGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGc                                                                                                                 <  7:376220/130‑1 (MQ=255)
                                               aGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGc                                                                                                                 >  8:376220/1‑130 (MQ=255)
                                                gggACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCtgtg                                                                                                                              <  7:164190/116‑1 (MQ=255)
                                                gggACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCtgtg                                                                                                                              >  8:164190/1‑116 (MQ=255)
                                                        ggggCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAg                                                                                     <  8:272772/149‑1 (MQ=255)
                                                          ggCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTGTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGa                                                                                    >  3:45099/1‑148 (MQ=255)
                                                           gCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGa                                                                                                                           <  8:136199/108‑1 (MQ=255)
                                                           gCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATc                                                                                  <  4:45099/149‑1 (MQ=255)
                                                                          gagaCAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGGGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTc                                                                   <  1:389267/149‑1 (MQ=255)
                                                                                 ttCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCa                                                                         <  1:266883/136‑1 (MQ=255)
                                                                                 ttCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCa                                                                         >  2:266883/1‑136 (MQ=255)
                                                                                    aaTTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACAccgcccgc                                                                                       >  4:57778/1‑119 (MQ=255)
                                                                                    aaTTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACAccgcccgc                                                                                       <  3:57778/119‑1 (MQ=255)
                                                                                               gCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTAcc                                              >  1:322170/1‑149 (MQ=255)
                                                                                                ccTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCcagc                                                               <  7:205708/131‑1 (MQ=255)
                                                                                                ccTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCcagc                                                               >  8:205708/1‑131 (MQ=255)
                                                                                                  tGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATAAGTCTTCCAGCAGCTCTTAGGCAGTAACGAg                                           >  4:312350/1‑149 (MQ=255)
                                                                                                     ggAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCg                                                      >  1:309520/1‑135 (MQ=255)
                                                                                                     ggAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCg                                                      <  2:309520/135‑1 (MQ=255)
                                                                                                         tttGCGGTGTGGCATTCGCGGGCAAGCCCCCACCCACATGGGTTACAGCAACCTCAGTGGGAGCGGGGTTGCCCGCCAATTCAAGCGCCACCACCCCCACATCAGTAAATAAGTCTTCAAGAACCTCTTAGGCCGATACGAGGATGTtc                                    >  8:355181/1‑149 (MQ=255)
                                                                                                                     cATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAg                        <  3:296921/149‑1 (MQ=255)
                                                                                                                                 aGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCg            <  7:355181/149‑1 (MQ=255)
                                                                                                                                    ccGCTCCCACAGAGGTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGaa          >  3:66083/1‑148 (MQ=37)
                                                                                                                                       cTCCCACAGAGGTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAAt      <  4:66083/149‑1 (MQ=38)
                                                                                                                                          ccACAGGGTTTCATGCAACCTATGTGGGAGCGCGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATAAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGAAGATCTCCAGGGTGAAGCAAAACTCATCGAACAATGcc   >  4:331729/1‑149 (MQ=255)
                                                                                                                                            cacagaggTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGcc   >  1:152224/7‑148 (MQ=38)
                                                                                                                                            acaGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCg  >  4:149701/1‑148 (MQ=255)
                                                                                                                                                  ggTTTCATGCAACCTCTGTGGGAGGGGGCTTGCCCGCGAATGAAGGCGACACCGCCCTCAGATCttt                                                                               <  6:131978/67‑1 (MQ=38)
                                                                                                                                                  ggTTTCATGAAACCTCTGTGGGAGCGGGCAAGCCCGCGAATGAAGGCGACACCGCCCGCAGATCttt                                                                               >  5:131978/1‑67 (MQ=37)
                                                                                                                                                      |                                                                                                                                           
CAGGACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCG  >  NC_002947/5655483‑5655770

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