Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,995,678 Δ1 bp 3.5% intergenic (+218/+2) RCAP_rcc02802 → / ← RCAP_rcc02803 membrane dipeptidase/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,995,6780G→. 3.5% 2.3 0.0 ‑0.186intergenic (+218/+2)RCAP_rcc02802/RCAP_rcc02803membrane dipeptidase/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (3/0):  ref base (65/18):  total (68/18)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.62e-01

GGCCCAGCGCCCTGCTTTTCCTTGGCAAAATACG.CCGGGGGGTGCGG..GGGGCGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAATTGCGGGAATTGCGCCTCGAA  >  NC_014034/2995581‑2995778
                                                                                                        |                                                                                                    
ggcCCAGCGCCCTGCTTTTCCTTGGCAAAATACG.CCGGGGGGTGCGG..GGGGCGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAAGCCCCCGCTc                                                                                                   >  1:686689/1‑101 (MQ=255)
       cccccTGCTTTTCCTTGGCAAAATCCG.CGGGGGGGTGGGG..GGGGCGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCa                                                                                            <  2:2924081/99‑1 (MQ=255)
           cTGCTTTTCCTTGGCAAAATACG.CCGGGGGGTGCGG..GGGGCGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAgaga                                                                                        >  1:632652/1‑101 (MQ=255)
            tGCTTTTCCTTGGCAAAATACG.CCGGGGGGTGCGG..GGGGCGG.AGCCCCCCG.CGGGTCTCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCagaggc                                                                                       >  2:3927275/1‑99 (MQ=255)
              tttttCCTTGGAAAAATACGCCCGGGGGGTGCGG..GGGGCGG.AGCCCCCCGCCGGGTCGCCTCGGCGTCAGCCGAGGCGAAATACCCC.CTCAGCccttcatc                                                                                        >  2:2611488/2‑95 (MQ=0)
                 ttCCTTGGAAAAATACGCCCGGGGGGTGCGG..GGGGCGG.AGCCCCCCGCCGGGTCGCCTCGGCGTCCGCCGAGGCGAAATACCCC.CTCAGCccttcatcttg                                                                                     >  1:789505/1‑92 (MQ=0)
                 ttCCTTGGAAAAATACGCCCGGGGGGTGCGG..GGGGCGG.AGCCCCCCGCCGGGTCGCCTCGGCGTCAGCCGAGGCGAAATACCCC.CTCAGCccttcatcttg                                                                                     >  2:4078231/1‑92 (MQ=0)
                       ggCAAAATACG.CCGGGGGGTGCGG..GGGGCGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGAccccc                                                                            >  1:1263423/1‑98 (MQ=255)
                           aaaTACG.CCGGGGGGTGCGG..GGGGCGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAAcc                                                                        >  1:1094053/1‑101 (MQ=255)
                           aaaTACG.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)
                                              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:311191/1‑101 (MQ=255)
                                                      cGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCg                                                >  2:4360965/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:3109657/1‑101 (MQ=255)
                                                                    cgGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCg                                    >  2:2708310/1‑101 (MQ=255)
                                                                    cgGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCg                                    >  1:3876056/1‑101 (MQ=255)
                                                                    cgGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCg                                    >  2:1835933/1‑101 (MQ=255)
                                                                    cgGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCg                                    >  1:1359693/1‑101 (MQ=255)
                                                                    cgGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCg                                    >  2:2508254/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                                >  1:2713617/1‑101 (MQ=255)
                                                                        tCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAg                                >  2:2216160/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:4689760/1‑101 (MQ=255)
                                                                                gCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGa                        <  2:1637647/101‑1 (MQ=255)
                                                                                gCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGa                        >  1:1780596/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)
                                                                                      gCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGGGAATTg                    >  2:1651572/1‑101 (MQ=255)
                                                                                      gCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAATTg                    <  1:2133291/101‑1 (MQ=255)
                                                                                       ccGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAATTGc                   >  1:146082/1‑101 (MQ=255)
                                                                                        cGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTTGCGAATTGCg                  >  2:513277/1‑101 (MQ=255)
                                                                                        cGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAATTGCg                  >  2:3603475/1‑101 (MQ=255)
                                                                                        cGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAATTGCg                  >  2:2800601/1‑101 (MQ=255)
                                                                                        cGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAATTGCg                  >  1:4702771/1‑101 (MQ=255)
                                                                                         ggggCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAATTGCgg                 >  2:4218708/1‑101 (MQ=255)
                                                                                         ggggCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAATTGCgg                 >  2:2539397/1‑101 (MQ=255)
                                                                                          gggCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAATTGCggg                >  2:769220/1‑101 (MQ=255)
                                                                                          gggCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAATTGCggg                >  1:1051651/1‑101 (MQ=255)
                                                                                            gCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAATTGCGGGaa              >  2:499780/1‑101 (MQ=255)
                                                                                              gAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAATTGCGGGAAtt            >  1:2618747/1‑101 (MQ=255)
                                                                                               aaaTCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAATTGCGGGAATTg           >  2:787752/1‑101 (MQ=255)
                                                                                               aaaTCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAATTGCGGGAATTg           >  1:3503856/1‑101 (MQ=255)
                                                                                                aaTCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGGGAATTGCGGGAATTgc          >  1:3215848/1‑101 (MQ=255)
                                                                                                aaTCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAATTGCGGGAATTgc          >  2:2325877/1‑101 (MQ=255)
                                                                                                 aTCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAATTGCGGGAATTgcg         <  1:2531338/101‑1 (MQ=255)
                                                                                                  tCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAATTGCGGGAATTgcgc        >  1:4403637/1‑101 (MQ=255)
                                                                                                   cccccGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAATTGCGGGAATTGCGcc       >  1:2859520/1‑101 (MQ=255)
                                                                                                   cccccGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAATTGCGGGAATTGCGcc       >  1:4740875/1‑101 (MQ=255)
                                                                                                    ccccGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAATTGCGGGAATTGCGCCt      >  2:3187393/1‑101 (MQ=255)
                                                                                                     cccGCTCAGTCCCAGAGTCTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAATTGCGGGAATTGCGCCTc     <  1:1635826/101‑1 (MQ=255)
                                                                                                     cccGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAATTGCGGGAATTGCGCCTc     >  1:2790299/1‑101 (MQ=255)
                                                                                                     cccGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAATTGCGGGAATTGCGCCTc     <  2:4403387/101‑1 (MQ=255)
                                                                                                     cccGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAATTGCGGGAATTGCGCCTc     >  1:82463/1‑101 (MQ=255)
                                                                                                       cGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAATTGCGGGAATTGCGCCTCGa   >  1:586749/1‑101 (MQ=255)
                                                                                                        gCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAATTGCGGGAATTGCGCCTCGaa  <  2:632588/101‑1 (MQ=255)
                                                                                                        gCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAATTGCGGGAATTGCGCCTCGaa  <  1:1651672/101‑1 (MQ=255)
                                                                                                        gCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAATTGCGGGAATTGCGCCTCGaa  <  2:686625/101‑1 (MQ=255)
                                                                                                        gCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAATTGCGGGAATTGCGCCTCGaa  >  1:762362/1‑101 (MQ=255)
                                                                                                        gCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAAGTGCGGGAATTGCGCCTCGaa  >  2:4633218/1‑101 (MQ=255)
                                                                                                        |                                                                                                    
GGCCCAGCGCCCTGCTTTTCCTTGGCAAAATACG.CCGGGGGGTGCGG..GGGGCGG.AGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCCGCAAGGCTGGCGAATTGCGGGAATTGCGCCTCGAA  >  NC_014034/2995581‑2995778

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 24 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG