Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,396,024 C→A 13.9% intergenic (+836/‑95) 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,396,0240C→A13.9% 103.0 ‑1.2 ‑0.0295intergenic (+836/‑95)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (16/25):  ref base (122/132):  total (138/157)
Fisher's exact test strand distribution p-value = 3.15e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.83e-01

Only 100 of 315 total aligned reads displayed.

TCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGAT  >  NC_014034/1395924‑1396124
                                                                                                    |                                                                                                       
tCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGc                                                                                                         <  1:1041999/101‑1 (MQ=255)
 cGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGcc                                                                                                        >  1:987302/1‑101 (MQ=255)
  gAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCt                                                                                                       <  2:4235324/101‑1 (MQ=255)
      cGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.cg                                                                                                  <  1:4381923/101‑1 (MQ=255)
       gATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.cgc                                                                                                 <  1:4312263/101‑1 (MQ=255)
        aTGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.cgc..g                                                                                              >  1:1973249/1‑101 (MQ=255)
           cGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCg                                                                                           >  2:544082/1‑101 (MQ=255)
           cGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCg                                                                                           >  1:3027118/1‑101 (MQ=255)
            gTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGc                                                                                          <  1:4495840/101‑1 (MQ=255)
              tCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTc                                                                                        >  1:637735/1‑101 (MQ=255)
              tCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTc                                                                                        >  2:3013636/1‑101 (MQ=255)
               cATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTcc                                                                                       >  1:1375475/1‑101 (MQ=255)
                aTGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTccc                                                                                      <  2:731998/101‑1 (MQ=255)
                 tGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTcccc                                                                                     >  1:1492346/1‑101 (MQ=255)
                 tGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTcccc                                                                                     <  1:1852969/101‑1 (MQ=255)
                 tGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTcccc                                                                                     <  1:2990033/101‑1 (MQ=255)
                  ggCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTccccc                                                                                    >  1:3504857/1‑101 (MQ=255)
                   gCCTACGTCGGGCAGGGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTcccccc                                                                                   <  2:3459970/101‑1 (MQ=255)
                    ccTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCcg                                                                                  <  1:4007528/101‑1 (MQ=255)
                    ccTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCcg                                                                                  >  2:3259117/1‑101 (MQ=255)
                    ccTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCcg                                                                                  >  2:2674060/1‑101 (MQ=255)
                      tACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCcgcg                                                                                >  2:4411338/1‑101 (MQ=255)
                        cGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCggg                                                                              >  2:2747184/1‑101 (MQ=255)
                            gggCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCg                                                                          <  2:2166966/101‑1 (MQ=255)
                            gggCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCg                                                                          >  2:388777/1‑101 (MQ=255)
                              gCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTg                                                                        <  1:38598/101‑1 (MQ=255)
                               cAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTgg                                                                       <  1:2065759/101‑1 (MQ=255)
                                aGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGa                                                                      >  1:363901/1‑101 (MQ=255)
                                  tGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                    >  1:3096508/1‑101 (MQ=255)
                                  tGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                    <  1:456814/101‑1 (MQ=255)
                                   tgcgcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTGCCTCGATCCGTTC.AGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                    <  1:5002561/96‑1 (MQ=2)
                                   ggACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                   >  1:3664543/1‑101 (MQ=255)
                                   ggACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCATCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                   >  2:4480835/1‑101 (MQ=255)
                                    gcgcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                   <  1:2924219/97‑1 (MQ=9)
                                    gcgcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                   <  2:819353/97‑1 (MQ=2)
                                    gcgcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGCGGATCg                                                                   >  1:1170307/5‑101 (MQ=2)
                                    gACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                  <  1:1595402/101‑1 (MQ=255)
                                      cTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                >  2:4675284/1‑101 (MQ=255)
                                      cTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                <  1:1775693/101‑1 (MQ=255)
                                       tGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                               <  1:893138/101‑1 (MQ=255)
                                        tgtccccccccccGGGGGGGCGTCGATCGAAACTCCCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                                              <  2:3558888/94‑1 (MQ=37)
                                        gtttgCTCCCCCCGCGTGGGCGTTGATCGAAACTCGCCCTCGATCCGCGCTTGCCCCACGCCCTC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                                              <  1:3373145/97‑1 (MQ=37)
                                        ggTCCCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                                              <  2:895575/101‑1 (MQ=255)
                                         gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                                              >  2:495815/1‑100 (MQ=2)
                                         gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAg                                                             >  2:179940/1‑101 (MQ=255)
                                           cGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGgc                                                           <  1:4210369/101‑1 (MQ=255)
                                           cGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGgc                                                           >  1:1905336/1‑101 (MQ=255)
                                            gCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtga                                                           <  2:4002224/101‑5 (MQ=2)
                                             cTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGgcgc                                                         >  2:1130162/1‑101 (MQ=255)
                                              tCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtgatg                                                         >  1:2409037/1‑95 (MQ=2)
                                               ccccccGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtgatga                                                        >  2:2586279/1‑94 (MQ=2)
                                               ccccccGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCAt                                                       >  2:2490563/1‑101 (MQ=255)
                                                 ccccGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtgatgacc                                                      <  1:4878672/101‑10 (MQ=2)
                                                 ccccGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCt                                                     <  2:4444974/101‑1 (MQ=255)
                                                   ccGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGa                                                   <  2:4332613/101‑1 (MQ=255)
                                                     gcgGGGGCGTGGATCGACACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATg                                                 <  2:4538307/101‑1 (MQ=255)
                                                       gggggCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATgga                                               >  2:1699548/1‑99 (MQ=255)
                                                       gggggCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGta                                               <  2:1060102/101‑1 (MQ=255)
                                                       gggggCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGta                                               <  1:1690687/101‑1 (MQ=255)
                                                       gggggCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGta                                               >  2:581701/1‑101 (MQ=255)
                                                        ggggCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGtat                                              >  1:1104785/1‑101 (MQ=255)
                                                        ggggCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGtat                                              <  1:4420703/101‑1 (MQ=255)
                                                        ggggCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGtat                                              <  1:662823/101‑1 (MQ=255)
                                                         gggCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGtata                                             <  1:565828/101‑1 (MQ=255)
                                                          ggCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTatat                                            <  1:2756562/101‑1 (MQ=255)
                                                          ggCGTGGATCGAAACGCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTatat                                            <  1:2184365/101‑1 (MQ=255)
                                                           gCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATg                                           <  2:3941850/101‑1 (MQ=255)
                                                            cGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGc                                          <  1:464256/101‑1 (MQ=255)
                                                             gTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGcc                                         >  2:281649/1‑101 (MQ=255)
                                                              tGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGccc                                        <  2:3130025/101‑1 (MQ=255)
                                                               ggATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCg                                       >  2:1765770/1‑101 (MQ=255)
                                                                gATCGAATCTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCgg                                      <  2:4152976/101‑1 (MQ=255)
                                                                  tCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTg                                    <  1:2679661/101‑1 (MQ=255)
                                                                   cGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTatatagatcggaa                                   >  1:3224743/1‑92 (MQ=255)
                                                                   cGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTga                                   >  2:118703/1‑101 (MQ=255)
                                                                   cGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTga                                   >  1:4805902/1‑101 (MQ=255)
                                                                   cGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTga                                   >  2:489527/1‑101 (MQ=255)
                                                                   cGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTga                                   >  1:738544/1‑101 (MQ=255)
                                                                   cGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTga                                   >  2:3712323/1‑101 (MQ=255)
                                                                    gAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTgag                                  <  1:1765523/101‑1 (MQ=255)
                                                                    gAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTgag                                  <  2:1620685/101‑1 (MQ=255)
                                                                      accTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCGCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGAc                                <  2:2194806/99‑1 (MQ=255)
                                                                      aaCTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGAc                                <  2:3319218/101‑1 (MQ=255)
                                                                          cGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCt                            >  2:2345658/1‑101 (MQ=255)
                                                                           gCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTg                           <  2:1413001/101‑1 (MQ=255)
                                                                             ccTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTc                         <  2:1289970/101‑1 (MQ=255)
                                                                               tCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGc                       <  2:1973105/101‑1 (MQ=255)
                                                                                cGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCt                      >  2:13508/1‑101 (MQ=255)
                                                                                 gATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTc                     <  1:121512/101‑1 (MQ=255)
                                                                                    ccGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTcccc                  >  1:2493150/1‑101 (MQ=255)
                                                                                    ccGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTcccc                  >  1:232131/1‑101 (MQ=255)
                                                                                      gcgcTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTcccccc                <  1:2239181/101‑1 (MQ=255)
                                                                                       cgcTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCcg               <  2:1492216/101‑1 (MQ=255)
                                                                                         ctcgCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCcgcg             <  2:1314725/98‑1 (MQ=255)
                                                                                           tGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCggg           >  2:754837/1‑101 (MQ=255)
                                                                                              cTCACGCCTGC.CGC..GTCGCTCCCCCCGCTGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGc        >  1:1270765/1‑101 (MQ=255)
                                                                                               tCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCg       >  1:3435647/1‑101 (MQ=255)
                                                                                                  cGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTgg    <  1:3161726/101‑1 (MQ=255)
                                                                                                    ccTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGAt  <  2:4805576/101‑1 (MQ=255)
                                                                                                    ccTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGAt  <  2:1043259/101‑1 (MQ=255)
                                                                                                    |                                                                                                       
TCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGAT  >  NC_014034/1395924‑1396124

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