Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,397 Δ1 bp 3.8% 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→. 3.8% 29.6 ‑0.5 ‑1.7364intergenic (+209/‑722)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (9/5):  ref base (186/164):  total (195/169)
Fisher's exact test strand distribution p-value = 5.86e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.05e-02

Only 100 of 386 total aligned reads displayed.

CTTGCGACAAGCACCTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC  >  NC_014034/1395297‑1395497
                                                                                                    |                                                                                                     
cTTGCGACAAGCACCTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                                                                                       >  2:1989426/1‑101 (MQ=255)
  tGCGACAAGCACCTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGc                                                                                                     <  1:72892/101‑1 (MQ=255)
  tGCGACAAGCACCTGCGCACCGCCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGc                                                                                                     <  1:815779/101‑1 (MQ=255)
    cGACAAGCACCTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCAc                                                                                                   >  2:2126560/1‑101 (MQ=255)
     gACAAGCACCTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCAcc                                                                                                  >  1:3305285/1‑101 (MQ=255)
      aCAAGCACCTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCAccc                                                                                                 >  2:1410004/1‑101 (MQ=255)
       cAAGCACCTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCg                                                                                                >  2:2784044/1‑101 (MQ=255)
       cAAGCACCTTCGCACCACCGGACCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCg                                                                                                >  2:1456526/1‑101 (MQ=255)
        aaGCACCTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGc                                                                                               >  2:608966/1‑101 (MQ=255)
            aCCTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGc                                                                                           <  2:4285465/101‑1 (MQ=255)
              cTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGccc                                                                                         <  2:6825137/101‑1 (MQ=255)
              cTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGccc                                                                                         <  2:3440060/101‑1 (MQ=255)
               ttCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGcccg                                                                                        >  2:6167726/1‑100 (MQ=255)
                 cGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGc                                                                                      >  2:853935/1‑101 (MQ=255)
                  gCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGcc                                                                                     >  2:3260226/1‑101 (MQ=255)
                    accaccGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcgg                                                                                   >  2:249691/1‑101 (MQ=255)
                    accaccGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcgg                                                                                   >  2:4812595/1‑101 (MQ=255)
                    accaccGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcgg                                                                                   <  1:1757156/101‑1 (MQ=255)
                      caccGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtg                                                                                 <  2:2913033/101‑1 (MQ=255)
                      caccGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtg                                                                                 >  1:6886452/1‑101 (MQ=255)
                        ccGGAGCAATGTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgat                                                                               <  2:3805712/101‑1 (MQ=255)
                          ggAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatcg                                                                             <  1:3435630/101‑1 (MQ=255)
                           gAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatcgg                                                                            >  2:6664969/1‑101 (MQ=255)
                             gCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatcggtg                                                                          <  2:1218923/101‑1 (MQ=255)
                                aTTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCggtgatcggtg.atc                                                                      >  2:5001777/1‑101 (MQ=255)
                                  ttACCCCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAg                                                                    <  2:4629020/101‑1 (MQ=255)
                                   tACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGc                                                                   <  2:6582824/101‑1 (MQ=255)
                                    acacCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCg                                                                  >  2:2685969/1‑101 (MQ=255)
                                    acacCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCg                                                                  >  2:3529469/1‑101 (MQ=255)
                                        cGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCg                                                              <  1:4225585/101‑1 (MQ=255)
                                         gTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGc                                                             <  1:5926775/101‑1 (MQ=255)
                                         gTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGCGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGGG.ATCAGCGGTCGc                                                             >  1:2023462/1‑101 (MQ=255)
                                             ggTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTccc                                                         >  2:1932141/1‑101 (MQ=255)
                                             ggTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTccc                                                         >  2:5438756/1‑101 (MQ=255)
                                               tCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTccccc                                                       <  1:6019582/101‑1 (MQ=255)
                                                cAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTcccccc                                                      >  2:3548268/1‑101 (MQ=255)
                                                cAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTcccccc                                                      <  1:349142/101‑1 (MQ=255)
                                                 aGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCcg                                                     <  1:2605823/101‑1 (MQ=255)
                                                 aGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCcg                                                     <  1:6121916/101‑1 (MQ=255)
                                                    gggTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCgg                                                  <  1:5003861/101‑1 (MQ=255)
                                                     ggTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCggg                                                 >  1:1813761/1‑101 (MQ=255)
                                                        tACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAACCTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGc                                              <  1:6665311/101‑1 (MQ=255)
                                                        tACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGc                                              <  1:2686065/101‑1 (MQ=255)
                                                        tACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGc                                              >  2:3143512/1‑101 (MQ=255)
                                                         aCGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCg                                             >  2:2202953/1‑101 (MQ=255)
                                                          cGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCgg                                            >  1:2127619/1‑100 (MQ=255)
                                                          cGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGt                                            <  1:3688529/101‑1 (MQ=255)
                                                           gAACCGGCAGTCCCGCCCCCTGGGGGGGCGTGGATCGAACCTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTg                                           <  1:4373671/101‑1 (MQ=38)
                                                            aaaCGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTgg                                          >  1:6589829/1‑101 (MQ=255)
                                                             caacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.GCTGGCGCGCGCCGCGCCG..GATCTTTGTCTGCGCGGTCGCTCCCCCCGCGGGGGCGTGGATc                                       >  1:5722560/8‑101 (MQ=2)
                                                             aaCGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGa                                         <  1:3106168/101‑1 (MQ=255)
                                                               cGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGGGGATc                                       >  1:211275/1‑101 (MQ=255)
                                                                ggCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                      >  1:4060283/1‑101 (MQ=255)
                                                                ggCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                      >  2:51109/1‑101 (MQ=255)
                                                                 gCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                     >  2:4148223/1‑101 (MQ=255)
                                                                 gCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                     >  2:7064180/1‑101 (MQ=255)
                                                                 gCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                     <  1:1007187/101‑1 (MQ=255)
                                                                 gCAGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                     <  1:2289806/101‑1 (MQ=255)
                                                                  ctgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.GCTGGCGCGCGCCGCGCCG..GATCTTTGTCTGCGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                  >  2:5630413/3‑101 (MQ=2)
                                                                   tgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.GCTGGGGCGGGCAGCGCCG..GATCTTTGTCTGCGCGGTCGCTCCCCCCGCGGGGGCGTGGATGGAAAct                                 >  1:2319677/2‑100 (MQ=2)
                                                                   tgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.GCTGGCGCGCGCCGCGCCG..GATCTTTGTCTGCGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                 >  1:2801650/2‑100 (MQ=2)
                                                                   tgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.GCTGGCGCGCGCCGCGCCG..GATCTTTGTCTGCGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                 >  1:2978594/2‑100 (MQ=2)
                                                                   aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                   <  2:4704209/101‑1 (MQ=255)
                                                                   aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                   <  2:6463978/101‑1 (MQ=255)
                                                                    gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                  <  2:440034/101‑1 (MQ=255)
                                                                    gTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.GCTGGCGCGCGCCGCGCCG..GATCTTTGTCTGCGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActc                                >  1:4448022/1‑99 (MQ=2)
                                                                    gTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                  <  1:2446846/101‑1 (MQ=255)
                                                                     tCGCCCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                 <  1:6798556/101‑1 (MQ=255)
                                                                      cGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACac                                >  2:2986513/1‑100 (MQ=255)
                                                                      cGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                >  2:5570832/1‑101 (MQ=255)
                                                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATc                               >  1:3156898/1‑101 (MQ=255)
                                                                        cTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATGAGCGGTCGCTCCGCCCGCGGGGGCGGGGGTGGAAACATca                              >  2:2275148/1‑100 (MQ=38)
                                                                          ccccccGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgcc                            >  1:1555381/1‑101 (MQ=255)
                                                                          ccccccGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgcc                            >  2:5748558/1‑101 (MQ=255)
                                                                          ccccccGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgcc                            >  2:524262/1‑101 (MQ=255)
                                                                           cccccGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgccg                           >  2:334981/1‑101 (MQ=255)
                                                                              ccGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCtt                        <  1:6197188/101‑1 (MQ=255)
                                                                               cgcgGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCttt                       <  1:1923486/101‑1 (MQ=255)
                                                                                  gggggCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGc                    >  1:5474860/1‑101 (MQ=255)
                                                                                  gggggCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGc                    >  2:721707/1‑101 (MQ=255)
                                                                                  gggggCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGc                    >  1:3951313/1‑101 (MQ=255)
                                                                                  gggggCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGc                    >  2:1316574/1‑101 (MQ=255)
                                                                                    gggCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCgg                  >  1:5672263/1‑101 (MQ=255)
                                                                                       cGTGGATCGAAACTGCCCCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCg               <  1:119886/101‑1 (MQ=255)
                                                                                        gTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGa              >  1:2847852/1‑101 (MQ=255)
                                                                                        gTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGa              >  2:6404143/1‑101 (MQ=255)
                                                                                         tGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGAc             >  2:174215/1‑101 (MQ=255)
                                                                                          ggATCGAAACTGCGCCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTGCTCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACa            <  1:4569773/101‑1 (MQ=255)
                                                                                             tCGAAACTGCCCCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGc         <  2:6330155/101‑1 (MQ=255)
                                                                                              cGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCt        >  1:3042537/1‑101 (MQ=255)
                                                                                              cGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGCCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCt        >  2:180493/1‑101 (MQ=255)
                                                                                                aaaCTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTgg      >  2:1096357/1‑101 (MQ=255)
                                                                                                aaaCTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTgg      <  2:3960147/101‑1 (MQ=255)
                                                                                                aaaCTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTgg      >  1:5349253/1‑101 (MQ=255)
                                                                                                 aaCTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTggg     >  1:507490/1‑100 (MQ=255)
                                                                                                  aCTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTa    <  1:94790/101‑1 (MQ=255)
                                                                                                  aCTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTa    >  1:4799972/1‑101 (MQ=255)
                                                                                                  aCTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTa    <  2:3050983/101‑1 (MQ=255)
                                                                                                   cTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTaa   <  1:6450452/101‑1 (MQ=255)
                                                                                                    tGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAc  >  2:7021530/1‑101 (MQ=255)
                                                                                                    |                                                                                                     
CTTGCGACAAGCACCTTCGCACCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC  >  NC_014034/1395297‑1395497

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