Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,406 T→G 4.0% intergenic (+218/‑713) 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,4060T→G 4.0% 20.3 ‑0.6 ‑0.4306intergenic (+218/‑713)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (6/7):  ref base (168/124):  total (174/132)
Fisher's exact test strand distribution p-value = 5.69e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.39e-01

Only 100 of 327 total aligned reads displayed.

AGCACCTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTC  >  NC_014034/1395306‑1395506
                                                                                                    |                                                                                                    
aGCACCTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCt                                                                                                      <  1:1527860/101‑1 (MQ=255)
 gCACCTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTg                                                                                                     <  2:3624990/101‑1 (MQ=255)
   aCCTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGc                                                                                                   >  1:4458463/1‑101 (MQ=255)
    ccTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGcc                                                                                                  <  1:1190543/101‑1 (MQ=255)
    ccTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGcc                                                                                                  >  1:3476683/1‑101 (MQ=255)
      ttCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCt                                                                                                <  1:1651224/101‑1 (MQ=255)
       tCGCACCACCGGAGCAAGTTACCCCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCCGCTGGCCCTg                                                                                               <  2:667190/101‑1 (MQ=255)
           accaccGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcgg                                                                                           <  1:5691619/101‑1 (MQ=255)
            ccaccGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggt                                                                                          <  2:4985008/101‑1 (MQ=255)
            ccaccGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggt                                                                                          <  2:5240162/101‑1 (MQ=255)
             caccGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtg                                                                                         >  2:3168307/1‑101 (MQ=255)
              agcGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATTTAAACTACACCCGCAGGCACGGCcgtcac                                                                                        >  2:3093789/3‑97 (MQ=37)
              accGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtga                                                                                        >  1:2043506/1‑101 (MQ=255)
                 ggAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatcg                                                                                     >  1:2249814/1‑101 (MQ=255)
                  gAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatcgg                                                                                    >  2:4095465/1‑101 (MQ=255)
                    gCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatcggtg                                                                                  <  1:3755202/101‑1 (MQ=255)
                    gCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatcggtg                                                                                  >  1:5241104/1‑101 (MQ=255)
                     cAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatcggtga                                                                                 >  2:1456616/1‑101 (MQ=255)
                      aaTTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggggatcggggat                                                                                >  1:2100420/1‑101 (MQ=255)
                          tACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATGGAAACTGCACCCGCTGGCCCTGCCCGGGATGGGGGAGCAGc                                                                            >  1:2649580/1‑101 (MQ=37)
                          tACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGc                                                                            >  2:2206604/1‑101 (MQ=255)
                            cacCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGGGATCAGCgg                                                                          >  2:4154398/1‑101 (MQ=255)
                               cGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCg                                                                       >  1:4366309/1‑101 (MQ=255)
                                 tCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCt                                                                     >  1:3249064/1‑101 (MQ=255)
                                  cTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGTGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTc                                                                    <  2:2280774/101‑1 (MQ=255)
                                  cTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTc                                                                    >  2:1646437/1‑101 (MQ=255)
                                  cTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTc                                                                    >  1:242816/1‑101 (MQ=255)
                                   tGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTcc                                                                   >  2:2233055/1‑101 (MQ=255)
                                   tGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGCGATCGGAGATCAGCGGTCGCTcc                                                                   >  1:3884149/1‑101 (MQ=255)
                                     gTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTcccc                                                                 >  2:2385594/1‑101 (MQ=255)
                                     gTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTcccc                                                                 >  1:5411197/1‑101 (MQ=255)
                                      tCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTccccc                                                                >  2:5537793/1‑101 (MQ=255)
                                      tCAGCGGGTTACGAAACGGCAGTCGCTCCCACCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCCGCGGTCGCTccccc                                                                >  1:5146193/1‑101 (MQ=255)
                                       cAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTcccccc                                                               >  1:4827706/1‑101 (MQ=255)
                                        aGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCcg                                                              >  1:768846/1‑101 (MQ=255)
                                        aGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCcg                                                              >  1:990406/1‑101 (MQ=255)
                                         gCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCcgc                                                             >  1:2778129/1‑101 (MQ=255)
                                           gggTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCgg                                                           <  2:4366043/101‑1 (MQ=255)
                                           gggTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCgg                                                           >  2:4709153/1‑101 (MQ=255)
                                             gTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCgggg                                                         >  1:1826115/1‑101 (MQ=255)
                                               tACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGc                                                       <  1:3168496/101‑1 (MQ=255)
                                                aCGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCg                                                      >  2:518280/1‑101 (MQ=255)
                                                  gAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCggg                                                    >  2:3797330/1‑99 (MQ=255)
                                                  gAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTg                                                    >  2:1951123/1‑101 (MQ=255)
                                                  gAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTg                                                    >  2:5295402/1‑101 (MQ=255)
                                                  gAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTg                                                    >  2:5053152/1‑101 (MQ=255)
                                                   aaaCGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGGgg                                                   >  2:4789620/1‑101 (MQ=255)
                                                    aaCGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGa                                                  >  1:3397763/1‑101 (MQ=255)
                                                    aaCGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGGGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGa                                                  >  1:1983220/1‑101 (MQ=255)
                                                     aCGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGAt                                                 <  2:730497/101‑1 (MQ=255)
                                                       gggaGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTGCAGCCACTGGTCCTGCCGGAGATCGGTGATCAGCGATCGCTCCCCACGCGAGGGCGGGGATca                                               >  2:2882346/4‑100 (MQ=25)
                                                          aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGgaa                                            >  2:445033/1‑101 (MQ=255)
                                                          aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCCCCGGGGCGGGGATCGaaa                                            >  2:2747744/1‑101 (MQ=255)
                                                           gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                           >  1:4164790/1‑101 (MQ=255)
                                                           gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                           >  1:4201485/1‑101 (MQ=255)
                                                             cGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                         >  1:1330323/1‑101 (MQ=255)
                                                              gCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATc                                        >  1:2834775/1‑101 (MQ=255)
                                                              gCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATc                                        <  1:353659/101‑1 (MQ=255)
                                                               cTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcg                                       >  1:2944526/1‑101 (MQ=255)
                                                                tCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgc                                      >  2:5127264/1‑101 (MQ=255)
                                                                cccccccGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgc                                      <  2:3476466/100‑1 (MQ=255)
                                                                 ccccccGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgcc                                     >  2:2561932/1‑101 (MQ=255)
                                                                      cgcgGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCttt                                >  1:3965018/1‑101 (MQ=255)
                                                                        cgGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTgg                              >  2:4576944/1‑101 (MQ=255)
                                                                         gggggCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGc                             >  1:5033205/1‑101 (MQ=255)
                                                                         gggggCGTGGACCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGc                             <  2:5603733/101‑1 (MQ=255)
                                                                          ggggCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCg                            >  2:423767/1‑101 (MQ=255)
                                                                           gggCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCgg                           >  1:4556070/1‑101 (MQ=255)
                                                                           gggCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCgg                           >  2:3384338/1‑101 (MQ=255)
                                                                           gggCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCgg                           >  1:861833/1‑101 (MQ=255)
                                                                             gCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCCCTTTGGCGGTc                         >  1:1277530/1‑101 (MQ=255)
                                                                              cgcgcagcGAGACTGCACCGGCTGGCCCTGCCGGTGAACGGTTATCAGCGGTCCCCCCCCCCGCGGGGGCGTGGATCGAATCATCGCCGCTTTGGAGATCg                        <  2:266407/94‑1 (MQ=25)
                                                                              cGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCg                        <  1:5335018/101‑1 (MQ=255)
                                                                                tGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGAc                      >  2:1554486/1‑101 (MQ=255)
                                                                                tGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGAc                      >  2:1087922/1‑101 (MQ=255)
                                                                                tGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGAc                      <  2:5371263/101‑1 (MQ=255)
                                                                                 ggATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACa                     >  1:5774280/1‑101 (MQ=255)
                                                                                 ggATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTCGCGGTCGACa                     >  2:3030113/1‑101 (MQ=255)
                                                                                  gATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAg                    >  2:2130315/1‑101 (MQ=255)
                                                                                  gATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAg                    >  1:3027284/1‑101 (MQ=255)
                                                                                   aTCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGGGGTCGCTCCCCCCGCGGGGACGTTGATCGAAACATCGCCGCTTTGGCGGTCGGCAgg                   >  1:4949447/1‑101 (MQ=38)
                                                                                   aTCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAgg                   >  2:129517/1‑101 (MQ=255)
                                                                                   aTCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAgg                   >  2:3062263/1‑101 (MQ=255)
                                                                                     cGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCt                 <  2:3248864/101‑1 (MQ=255)
                                                                                      gAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTGGACAGGcgg                >  2:798409/1‑99 (MQ=255)
                                                                                      gAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTg                <  1:4074273/101‑1 (MQ=255)
                                                                                      gAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTg                >  2:5676491/1‑101 (MQ=255)
                                                                                       aaaCTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGCGGACAGGCTgg               >  1:2584412/1‑101 (MQ=255)
                                                                                        aaCTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGt              >  2:570184/1‑101 (MQ=255)
                                                                                         aCTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTa             <  1:4667175/101‑1 (MQ=255)
                                                                                         aCTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGGCGACAGGCTGGTa             >  1:5560209/1‑101 (MQ=255)
                                                                                          cTGCACCCGCTGGCCCTGCCGGTGGTCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGGGTGGACCGAAACCTCGCCGCTTTGGCGGTCGCCAGGCCCGTaa            >  2:2300974/1‑101 (MQ=37)
                                                                                          cTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTaa            >  1:637597/1‑101 (MQ=255)
                                                                                             cACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACgt         <  1:1008518/101‑1 (MQ=255)
                                                                                                 cGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGGAACGTGTCg     >  2:3961971/1‑101 (MQ=255)
                                                                                                  gCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGc    <  1:1075538/101‑1 (MQ=255)
                                                                                                  gCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGAGGTCGACAGGCTGGTAACGTGTCGc    <  1:5627598/101‑1 (MQ=255)
                                                                                                   cTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCt   >  1:1746312/1‑101 (MQ=255)
                                                                                                   cTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCt   <  1:144944/101‑1 (MQ=255)
                                                                                                   cTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGGAAAGTGTCGCt   >  2:1030972/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
AGCACCTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTC  >  NC_014034/1395306‑1395506

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