Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,185,930 +C 89.9% intergenic (+203/‑302) mtlK → / → RCAP_rcc02028 mannitol 2‑dehydrogenase/FeoA family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,185,9301.→C89.9% 16.9 ‑2.9 ‑0.369intergenic (+203/‑302)mtlK/RCAP_rcc02028mannitol 2‑dehydrogenase/FeoA family protein
Reads supporting (aligned to +/- strand):  new base (34/28):  ref base (1/6):  total (35/34)
Fisher's exact test strand distribution p-value = 5.48e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.27e-01

TGCGGGGGCTCCG.CCCCGCACCCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTGTTTTTGATGGCACAAGAAACGCCGA.ACCG.CCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAGCCCGCCGGGTCTC  >  NC_014034/2185829‑2185988
                                                                                                        |                                                          
tGCGGGGGCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCg                                                             <  2:4204848/101‑1 (MQ=37)
  cGGGGGCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCG.cc                                                          <  2:3452597/101‑1 (MQ=37)
   gggggCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCCA.ACCG.ccc                                                         >  2:1462133/1‑101 (MQ=37)
       gCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCG.CCCCctt                                                     <  1:1979636/101‑1 (MQ=37)
       gCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCG.CCCCctt                                                     <  2:3493932/101‑1 (MQ=37)
       gCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCG.CCCCctt                                                     <  2:2433821/101‑1 (MQ=37)
        cTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCG.CCCCcttt                                                    <  2:4077923/101‑2 (MQ=37)
        cTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCG.CCCCcttt                                                    <  1:1960695/101‑2 (MQ=37)
         tccgcccccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCG.CCCCctttc                                                   >  1:3418966/6‑99 (MQ=37)
         tccgcccccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCG.CCCCctttc                                                   <  2:274819/96‑3 (MQ=37)
           cgcccccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCAc                                                  <  1:3691128/98‑1 (MQ=38)
           cgcccccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCAc                                                  >  1:2831382/4‑101 (MQ=38)
            gcccccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCCCAAGAAACGCCGA.ACCGCCCCCTTTCAcc                                                 >  1:3387795/3‑101 (MQ=38)
             cccccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCt                                                <  1:3415297/100‑1 (MQ=38)
             cccccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCt                                                <  2:3250880/100‑1 (MQ=38)
              ccccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCtt                                               <  2:3947075/101‑1 (MQ=38)
               cccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGGACAAGAAACGCCGA.ACCGCCCCCTTTCCCCTTg                                              >  1:4051642/1‑101 (MQ=255)
               cccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTg                                              >  1:3712301/1‑101 (MQ=38)
                cccgcgcccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTg                                              <  1:1602509/95‑1 (MQ=38)
                  cgcacccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTGTTTTT.ATGGCACAAAAAAAGCCGA.ACCGCCCCCTTTCACCCTGGc                                            >  2:3666089/5‑101 (MQ=255)
                  cgcacccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGc                                            >  2:3370271/5‑101 (MQ=255)
                   gcacccccGGCGTATTTATCCCAAGAAAAAGCATCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCt                                           >  1:2030297/4‑101 (MQ=255)
                   gcacccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCt                                           >  1:3013979/4‑101 (MQ=255)
                   gcacccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCt                                           >  1:1829755/4‑101 (MQ=255)
                    cacccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTTCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGTCTc                                          >  1:1387385/3‑101 (MQ=38)
                    cacccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTc                                          >  1:3602962/3‑101 (MQ=255)
                      cccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCaa                                        >  1:3731836/1‑101 (MQ=255)
                      cccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCaa                                        <  1:3957476/101‑1 (MQ=255)
                       ccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCaaa                                       >  2:1173591/1‑101 (MQ=255)
                       ccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCaaa                                       <  1:3370477/101‑1 (MQ=255)
                        cccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAat                                      <  2:3087042/101‑1 (MQ=255)
                        cccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAat                                      <  1:2398340/101‑1 (MQ=255)
                          cGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAatat                                    <  1:2957452/101‑1 (MQ=255)
                           ggCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATc                                   >  2:1590716/1‑101 (MQ=255)
                           ggCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATc                                   >  2:1579731/1‑101 (MQ=255)
                           ggCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATc                                   <  2:1314124/101‑1 (MQ=255)
                            gCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATcc                                  >  1:4748640/1‑101 (MQ=255)
                             cGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATccc                                 <  2:546385/101‑1 (MQ=255)
                             cGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATccc                                 <  1:3319079/101‑1 (MQ=255)
                              gTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATCCCg                                <  2:4506427/101‑1 (MQ=255)
                              gTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATCCCg                                >  1:942191/1‑101 (MQ=255)
                              gTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATCCCg                                >  1:420665/1‑101 (MQ=255)
                               tATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATCCCgg                               >  2:3443214/1‑101 (MQ=255)
                               tATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATCCCgg                               >  2:638831/1‑101 (MQ=255)
                                aTTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATCCCggg                              <  2:2484188/101‑1 (MQ=255)
                                aTTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATCCCggg                              >  2:255395/1‑101 (MQ=255)
                                aTTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATCCCggg                              <  1:1579831/101‑1 (MQ=255)
                                  ctATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATCCCggggg                            <  2:1526171/100‑1 (MQ=39)
                                   tATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGt                           >  1:4350656/1‑101 (MQ=255)
                                   tATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGt                           >  1:3030817/1‑101 (MQ=255)
                                   tATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGt                           >  1:2390317/1‑101 (MQ=255)
                                    aTCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCCCCTTGGCTCAAATATCCCggggggg                          >  1:1526272/1‑99 (MQ=255)
                                    aTCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTc                          <  1:3817493/101‑1 (MQ=255)
                                     tCCCAAGAAAAAGCATCGGCATGCCCCGAAGGGCTTGTTTTT.ATGGGACAAGAAAAGCCAA.ACAGCCCCCCCTCACCTTGGCTCAAATATCCCGGGGGTcc                         >  2:330770/1‑101 (MQ=19)
                                     tCCCAAGAAAAAGCAGCGGCTTTCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCGGA.ACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTcc                         <  2:3679886/101‑1 (MQ=38)
                                     tCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTcc                         <  2:2831217/101‑1 (MQ=255)
                                     tCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTcc                         >  1:355939/1‑101 (MQ=255)
                                     tCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTcc                         <  2:3598873/101‑1 (MQ=255)
                                      cccAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCg                        >  2:2994983/1‑101 (MQ=255)
                                         aaGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCgggg                     <  1:2465401/101‑1 (MQ=255)
                                         aaGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCgggg                     <  2:3395402/101‑1 (MQ=255)
                                         aaGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCgggg                     >  1:199110/1‑101 (MQ=255)
                                            aaaaaGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCg                  >  2:2921070/1‑101 (MQ=255)
                                            aaaaaGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGGCCGGGGGCg                  >  2:124648/1‑101 (MQ=38)
                                             aaaaGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCgg                 >  2:4312070/1‑101 (MQ=255)
                                              aaaGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCggc                >  2:2003376/1‑100 (MQ=38)
                                                   agcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAGcccc           <  1:664908/101‑2 (MQ=37)
                                                   agcGGCTTGCCCCGAAAGACTTG.TTTTGATGGCACAAGAAATACCCA.AACGCCCCCTTTCACCTTTGCTCAAAAATGCCGGGGGGACGGGGGCGGAgggcc           >  2:1992097/1‑97 (MQ=255)
                                                    tcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAGccccc          <  2:2030180/100‑3 (MQ=37)
                                                    gcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAGccccc          <  1:4261302/101‑3 (MQ=37)
                                                     cgtcgtGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAGcccccg         <  2:706048/96‑4 (MQ=25)
                                                     cGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAGcccccg         <  1:3912014/101‑4 (MQ=37)
                                                         ttGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAGcccccggcgg     >  1:4522225/1‑94 (MQ=25)
                                                          tGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAGcccccggcggg    >  1:2443020/1‑93 (MQ=21)
                                                           cccccGAATGCCTTG.TTTTGATGGCACAAGAAACCCCGACCCCC.CCCCTTTCCCCTGGGCTAAAATTCCCGGGGGGTCGGGGGGCGGAGcccccggcgggt   <  2:3418760/100‑10 (MQ=12)
                                                            ccccGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGA.ACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAGcccccggcggggg  >  2:1550205/1‑91 (MQ=21)
                                                                                                        |                                                          
TGCGGGGGCTCCG.CCCCGCACCCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTGTTTTTGATGGCACAAGAAACGCCGA.ACCG.CCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAGCCCGCCGGGTCTC  >  NC_014034/2185829‑2185988

Base quality scores:  ATCG  < 3 ≤  ATCG  < 19 ≤  ATCG  < 28 ≤  ATCG  < 34 ≤  ATCG  < 40 ≤  ATCG