Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,995,659 +A 5.0% intergenic (+199/+21) RCAP_rcc02802 → / ← RCAP_rcc02803 membrane dipeptidase/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,995,6592.→A 5.0% 2.8 0.0 ‑0.558intergenic (+199/+21)RCAP_rcc02802/RCAP_rcc02803membrane dipeptidase/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (42/13):  total (45/13)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.42e-01

GCCGCCCCGCAACCCGTCGGCGGCCCAGCGCCCTGCTTTTCCTTGGCAAAAT.ACG.CCGGGGGGTGCGG..GGGGCGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAATT  >  NC_014034/2995560‑2995759
                                                                                                          |                                                                                                     
ccccccccCCACCCCGTGGGCGCCCCAGCCCCCTTTTTTCCCGGGGCAAAAT.CCG.CGGGGGGGTGGGG..GGGGGGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..g                                                                                                     <  1:1419898/97‑1 (MQ=12)
             ccGTCGGCGGCCCAGCGCCCTGCTTTTCCTTGGCAAAAT.ACG.CCGGGGGGTGCGG..GGGGCGG.AGCCCCCCG.CGGGGCGCCCCGGCTT..GCCGGGGCGGAATc                                                                                        >  2:4057193/1‑101 (MQ=255)
             ccGCCGGCGCCCCAGCCCCCTGCTTTTCCTGGGAAAAAT.GGG.CGGGGGGGTGGGG..GGGGGGG.ACCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATc                                                                                        <  1:4601521/101‑1 (MQ=18)
                gccgcggcCCCGGCCCCTTTTTTTCCTTGGGAAAATCGCG.CGGGGGGGGGGGG..GGGGGGG.CCCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATccc                                                                                      <  2:2570271/98‑1 (MQ=16)
                     ggcCCAGCGCCCTGCTTTTCCTTGGCAAAAT.ACG.CCGGGGGGTGCGG..GGGGCGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAAGCCCCCGCTc                                                                                >  1:686689/1‑101 (MQ=255)
                            cccccTGCTTTTCCTTGGCAAAAT.CCG.CGGGGGGGTGGGG..GGGGCGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCa                                                                         <  2:2924081/99‑1 (MQ=255)
                                cTGCTTTTCCTTGGCAAAAT.ACG.CCGGGGGGTGCGG..GGGGCGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAgaga                                                                     >  1:632652/1‑101 (MQ=255)
                                 tGCTTTTCCTTGGCAAAAT.ACG.CCGGGGGGTGCGG..GGGGCGG.AGCCCCCCG.CGGGTCTCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCagaggc                                                                    >  2:3927275/1‑99 (MQ=255)
                                   tttttCCTTGGAAAAAT.ACGCCCGGGGGGTGCGG..GGGGCGG.AGCCCCCCGCCGGGTCGCCTCGGCGTCAGCCGAGGCGAAATACCCC.CTCAGCccttcatc                                                                     >  2:2611488/2‑95 (MQ=0)
                                      ttCCTTGGAAAAAT.ACGCCCGGGGGGTGCGG..GGGGCGG.AGCCCCCCGCCGGGTCGCCTCGGCGTCCGCCGAGGCGAAATACCCC.CTCAGCccttcatcttg                                                                  >  1:789505/1‑92 (MQ=0)
                                      ttCCTTGGAAAAAT.ACGCCCGGGGGGTGCGG..GGGGCGG.AGCCCCCCGCCGGGTCGCCTCGGCGTCAGCCGAGGCGAAATACCCC.CTCAGCccttcatcttg                                                                  >  2:4078231/1‑92 (MQ=0)
                                            ggCAAAAT.ACG.CCGGGGGGTGCGG..GGGGCGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGAccccc                                                         >  1:1263423/1‑98 (MQ=255)
                                                aaaT.ACG.CCGGGGGGTGCGG..GGGGCGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAAcc                                                     >  1:1094053/1‑101 (MQ=255)
                                                aaaT.ACG.CCGGGGGGTGCGG..GGGGCGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAAGCCACCCCTCCGTCCCAGAGCCTTTCGACCTCGAAcc                                                     >  1:4107436/1‑101 (MQ=25)
                                                     aCG.CCGGGGGGTGCGG..GGGGCGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCCGAGACTTTCGACCTCGAACCGGGt                                                 >  2:675890/1‑101 (MQ=255)
                                                      cG.CCGGGGGGTGGGG..GGGGCGG.CGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTc                                                <  1:3787334/101‑1 (MQ=255)
                                                       cgctgcggggTGCGGGCGGGGCGGTCGCCCCCCGCCGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCg                                                    <  1:19759/95‑1 (MQ=255)
                                                           ggggggGGGGG..GGGGGGC.CCCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGAc                                            <  1:1652402/101‑1 (MQ=38)
                                                           ggggggCGGGG..GGGGGGC.CCCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGAc                                            <  1:3013837/101‑1 (MQ=38)
                                                              gggTGCGG..GGGGCGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCCGAcgcg                                         >  2:4639260/1‑101 (MQ=255)
                                                              gggTGCGG..GGGGCGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGAcgcg                                         >  2:1120709/1‑101 (MQ=255)
                                                               ggTGCGG..GGGGCGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCgg                                        >  1:1332710/1‑101 (MQ=255)
                                                                gTGCGG..GGGGAGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GGCGGGGCGAAACCCCCGGCTCCGTCCCAGAGGCTTTGGACCTCGGACCGGGCCACAccccga                                       >  2:1069708/1‑96 (MQ=21)
                                                                   cGG..GGGGCGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATgc                                    >  1:2073503/1‑101 (MQ=255)
                                                                   cGG..GGGGCGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATgc                                    >  1:3848689/1‑101 (MQ=255)
                                                                    gg..ggggCGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGAGTTTCGACCTCGAACCGGGTCAGACGCGGATgcg                                   >  1:2287718/1‑101 (MQ=255)
                                                                    gg..ggggCGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATgcg                                   >  1:2541229/1‑101 (MQ=255)
                                                                         gggcagggcccccccG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCa                                 <  1:4400751/92‑1 (MQ=255)
                                                                          ggCGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAgg                               >  1:3595288/1‑101 (MQ=255)
                                                                           gCGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTGTCGACCTCGAACCCGGGCACACGCGGATGGGCGGGc                              >  2:3235161/1‑101 (MQ=38)
                                                                            cGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCg                             >  2:4360965/1‑101 (MQ=255)
                                                                            cGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCg                             >  2:311191/1‑101 (MQ=255)
                                                                            cGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCg                             >  1:709107/1‑101 (MQ=255)
                                                                            cGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCg                             >  2:956431/1‑101 (MQ=255)
                                                                                 gCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGcc                         >  2:278433/1‑101 (MQ=255)
                                                                                 gCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGGCAGACGCGGATGCGCAGGCGGGcc                         >  2:3955297/1‑101 (MQ=255)
                                                                                  ccccccG.CGGGTCGCACCGGCTT..GCCGAGGCGAAATCCCCCGCTCAGTCCCATAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCg                        >  2:3157541/1‑101 (MQ=255)
                                                                                   cccccG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGa                       <  1:311225/101‑1 (MQ=255)
                                                                                    ccccG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGCCAGACGGGGAGGCGCAGGGGGGcgcaa                      >  1:928986/1‑97 (MQ=38)
                                                                                     cccG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGaaa                     <  1:1957842/101‑1 (MQ=255)
                                                                                        g.cgGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGcc                  <  1:2959345/101‑1 (MQ=255)
                                                                                         gggggTCGCCCCGGCTT..GCCGGGGCGAAATCCCGAGCCCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGcc                  <  1:4057435/99‑1 (MQ=255)
                                                                                          cgGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCg                 >  2:1835933/1‑101 (MQ=255)
                                                                                          cgGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCg                 >  2:3109657/1‑101 (MQ=255)
                                                                                          cgGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCg                 >  2:2508254/1‑101 (MQ=255)
                                                                                          cgGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCg                 >  1:1359693/1‑101 (MQ=255)
                                                                                          cgGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCg                 >  1:3876056/1‑101 (MQ=255)
                                                                                          cgGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCg                 >  2:2708310/1‑101 (MQ=255)
                                                                                           gggTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGc                >  1:3070971/1‑101 (MQ=255)
                                                                                           gggTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGAGCAGGCGGGCCGAAAGCCGc                >  2:1635726/1‑101 (MQ=255)
                                                                                            ggTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCa               >  2:4623684/1‑101 (MQ=255)
                                                                                              tCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAg             >  2:2703573/1‑101 (MQ=255)
                                                                                              tCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAg             >  2:2216160/1‑101 (MQ=255)
                                                                                              tCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAg             >  1:2713617/1‑101 (MQ=255)
                                                                                                   ccGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTgg        >  1:508535/1‑101 (MQ=255)
                                                                                                    cGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTggg       >  1:3754833/1‑100 (MQ=255)
                                                                                                    cGGCTT..GCCGGGGCGAAACCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGc       <  2:1186776/101‑1 (MQ=255)
                                                                                                    cGGCTT..CCCGGGGCGAAATCCCCCGCTCAGTCCCGGAAACTTTCGACCTCGCACCGAGTCAGTCGCGGATAGGCAGGCGGGCCCAAAGCCGCAAGGCTggt       >  2:726118/1‑100 (MQ=255)
                                                                                                     ggCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCg      <  2:3848447/101‑1 (MQ=255)
                                                                                                     ggCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCg      >  1:1946509/1‑101 (MQ=255)
                                                                                                      gCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGa     >  1:1780596/1‑101 (MQ=255)
                                                                                                      gCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGa     <  2:1637647/101‑1 (MQ=255)
                                                                                                      gCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGa     >  1:4689760/1‑101 (MQ=255)
                                                                                                       cTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGaa    >  2:744517/1‑101 (MQ=255)
                                                                                                        tt..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAAt   <  2:964186/101‑1 (MQ=255)
                                                                                                        tt..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAAt   >  1:1217248/1‑101 (MQ=255)
                                                                                                         t..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAAtt  >  1:4257252/1‑101 (MQ=255)
                                                                                                         t..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGACGAAAGCCGCAAGGCTGGCGAAtt  >  2:1005930/1‑101 (MQ=255)
                                                                                                          |                                                                                                     
GCCGCCCCGCAACCCGTCGGCGGCCCAGCGCCCTGCTTTTCCTTGGCAAAAT.ACG.CCGGGGGGTGCGG..GGGGCGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAATT  >  NC_014034/2995560‑2995759

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 23 ≤  ATCG  < 32 ≤  ATCG  < 39 ≤  ATCG