Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,995,625 C→G 14.3% intergenic (+165/+55) RCAP_rcc02802 → / ← RCAP_rcc02803 membrane dipeptidase/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,995,6250C→G14.3% 2.9 ‑3.8 0.039intergenic (+165/+55)RCAP_rcc02802/RCAP_rcc02803membrane dipeptidase/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (0/6):  ref base (18/15):  total (18/21)
Fisher's exact test strand distribution p-value = 2.23e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

AGAGGCGTATTTTGCCAAGGAAAACCCCTTCGGTTG.CCGCCCCGCAACCCGTCGGCGGCCCAGCGCCCTGCTTTTCCTTGG.CAAAAT.ACG.CCGGGGGGTGCGG..GGGGCGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGC  >  NC_014034/2995525‑2995725
                                                                                                        |                                                                                                          
aggggCTTATTTGGCCAAGGAAACCCCCTTCGGTTG.CCGCCCCGCACCCCGTCGGCGGCCCAGCGCCCTGCTTTTCCTTGG.CAAAAT.ACG.CCGGGGGGTGc                                                                                                            <  2:1711956/98‑1 (MQ=37)
agagGCGTTTTTTGCCAAGAAAACCCCCTTGGGTTC.CCGCCCCGCACCCCGTCGGCGCCCCAGCGCCCTGCTTTTCCTTGG.AAAAAT.ACG.CCGGGGGGTGc                                                                                                            <  2:3587280/101‑1 (MQ=25)
agagGCGTATTTTGCCAAGGAAAACCCCTTCGGTTG.CCGCCCCGCAACCCGTCGGCGGCCCAGCGCCCTGCTTTTCCTTGG.CAAAAT.ACG.CCGGGGGGTGc                                                                                                            <  2:4729626/101‑1 (MQ=255)
 gagGCGTATTTTGCCAAGGAAAACCCCTTCGGTTG.CCGCCCCGCAACCCGTCGGCGGCCCAGCGCCCTGCTTTTCCTTGG.CAAAAT.ACG.CCGGGGGGTGCg                                                                                                           <  2:2875893/101‑1 (MQ=255)
  agGCGTATTTTGCCAAGGAAACCCCCTTCGGCTG.CCGCCCCGCAGCCCGTCGGCGGCCCGGCGCCCTGCTTTTCCTTCG.CAAACT.ACG.CCGGGGGGTGCgg                                                                                                          <  2:4735559/101‑1 (MQ=37)
   ggCGTATTTTGCCAAGGAAAACCCCTTCGGTTG.CCGCCCCGCAACCCGTCGGCGGCCCAGCGCCCTGCTTTTCCTTGG.CAAAAT.ACG.CCGGGGGGTGCgg..g                                                                                                       <  2:702095/101‑1 (MQ=255)
       ttgtttGCCAAGGAAAACCCCCTCGGGTGCCCCCCCCGCACCCGGCGGGCGGCCCAGCCCCTTGTTTTTCCTTGG.AAAAAT.ACG.CCGGGGGGTGCgg..gggg                                                                                                    <  2:3437670/98‑1 (MQ=16)
        aTTTTGCCAAGGAAACCCCCTTCGGTTG.CCGCCCCGCAACCCGTCGGCGGCCCAGCGCCCTGCTTTTCCTTGG.CAAAAT.ACG.CCGGGGGGTGCGG..GGGGCg                                                                                                  <  2:749327/101‑1 (MQ=255)
        aTTTTGCCAAGGAAAACCCCTTCGGTTG.CCGCCCCGCAACCCGTCGGCGGCCCAGCGCCCTGCTTTTCCTTGG.CAAAAT.ACG.CCGGGGGGTGCGG..GGGGCg                                                                                                  >  2:19752/1‑101 (MQ=255)
          tttGCCAAGGAAAACCCCTTCGGTTG.CCGCCCCGCAACCCGTCGGCGGCCCAGCGCCCTGCTTTTCCTTGG.CAAAAT.ACG.CCGGGGGGTGCGG..GGGGCGG.a                                                                                               <  1:3844525/101‑1 (MQ=255)
          tttGCCAAGAAACCCCCCTTCGTTTC.CCCCCCCGCACCCCGTCGGCGGCCCAGCCCCTTGTTTTTCCTTGG.CAAAAT.ACG.CCGGGGGGTGCGG..GGGGCGG.a                                                                                               <  2:3212058/101‑1 (MQ=18)
           ttGCCAAGGAAACCCCCTTCGGTTG.CCGCCCCGCAACCCGTCGGCGGCCCAGCGCCCTGCTTTTCCTTGG.CAAAAT.ACG.CCGGGGGGTGCGG..GGGGCGG.Ag                                                                                              <  2:4203650/101‑1 (MQ=255)
            tGCCAAGGAAAACCCCTTCGGTTG.CCGCCCCGCAACCCGTCGGCGGCCCAGCGCCCTGCTTTTCCTTGG.CAAAAT.ACG.CCGGGGGGTGCGG..GGGGCGG.AGc                                                                                             >  2:362088/1‑101 (MQ=255)
             gCCAAGGAAACCCCCTTCGGTTG.CCGCCCCGCAGCCCGTCGGCGGCCCAGCGCCCTGCTTTTCCTTGG.CAAAAT.ACG.CCGGGGGGTGCGG..GGGGCGG.AGcc                                                                                            <  1:362122/101‑1 (MQ=255)
             gCCAAGGAAAACCCCTTCGGTTG.CCGCCCCGCAACCCGTCGGCGGCCCAGCGCCCTGCTTTTCCTTGG.CAAAAT.ACG.CCGGGGGGTGCGG..GGGGCGG.AGcc                                                                                            <  2:3173967/101‑1 (MQ=255)
              caaagaaaccccccTTCGGTCC.CCCCCCCGCCCCCCGGGGGC.CCCCCCCCCCCTTTTTTTTCTTGGAAAAAAG.CCC.CCGGGGGGTGCGG..GGGGCGG.AGccc                                                                                           <  1:1623871/91‑1 (MQ=2)
                aaGGAAACCCCCTTCGGTTG.CCGCCCCGCAACCCGTCGGCGGCCCAGCGCCCTGCTTTTCCTTGG.CAAAAT.ACG.CCGGGGGGTGCGG..GGGGCGG.AGccccc                                                                                         <  2:1944041/101‑1 (MQ=255)
                  ggAACACCCCTTCGGTTC.CCGCCCCGCACCCCGTCGGCGGCCCAGCGCCCTGCTTTTCCTTGG.CAAAAT.ACG.CCGGGGGGTGCGG..GGGGCGG.AGCCCCCcg                                                                                       <  1:4384569/101‑1 (MQ=255)
                    gaaaCCCCTTCGGTTG.CCGCCCCGCAACCCGTCGGCGGCCCAGCGCCCTGCTTTTCCTTGG.CACAAT.ACG.CCGGGGGGTGCGG..GGGGCGG.AGCCCCCcc.cg                                                                                    >  1:2061467/2‑101 (MQ=255)
                        ccccTTCGTTTC.CCGCCCCGCACCCCGTGGGCGGCCCGGCGCCCTGTTTTTCTTGGG.AAAATT.ACG.CGGGGGGGTGCGG..GGGGGGG.GGCCCCCCG.CGGGTc                                                                                <  2:1544220/101‑1 (MQ=16)
                         cccTGCGGGTG.CGGCCCCGAAGCCCGTGGGGGGCCCCGCCCCTTGCTTTGCCTGGC.CAAAAG.GGG.GGGGGGGGGGGGG..GGGGCGG.AGCCCCCCG.CGGGTCg                                                                               <  2:1929713/101‑1 (MQ=9)
                            ttCGGTTG.CCCCCGCGCACCCCGTCGGCGGCCCAGCGCCCTGCTTTCCCTTGG.CAAAAT.AGG.CCGGGGGGTCCGG..GGGGCGG.AGCCCCCCG.CGGGTCGccc                                                                            <  2:3893328/101‑1 (MQ=18)
                                    ccccccccCCACCCCGTGGGCGCCCCAGCCCCCTTTTTTCCCGGGG.CAAAAT.CCG.CGGGGGGGTGGGG..GGGGGGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..g                                                                   <  1:1419898/97‑1 (MQ=12)
                                                 ccGTCGGCGGCCCAGCGCCCTGCTTTTCCTTGG.CAAAAT.ACG.CCGGGGGGTGCGG..GGGGCGG.AGCCCCCCG.CGGGGCGCCCCGGCTT..GCCGGGGCGGAATc                                                      >  2:4057193/1‑101 (MQ=255)
                                                 ccGCCGGCGCCCCAGCCCCCTGCTTTTCCTGGG.AAAAAT.GGG.CGGGGGGGTGGGG..GGGGGGG.ACCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATc                                                      <  1:4601521/101‑1 (MQ=18)
                                                    gccgcggcCCCGGCCCCTTTTTTTCCTTGG.GAAAATCGCG.CGGGGGGGGGGGG..GGGGGGG.CCCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATccc                                                    <  2:2570271/98‑1 (MQ=16)
                                                         ggcCCAGCGCCCTGCTTTTCCTTGG.CAAAAT.ACG.CCGGGGGGTGCGG..GGGGCGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAAGCCCCCGCTc                                              >  1:686689/1‑101 (MQ=255)
                                                                cccccTGCTTTTCCTTGG.CAAAAT.CCG.CGGGGGGGTGGGG..GGGGCGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCa                                       <  2:2924081/99‑1 (MQ=255)
                                                                    cTGCTTTTCCTTGG.CAAAAT.ACG.CCGGGGGGTGCGG..GGGGCGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAgaga                                   >  1:632652/1‑101 (MQ=255)
                                                                     tGCTTTTCCTTGG.CAAAAT.ACG.CCGGGGGGTGCGG..GGGGCGG.AGCCCCCCG.CGGGTCTCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCagaggc                                  >  2:3927275/1‑99 (MQ=255)
                                                                       tttttCCTTGG.AAAAAT.ACGCCCGGGGGGTGCGG..GGGGCGG.AGCCCCCCGCCGGGTCGCCTCGGCGTCAGCCGAGGCGAAATACCCC.CTCAGCccttcatc                                   >  2:2611488/2‑95 (MQ=0)
                                                                          ttCCTTGG.AAAAAT.ACGCCCGGGGGGTGCGG..GGGGCGG.AGCCCCCCGCCGGGTCGCCTCGGCGTCCGCCGAGGCGAAATACCCC.CTCAGCccttcatcttg                                >  1:789505/1‑92 (MQ=0)
                                                                          ttCCTTGG.AAAAAT.ACGCCCGGGGGGTGCGG..GGGGCGG.AGCCCCCCGCCGGGTCGCCTCGGCGTCAGCCGAGGCGAAATACCCC.CTCAGCccttcatcttg                                >  2:4078231/1‑92 (MQ=0)
                                                                                gg.CAAAAT.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:3848689/1‑101 (MQ=255)
                                                                                                        cGG..GGGGCGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATgc  >  1:2073503/1‑101 (MQ=255)
                                                                                                        |                                                                                                          
AGAGGCGTATTTTGCCAAGGAAAACCCCTTCGGTTG.CCGCCCCGCAACCCGTCGGCGGCCCAGCGCCCTGCTTTTCCTTGG.CAAAAT.ACG.CCGGGGGGTGCGG..GGGGCGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGC  >  NC_014034/2995525‑2995725

Base quality scores:  ATCG  < 3 ≤  ATCG  < 12 ≤  ATCG  < 20 ≤  ATCG  < 31 ≤  ATCG  < 39 ≤  ATCG