Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,411 C→G 5.0% intergenic (+223/‑708) 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,4110C→G 5.0% 20.3 ‑1.4 ‑0.0256intergenic (+223/‑708)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (9/4):  ref base (121/122):  total (130/126)
Fisher's exact test strand distribution p-value = 2.55e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.57e-01

Only 100 of 281 total aligned reads displayed.

CTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCC  >  NC_014034/1395311‑1395511
                                                                                                    |                                                                                                     
cTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGccc                                                                                                       <  1:1745609/101‑1 (MQ=255)
cTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGccc                                                                                                       <  1:1714047/101‑1 (MQ=255)
 ttCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCt                                                                                                      <  1:3969469/101‑1 (MQ=255)
   cGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGc                                                                                                    >  1:3032383/1‑101 (MQ=255)
    gCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGcc                                                                                                   <  1:3376779/101‑1 (MQ=255)
     caccacCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcg                                                                                                  >  2:4755104/1‑101 (MQ=255)
       ccaccGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCACCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGAcggg                                                                                                >  2:1944473/1‑100 (MQ=255)
        caccGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtg                                                                                               <  2:1842893/101‑1 (MQ=255)
         accGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtga                                                                                              <  1:2062073/101‑1 (MQ=255)
             gAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCCCCGCTGGCCCTGCcggtgatcgg                                                                                          <  1:1466262/101‑1 (MQ=255)
             gAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatcgg                                                                                          >  2:4795798/1‑101 (MQ=255)
             gAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatcgg                                                                                          >  2:2262903/1‑101 (MQ=255)
               gCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatcggtg                                                                                        >  2:1600595/1‑101 (MQ=255)
                cAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatcggtg.a                                                                                      >  2:3137152/1‑101 (MQ=255)
                 aaTTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatcggtg.at                                                                                     <  1:2448610/101‑1 (MQ=255)
                 aaTTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggcgatcggtg.at                                                                                     >  2:522397/1‑101 (MQ=255)
                  aTTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCCCCGCTGGCCCTGCCggtgatcggtg.atc                                                                                    <  1:4950614/101‑1 (MQ=255)
                  aTTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCggtgatcggtg.atc                                                                                    <  1:4406872/101‑1 (MQ=255)
                  aTTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCggtgatcggtg.atc                                                                                    <  1:851536/101‑1 (MQ=255)
                   tttACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCGCTGCCGGTGATGGGGG.ATCa                                                                                   >  1:2532430/1‑101 (MQ=255)
                   tttACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCa                                                                                   <  1:2230667/101‑1 (MQ=255)
                   tttACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCa                                                                                   <  1:656725/101‑1 (MQ=255)
                   tttACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCa                                                                                   <  1:3605008/101‑1 (MQ=255)
                   tttAAACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCa                                                                                   >  2:1226659/1‑101 (MQ=255)
                    ttACCCCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAg                                                                                  <  2:270826/101‑1 (MQ=255)
                    ttACCCCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAg                                                                                  <  1:3726104/101‑1 (MQ=255)
                           gTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGc                                                                           >  2:2898794/1‑101 (MQ=255)
                             cTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGGG.ATCAGCGGTCGCTc                                                                         >  2:3444396/1‑101 (MQ=255)
                              tGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTcc                                                                        <  2:4531187/101‑1 (MQ=255)
                              tGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGGGATCGGTG.ATCAGCGGGCGCTcc                                                                        >  2:1319656/1‑101 (MQ=255)
                               ggTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCCGCGGTCGCTccc                                                                       >  2:1008494/1‑101 (MQ=255)
                                  cAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTcccccc                                                                    <  2:4301920/101‑1 (MQ=255)
                                  cAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTcccccc                                                                    >  2:3305349/1‑101 (MQ=255)
                                  cAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTcccccc                                                                    >  2:638569/1‑101 (MQ=255)
                                  cAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGGGATCGGTG.ATCAGCGGTCGCTcccccc                                                                    >  1:444908/1‑101 (MQ=255)
                                    gCGGGTTACGAAACGGCAGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACTGCTCCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCcgc                                                                  <  1:1944610/101‑1 (MQ=255)
                                       ggTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCggg                                                               >  2:1717411/1‑101 (MQ=255)
                                        gtgaCGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCTCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCgggg                                                              <  2:826987/98‑1 (MQ=255)
                                          tACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGc                                                            >  1:2124701/1‑101 (MQ=255)
                                          tACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGc                                                            >  2:3228192/1‑101 (MQ=255)
                                          tACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCCCGGGGGc                                                            >  2:4819848/1‑101 (MQ=255)
                                            cGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGt                                                          <  1:1388938/101‑1 (MQ=255)
                                             gAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTg                                                         <  1:4117460/101‑1 (MQ=255)
                                             gAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTg                                                         <  2:21498/101‑1 (MQ=255)
                                             gAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTg                                                         <  2:3787406/101‑1 (MQ=255)
                                             gAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGGCGCTCCCCCCGCGGGGGCGTg                                                         >  1:1932015/1‑101 (MQ=255)
                                              tcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.GCTGGCGCGCGCCGCGCCG..GATCTTTGTCTGCGCGGTCGCTCCCCCCGCGGGGGCGTGGAt                                                      >  2:3656624/9‑101 (MQ=2)
                                               caacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.GCTGGCGCGCGCCGCGCCG..GATCTTTGTCTGCGCGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                     >  2:2846639/8‑101 (MQ=2)
                                               caacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.GCTGGCGCGCGCCGCGCCG..GATCTTTGTCTGCGCGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                     >  1:10638/8‑101 (MQ=2)
                                               aaCGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGa                                                       <  1:4357941/101‑1 (MQ=255)
                                               aaCGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGAG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGa                                                       >  1:4311281/1‑101 (MQ=255)
                                                aCGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTACCGGTGATTGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGAt                                                      >  2:4017572/1‑101 (MQ=255)
                                                aCGGCAGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTGCCCCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGAt                                                      <  1:486740/101‑1 (MQ=255)
                                                  ggCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                    >  1:3854686/1‑101 (MQ=255)
                                                  ggCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                    <  2:5050779/101‑1 (MQ=255)
                                                  ggCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                    <  1:1600727/101‑1 (MQ=255)
                                                  cactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.GCTGGCGCGCGCCGCGCCG..GATCTTTGTCTGCGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                  >  1:3271127/5‑101 (MQ=2)
                                                   gCAGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                   <  2:4376073/101‑1 (MQ=255)
                                                    cAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGGGGATCGaa                                                  >  1:2983161/1‑101 (MQ=255)
                                                     tgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.GCTGGCGCGCGCCGCGCCG..GATCTTTGTCTGCGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                               >  1:3547550/2‑100 (MQ=2)
                                                     aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                 >  1:247630/1‑101 (MQ=255)
                                                      gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                <  1:4546390/101‑1 (MQ=255)
                                                       tCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                               <  2:3000025/101‑1 (MQ=255)
                                                         gCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATc                                             <  2:3242635/101‑1 (MQ=255)
                                                         gCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATc                                             >  1:1292418/1‑101 (MQ=255)
                                                          ccccccccGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcg                                            <  2:2330225/99‑1 (MQ=255)
                                                          ccccccccGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcg                                            <  1:518709/99‑1 (MQ=255)
                                                           tCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgc                                           >  2:119696/1‑101 (MQ=255)
                                                             cccccGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgccg                                         <  1:625598/101‑1 (MQ=255)
                                                                  gcgGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTg                                    <  1:3137390/101‑1 (MQ=255)
                                                                   cgGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTgg                                   >  2:4940518/1‑101 (MQ=255)
                                                                     ggggCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCg                                 <  2:4220938/101‑1 (MQ=255)
                                                                       ggCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCggg                               >  2:3548410/1‑100 (MQ=255)
                                                                             gATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAg                         >  1:90282/1‑101 (MQ=255)
                                                                             gATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAg                         <  2:3721125/101‑1 (MQ=255)
                                                                              aTCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAgg                        <  1:1157456/101‑1 (MQ=255)
                                                                                 gAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTg                     <  1:2880644/101‑1 (MQ=255)
                                                                                    aCTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTa                  >  2:1690843/1‑101 (MQ=255)
                                                                                     cTGCACCCGCTGGCCCCGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTaa                 >  1:149226/1‑101 (MQ=255)
                                                                                       gcccccGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACg               <  2:4082931/98‑1 (MQ=255)
                                                                                       gCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACg               >  1:3891483/1‑101 (MQ=255)
                                                                                       gCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACg               <  1:2589447/101‑1 (MQ=255)
                                                                                       gCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCCACAGGCTGGTAACg               >  2:1554630/1‑101 (MQ=255)
                                                                                        cACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTCACgt              <  2:167720/101‑1 (MQ=255)
                                                                                           ccGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTc           >  1:1097574/1‑101 (MQ=255)
                                                                                            cccTGTCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCg          <  2:2532249/99‑1 (MQ=255)
                                                                                              cTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTcccc        >  2:1791523/1‑98 (MQ=255)
                                                                                               tGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTc       >  2:3934627/1‑101 (MQ=255)
                                                                                                ggCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTcc      >  1:3254709/1‑101 (MQ=255)
                                                                                                ggCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTcc      >  2:4438311/1‑101 (MQ=255)
                                                                                                ggCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTcc      <  1:275460/101‑1 (MQ=255)
                                                                                                ggCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGGGTCGCTcc      >  2:1100560/1‑101 (MQ=255)
                                                                                                ggCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGATGGTAACGTGTCGCTcc      >  2:719928/1‑101 (MQ=255)
                                                                                                ggCCCTGCCGGTCATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTcc      >  1:725385/1‑101 (MQ=255)
                                                                                                 gCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTccc     <  1:2291375/101‑1 (MQ=255)
                                                                                                 gCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTccc     <  1:4778921/101‑1 (MQ=255)
                                                                                                  cccTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTcccc    >  2:4610795/1‑101 (MQ=255)
                                                                                                  cccTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTcccc    >  2:1893571/1‑101 (MQ=255)
                                                                                                   ccTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTccccc   >  1:2344285/1‑101 (MQ=255)
                                                                                                   ccTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTccccc   <  2:90272/101‑1 (MQ=255)
                                                                                                    |                                                                                                     
CTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCC  >  NC_014034/1395311‑1395511

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