Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,962 C→A 2.0% intergenic (+774/‑157) 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,9620C→A 2.0% 5.0 ‑2.2 ‑0.0261intergenic (+774/‑157)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (1/5):  ref base (132/122):  total (134/127)
Fisher's exact test strand distribution p-value = 1.13e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.81e-01

Only 100 of 282 total aligned reads displayed.

ACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAA  >  NC_014034/1395862‑1396062
                                                                                                    |                                                                                                    
aCTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGAc                                                                                                      <  1:4702660/101‑1 (MQ=255)
   gcgcCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGGGGACTgg                                                                                                   >  1:5558568/1‑101 (MQ=255)
    cgccgcCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGt                                                                                                  >  2:4410877/1‑101 (MQ=255)
     gccgccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTc                                                                                                 >  1:2850069/1‑101 (MQ=255)
      tccgatctCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCg                                                                                                <  2:838141/94‑1 (MQ=255)
      ccgccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCg                                                                                                >  1:5634770/1‑101 (MQ=255)
      ccgccTATCGCATCCGGTCGAGCGCCGCGCAGTCCCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCg                                                                                                <  1:4972221/101‑1 (MQ=255)
       cgccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGc                                                                                               >  2:1938377/1‑101 (MQ=255)
          cagccGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTcc                                                                                            >  2:580664/5‑101 (MQ=255)
          cTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTcc                                                                                            >  1:1582974/1‑101 (MQ=255)
          cTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTcc                                                                                            <  2:1157876/101‑1 (MQ=255)
          cTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTcc                                                                                            >  1:2432327/1‑101 (MQ=255)
           tATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTccc                                                                                           >  1:59089/1‑101 (MQ=255)
            aTCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTcccc                                                                                          >  1:5357499/1‑101 (MQ=255)
             tcccATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTccccc                                                                                         <  1:1696145/98‑1 (MQ=255)
             tCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTccccc                                                                                         <  1:2707817/101‑1 (MQ=255)
              cacATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTcccccc                                                                                        <  1:258944/99‑1 (MQ=255)
              cGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCgagatcgga                                                                                        >  1:838209/1‑93 (MQ=255)
              cGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTcccccc                                                                                        >  2:3498968/1‑101 (MQ=255)
              cGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTcccccc                                                                                        <  2:471978/101‑1 (MQ=255)
              cGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTcccccc                                                                                        <  2:2838893/101‑1 (MQ=255)
               gCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCcg                                                                                       <  2:3121597/101‑1 (MQ=255)
               gCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCcg                                                                                       >  1:915915/1‑101 (MQ=255)
               gCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCcg                                                                                       >  1:3256612/1‑101 (MQ=255)
               gCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCGCCCCCcg                                                                                       >  1:2649059/1‑101 (MQ=255)
               gCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGACGGGCAGTGGACTGGTCGCTCCCCCcg                                                                                       >  1:3327905/1‑101 (MQ=255)
                cATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTcccccccc                                                                                      >  2:2613685/1‑99 (MQ=255)
                 aTCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCcgcg                                                                                     >  2:5768935/1‑101 (MQ=255)
                 aTCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCcgcg                                                                                     <  2:172947/101‑1 (MQ=255)
                 aTCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCcgcg                                                                                     <  2:4071010/101‑1 (MQ=255)
                      gTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGc                                                                                >  2:3759963/1‑101 (MQ=255)
                         gAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTg                                                                             <  2:431219/101‑1 (MQ=255)
                           gcgcCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGa                                                                           <  1:3515086/101‑1 (MQ=255)
                            gaacaggatgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGCTCGACCGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                          <  1:4820825/92‑1 (MQ=2)
                            cgccgcGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                          >  2:3853994/1‑101 (MQ=255)
                            cgccgcGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                          >  2:5335564/1‑101 (MQ=255)
                            cgccgcGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                          >  1:3077662/1‑101 (MQ=255)
                             gccgcGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                         >  1:2057085/1‑101 (MQ=255)
                                 cgcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                     <  2:4206024/101‑1 (MQ=255)
                                  gcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                    <  1:689988/101‑1 (MQ=255)
                                     gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCg                                                                 <  2:772034/101‑1 (MQ=255)
                                      gcGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGc                                                                <  1:115015/100‑1 (MQ=255)
                                       cGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAATCGcc                                                               >  2:1117689/1‑101 (MQ=255)
                                         cTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCt                                                             >  2:1312031/1‑101 (MQ=255)
                                         cTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGAGGGGGCGTGGATCGAAACTCGcccc                                                             >  2:1025410/1‑100 (MQ=255)
                                          tCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTc                                                            >  1:4583139/1‑101 (MQ=255)
                                          tCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTc                                                            <  2:3967558/101‑1 (MQ=255)
                                          tCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTc                                                            >  2:457168/1‑101 (MQ=255)
                                          tCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTc                                                            <  1:794596/101‑1 (MQ=255)
                                           tccgctgtgggggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCg                                                           <  1:2888963/93‑1 (MQ=255)
                                           ccccccGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCg                                                           >  2:2132558/1‑101 (MQ=255)
                                             ccccGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGAt                                                         <  1:3984056/101‑1 (MQ=255)
                                             ccccGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGAt                                                         >  1:3710649/1‑101 (MQ=255)
                                               ccGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATcc                                                       >  1:5299046/1‑101 (MQ=255)
                                               ccGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATcc                                                       >  1:5480695/1‑101 (MQ=255)
                                                     gggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCtt                                                 >  1:4342230/1‑101 (MQ=255)
                                                     gggCGTGGATCGAAACGATGCGTTCATGGCCCACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCtt                                                 >  1:5163584/1‑101 (MQ=255)
                                                      ggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTg                                                <  1:1277307/101‑1 (MQ=255)
                                                       gCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGc                                               >  1:4908667/1‑101 (MQ=255)
                                                          tGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTc                                            >  2:1251179/1‑101 (MQ=255)
                                                          tGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTc                                            >  1:905913/1‑101 (MQ=255)
                                                          tGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTc                                            >  1:4285857/1‑101 (MQ=255)
                                                             aTCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACg                                         >  2:4832365/1‑101 (MQ=255)
                                                             aTCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACg                                         <  1:3450828/101‑1 (MQ=255)
                                                             aTCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACg                                         >  2:4445474/1‑101 (MQ=255)
                                                              tCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGc                                        <  2:3077488/101‑1 (MQ=255)
                                                               cGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGcc                                       <  2:2056968/101‑1 (MQ=255)
                                                                 aaaCGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTg                                     >  1:4184517/1‑101 (MQ=255)
                                                                 aaaCGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTg                                     <  2:5586936/101‑1 (MQ=255)
                                                                  aaCGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCCTGCCTCACGCCTGc                                    >  2:2493019/1‑101 (MQ=255)
                                                                    cGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCcg                                  <  1:238145/101‑1 (MQ=255)
                                                                    cGAAGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCcg                                  <  1:1938485/101‑1 (MQ=255)
                                                                     gATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCcgc                                 <  1:2132675/101‑1 (MQ=255)
                                                                      aTGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCcgcg                                <  1:2963119/101‑1 (MQ=255)
                                                                      aTGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCCCGGGGGCGTGGATCGAAACCCGCCCTCGATCCGCGCTTGCCTCACGCCTGCcgcg                                >  1:2249449/1‑101 (MQ=255)
                                                                       tGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGt                               <  1:1326572/101‑1 (MQ=255)
                                                                          gTTCATGGCCTACGGCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGc                            <  2:629268/101‑1 (MQ=255)
                                                                           ttCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCt                           >  2:953081/1‑101 (MQ=255)
                                                                           ttCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCt                           >  2:5294532/1‑101 (MQ=255)
                                                                              aTGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTccc                        >  2:1637258/1‑101 (MQ=255)
                                                                              aTGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTccc                        <  2:5163298/101‑1 (MQ=255)
                                                                                 gCCTCCGTGGGGCAGTGGGCTGGTCGCGCCCCCCGCGGGGGCGTGGATCGAATCCCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTcccccc                     <  2:3206465/101‑1 (MQ=37)
                                                                                   cTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCcgc                   >  2:3286711/1‑101 (MQ=255)
                                                                                    tACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCcgcg                  <  2:5213403/101‑1 (MQ=255)
                                                                                       gTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCgggg               <  2:3686590/101‑1 (MQ=255)
                                                                                       gTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCgggg               <  1:3242463/101‑1 (MQ=255)
                                                                                         cGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCGGGGGc             >  1:3645473/1‑101 (MQ=255)
                                                                                           atctgTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCCCGCCTGCCGCGTCGCTCCCCCCGCGGGGGCGt           <  1:1444162/97‑1 (MQ=255)
                                                                                            gCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCGGGGGCGTg          <  1:3854236/101‑1 (MQ=255)
                                                                                             cAGTCGAGTTGTCGCTCCCCCCGCGGGGGCGTGGACCGAAACTCGCCCCCGCTCCGCGCTGGCCTCACCCCTGCCGCGTCTCTCCCCCCGCGGGGGCGTgg         <  2:3723268/101‑1 (MQ=25)
                                                                                              aggggccTGGTCGCCCCCCCCGCGGGGGCGTGGATCGAACCTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCGGGGGCGTGGa        <  1:1007181/95‑1 (MQ=38)
                                                                                               gTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCGGGGGCGTAGAt       >  2:1444069/1‑101 (MQ=255)
                                                                                                tGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCGGGGGCGTGGATc      <  2:3341375/101‑1 (MQ=255)
                                                                                                gggACTGGGCCCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCGGGGGCGTGGATc      <  2:4938782/100‑1 (MQ=255)
                                                                                                 ggACTGGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCGGGGGCGTGGATCg     <  2:2395372/101‑1 (MQ=255)
                                                                                                  gACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGa    <  2:2857706/101‑1 (MQ=255)
                                                                                                   aCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa   <  1:1526755/101‑1 (MQ=255)
                                                                                                   aCTGGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa   <  1:2493161/101‑1 (MQ=255)
                                                                                                   aCTGGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa   <  1:4371265/101‑1 (MQ=255)
                                                                                                    cTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa  <  2:2434526/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
ACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAA  >  NC_014034/1395862‑1396062

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