Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,995,643 +C 7.7% intergenic (+183/+37) RCAP_rcc02802 → / ← RCAP_rcc02803 membrane dipeptidase/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,995,6431.→C 7.7% 3.7 ‑0.6 ‑0.039intergenic (+183/+37)RCAP_rcc02802/RCAP_rcc02803membrane dipeptidase/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (3/0):  ref base (24/12):  total (27/12)
Fisher's exact test strand distribution p-value = 5.39e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.34e-01

GGAAAACCCCTTCGGTTGCCGCCCCGCAACCCGTCGGCGGCCCAGCGCCCTGCTTTTCCTTGGCAAAAT.ACG.CCGGGGGGTGCGG..GGGGCGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCC  >  NC_014034/2995543‑2995743
                                                                                                          |                                                                                                      
ggAACACCCCTTCGGTTCCCGCCCCGCACCCCGTCGGCGGCCCAGCGCCCTGCTTTTCCTTGGCAAAAT.ACG.CCGGGGGGTGCGG..GGGGCGG.AGCCCCCcg                                                                                                         <  1:4384569/101‑1 (MQ=255)
  gaaaCCCCTTCGGTTGCCGCCCCGCAACCCGTCGGCGGCCCAGCGCCCTGCTTTTCCTTGGCACAAT.ACG.CCGGGGGGTGCGG..GGGGCGG.AGCCCCCcc.cg                                                                                                      >  1:2061467/2‑101 (MQ=255)
      ccccTTCGTTTCCCGCCCCGCACCCCGTGGGCGGCCCGGCGCCCTGTTTTTCTTGGGAAAATT.ACG.CGGGGGGGTGCGG..GGGGGGG.GGCCCCCCG.CGGGTc                                                                                                  <  2:1544220/101‑1 (MQ=16)
       cccTGCGGGTGCGGCCCCGAAGCCCGTGGGGGGCCCCGCCCCTTGCTTTGCCTGGCCAAAAG.GGG.GGGGGGGGGGGGG..GGGGCGG.AGCCCCCCG.CGGGTCg                                                                                                 <  2:1929713/101‑1 (MQ=9)
          ttCGGTTGCCCCCGCGCACCCCGTCGGCGGCCCAGCGCCCTGCTTTCCCTTGGCAAAAT.AGG.CCGGGGGGTCCGG..GGGGCGG.AGCCCCCCG.CGGGTCGccc                                                                                              <  2:3893328/101‑1 (MQ=18)
                 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:956431/1‑101 (MQ=255)
                                                                                             cGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCg             >  2:311191/1‑101 (MQ=255)
                                                                                             cGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCg             >  1:709107/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)
                                                                                                          |                                                                                                      
GGAAAACCCCTTCGGTTGCCGCCCCGCAACCCGTCGGCGGCCCAGCGCCCTGCTTTTCCTTGGCAAAAT.ACG.CCGGGGGGTGCGG..GGGGCGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCC  >  NC_014034/2995543‑2995743

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