Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,397,902 A→G 2.6% intergenic (+206/+1692) RCAP_rcc01307 → / ← RCAP_rcc01309 hypothetical protein/IS66 family transposase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,397,9020A→G 2.6% 2.5 ‑0.9 ‑0.0216intergenic (+206/+1692)RCAP_rcc01307/RCAP_rcc01309hypothetical protein/IS66 family transposase
Reads supporting (aligned to +/- strand):  new base (4/2):  ref base (94/116):  total (98/118)
Fisher's exact test strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 234 total aligned reads displayed.

TCCAGCGTCACATGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC  >  NC_014034/1397802‑1398002
                                                                                                       |                                                                                                    
tCCAGCGTCACATGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACa                                                                                                      <  1:2717283/101‑1 (MQ=255)
tCCAGCGTCACATGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACa                                                                                                      <  2:2958926/101‑1 (MQ=255)
  cAGCGTCACATGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGt                                                                                                    >  1:2296803/1‑101 (MQ=255)
  cAGCGTCACATGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGt                                                                                                    <  1:3098531/101‑1 (MQ=255)
   aGCGTCACATGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTc                                                                                                   <  1:77548/101‑1 (MQ=255)
    gCGTCACATGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCg                                                                                                  <  2:2792113/101‑1 (MQ=255)
       tCACATGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTc                                                                                               <  1:2123664/101‑1 (MQ=255)
       tCACATGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTc                                                                                               <  2:496797/101‑1 (MQ=255)
        cacaTGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTcc                                                                                              <  2:3173277/101‑1 (MQ=255)
           aTGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTccccc                                                                                           <  2:2301919/101‑1 (MQ=255)
            tGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTcccccc                                                                                          <  1:4411794/101‑1 (MQ=255)
            tGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTcccccc                                                                                          <  1:4506552/101‑1 (MQ=255)
             gCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTccccccc                                                                                         >  1:4106671/1‑100 (MQ=255)
             gCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCcg                                                                                         <  1:3502049/101‑1 (MQ=255)
               aGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCcgcg                                                                                       >  2:4329540/1‑101 (MQ=255)
                 cggcgTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCggg                                                                                     >  2:3390576/1‑101 (MQ=255)
                  ggcgTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCgggg                                                                                    <  1:386495/101‑1 (MQ=255)
                  ggcgTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCgggg                                                                                    >  1:1592604/1‑101 (MQ=255)
                  ggcgTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCgggg                                                                                    >  1:1195844/1‑101 (MQ=255)
                     gTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCg                                                                                 <  1:1082764/101‑1 (MQ=255)
                       cGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTg                                                                               <  1:1051275/101‑1 (MQ=255)
                         cATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGa                                                                             >  1:2923185/1‑101 (MQ=255)
                          aTCGCTATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                            <  1:2969641/101‑1 (MQ=255)
                          aTCGCCATGTCGCTCCCCCCGCGGGGGCGTGGAGCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                            <  1:4032821/101‑1 (MQ=255)
                          aTCGCCATGTCGCTCCCACCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                            >  2:2815646/1‑101 (MQ=255)
                           tCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                           <  1:3766457/101‑1 (MQ=255)
                           tCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                           <  2:1731624/101‑1 (MQ=255)
                           tCGCCATGTCGCTCCCCCCGCGGGGGCGTCGATCAAACC.CCGGCAC..CTCAGCCCCCCCCCCTGTCGCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                           <  2:3665611/101‑1 (MQ=25)
                           cgatcttGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                           <  2:839584/95‑1 (MQ=255)
                            cGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                          >  2:4550602/1‑101 (MQ=255)
                               ggcctCGCCCCCCCCGCGGGGGCGTGGATCGAACC.CCGCCACGTCTTGGAACAGCGTGCTAT.TGGCC.CCGGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                       <  2:1082909/97‑1 (MQ=2)
                                aTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                      >  2:3483505/1‑101 (MQ=255)
                                aTGTCGCCCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                      <  1:2438615/101‑1 (MQ=255)
                                 tGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGGGGGTcagatcg                                                                     >  1:839660/1‑95 (MQ=37)
                                  gTCGCTCCCCCCGCGGGGGCGTGGATCGAAGC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCCCCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActc                                                                    <  2:3562664/101‑1 (MQ=255)
                                  gTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGGTCGAAActc                                                                    >  1:859420/1‑101 (MQ=255)
                                   tCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActct                                                                   <  2:219608/101‑1 (MQ=255)
                                    cGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCCCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTa                                                                  <  1:4619721/101‑1 (MQ=255)
                                       tCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATgg                                                               >  1:2785321/1‑101 (MQ=255)
                                       tCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATgg                                                               >  1:290704/1‑101 (MQ=255)
                                         cccccGCGGGGGCGTGGATCGAAACACCAATGC..CGAAGCCAAGGCCGCAGCAGCGTCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcaaaaagcc                                                              >  1:3924347/1‑93 (MQ=0)
                                         cccccGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGTATCGAAACTCTATggga                                                             >  2:4193135/1‑101 (MQ=255)
                                         cccccGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATggga                                                             <  1:2836061/101‑1 (MQ=255)
                                         cccccGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATggga                                                             <  2:678333/101‑1 (MQ=255)
                                          ccccGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATgggag                                                            <  1:4329832/101‑1 (MQ=255)
                                          ccccGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATgggag                                                            <  1:4176870/101‑1 (MQ=255)
                                          ccccGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATgggag                                                            <  1:3793499/101‑1 (MQ=255)
                                              gcgGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATgggaggcta                                                        >  2:4114109/1‑101 (MQ=255)
                                               cgGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATgggaggctag                                                       >  2:623629/1‑101 (MQ=255)
                                                gggggCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATgggaggctagg                                                      >  2:3688212/1‑101 (MQ=255)
                                                  gggCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATgggaggctaggga                                                    <  2:2532301/101‑1 (MQ=255)
                                                    gCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATgggaggctagggagg                                                  <  2:4118583/101‑1 (MQ=255)
                                                    gCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATgggaggctagggagg                                                  <  1:3633921/101‑1 (MQ=255)
                                                     cGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATgggcggctagggaggc                                                 >  2:1153129/1‑101 (MQ=255)
                                                     cGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATgggaggctagggaggc                                                 >  1:4278602/1‑101 (MQ=255)
                                                     cGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATgggaggctagggaggc                                                 >  2:1444525/1‑101 (MQ=255)
                                                      gTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATgggaggctagggaggct                                                <  1:4462904/101‑1 (MQ=255)
                                                       tGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATgggaggctagggaggcta                                               >  2:4452062/1‑101 (MQ=255)
                                                       tGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATgggaggctagggaggcta                                               <  1:47978/101‑1 (MQ=255)
                                                        ggATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTAc                                              <  1:1966589/101‑1 (MQ=255)
                                                         gATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACg                                             <  1:1475859/101‑1 (MQ=255)
                                                          aTCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGc                                            >  1:520079/1‑101 (MQ=255)
                                                            cGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCAc                                          >  2:1295068/1‑101 (MQ=255)
                                                             gAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCAcc                                         >  2:899211/1‑101 (MQ=255)
                                                               aaC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCtt                                       >  1:3054151/1‑101 (MQ=255)
                                                                aC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTc                                      <  2:764942/101‑1 (MQ=255)
                                                                aC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTc                                      <  1:1905722/101‑1 (MQ=255)
                                                                aC.CCGGCAC..CTCAGCCCCGCCCGCTGTATCACCAACAGCCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCtcc                                      <  2:4081661/101‑3 (MQ=255)
                                                                 c.ccGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCa                                     <  1:3972157/101‑1 (MQ=255)
                                                                   ccGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCaa                                    <  2:1939937/101‑1 (MQ=255)
                                                                   ccGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCaa                                    >  2:1917612/1‑101 (MQ=255)
                                                                    cGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCaaa                                   >  2:3588667/1‑101 (MQ=255)
                                                                     ggCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGCGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCaaaa                                  <  2:4030613/101‑1 (MQ=255)
                                                                      gCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAg                                 >  1:1536681/1‑101 (MQ=255)
                                                                      gCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAg                                 <  2:228437/101‑1 (MQ=255)
                                                                       cAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGt                                <  2:1547175/101‑1 (MQ=255)
                                                                       cAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGt                                >  1:2221695/1‑101 (MQ=255)
                                                                         c..cTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCg                              >  1:2637292/1‑101 (MQ=255)
                                                                              cAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTc                           <  2:1592499/101‑1 (MQ=255)
                                                                               aGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTcc                          <  1:1827272/101‑1 (MQ=255)
                                                                                gCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTccc                         >  1:2015765/1‑101 (MQ=255)
                                                                                 ccccGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTcccc                        >  1:970395/1‑101 (MQ=255)
                                                                                  cccgcccgCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTccccc                       <  1:720301/101‑1 (MQ=255)
                                                                                  cccgcccgCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTccccc                       <  1:623690/101‑1 (MQ=255)
                                                                                   ccgcccgcTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTcccccc                      <  2:3243515/101‑1 (MQ=255)
                                                                                    cgcccgcTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCcg                     <  1:419774/101‑1 (MQ=255)
                                                                                       ccgcTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCCCgg                  <  2:3951668/101‑1 (MQ=255)
                                                                                         gcTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCgggg                <  2:4277920/101‑1 (MQ=255)
                                                                                             tctcACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCgg            >  2:538225/1‑100 (MQ=255)
                                                                                             tctcACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGt            >  1:1000050/1‑101 (MQ=255)
                                                                                              ctcACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTg           <  1:1743379/101‑1 (MQ=255)
                                                                                               tcACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTgg          >  1:235897/1‑101 (MQ=255)
                                                                                                 aCCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGAt        <  2:4619011/101‑1 (MQ=255)
                                                                                                 aCCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGAt        >  1:3383482/1‑101 (MQ=255)
                                                                                                    aaCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGa     >  1:2527/1‑101 (MQ=255)
                                                                                                      cAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa   >  2:4506025/1‑101 (MQ=255)
                                                                                                      cAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa   <  2:251827/101‑1 (MQ=255)
                                                                                                       aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGGGTGGATCGAAAc  >  2:3236093/1‑101 (MQ=255)
                                                                                                       aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc  <  2:4348779/101‑1 (MQ=255)
                                                                                                       aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAAGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc  >  2:2383897/1‑101 (MQ=255)
                                                                                                       |                                                                                                    
TCCAGCGTCACATGCAGCGGCGTCGCATCGCCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGGCAC..CTCAGCCCCGCCCGCTGTCTCACCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC  >  NC_014034/1397802‑1398002

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG