Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,731 +C 15.7% intergenic (+543/‑388) RCAP_rcc01306 → / → RCAP_rcc01307 Cas2 family CRISPR‑associated protein/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,395,7311.→C15.7% 150.7 ‑0.3 ‑0.0293intergenic (+543/‑388)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (20/26):  ref base (116/131):  total (136/157)
Fisher's exact test strand distribution p-value = 7.48e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.99e-01

Only 100 of 308 total aligned reads displayed.

GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGAC  >  NC_014034/1395631‑1395831
                                                                                                     |                                                                                                            
gCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                                                                                               <  1:1695066/101‑1 (MQ=255)
  gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.gg                                                                                                            >  2:2690386/1‑101 (MQ=255)
   gcGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGt                                                                                                           <  2:2724383/100‑1 (MQ=255)
      cTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATg                                                                                                        <  2:3657155/101‑1 (MQ=255)
      cTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATg                                                                                                        >  1:4599649/1‑101 (MQ=255)
       cccccccGGGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGGGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATgg                                                                                                       <  1:1760071/100‑1 (MQ=255)
       cccccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATgg                                                                                                       <  1:2108345/100‑1 (MQ=255)
         cccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGAt                                                                                                     >  2:1041889/1‑101 (MQ=255)
         cccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGAt                                                                                                     <  1:2087672/101‑1 (MQ=255)
         cccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGAt                                                                                                     <  1:3621762/101‑1 (MQ=255)
          ccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATc                                                                                                    <  2:1935384/101‑1 (MQ=255)
          ccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATc                                                                                                    <  1:340998/101‑1 (MQ=255)
           cccGGGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCt                                                                                                   <  1:2611943/101‑1 (MQ=255)
              gcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..gc                                                                                             >  1:1235156/1‑101 (MQ=255)
              gcaggggCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..gc                                                                                             >  1:2892064/4‑101 (MQ=255)
                gggggGCTGGACCGAAATATCTGCGTGTAGTCGGGGGGCTGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..gcgc                                                                                           <  2:3835088/101‑1 (MQ=25)
                 ggggCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCa                                                                                          >  2:3915629/1‑101 (MQ=255)
                  gggCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCaa                                                                                         <  1:3254446/101‑1 (MQ=255)
                   ggggTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCaaa                                                                                        <  1:3068081/98‑1 (MQ=255)
                   ggCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCaaa                                                                                        >  2:521171/1‑101 (MQ=255)
                   ggCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCaaa                                                                                        >  1:728828/1‑101 (MQ=255)
                      gTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGAGCGAAACT.GGTATGGATCT.G..GCGCAAACgg                                                                                     >  2:479927/1‑101 (MQ=255)
                       tGGATCGAAACATCTGCATGTAGTGGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..a                                                                                  <  2:1860631/101‑1 (MQ=255)
                        ggATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..aa                                                                                 >  2:2285565/1‑101 (MQ=255)
                         gATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..aaa                                                                                >  1:138248/1‑101 (MQ=255)
                            cGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTc                                                                             >  1:516782/1‑101 (MQ=255)
                            cGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTc                                                                             >  2:2767760/1‑101 (MQ=255)
                             gAAACATCTGCATGTAGTGGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGCCt                                                                            >  2:2138198/1‑101 (MQ=255)
                                  ctCTGCATGTGGCGGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.g.tcgt                                                                     <  2:1548907/100‑1 (MQ=255)
                                  aTCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.g.tcgt                                                                     <  2:1616498/101‑1 (MQ=255)
                                    cTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.tcgtcg                                                                   >  1:2398007/1‑101 (MQ=255)
                                        aTGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.cc                                                              >  1:3816673/1‑101 (MQ=255)
                                          gTAGTCGGGGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.cccc                                                            <  1:1161287/101‑1 (MQ=255)
                                          gTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.cccc                                                            <  2:4599327/101‑1 (MQ=255)
                                              tgggggggCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCcgc                                                        <  2:4502216/96‑1 (MQ=255)
                                                   gggCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCggggg                                                   >  1:4037085/1‑101 (MQ=255)
                                                     gCAGTGCGAGCCGCAGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCCGCGGGGGc                                                  <  2:2617497/101‑1 (MQ=38)
                                                      cAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGt                                                <  1:2690571/101‑1 (MQ=255)
                                                       aGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTg                                               >  2:3345035/1‑101 (MQ=255)
                                                       aGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCN.CCCCCCGCGGGGGCGTg                                               >  2:1196052/1‑101 (MQ=255)
                                                        gTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTgg                                              >  2:3049997/1‑101 (MQ=255)
                                                         ttCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGa                                             <  1:4930463/101‑1 (MQ=255)
                                                          tcaagattcAGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGTg                                               <  2:3251915/93‑1 (MQ=2)
                                                          tCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGAt                                            <  1:3604172/101‑1 (MQ=255)
                                                           caagattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGTgg                                              <  1:4882990/94‑1 (MQ=2)
                                                           cGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATc                                           <  2:721739/101‑1 (MQ=255)
                                                            gAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCg                                          <  1:2328698/101‑1 (MQ=255)
                                                            aagattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGTGGa                                             <  1:2524565/95‑1 (MQ=2)
                                                            aagattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGTGGa                                             <  1:3662728/95‑1 (MQ=2)
                                                              gcgctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCT.CCCCCCGCGGGGGCGTGGATCg                                          <  1:851139/95‑1 (MQ=0)
                                                              gcgctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCT.CCCCCCGCGGGGGCGTGGATCg                                          <  1:3977772/95‑1 (MQ=0)
                                                              gattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.GGGTCGCT.CCCCCCGCGGGGGCGTGGATc                                           >  2:506593/5‑101 (MQ=2)
                                                               ccGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TGGTCGCT.CCCCCCGCGGGGGCGTGGATCGaaa                                       <  1:4527538/101‑1 (MQ=255)
                                                               attcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCG..GCGATCTGGGT.AAGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGTGGATCg                                          >  1:2002408/4‑101 (MQ=2)
                                                               attcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGTGGATCg                                          <  1:4809219/98‑1 (MQ=2)
                                                                ttcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTACGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGCGGATCGa                                         >  1:4263048/3‑101 (MQ=0)
                                                                ttcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGTGGATCGa                                         <  2:4752545/99‑1 (MQ=2)
                                                                 tcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGTGGATCGaa                                        >  2:3575117/2‑101 (MQ=2)
                                                                 gCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAAcc                                     >  2:4111666/1‑101 (MQ=255)
                                                                 ctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCT.CCCCCCGCGGGGGCGTGGATCGaaa                                       <  2:4923984/98‑1 (MQ=2)
                                                                  tggTCGCCCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAAc                                      <  2:1387470/99‑1 (MQ=2)
                                                                  cAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCt                                    >  1:3220649/1‑101 (MQ=255)
                                                                   ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAAcc                                     >  2:3711699/2‑101 (MQ=2)
                                                                   aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAAc                                      <  1:955020/101‑1 (MQ=2)
                                                                   aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTg                                   <  2:1488763/101‑1 (MQ=255)
                                                                    gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTgc                                  >  1:1314846/1‑101 (MQ=255)
                                                                    gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTgc                                  >  2:4007226/1‑101 (MQ=255)
                                                                    gTCGCTCCCCCCGCGGGGGAGTGGATCGAAACT.GGTATGGATCT.G..GCGTAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTgc                                  <  1:2830236/101‑1 (MQ=255)
                                                                    gTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGTGGATCGAGAca                                     <  1:1896147/101‑2 (MQ=2)
                                                                      cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCtat                                  <  2:4289950/101‑3 (MQ=2)
                                                                      cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTgcgg                                >  2:4903140/1‑101 (MQ=255)
                                                                      cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCCgcgg                                >  1:2222896/1‑101 (MQ=255)
                                                                      cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGGCT.G.TCCTCGCT.CCCCCCGCGGGGGCGGGGCTCGAAAGCTgtgg                                >  2:4184789/1‑101 (MQ=25)
                                                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAAcaccc                                  >  2:4984721/1‑97 (MQ=2)
                                                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAAcaccc                                  >  1:2790962/1‑97 (MQ=2)
                                                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCtatg                                 <  1:155601/101‑4 (MQ=2)
                                                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCtatg                                 >  1:3715007/1‑98 (MQ=2)
                                                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCtatg                                 <  1:1057341/101‑4 (MQ=2)
                                                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTgcggt                               <  2:4565843/101‑2 (MQ=255)
                                                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTgcggc                               <  1:618059/101‑1 (MQ=255)
                                                                        cTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAAcacccg                                 >  2:2897307/1‑96 (MQ=2)
                                                                        cTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCtatga                                <  1:4987852/101‑5 (MQ=2)
                                                                         cccccccGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAAcacccgg                                <  2:2086059/100‑7 (MQ=0)
                                                                          ccccccGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAAcacccggc                               >  1:2159312/1‑94 (MQ=0)
                                                                          ccccccGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAAcacccggc                               <  2:4338879/101‑8 (MQ=0)
                                                                           cccccGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAAcacccggcg                              >  2:3771739/1‑93 (MQ=2)
                                                                           cccccGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT..CCCCCGCGGGGGCGTGGATCGAAACCCGCGGCGGGAg                          >  2:3491625/1‑101 (MQ=255)
                                                                             cccGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGGCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGt                         <  1:4623008/101‑1 (MQ=255)
                                                                               cgcgGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAgatc                       >  2:5017295/1‑98 (MQ=255)
                                                                                 cgGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGcc                     <  1:2679592/101‑1 (MQ=255)
                                                                                    gggCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAgg                  >  2:2848348/1‑100 (MQ=255)
                                                                                       cGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGAt               <  1:4326667/101‑1 (MQ=255)
                                                                                        gTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATa              >  2:1974341/1‑101 (MQ=255)
                                                                                         gggATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT..CCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGa            <  1:3491890/100‑1 (MQ=255)
                                                                                           gATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGACAGAc           <  2:4273559/101‑1 (MQ=255)
                                                                                            aTCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACg          >  1:2985520/1‑101 (MQ=255)
                                                                                            aTCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACg          <  1:444478/101‑1 (MQ=255)
                                                                                               gAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCg       >  1:3327564/1‑101 (MQ=255)
                                                                                               gAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCg       >  2:2298444/1‑101 (MQ=255)
                                                                                                   cT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGa   >  2:2219080/1‑101 (MQ=255)
                                                                                                     |                                                                                                            
GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGAC  >  NC_014034/1395631‑1395831

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG