Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,397 Δ1 bp 4.6% intergenic (+209/‑722) 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,3970T→. 4.6% 23.8 ‑1.4 ‑0.5241intergenic (+209/‑722)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (4/7):  ref base (127/103):  total (131/110)
Fisher's exact test strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.83e-01

Only 100 of 258 total aligned reads displayed.

CTTGCGACAAGCACCTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC  >  NC_014034/1395297‑1395497
                                                                                                    |                                                                                                     
cTTGCGACAAGCACCTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                                                                                       <  1:1439316/101‑1 (MQ=255)
cTTGCGACAAGCACCTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                                                                                       <  2:327006/101‑1 (MQ=255)
 ttGCGACAAGCACCTTCGCACCACCGGTGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCCCGGGGGCGTGGATCGAAACTg                                                                                                      >  1:1973563/1‑101 (MQ=255)
 ttGCGACAAGCACCTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTg                                                                                                      >  1:2940865/1‑101 (MQ=255)
 ttGCGACAAGCACCTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTg                                                                                                      <  1:4619648/101‑1 (MQ=255)
 ttGCGACAAGCACCTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTg                                                                                                      >  2:3023746/1‑101 (MQ=255)
 ttGCGACAAGCACCTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTg                                                                                                      >  1:3441418/1‑101 (MQ=255)
  tGCGACAAGCACCTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCAGGGGCGTGGATCGAAACTGc                                                                                                     <  1:1664328/101‑1 (MQ=255)
    cGACAAGCACCTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCAc                                                                                                   <  2:1844203/101‑1 (MQ=255)
     gACAAGCACCTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCAcc                                                                                                  <  1:3134873/101‑1 (MQ=255)
      aCAAGCACCTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCAccc                                                                                                 <  1:3717022/101‑1 (MQ=255)
        aaGCACCTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGc                                                                                               <  1:754359/101‑1 (MQ=255)
            aCCTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGc                                                                                           >  1:3895744/1‑101 (MQ=255)
            aCCTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGc                                                                                           >  1:1169551/1‑101 (MQ=255)
            aCCTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGc                                                                                           <  2:2908269/101‑1 (MQ=255)
             ccTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGcc                                                                                          <  2:1066353/101‑1 (MQ=255)
             ccTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCCCCGGcc                                                                                          >  1:903505/1‑101 (MQ=255)
              cTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGccc                                                                                         <  2:4518306/101‑1 (MQ=255)
               ttCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCt                                                                                        <  2:380919/101‑1 (MQ=255)
                 cGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGc                                                                                      <  1:2102787/101‑1 (MQ=255)
                 cGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGc                                                                                      <  2:2636976/101‑1 (MQ=255)
                  gCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGcc                                                                                     >  1:4085139/1‑101 (MQ=255)
                  gCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGcc                                                                                     <  2:2145914/101‑1 (MQ=255)
                  gCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGcc                                                                                     <  2:2947890/101‑1 (MQ=255)
                   caccacCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcg                                                                                    <  2:1719706/101‑1 (MQ=255)
                      caccGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcgggg                                                                                 >  1:1534327/1‑99 (MQ=255)
                        ccGGAGCAATTTACACCGTCTGGTCAGCGGGCTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgat                                                                               >  1:4232386/1‑101 (MQ=255)
                         cGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatc                                                                              <  1:2326153/101‑1 (MQ=255)
                          ggAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatcg                                                                             >  2:141892/1‑101 (MQ=255)
                           gAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatcgg                                                                            <  1:4027457/101‑1 (MQ=255)
                             gCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatcggtg                                                                          <  1:3023934/101‑1 (MQ=255)
                              cAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatcgggg.a                                                                        >  1:1240469/1‑101 (MQ=255)
                               aaTTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGCGGGGGGTCGAAACTGCACCGCTTGGCCCCGCcggtggccggtgtt                                                                        >  1:2891482/1‑99 (MQ=255)
                                aTTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCggtgatcggtg.atc                                                                      >  2:4142524/1‑101 (MQ=255)
                                aTTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCggtgatcggtg.atc                                                                      <  2:1541292/101‑1 (MQ=255)
                                aTTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCggtgatcggtg.atc                                                                      >  2:54218/1‑101 (MQ=255)
                                aTTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCggggatcgggg.atc                                                                      >  2:851077/1‑101 (MQ=255)
                                 tttACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCa                                                                     <  2:4181581/101‑1 (MQ=255)
                                 tttACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCa                                                                     <  1:3251464/101‑1 (MQ=255)
                                  ttACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAg                                                                    >  1:3446605/1‑101 (MQ=255)
                                    acacCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCg                                                                  >  1:2781062/1‑101 (MQ=255)
                                    acacCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGGG.ATCAGCg                                                                  >  2:2692648/1‑101 (MQ=255)
                                       ccGTCTGGTCAGCGTGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTc                                                               <  2:3118557/101‑1 (MQ=255)
                                       ccGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTc                                                               >  1:805163/1‑101 (MQ=255)
                                         gTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGc                                                             >  2:4690501/1‑101 (MQ=255)
                                         gTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGc                                                             <  2:901034/101‑1 (MQ=255)
                                         gTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGc                                                             >  1:1542723/1‑101 (MQ=255)
                                         gTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGc                                                             >  1:1128238/1‑101 (MQ=255)
                                           cTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTc                                                           <  2:2509369/101‑1 (MQ=255)
                                            tGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTcc                                                          >  2:3478999/1‑101 (MQ=255)
                                            tGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTcc                                                          >  2:241564/1‑101 (MQ=255)
                                             ggTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTccc                                                         <  1:3285908/101‑1 (MQ=255)
                                              gTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTcccc                                                        >  2:859224/1‑101 (MQ=255)
                                                cAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGCTG.ATCAGCGGTCGCTcccccc                                                      <  1:2690938/101‑1 (MQ=255)
                                                 aGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCcg                                                     >  1:3334749/1‑101 (MQ=255)
                                                   cGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCcgcg                                                   <  1:1591842/101‑1 (MQ=255)
                                                      gTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCgggg                                                >  2:361635/1‑101 (MQ=255)
                                                       ttACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCggggg                                               >  2:3899445/1‑101 (MQ=255)
                                                        tACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGc                                              <  1:1773351/101‑1 (MQ=255)
                                                           caaccGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTg                                           <  2:2555237/97‑1 (MQ=255)
                                                            aaaCGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTGCCCCCGCGGGGGCGTgg                                          >  1:119442/1‑101 (MQ=255)
                                                               acactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.GCTGGCGCGCGCCGCGCCG..GATCTTTGTCTGCGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                     >  2:2395214/6‑101 (MQ=2)
                                                                cactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.GCTGGCGCGCGCCGCGCCG..GATCTTTGTCTGCGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                    <  1:550018/97‑1 (MQ=2)
                                                                   tgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.GCTGGCGCGCGCCGCGCCG..GATCTTTGTCTGCGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                 <  2:4389553/100‑2 (MQ=2)
                                                                    gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                  <  2:1576505/101‑1 (MQ=255)
                                                                     tCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                 >  1:1952128/1‑101 (MQ=255)
                                                                      cGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCCCCCCCCCCGCGGGGGCGTGGATCGAACCAt                                >  1:1897572/1‑101 (MQ=255)
                                                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATc                               >  2:4056848/1‑101 (MQ=255)
                                                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCGGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATc                               >  1:4719794/1‑101 (MQ=255)
                                                                        cTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcg                              >  1:1567663/1‑101 (MQ=255)
                                                                        cTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGGGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTcg                              >  2:3964064/1‑101 (MQ=255)
                                                                         tCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgc                             >  2:4343137/1‑101 (MQ=255)
                                                                          ccccccGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGCTCGAAACATcgcc                            >  2:1314802/1‑101 (MQ=255)
                                                                           cccccGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgccg                           <  1:2376186/101‑1 (MQ=255)
                                                                           cccccGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGGCGCTCCCCCCCCGGGGGCGAGGAGCGAAACATcgccg                           >  2:4262527/1‑101 (MQ=39)
                                                                            ccccGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgccgc                          <  1:1042531/101‑1 (MQ=255)
                                                                               cgcgGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCttt                       <  2:4602113/101‑1 (MQ=255)
                                                                               cgcgGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCttt                       <  1:4162914/101‑1 (MQ=255)
                                                                                 ggggggCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTgg                     <  2:3673472/100‑1 (MQ=255)
                                                                                 cgGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTgg                     <  2:4232139/101‑1 (MQ=255)
                                                                                  gggggCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGGCGCTCCCCCCGCGGGGGCGTCGATCGAAACATCGCCTCCTTggg                    >  2:1805213/1‑100 (MQ=38)
                                                                                    gggCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCgg                  <  2:3446393/101‑1 (MQ=255)
                                                                                     ggCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGt                 >  1:3108544/1‑101 (MQ=255)
                                                                                      gCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTc                >  2:292233/1‑101 (MQ=255)
                                                                                        gTGGATCGAAACTGCCCCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGa              <  1:1354465/101‑1 (MQ=255)
                                                                                         tGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGAc             <  2:1169461/101‑1 (MQ=255)
                                                                                           gATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAg           >  1:487636/1‑101 (MQ=255)
                                                                                           gATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAg           >  1:4496414/1‑101 (MQ=255)
                                                                                            aTCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAgg          <  1:1705247/101‑1 (MQ=255)
                                                                                             tCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGc         <  1:377253/101‑1 (MQ=255)
                                                                                              cGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCt        >  1:99373/1‑101 (MQ=255)
                                                                                               gAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTg       >  1:2466069/1‑101 (MQ=255)
                                                                                               gAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTg       >  1:1802080/1‑101 (MQ=255)
                                                                                                aaaCTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTgg      >  2:1999207/1‑101 (MQ=255)
                                                                                                aaaCTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTgg      <  2:805090/101‑1 (MQ=255)
                                                                                                 aaCTGCCCCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGt     <  1:2263163/101‑1 (MQ=255)
                                                                                                 aaCTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGt     >  2:2451969/1‑101 (MQ=255)
                                                                                                  aCTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTa    <  2:2031517/101‑1 (MQ=255)
                                                                                                  aCTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTa    >  1:87160/1‑101 (MQ=255)
                                                                                                   cTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTaa   <  2:2835698/101‑1 (MQ=255)
                                                                                                    |                                                                                                     
CTTGCGACAAGCACCTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC  >  NC_014034/1395297‑1395497

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