Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,590,199 +C 100% intergenic (‑98/‑158) glnE ← / → RCAP_rcc02415 glutamate‑ammonia‑ligase adenylyltransferase/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,590,1991.→C100.0% 42.8 208intergenic (‑98/‑158)glnE/RCAP_rcc02415glutamate‑ammonia‑ligase adenylyltransferase/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (187/5):  ref base (0/16):  total (187/21)
Fisher's exact test strand distribution p-value = 6.27e-20
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.96e-01

Only 100 of 280 total aligned reads displayed.

CATGGCGCCAGACTAGCCGCCGCAAGCCCCGGGTCAAGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CG.CCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTAAAAATCTTTCACATCCCCTCTTGACCCG  >  NC_014034/2590099‑2590294
                                                                                                           |                                                                                                    
cATGGCGCCAGACTAGCCGCCGCAAGCCCCGGGTCAAGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.g.cg                                                                                                       >  2:671400/1‑101 (MQ=255)
cATGGCGCCAGACTAGCCGCCGCAAGCCCCGGGTCAAGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.g.cg                                                                                                       >  2:3347350/1‑101 (MQ=255)
 aTGGCGCCAGACTAGCCGCCGCAAGCCCCGGGTCAAGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.g.cg.c                                                                                                     >  2:1729160/1‑101 (MQ=255)
  gggCCCCAGACTAGCCGCCGCAAGCCCCGGGTCAAGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CG.cc                                                                                                    <  2:5616193/100‑1 (MQ=255)
   ggCGCCAGACTAGCCGCCGCAAGCCCCGGGTCAAGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CG.ccc                                                                                                   <  1:3138742/101‑1 (MQ=255)
    gcgcCAGCCTAGCCGCCGCAAGCCCCGGGTCAAGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CG.ccc.c                                                                                                 <  2:3569481/101‑1 (MQ=255)
     cgcCAGACTAGCCGCCGCAAGCCCCGGGTCAAGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CG.ccc.cc                                                                                                >  1:6266505/1‑101 (MQ=255)
      gcCAGACTAGCCGCCGCAAGCCCCGGGTCAAGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CG.ccc.ccc                                                                                               >  1:3336921/1‑100 (MQ=255)
      ccaagccTAGCCGCCGCAAGCCCCGGGTCAAGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CG.ccc.ccc                                                                                               <  1:1198182/95‑2 (MQ=255)
        cAGACTAGCCGCCGCAAGCCCCGGGTCAAGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CG.ccc.cccgc                                                                                             >  2:6830138/1‑98 (MQ=255)
            cTAGCCGCCGCAAGCCCCGGGTCAAGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCA..GGGGGGTGGGGGGGGC.A.G.CGCCCC.CC..G..CCggg                                                                                      >  2:1122323/1‑101 (MQ=38)
             tAGCCGCCGCAAGCCCCGGGTCAAGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.G.G.GGCCCC.CC..G..CCgggg                                                                                     >  2:3887293/1‑100 (MQ=38)
             tAGCCGCCGCAAGCCCCGGGTCAAGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGt                                                                                     >  2:8486488/1‑101 (MQ=255)
             tAGCCGCCGCAAGCCCCGGGTCAAGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGt                                                                                     >  2:7566196/1‑101 (MQ=255)
              aGCCGCCGCAAGCCCCGGGTCAAGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTc                                                                                    >  2:6426764/1‑101 (MQ=255)
                     gcAAGCCCCGGGTCAAGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCgg                                                                             >  2:4902762/1‑101 (MQ=255)
                      cAAGCCCCGGGTCAAGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGgc                                                                            >  2:8347909/1‑101 (MQ=255)
                      cAAGCCCCGGGTCAAGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGgc                                                                            >  2:5445080/1‑101 (MQ=255)
                      cAAGCCCCGGGTCAAGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CG.CCCGCC..G..CCGGGTCGCCTCGgc                                                                            >  2:5844013/1‑101 (MQ=255)
                       aaGCCCCGGGTCAAGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCGTCGgcg                                                                           >  2:4008841/1‑101 (MQ=255)
                        aGCCCCGGGTCAAGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGgcgc                                                                          >  1:6596420/1‑101 (MQ=255)
                        aGCCCCGGGTCAAGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGgcgc                                                                          >  1:1139485/1‑101 (MQ=255)
                         gCCCCGGGTCAAGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGgcttg                                                                         >  1:2292980/1‑98 (MQ=255)
                             cGGGTCAAGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGcc                                                                     >  2:1407603/1‑101 (MQ=255)
                             cGGGTCAAGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGcc                                                                     >  1:727394/1‑101 (MQ=255)
                              gggTCAAGGCGCACAAATGCCCCGCCCCGGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGAGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCg                                                                    >  1:6520340/1‑101 (MQ=38)
                              gagTCAAGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGGCGCCTCGGCGCAAGCCg                                                                    >  2:5045244/3‑101 (MQ=255)
                               ggTCAAGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGa                                                                   >  1:7494311/1‑101 (MQ=255)
                                gTCAAGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAg                                                                  >  1:5033124/1‑101 (MQ=255)
                                gTCAAGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAg                                                                  >  1:1676723/1‑101 (MQ=255)
                                 tCAAGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAgg                                                                 >  2:4554022/1‑101 (MQ=255)
                                 tCAAGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAgg                                                                 >  1:5978192/1‑101 (MQ=255)
                                 tCAAGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAgg                                                                 >  2:4337441/1‑101 (MQ=255)
                                   aaGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAgcga                                                               >  2:8045881/1‑98 (MQ=255)
                                   aaGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCg                                                               >  2:4612483/1‑101 (MQ=255)
                                   aaGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATACC..GGGGGGTGCGGGGGGA.A.G.CGCCCC.CC..G..CCGGGGAGCCTCGGCGCAAACGGAggga                                                               >  2:4240888/1‑99 (MQ=255)
                                     ggCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGaa                                                             >  2:1941703/1‑101 (MQ=255)
                                        gccgaaacccccccccccTGTCCGTCTGGGCGCGAATCCCGG..GGGGGGGGGGGGGGGC.G.GCCC.CCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAAt                                                           <  2:2632656/92‑1 (MQ=12)
                                          gcAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAAtctc                                                        >  1:445744/1‑101 (MQ=255)
                                           caagccccccgccctTGT.CCTCTGGGCGCAAATCCCCG..GGGGGGGGGGGGGGGC.G.GCCC.CCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAAtctct                                                       <  2:2160030/95‑1 (MQ=16)
                                           cAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAAtctct                                                       >  2:6320970/1‑101 (MQ=255)
                                            aaaTGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAAtctctc                                                      >  2:296735/1‑101 (MQ=255)
                                            aaaTGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAAtctctc                                                      >  2:1614950/1‑101 (MQ=255)
                                             aaTGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTcc                                                     >  1:6844545/1‑101 (MQ=255)
                                              ccccccccccccTGT.CCCTTGGGGGCGAA.ATCCCGGGGGGGGGGGGGGGGGC.C.GCCC.CCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTcc                                                     <  2:4665892/98‑1 (MQ=9)
                                              atccccccccccTGT.CTTTGTGGCGCGAAGCCCGG..GGGGGGTGCGGGGGGC.G.GCCC.CCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTccc                                                    <  1:5682817/98‑1 (MQ=16)
                                              aTGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTccc                                                    >  1:8042635/1‑101 (MQ=255)
                                              aTGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTccc                                                    >  1:6703641/1‑101 (MQ=255)
                                                gccccgccccTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCg                                                  >  1:9163846/1‑101 (MQ=255)
                                                   ccgccccTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGAcc                                               >  1:1329273/1‑101 (MQ=255)
                                                     gccccTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGAcccc                                             >  1:5336377/1‑101 (MQ=255)
                                                       cccTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCtt                                           >  1:4781477/1‑101 (MQ=255)
                                                       cccTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCtt                                           >  2:5877105/1‑101 (MQ=255)
                                                         cgtttcctCTGGGCGCAAATCCCCG..GGGGGGGGGGGGGGGC.A.GCCC.CCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTg                                          <  2:7191260/94‑1 (MQ=21)
                                                          tGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTgcg                                        >  2:5155344/1‑101 (MQ=255)
                                                           tgcctgctcggCGCGAATCCCGG..GGGGGGGGGGGGGGGCGC.C.CC.CCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTgcg                                        <  2:568904/92‑1 (MQ=16)
                                                           gT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTgcgg                                       >  2:2717364/1‑101 (MQ=255)
                                                           gT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTgcgg                                       >  1:2113359/1‑101 (MQ=255)
                                                            t.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGCAATCTCTCCCCCAACCCTTgcgtc                                      >  2:7339581/1‑99 (MQ=25)
                                                              cttttgggCGCAAATCCCCG..GGGGGGGGGGGGGGGC.A.GCCC.CCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTgcggcg                                     <  2:6629659/95‑1 (MQ=21)
                                                               aTCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGgcgc                                    >  2:7126104/1‑101 (MQ=255)
                                                                 cTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCaa                                  >  2:8722580/1‑101 (MQ=255)
                                                                  ttGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAAt                                 >  2:9104203/1‑101 (MQ=255)
                                                                  ttGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAAt                                 >  1:404151/1‑101 (MQ=255)
                                                                   tGGCGCGAATCTCCG..GGGGGGGGGGGGGGCC.G.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATg                                <  2:8310110/101‑1 (MQ=25)
                                                                     gggggatggtCCG..GGGGGGGGGCGGGGGC.A.C.CC.CCC.CC..GGCCCGGGTCGCCACGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGt                               <  1:5793140/92‑1 (MQ=19)
                                                                       gcgAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTaaa                            >  1:3370824/1‑101 (MQ=255)
                                                                        gcaatgcccc..GGGGGGGGGGGGGGGC.C.C.CC.CCC.CCCGG..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTaaa                            <  1:8474097/94‑1 (MQ=25)
                                                                         gAATATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTaaaaa                          >  1:1874012/1‑101 (MQ=255)
                                                                          gaaccccc..GGGGGGGGGGGGGGGCGAGG.CC.CCC.CG..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTaaaaa                          <  1:7566520/96‑1 (MQ=25)
                                                                          aaTATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTAAAAAt                         >  2:7797166/1‑101 (MQ=255)
                                                                          aaTATCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTAAAAAt                         >  1:9187493/1‑101 (MQ=255)
                                                                           atatCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTAAAAATc                        >  1:3664511/1‑101 (MQ=255)
                                                                            tatCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTAAAAATCt                       >  1:696435/1‑101 (MQ=255)
                                                                             atCCC..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTAAAAATCtt                      >  1:4501996/1‑101 (MQ=255)
                                                                               ccc..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCGCGACCCCTTGCGGCGCAATGTAAAAATCTTTc                    >  1:4669667/1‑101 (MQ=255)
                                                                               ccc..GGGGGGTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTAAAAATCTTTc                    >  2:8123320/1‑101 (MQ=255)
                                                                               acccgggggggGGGGGGGGGC.G.C.CCCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTAAAAATCtt                      <  2:3696979/96‑1 (MQ=25)
                                                                                cgggggggggGGGGGGGCGC.C.C.CC.CCC.CC..C..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTAAAAATCTTTc                    <  2:5992951/97‑1 (MQ=37)
                                                                                cgcgggggggGGCGGGGGGC.C.GCCC.CCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTAAAAATCttt                     <  2:2863103/97‑1 (MQ=37)
                                                                                  ggggggggGGGGGGGGCC.C.C.CC.CCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTAAAAATCTTTcac                  <  1:3855317/99‑1 (MQ=37)
                                                                                    ggggggTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTAAAAATCTTTcaca                 >  2:448082/1‑101 (MQ=255)
                                                                                    ggggggGGGGGGGGGC.G.GCCC.CCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTAAAAATCTTTcaca                 <  2:7397945/101‑1 (MQ=37)
                                                                                    ggggggGGGGGGGGCC.C.G.CCCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTAAAAATCTTTcaca                 <  2:4072036/101‑1 (MQ=37)
                                                                                    ggggggGGGGGGGGCC.C.G.CCCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTAAAAATCTTTcaca                 <  1:4821421/101‑1 (MQ=37)
                                                                                     gggggTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTAAAAATCTTTCACAt                >  1:8362537/1‑101 (MQ=255)
                                                                                      ggggTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTAAAAATCTTTCACATc               >  2:7005898/1‑101 (MQ=255)
                                                                                      ggggTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTAAAAATCTTTCACATc               >  2:1563711/1‑101 (MQ=255)
                                                                                        ggTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTAAAAATCTTTCACATccc             >  2:329892/1‑101 (MQ=255)
                                                                                         gTGCGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCATTGTAAAAATCTTTCACATcccc            >  2:3813538/1‑101 (MQ=255)
                                                                                            gggggggC.G.G.CCCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTAAAAATCTTTCACATCCCctct         <  1:4378896/100‑1 (MQ=38)
                                                                                            cGGGGGGC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAACGTAAAAATCTTTCACATCCCctct         >  2:9000831/1‑101 (MQ=255)
                                                                                             ggggggC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTAAAAATCTTTCACATCCCCTCtt        >  2:7716798/1‑101 (MQ=255)
                                                                                             ggggggC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTAAAAATCTTTCACATCCCCTCtt        >  1:2585211/1‑101 (MQ=255)
                                                                                              gggggC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTAAAAATCTTTCACATCCCCTCTTg       >  2:6523669/1‑101 (MQ=255)
                                                                                               ggggC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTAAAAATCTTTCACATCCCCTCTTga      >  1:9068031/1‑101 (MQ=255)
                                                                                                gggC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTAAAAATCTTTCACATCCCCTCTTgac     >  1:5486638/1‑101 (MQ=255)
                                                                                                 ggC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTAAAAATCTTTCACATCCCCTCTTgacc    >  2:2457419/1‑101 (MQ=255)
                                                                                                 ggC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTAAAAATCTTTCACATCCCCTCTTgacc    >  1:1984571/1‑101 (MQ=255)
                                                                                                  gC.A.G.CGCCCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTAAAAATCTTTCACATCCCCTCTTgaccc   >  1:3530437/1‑101 (MQ=255)
                                                                                                           |                                                                                                    
CATGGCGCCAGACTAGCCGCCGCAAGCCCCGGGTCAAGGCGCGCAAATGCCCCGCCCCTGT.CATCTTGGCGCGAATATCCC..GGGGGGTGCGGGGGGC.A.G.CG.CCC.CC..G..CCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTAAAAATCTTTCACATCCCCTCTTGACCCG  >  NC_014034/2590099‑2590294

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