Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 395,848 T→G 11.4% intergenic (+30/+174) 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,8480TG11.4% 141.0 / 9.0 51intergenic (+30/+174)fdhA/PP_0329formaldehyde dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (22/23);  new base G (2/4);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

TGGTAGGTGTGCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTGGAGCGAAGGGCCACTGCG  >  NC_002947/395703‑395996
                                                                                                                                                 |                                                                                                                                                    
tGGTAGGTGTGCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGa                                                                                                                                                   <  1:230889/149‑1 (MQ=255)
  gTAGGTGTGCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAgg                                                                                                                                                 >  8:103784/1‑149 (MQ=255)
  gTAGGTGTGCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAg                                                                                                                                                  <  4:228221/148‑1 (MQ=255)
     ggTGTGCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAgggg                                                                                                                                               <  3:12495/148‑1 (MQ=255)
       tgtgCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGATGCGAGGGGCtc                                                                                                                                            <  5:322385/149‑2 (MQ=255)
         tgCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGTTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCttt                                                                                                                                           <  4:99052/148‑1 (MQ=255)
             ggTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGATGCGAGGGGCtcttttta                                                                                                                                      <  3:83306/149‑4 (MQ=255)
                gATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTTTCCAt                                                                                                                                   >  4:272256/1‑149 (MQ=255)
                   cAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCt                                                                                                                                >  4:139728/1‑149 (MQ=255)
                                  gcCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCgg                                                                                                                 <  3:311032/149‑1 (MQ=255)
                                       aaGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGa                                                                                                            <  7:143277/149‑1 (MQ=255)
                                        aGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAg                                                                                                           <  2:37472/149‑1 (MQ=255)
                                              tGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCtt                                                                                                     >  1:272916/1‑149 (MQ=255)
                                                gCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTg                                                                                                    >  6:151780/1‑148 (MQ=255)
                                                gCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGt                                                                                                   <  7:103784/149‑1 (MQ=255)
                                                          tGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCtt                                                                                         <  6:90299/149‑1 (MQ=255)
                                                                tgtgCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCg                                                                                   <  6:46021/149‑1 (MQ=255)
                                                                      gaagaaATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTg                                                                                                    <  2:71626/126‑1 (MQ=255)
                                                                      gaagaaATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTg                                                                                                    >  1:71626/1‑126 (MQ=255)
                                                                        agaaATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGaa                                                                           <  5:151780/149‑1 (MQ=255)
                                                                             ttCGTTATCGACCCGCACAGAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTc                                                                        >  5:86402/1‑147 (MQ=255)
                                                                                       aCCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGccgcccg                                                            >  4:297756/1‑149 (MQ=255)
                                                                                        cccGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCCAGGGGGTTGTTGTCCATTCTGATAAGACAACTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGATCGccgcccg                                                            >  7:198905/1‑148 (MQ=255)
                                                                                                aaaTGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCt                                                   >  5:105844/1‑149 (MQ=255)
                                                                                                  aTGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCtca                                                 >  7:70993/1‑149 (MQ=255)
                                                                                                      ggggTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTATCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAg                                             <  6:316157/149‑1 (MQ=255)
                                                                                                          tGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGccgcccg                                                            >  5:222015/1‑130 (MQ=255)
                                                                                                          tGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGccgcccg                                                            <  6:222015/130‑1 (MQ=255)
                                                                                                          tGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCACTCGCATCGAGGGGCTTTTTGTCGATTCTGATAAGACATCTCGGCCGGAGGTGGTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTcc                                         >  1:337096/1‑149 (MQ=255)
                                                                                                             ggcgTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCTGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGaa                                      <  1:136632/149‑1 (MQ=255)
                                                                                                              gcgTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAg                                     >  4:62167/1‑149 (MQ=255)
                                                                                                               cgTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGc                                    <  3:297756/149‑1 (MQ=255)
                                                                                                               cgTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGc                                    <  3:62167/149‑1 (MQ=255)
                                                                                                                   aTCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTgg                                <  8:70993/149‑1 (MQ=255)
                                                                                                                    tCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACTTACCTCATCAg                                             <  7:147286/135‑1 (MQ=255)
                                                                                                                    tCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACTTACCTCATCAg                                             >  8:147286/1‑135 (MQ=255)
                                                                                                                    tCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGt                               <  6:311282/149‑1 (MQ=255)
                                                                                                                    tCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGt                               <  6:105844/149‑1 (MQ=255)
                                                                                                                        cGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCTTCAGTTCCGAAGCCTGGTGAGt                           >  8:78170/1‑149 (MQ=255)
                                                                                                                         gTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTg                          <  3:139728/149‑1 (MQ=255)
                                                                                                                         gTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTg                          >  7:102731/1‑149 (MQ=255)
                                                                                                                                 caaaaAGCCCCTCGATGCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCtt                                                                                         <  6:35999/77‑1 (MQ=37)
                                                                                                                                 caaaaAGCCCCTCGATGCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCtt                                                                                         >  5:35999/2‑78 (MQ=37)
                                                                                                                                   aaGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTGGAGc                >  3:288352/1‑149 (MQ=255)
                                                                                                                                    aaagCCCCTCGATGCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCtttt                                                                              <  3:30179/84‑1 (MQ=39)
                                                                                                                                    aaagCCCCTCGATGCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCtttt                                                                              >  4:30179/3‑86 (MQ=39)
                                                                                                                                        ccccTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGtt                                           >  6:91038/1‑117 (MQ=255)
                                                                                                                                        ccccTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGtt                                           <  5:91038/117‑1 (MQ=255)
                                                                                                                                        ccccTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTGGAGCGAAgg           >  6:240288/1‑149 (MQ=255)
                                                                                                                                          ccTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTGGAGCGAAGGGc         >  7:67267/1‑149 (MQ=255)
                                                                                                                                              gCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTGGAGCGAAGGGCCACt     >  7:179840/1‑149 (MQ=255)
                                                                                                                                                 tCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTGGAGCGAAGGGCCACTgcg  <  3:272256/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TGGTAGGTGTGCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTGGAGCGAAGGGCCACTGCG  >  NC_002947/395703‑395996

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