Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,687 G→C 5.7% intergenic (+499/‑432) 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,6870G→C 5.7% 37.2 0.0 ‑0.5375intergenic (+499/‑432)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (12/9):  ref base (196/155):  total (209/166)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.52e-01

Only 100 of 429 total aligned reads displayed.

GTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGG  >  NC_014034/1395587‑1395787
                                                                                                      |                                                                                                     
gTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTGGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGACGGCGG..GCAg                                                                                                       >  2:3326100/1‑101 (MQ=255)
gTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAg                                                                                                       >  1:116810/1‑101 (MQ=255)
 tGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGt                                                                                                      <  2:3791909/101‑1 (MQ=255)
 tGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGt                                                                                                      >  1:6544004/1‑101 (MQ=255)
   gATCGAAACCGCCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTc                                                                                                    <  1:5544895/101‑1 (MQ=255)
   gATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTc                                                                                                    <  1:1351721/101‑1 (MQ=255)
   gATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTc                                                                                                    <  2:1693857/101‑1 (MQ=255)
   gATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTc                                                                                                    <  1:6798527/101‑1 (MQ=255)
   gATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTc                                                                                                    <  1:2178424/101‑1 (MQ=255)
    aTCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCg                                                                                                   <  1:4351392/101‑1 (MQ=255)
    aTCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCg                                                                                                   <  2:2671088/101‑1 (MQ=255)
    aTCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCg                                                                                                   <  1:968453/101‑1 (MQ=255)
    aTCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCg                                                                                                   <  1:1870073/101‑1 (MQ=255)
    aTCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCg                                                                                                   >  1:6719439/1‑101 (MQ=255)
    aTCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCg                                                                                                   <  1:5132367/101‑1 (MQ=255)
     tCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGCAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTTCATGTAGTCGGCGG..GCAGTTCGa                                                                                                  >  2:832486/1‑101 (MQ=255)
       gAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.Gc                                                                                               >  2:2596562/1‑101 (MQ=255)
       gAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.Gc                                                                                               >  2:441068/1‑101 (MQ=255)
          accaccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGc                                                                                            <  1:6141863/101‑1 (MQ=255)
           ccaccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCa                                                                                           >  2:6783074/1‑101 (MQ=255)
            caccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAATTCGA.GCCGCAg                                                                                          <  1:5076069/101‑1 (MQ=255)
             cccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGt                                                                                         <  1:4308332/100‑1 (MQ=255)
                ttGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGc                                                                                      >  1:240367/1‑101 (MQ=255)
                ttGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGc                                                                                      >  2:6102644/1‑101 (MQ=255)
                  gCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTc                                                                                    >  2:5240648/1‑101 (MQ=255)
                   cGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTcc                                                                                   >  1:4086726/1‑101 (MQ=255)
                    ggATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTccc                                                                                  >  1:6875185/1‑101 (MQ=255)
                      gttctccctctCGCGGGGCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTccccc                                                                                <  2:2277073/91‑1 (MQ=38)
                       ttGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTcccccc                                                                               <  2:476159/101‑1 (MQ=255)
                       ttGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTcccccc                                                                               >  1:6412029/1‑101 (MQ=255)
                        tGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGTATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCcg                                                                              >  1:4609438/1‑101 (MQ=255)
                        tGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCACTCGCTCCCCCcg                                                                              <  1:6268872/101‑1 (MQ=255)
                         gTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCcgc                                                                             >  2:6723088/1‑101 (MQ=255)
                             cTTTCGCGGGACGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCgggg                                                                         >  2:356006/1‑101 (MQ=255)
                              tttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCggggg                                                                        >  1:4673702/1‑101 (MQ=255)
                               ttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGc                                                                       >  1:1844897/1‑101 (MQ=255)
                               ttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGc                                                                       >  2:344877/1‑101 (MQ=255)
                                  tggtgccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGC.GAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGt                                                                     >  1:2251706/11‑101 (MQ=2)
                                    ggggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGa                                                                  >  1:2488067/1‑101 (MQ=255)
                                    gggccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGC.GAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTgg                                                                   <  1:1446459/93‑1 (MQ=2)
                                      tctgAAGCGGCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                <  2:1805845/98‑1 (MQ=255)
                                      gccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGC.GAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                                 <  2:4224277/95‑1 (MQ=2)
                                      gccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGC.GAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                                 >  2:5553273/7‑101 (MQ=2)
                                        cGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                              >  1:1221912/1‑101 (MQ=255)
                                         gctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGC.GAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                              <  1:4546228/98‑1 (MQ=2)
                                         gctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGC.GAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                              >  1:3129202/4‑101 (MQ=2)
                                         gAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                             >  1:105771/1‑101 (MQ=255)
                                           tgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGC.GAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                            >  1:6690887/2‑101 (MQ=2)
                                           tgCGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACACCCCGAACCTGTCGGCGGC.GAACACCGATACCGC.GTCGCTCCCGCCGCGGGGGCGTGGATCGagc                                                            >  2:5848258/2‑99 (MQ=2)
                                             cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCGGCAGTGGCAGGTTTTCCA.GCGGC.GTCGCTCCCCCCGCGGGGGCGTGGAGCGAAAcc                                                          >  2:6877772/1‑100 (MQ=2)
                                             cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTg                                                         >  2:2073894/1‑101 (MQ=255)
                                              gTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTgg                                                        <  1:2600059/101‑1 (MQ=255)
                                                 gCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGCGGATCGGAACTGGTAt                                                     >  2:5088065/1‑101 (MQ=255)
                                                 ccTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTAt                                                     <  2:1867012/100‑1 (MQ=255)
                                                     cccccGGGGGGGCGGGGATCGAAACATCTGCATGTAGTCGGGGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGAt                                                 <  1:783163/101‑1 (MQ=255)
                                                        ccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTg                                              <  1:6005930/101‑1 (MQ=255)
                                                         cgcgGGGGCGTGGATCGAACCATCTGAATGGAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTgg                                             <  2:1308168/101‑1 (MQ=255)
                                                         cgcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTgg                                             <  2:4579664/101‑1 (MQ=255)
                                                           cgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgcg                                           >  2:1231782/1‑101 (MQ=255)
                                                            gggggCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgcgc                                          >  2:4899301/1‑101 (MQ=255)
                                                              gggCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGGATGGATCTGGCGCaa                                        >  1:3304920/1‑101 (MQ=255)
                                                               ggCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCaaa                                       >  2:3937816/1‑101 (MQ=255)
                                                                gCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAAc                                      <  2:2900898/101‑1 (MQ=255)
                                                                gCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAAc                                      <  2:3035281/101‑1 (MQ=255)
                                                                gCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAAc                                      <  1:3238987/101‑1 (MQ=255)
                                                                  gTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACgg                                    >  2:1461481/1‑101 (MQ=255)
                                                                  gTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACgg                                    >  1:2791315/1‑101 (MQ=255)
                                                                   gggATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTGGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGa                                   <  2:6298452/100‑1 (MQ=255)
                                                                    ggATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaa                                  <  2:2487973/101‑1 (MQ=255)
                                                                    ggATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaa                                  >  2:4091292/1‑101 (MQ=255)
                                                                       tCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGt                               <  2:2811739/101‑1 (MQ=255)
                                                                       tCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGt                               >  2:6261544/1‑101 (MQ=255)
                                                                        cGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTc                              <  1:3937991/101‑1 (MQ=255)
                                                                        cGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTc                              >  2:6994369/1‑101 (MQ=255)
                                                                          aaaCATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTg                            >  1:4041107/1‑101 (MQ=255)
                                                                           ccgaTCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgt                           <  2:1195826/98‑1 (MQ=255)
                                                                           aaCATCTGCATGTAGGCGGGGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgt                           <  1:1107614/101‑1 (MQ=255)
                                                                              aTCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTNTGGATCTGGCGCAAACGGAAAGTCTgtcgt                        >  1:1984913/1‑101 (MQ=255)
                                                                              aTCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgt                        >  1:6627045/1‑101 (MQ=255)
                                                                              aTCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgt                        >  1:3850826/1‑101 (MQ=255)
                                                                                cTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgtcc                      >  1:1668764/1‑100 (MQ=255)
                                                                                 tGCATGTGGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGc                     <  2:1619603/101‑1 (MQ=255)
                                                                                 tGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGc                     >  2:5899082/1‑101 (MQ=255)
                                                                                  gCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCt                    >  2:1930570/1‑101 (MQ=255)
                                                                                   cATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTc                   <  2:3992258/101‑1 (MQ=255)
                                                                                     tattgtgggcgg..ACAGTCCGA.GCCGCAGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTccc                 <  2:3199969/94‑1 (MQ=37)
                                                                                      gTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTcccc                >  1:5689069/1‑101 (MQ=255)
                                                                                       tAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTccccc               <  2:6461011/101‑1 (MQ=255)
                                                                                        agccggggg..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTcccccc              <  2:6594548/98‑1 (MQ=255)
                                                                                         ggcggcgg..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCcg             <  2:4673463/99‑1 (MQ=255)
                                                                                         gTCGGCGG..GCAGTTCGA.GCCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCcg             <  1:3553332/101‑1 (MQ=255)
                                                                                            ggcgg..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCgg          >  2:5031797/1‑101 (MQ=255)
                                                                                              cgg..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCgggg        <  1:4921565/101‑1 (MQ=255)
                                                                                                g..gCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGc      <  1:1536426/101‑1 (MQ=255)
                                                                                                g..gCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGc      <  1:995056/101‑1 (MQ=255)
                                                                                                    cAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCgg    >  2:566021/1‑100 (MQ=255)
                                                                                                    cAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGt    >  1:4569271/1‑101 (MQ=255)
                                                                                                     aGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTg   >  1:5938016/1‑101 (MQ=255)
                                                                                                     aGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTg   >  1:5317281/1‑101 (MQ=255)
                                                                                                     aGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTg   >  1:2979459/1‑101 (MQ=255)
                                                                                                      |                                                                                                     
GTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGG  >  NC_014034/1395587‑1395787

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG