Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 395,847 A→T 13.4% intergenic (+29/+175) fdhA → / ← PP_0329 formaldehyde dehydrogenase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947395,8470AT13.4% 129.3 / 10.0 51intergenic (+29/+175)fdhA/PP_0329formaldehyde dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (24/20);  new base T (3/4);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

GGTGTGCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTGGAGCGAAGGGCCACTG  >  NC_002947/395708‑395994
                                                                                                                                           |                                                                                                                                                   
ggTGTGCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGc                                                                                                                                            <  2:3367/149‑1 (MQ=255)
ggTGTGCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGATGCGAGGGGc                                                                                                                                            <  2:80319/149‑1 (MQ=255)
    tgCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCttt                                                                                                                                         <  4:138259/148‑1 (MQ=255)
     gCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCttttt                                                                                                                                       >  4:32858/1‑149 (MQ=255)
      cAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGATGCGAGGGGCtctttt                                                                                                                                      <  3:31901/149‑2 (MQ=255)
       aGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGt                                                                                                                                     <  4:49812/149‑1 (MQ=255)
        ggTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTc                                                                                                                                    <  2:302645/149‑1 (MQ=255)
        ggTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTc                                                                                                                                    >  2:84941/1‑149 (MQ=255)
        ggTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTAGCATCGAGGGGCTTTTTGTc                                                                                                                                    >  6:30775/1‑149 (MQ=255)
              cAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCAGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTc                                                                                                                               >  7:70254/1‑148 (MQ=255)
                aCCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGa                                                                                                                            >  1:245799/1‑149 (MQ=255)
                 ccTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGAt                                                                                                                           <  5:30775/149‑1 (MQ=255)
                 ccTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGCGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGa                                                                                                                            >  2:113745/1‑148 (MQ=255)
                               cGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAATAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGcc                                                                                                             <  1:84941/149‑1 (MQ=255)
                               cGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAgggg                                                                                                                                             >  6:59728/1‑117 (MQ=255)
                                       tATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTgt                                                                                                     <  3:32858/149‑2 (MQ=255)
                                            cGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTc                                                                                                                                    <  2:246576/113‑1 (MQ=255)
                                            cGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTc                                                                                                                                    >  1:246576/1‑113 (MQ=255)
                                                tttGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCAt                                                                                              >  8:243904/1‑147 (MQ=255)
                                                 ttGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATaa                                                                                            <  8:115359/148‑1 (MQ=255)
                                                    aTGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCt                                                                                        <  7:243904/149‑1 (MQ=255)
                                                                cGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCtttt                                                                            >  8:44370/1‑149 (MQ=255)
                                                                 gaagaaATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCtttt                                                                            >  6:180275/1‑148 (MQ=255)
                                                                 gaagaaATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCtttt                                                                            >  5:316509/1‑148 (MQ=255)
                                                                 gaagaaATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCttga                                                                            >  1:114213/1‑146 (MQ=255)
                                                                 gaagaaATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTg                                                                           <  8:70254/149‑1 (MQ=255)
                                                                    gaaATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACAtc                                                                                                                   >  8:110996/1‑106 (MQ=255)
                                                                    gaaATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACAtc                                                                                                                   <  7:110996/106‑1 (MQ=255)
                                                                           gTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAgcg                                                                 <  3:151952/149‑1 (MQ=255)
                                                                            ttATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAgcgc                                                                >  7:108875/1‑149 (MQ=255)
                                                                               tCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGcgccgc                                                             >  2:122385/1‑149 (MQ=255)
                                                                                    ccGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTc                                                                                                <  2:278843/109‑1 (MQ=255)
                                                                                    ccGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTc                                                                                                >  1:278843/1‑109 (MQ=255)
                                                                                      gCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTAATCGCTTTTGAAGTCGAGCGCCCCCCGcaca                                                      >  4:75972/1‑149 (MQ=255)
                                                                                         caAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCTGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATAc                                                   <  6:259699/149‑1 (MQ=255)
                                                                                                tgGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCtcatca                                            <  6:316509/149‑1 (MQ=255)
                                                                                                    gTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTc                                        >  4:270014/1‑149 (MQ=255)
                                                                                                             aaTCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTg                               >  6:249454/1‑149 (MQ=255)
                                                                                                              aTCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTgg                              <  3:270014/149‑1 (MQ=255)
                                                                                                               tCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGt                             <  8:42743/149‑1 (MQ=255)
                                                                                                                 aGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGa                           <  2:245799/149‑1 (MQ=255)
                                                                                                                  gCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAg                          >  5:67644/1‑149 (MQ=255)
                                                                                                                   cGTCGAGTAAAAAGAGCCCCTCGCATCGTGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGt                         >  1:180984/1‑149 (MQ=255)
                                                                                                                       gAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGcc                     >  3:221321/1‑149 (MQ=255)
                                                                                                                        tggacaaaaAGCCCCTCGATGCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGccc                    >  7:193131/6‑149 (MQ=255)
                                                                                                                         gTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCt                   >  6:243209/1‑149 (MQ=255)
                                                                                                                             aaaaatccccTCGATGCGATGGGCTTTTTGTCCTTTCTGATAAGACATCTCGTCCTGATGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGAGCCGCCCGCACATACCTCATCAGtt                                         <  4:227012/117‑1 (MQ=255)
                                                                                                                             aaaaAGCCCCTCGATGCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAttt                                         >  3:227012/1‑123 (MQ=255)
                                                                                                                             aaaaAGCCCCTCGATGCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTGGAg               <  8:193131/149‑1 (MQ=255)
                                                                                                                               aaagCCCCTCGATGCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTGGAGc              >  8:126186/3‑148 (MQ=255)
                                                                                                                                          cATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTGGAGCGAAGGGCCACTg  <  5:180275/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
GGTGTGCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTGGAGCGAAGGGCCACTG  >  NC_002947/395708‑395994

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