Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,406 T→G 8.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 8.0% 38.8 ‑0.4 ‑0.3243intergenic (+218/‑713)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (12/8):  ref base (120/103):  total (132/111)
Fisher's exact test strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.90e-01

Only 100 of 266 total aligned reads displayed.

AGCACCTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTC  >  NC_014034/1395306‑1395506
                                                                                                    |                                                                                                     
      ttCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCt                                                                                                 >  1:2979026/1‑101 (MQ=255)
       tCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTg                                                                                                <  2:1314627/101‑1 (MQ=255)
       tCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGATACGAAACGGCAGTCGCTCCCCCCTCGGGGGCGTGGATCGGAAATGCACCCGCCAGCCCTg                                                                                                >  2:3978216/1‑101 (MQ=255)
        cGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGc                                                                                               >  2:2965934/1‑101 (MQ=255)
         gCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGcc                                                                                              >  2:3750790/1‑101 (MQ=255)
            ccaccGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggt                                                                                           <  2:183921/101‑1 (MQ=255)
             caccGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGAGGATCGAAACTGCACCCGCTGGCCCTGCcggtg                                                                                          >  1:3998650/1‑101 (MQ=255)
              cccGGAGCAATTTCCCCCGTCTGGTCAGCGGGTTACGAACCGGCAGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtga                                                                                         <  1:3883878/100‑1 (MQ=38)
              accGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtga                                                                                         >  2:1696756/1‑101 (MQ=255)
               ccGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgat                                                                                        <  1:98353/101‑1 (MQ=255)
                cGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatc                                                                                       <  1:1524043/101‑1 (MQ=255)
                 ggAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatcg                                                                                      >  1:1986137/1‑101 (MQ=255)
                  gAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatcgg                                                                                     >  2:1164164/1‑101 (MQ=255)
                  gAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatcgg                                                                                     <  2:396676/101‑1 (MQ=255)
                   aGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatcggt                                                                                    >  1:2986407/1‑101 (MQ=255)
                     cAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatcggtg.a                                                                                 >  1:531958/1‑101 (MQ=255)
                     cAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggcgatcggtg.a                                                                                 >  1:3469367/1‑101 (MQ=255)
                        tttACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCa                                                                              <  1:1468023/101‑1 (MQ=255)
                        tttACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGGGATCGGTG.ATCa                                                                              >  2:4031935/1‑101 (MQ=255)
                         ttACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGGCCTGCCGGTGATCGGTG.ATCAg                                                                             >  1:4574207/1‑101 (MQ=255)
                         ttACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAg                                                                             <  1:2484104/101‑1 (MQ=255)
                         ttACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAg                                                                             >  1:2888223/1‑101 (MQ=255)
                         ttACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAg                                                                             >  1:4300910/1‑101 (MQ=255)
                         ttACACCGGCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAg                                                                             <  2:20289/101‑1 (MQ=255)
                           acacCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCg                                                                           >  2:1791397/1‑101 (MQ=255)
                           acacCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCg                                                                           >  2:3457547/1‑101 (MQ=255)
                             cccGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGt                                                                         <  1:1106706/100‑1 (MQ=255)
                             acCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCggc                                                                         >  2:3178039/1‑100 (MQ=255)
                             acCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGt                                                                         >  1:394741/1‑101 (MQ=255)
                               cGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCg                                                                       <  1:3351839/101‑1 (MQ=255)
                               cGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCg                                                                       <  1:823746/101‑1 (MQ=255)
                                gTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGc                                                                      <  1:1874865/101‑1 (MQ=255)
                                gTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGc                                                                      >  2:2300353/1‑101 (MQ=255)
                                 tCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCt                                                                     >  1:182349/1‑101 (MQ=255)
                                      tCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTccccc                                                                >  1:3200532/1‑101 (MQ=255)
                                      tCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTccccc                                                                <  1:2381168/101‑1 (MQ=255)
                                         gCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCcgc                                                             <  1:3417436/101‑1 (MQ=255)
                                          cGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCcgcg                                                            <  2:3143124/101‑1 (MQ=255)
                                          cGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCcgcg                                                            >  2:4319255/1‑101 (MQ=255)
                                             gTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCgggg                                                         <  2:1986002/101‑1 (MQ=255)
                                             gTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCgggg                                                         <  1:2338731/101‑1 (MQ=255)
                                                  gAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTg                                                    >  2:30664/1‑101 (MQ=255)
                                                  gAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTg                                                    <  1:3873749/101‑1 (MQ=255)
                                                   tcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.GCTGGCGCGCGCCGCGCCG..GATCTTTGTCTGCGCGGTCGCTCCCCCCGCGGGGGCGTGGAt                                                 >  1:3692723/9‑101 (MQ=2)
                                                   tcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.GCTGGCGCGCGCCGCGCCG..GATCTTTGTCTGCGCGGTCGCTCCCCCCGCGGGGGCGTGGAt                                                 >  1:4524058/9‑101 (MQ=2)
                                                   aaaCGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.CTCAGCGGGCGCTCCCCCCCCGGGGGCGTgg                                                   >  2:4519730/1‑101 (MQ=255)
                                                    aaCGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGa                                                  >  1:4494640/1‑101 (MQ=255)
                                                    aaCGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGa                                                  <  2:4180447/101‑1 (MQ=255)
                                                      cGGCAGTCGCTCGCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                >  1:2531988/1‑101 (MQ=255)
                                                       ggCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                               >  1:413753/1‑101 (MQ=255)
                                                       cactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.GCTGGCGCGCGCCGCGCCG..GATCTTTGTCTGCGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                             >  2:1610708/5‑101 (MQ=2)
                                                        gCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCCCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                              <  1:749198/101‑1 (MQ=255)
                                                         cAGACGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGATGGACCTGCCGGTGATCGGTG.ATCAGCGGTCGCACCCCCCGCGGGGGCGTGGATCGaa                                             >  2:438207/1‑101 (MQ=38)
                                                          tgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.GCTGGCGCGCGCCGCGCCG..GATCTTTGTCTGCGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                          <  2:2221824/100‑2 (MQ=2)
                                                          tgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.GCTGGCGCGCGCCGCGCCG..GATCTTTGTCTGCGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                          >  2:1752721/2‑100 (MQ=2)
                                                          aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                            >  2:3251534/1‑101 (MQ=255)
                                                          aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                            >  2:2555926/1‑101 (MQ=255)
                                                          aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                            <  1:4347029/101‑1 (MQ=255)
                                                          aGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                            <  1:4633246/101‑1 (MQ=255)
                                                            tCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                          <  2:182325/101‑1 (MQ=255)
                                                            tCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                          <  1:1755816/101‑1 (MQ=255)
                                                             cGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                         >  1:592815/1‑101 (MQ=255)
                                                             accccaccccGTGGGGGACTGGATCGGAACGGCACCCCCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                         <  2:711214/95‑1 (MQ=25)
                                                                tCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgc                                      <  2:1389185/101‑1 (MQ=255)
                                                                tCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgc                                      >  2:234818/1‑101 (MQ=255)
                                                                 ccccccGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgcc                                     <  1:1019681/101‑1 (MQ=255)
                                                                 ccccccGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgcc                                     <  2:3002593/101‑1 (MQ=255)
                                                                       gcgGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTg                               >  1:2463130/1‑101 (MQ=255)
                                                                        cgGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTgg                              >  1:1191417/1‑101 (MQ=255)
                                                                         gggggCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGc                             <  1:359693/101‑1 (MQ=255)
                                                                          ggggCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCg                            >  1:2100205/1‑101 (MQ=255)
                                                                             gCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTc                         <  1:3508257/101‑1 (MQ=255)
                                                                              cGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCg                        <  2:4435355/101‑1 (MQ=255)
                                                                              cGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCg                        >  1:1789802/1‑101 (MQ=255)
                                                                               gtcgATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGa                       <  2:3931072/98‑1 (MQ=255)
                                                                               gTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGa                       >  1:2759213/1‑101 (MQ=255)
                                                                               gTGGATCGAAACCGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGa                       <  2:2742043/101‑1 (MQ=255)
                                                                                tGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGGGTGGATCGAAACATCGCCGCTTTGGCGGTAGAc                      >  1:1643886/1‑101 (MQ=255)
                                                                                tGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTTGGGGTCggc                      >  1:4392100/1‑99 (MQ=255)
                                                                                tGGATCGAAACTGCACCCGCTGGCCCTGCCGGCGATCGGTG.ATCAGCGGTCCCTCCCCCCGCGGGGGGGTGGATCGAAACATCGCCGCTTTGGCGGGCGAc                      >  1:623051/1‑101 (MQ=255)
                                                                                 ggATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACa                     <  1:2348909/101‑1 (MQ=255)
                                                                                 ggATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACa                     >  2:3524531/1‑101 (MQ=255)
                                                                                  gATCGAAACTGCCCCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAg                    <  1:3081036/101‑1 (MQ=255)
                                                                                  gATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAg                    >  2:2815872/1‑101 (MQ=255)
                                                                                    tCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGc                  <  2:2984023/101‑1 (MQ=255)
                                                                                     cGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCt                 >  1:1286704/1‑101 (MQ=255)
                                                                                        aaCTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGt              <  2:4587914/101‑1 (MQ=255)
                                                                                        aaCTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGt              <  2:1450177/101‑1 (MQ=255)
                                                                                         aCTGCCCCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTa             <  1:1164256/101‑1 (MQ=255)
                                                                                         aCTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTa             >  2:3229170/1‑101 (MQ=255)
                                                                                           tGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAc           >  2:3542194/1‑101 (MQ=255)
                                                                                           tGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAc           <  1:4260757/101‑1 (MQ=255)
                                                                                           tGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAc           <  1:4413331/101‑1 (MQ=255)
                                                                                             cACCCGCTGGCCCTGCCGGTGATCGGGG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGGGGATCGAAAGATCGCCGCTTTGGCGGCTTCCAGGCTggataagc         >  2:2461677/1‑95 (MQ=25)
                                                                                               cccGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACgtgt       <  2:2888007/101‑1 (MQ=255)
                                                                                                ccGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTc      <  2:60773/101‑1 (MQ=255)
                                                                                                  gCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGc    <  2:3362020/101‑1 (MQ=255)
                                                                                                    tGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTc  >  1:402515/1‑101 (MQ=255)
                                                                                                    tGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTc  >  2:20757/1‑101 (MQ=255)
                                                                                                    tGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTc  >  1:924648/1‑101 (MQ=255)
                                                                                                    |                                                                                                     
AGCACCTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTC  >  NC_014034/1395306‑1395506

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