Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,185,930 +C 93.5% 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.→C93.5% 2.9 ‑0.5 ‑1.146intergenic (+203/‑302)mtlK/RCAP_rcc02028mannitol 2‑dehydrogenase/FeoA family protein
Reads supporting (aligned to +/- strand):  new base (20/23):  ref base (2/1):  total (22/24)
Fisher's exact test strand distribution p-value = 6.00e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.80e-02

TGCGGGGGCTCCG.CCCCGCACCCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTGTTTTTGATGGCACAAGAAACGCCGAACCG.CCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..CCCGCCGGGTCTCCTCCGGGAGCAAGGCCGAGGCGAGACACCC  >  NC_014034/2185829‑2186018
                                                                                                       |                                                                                          
tGCGGGGGCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCTGCTTGCCCCGAAGGCCTTG.GTTTGATGGCACAAGAAACGCCGAACCg                                                                                             <  2:2552738/101‑1 (MQ=25)
tGCGGGGGCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCg                                                                                             <  1:1650831/101‑1 (MQ=37)
  cGGGGGCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCG.cc                                                                                          <  2:4176632/101‑1 (MQ=37)
  cGGGGGCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCG.cc                                                                                          >  2:2101707/1‑101 (MQ=37)
     gggCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCG.ccccc                                                                                       <  1:3041073/101‑2 (MQ=37)
      ggCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCG.ccccct                                                                                      <  1:5353620/101‑3 (MQ=37)
        cgatctccccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCG.CCCCcttt                                                                                    <  1:5273938/95‑2 (MQ=38)
         tccgcccccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCG.CCCCctttc                                                                                   >  1:3924062/6‑99 (MQ=37)
         tccgcccccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCCTGTTTTTG.TGGCAAAAGAAACGCCGAACCG.CCCCctttc                                                                                   >  1:4397715/6‑99 (MQ=25)
          ccgcccccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCa                                                                                   >  2:3966850/5‑101 (MQ=38)
              ccccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCtt                                                                               <  1:5422582/101‑1 (MQ=38)
              ccccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCG.CCCCctttagatcg                                                                              >  2:5273648/1‑94 (MQ=255)
              ccccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAAAGCCGAACCGCCCCCTTTCACCtt                                                                               >  1:3989617/1‑101 (MQ=38)
                ccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCCAACCGCCCCCTTTCACCATgg                                                                             >  2:5709316/1‑101 (MQ=37)
                  cgcacccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAAAAACGCCGAACCGCCCCCTTTCACCTTGGc                                                                            >  1:859876/5‑101 (MQ=255)
                    cacccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCTTGGCTc                                                                          >  2:2041202/3‑101 (MQ=255)
                    cacccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCTTGGCTc                                                                          <  1:962342/99‑1 (MQ=39)
                    cacccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCTTGGCTc                                                                          >  1:5270822/3‑101 (MQ=39)
                     gcatccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCTTGGCTCa                                                                         <  1:2258579/97‑1 (MQ=255)
                     acccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGGCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCTTGGCTCa                                                                         >  2:4268655/2‑101 (MQ=38)
                       ccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCTTGGCTCaaa                                                                       >  1:5002264/1‑101 (MQ=255)
                         tctgCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCTTGGCTCAAata                                                                     <  2:3572270/98‑1 (MQ=255)
                            gCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCTTGGCTCAAataaga                                                                  >  1:3572487/1‑98 (MQ=255)
                             cGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCTTGGCTCAAATATccc                                                                 <  2:4397448/101‑1 (MQ=255)
                             cGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCTTGGCTCAAATATccc                                                                 >  2:5598355/1‑101 (MQ=255)
                              gTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCTTGGCTCAAATATCCCg                                                                <  2:1573167/101‑1 (MQ=255)
                               tagttATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCTTGGCTCAAATATCCCgg                                                               >  1:3297020/4‑101 (MQ=255)
                               tATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCTTGGCTCCAATATCCCgg                                                               >  2:5488185/1‑101 (MQ=255)
                                aTTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCTTGGCTCAAATATCCCggg                                                              <  2:4008703/101‑1 (MQ=255)
                                   tATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCCCCTTGGCTCAAATATCCCgggggg                                                           >  2:4169457/1‑100 (MQ=255)
                                   tATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGt                                                           >  2:324286/1‑101 (MQ=255)
                                     tCCGAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTcc                                                         <  1:4268915/101‑1 (MQ=38)
                                     tCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTcc                                                         <  1:5488475/101‑1 (MQ=255)
                                     tCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTcc                                                         <  2:2336895/101‑1 (MQ=255)
                                     tCCCAAGAAAAAGCAGCGGCTGGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTcc                                                         <  1:5415304/101‑1 (MQ=38)
                                      cccAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCg                                                        >  1:5465545/1‑101 (MQ=255)
                                       ccAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCgg                                                       <  2:859807/101‑1 (MQ=255)
                                       ccAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCgg                                                       <  1:1318127/101‑1 (MQ=255)
                                        cAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCggg                                                      <  2:2412847/101‑1 (MQ=255)
                                         aaGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.GTTTGATGGGACAAGAAACGCCGAACCGCCCCCTTTCACCTTGGCTCAAATATCCCCGGGGTCCgggg                                                     >  2:4290930/1‑101 (MQ=255)
                                          aGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCggggg                                                    <  1:2041319/101‑1 (MQ=255)
                                            aaaaaGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCg                                                  <  2:4387174/101‑1 (MQ=255)
                                            aaaaaGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCg                                                  >  1:4561913/1‑101 (MQ=255)
                                             aaaaGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCgg                                                 >  2:1065633/1‑101 (MQ=255)
                                             aaaaGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGGCCGGGGGCgg                                                 >  2:2196355/1‑101 (MQ=38)
                                              aaaGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCCCCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGa                                                <  1:1976129/101‑1 (MQ=38)
                                              aaaGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGa                                                >  1:2400417/1‑101 (MQ=255)
                                              aaaGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGa                                                <  2:756933/101‑1 (MQ=255)
                                                 gcagcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..cc                                           <  1:4696570/101‑1 (MQ=255)
                                                    gcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..ccccc                                        <  1:812045/101‑3 (MQ=37)
                                                         ttGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCCCCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..cccccggcgg                                   <  1:2101824/101‑8 (MQ=25)
                                                         ttGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..cccccggcgg                                   >  2:3585605/1‑94 (MQ=25)
                                                            ccccGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCCCCTTGGCTCAAATACCCCGGGGGTCGGGGGGCGGAG..cccccggcgggtc                                <  2:2361569/101‑11 (MQ=18)
                                                                                           accgcgcaaCCGCCCCCTTTCCCCTTGGCTCAAATATCCGGGGGGTCCGGGGGCCGAGCCCCCGGCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACAccc  <  2:696448/94‑1 (MQ=2)
                                                                                                       |                                                                                          
TGCGGGGGCTCCG.CCCCGCACCCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTGTTTTTGATGGCACAAGAAACGCCGAACCG.CCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..CCCGCCGGGTCTCCTCCGGGAGCAAGGCCGAGGCGAGACACCC  >  NC_014034/2185829‑2186018

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 24 ≤  ATCG  < 32 ≤  ATCG  < 40 ≤  ATCG