Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,417 Δ1 bp 6.0% intergenic (+229/‑702) 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,4170G→. 6.0% 31.1 ‑1.8 ‑0.0217intergenic (+229/‑702)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (9/4):  ref base (99/105):  total (108/109)
Fisher's exact test strand distribution p-value = 1.65e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.99e-01

Only 100 of 269 total aligned reads displayed.

CCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGG  >  NC_014034/1395318‑1395517
                                                                                                   |                                                                                                     
ccaccGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCACCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGAcggg                                                                                                      >  2:1944473/1‑100 (MQ=255)
 caccGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtg                                                                                                     <  2:1842893/101‑1 (MQ=255)
  accGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCCCGGGGGCGTCGATCGAAACTGCACCCCCCGGCCCCGCcgggga                                                                                                    >  1:4376390/1‑101 (MQ=37)
   ccGGAGCTATTTACTCTATCTGGTCAGCGGGTTACGAAACGCCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgat                                                                                                   <  2:4543096/101‑1 (MQ=37)
    cGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatc                                                                                                  <  2:1892528/101‑1 (MQ=255)
     ggAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatcg                                                                                                 >  2:2186697/1‑101 (MQ=255)
     ggAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatcg                                                                                                 >  2:3254195/1‑101 (MQ=255)
      gAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCCCCGCTGGCCCTGCcggtgatcgg                                                                                                <  1:1466262/101‑1 (MQ=255)
      gAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatcgg                                                                                                >  2:2262903/1‑101 (MQ=255)
      gAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatcgg                                                                                                >  2:4795798/1‑101 (MQ=255)
      gAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCCGCTGGCCCTGCcggtgatcgg                                                                                                >  1:4244214/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)
           aTTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCCCCGCTGGCCCTGCCggtgatcggtg.atc                                                                                          <  1:4950614/101‑1 (MQ=255)
           aTTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCggtgatcggtg.atc                                                                                          <  1:851536/101‑1 (MQ=255)
            tttACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCGCTGCCGGTGATGGGGG.ATCa                                                                                         >  1:2532430/1‑101 (MQ=255)
            tttACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCa                                                                                         <  1:656725/101‑1 (MQ=255)
            tttACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCa                                                                                         <  1:3605008/101‑1 (MQ=255)
             ttACCCCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAg                                                                                        <  1:3726104/101‑1 (MQ=255)
              tACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGc                                                                                       <  2:1296915/101‑1 (MQ=255)
              tACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGc                                                                                       >  1:3092451/1‑101 (MQ=255)
               acacCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCg                                                                                      >  1:182129/1‑101 (MQ=255)
                cacCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCgg                                                                                     >  1:1376203/1‑101 (MQ=255)
                 acCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGt                                                                                    <  1:4832032/101‑1 (MQ=255)
                    gTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGc                                                                                 >  2:2898794/1‑101 (MQ=255)
                      cTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGGG.ATCAGCGGTCGCTc                                                                               >  2:3444396/1‑101 (MQ=255)
                       tGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTcc                                                                              >  1:3247124/1‑101 (MQ=255)
                       tGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTcc                                                                              <  2:263749/101‑1 (MQ=255)
                         gTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTcccc                                                                            >  1:3336423/1‑101 (MQ=255)
                           cAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTcccccc                                                                          >  2:3305349/1‑101 (MQ=255)
                             gCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACTGCACCCGCTGGCCCTGCCGGGGATCGGGG.ATCAGCGGGCGCTcccccccc                                                                        >  1:951529/1‑99 (MQ=38)
                                ggTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCggg                                                                     >  2:1717411/1‑101 (MQ=255)
                                 gtgaCGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCTCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCgggg                                                                    <  2:826987/98‑1 (MQ=255)
                                  ttACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCggggg                                                                   >  2:1157339/1‑101 (MQ=255)
                                  ttACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCggggg                                                                   >  1:4444235/1‑101 (MQ=255)
                                   tACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGc                                                                  <  1:4064809/101‑1 (MQ=255)
                                   tACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGc                                                                  >  2:3228192/1‑101 (MQ=255)
                                   tACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGc                                                                  >  1:2124701/1‑101 (MQ=255)
                                   tACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCCCGGGGGc                                                                  >  2:4819848/1‑101 (MQ=255)
                                   tACGAAACGGCAGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGc                                                                  <  2:4018521/101‑1 (MQ=255)
                                      gAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTg                                                               <  2:3787406/101‑1 (MQ=255)
                                      gAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTg                                                               >  1:3955551/1‑101 (MQ=255)
                                       tcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.GCTGGCGCGCGCCGCGCCG..GATCTTTGTCTGCGCGGTCGCTCCCCCCGCGGGGGCGTGGAt                                                            >  1:4539766/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)
                                         aCGGCAGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTGCCCCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGAt                                                            <  1:486740/101‑1 (MQ=255)
                                          cGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGACCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                           >  2:3981342/1‑101 (MQ=255)
                                           ggCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                          <  2:5050779/101‑1 (MQ=255)
                                           ggCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                          >  1:1932425/1‑101 (MQ=255)
                                           ggCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                          <  1:1600727/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)
                                              aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.AGCAGCGGTCGCCCCCCCCGCGGGGGCGTGGATCGaaa                                                       <  2:4389555/101‑1 (MQ=255)
                                               gTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.GCTGGCGCGCGCCGCGCCG..GATCTTTGTCTGCGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActc                                                    <  1:4575325/101‑3 (MQ=2)
                                                tCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                                     <  1:1717541/101‑1 (MQ=255)
                                                tCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                                     <  2:3000025/101‑1 (MQ=255)
                                                  gCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATc                                                   >  1:1292418/1‑101 (MQ=255)
                                                   ccccccccGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcg                                                  <  1:518709/99‑1 (MQ=255)
                                                   ccccccccGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcg                                                  <  2:2330225/99‑1 (MQ=255)
                                                     ccccccGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgcc                                                >  2:1693212/1‑101 (MQ=255)
                                                     ccccccGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgcc                                                >  1:2082614/1‑101 (MQ=255)
                                                      cccccGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgccg                                               <  1:625598/101‑1 (MQ=255)
                                                          cgcgGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCttt                                           <  1:1497987/101‑1 (MQ=255)
                                                            cgGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTgg                                         >  2:4940518/1‑101 (MQ=255)
                                                              ggggCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCg                                       <  2:4220938/101‑1 (MQ=255)
                                                               gggCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCgg                                      >  1:3857868/1‑101 (MQ=255)
                                                                ggCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCggg                                     >  2:3548410/1‑100 (MQ=255)
                                                                ggCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGt                                     >  2:4768201/1‑101 (MQ=255)
                                                                 gCGTGGATCGAAACTGCCCCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTc                                    <  1:2286798/101‑1 (MQ=255)
                                                                    tGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGAc                                 <  2:537328/101‑1 (MQ=255)
                                                                      gATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAg                               <  1:2928791/101‑1 (MQ=255)
                                                                      gATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAg                               <  2:3721125/101‑1 (MQ=255)
                                                                       aTCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAgg                              <  1:1157456/101‑1 (MQ=255)
                                                                          gAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTg                           >  2:1218574/1‑101 (MQ=255)
                                                                           aaaCTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTgg                          >  1:2146742/1‑101 (MQ=255)
                                                                           aaaCTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTgg                          >  2:3597497/1‑101 (MQ=255)
                                                                             aCTGCCCCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTa                        <  1:3270020/101‑1 (MQ=255)
                                                                                gcccccGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACg                     <  2:4082931/98‑1 (MQ=255)
                                                                                gcccccGCGGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACg                     <  2:639599/98‑1 (MQ=255)
                                                                                gCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTaagg                     >  1:679051/1‑99 (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)
                                                                                     cGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCg                <  2:1097460/101‑1 (MQ=255)
                                                                                       cTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTcccc              >  2:1791523/1‑98 (MQ=255)
                                                                                        tGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTc             <  2:4310964/101‑1 (MQ=255)
                                                                                         ggCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTcc            >  2:4438311/1‑101 (MQ=255)
                                                                                         ggCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTcc            <  1:275460/101‑1 (MQ=255)
                                                                                         ggCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGATGGTAACGTGTCGCTcc            >  2:719928/1‑101 (MQ=255)
                                                                                            ccTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTccccc         >  1:2344285/1‑101 (MQ=255)
                                                                                            ccTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTccccc         <  2:90272/101‑1 (MQ=255)
                                                                                             cTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTcccccc        <  2:3854401/101‑1 (MQ=255)
                                                                                                ccGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCcgcg     >  1:993649/1‑101 (MQ=255)
                                                                                                 cggtgatcggtg.atCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCgg    >  2:1387923/1‑101 (MQ=255)
                                                                                                 cggtgatcggtg.atCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCgg    >  1:202781/1‑101 (MQ=255)
                                                                                                   |                                                                                                     
CCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGG  >  NC_014034/1395318‑1395517

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