Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,418 Δ1 bp 6.0% intergenic (+230/‑701) 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,4180T→. 6.0% 30.4 ‑1.8 ‑0.2217intergenic (+230/‑701)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (9/4):  ref base (96/105):  total (108/109)
Fisher's exact test strand distribution p-value = 1.60e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.01e-01

Only 100 of 270 total aligned reads displayed.

CCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGG  >  NC_014034/1395318‑1395518
                                                                                                    |                                                                                                     
 caccGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtg                                                                                                      <  2:284773/101‑1 (MQ=255)
  accGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCCCGGGGGCGTCGATCGAAACTGCACCCCCCGGCCCCGCcgggga                                                                                                     >  1:4376390/1‑101 (MQ=37)
   ccGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgat                                                                                                    >  2:183147/1‑101 (MQ=255)
    cGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatc                                                                                                   <  2:1892528/101‑1 (MQ=255)
     ggAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatcg                                                                                                  >  2:2186697/1‑101 (MQ=255)
     ggAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatcg                                                                                                  >  2:3254195/1‑101 (MQ=255)
      gAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCCCCGCTGGCCCTGCcggtgatcgg                                                                                                 <  1:1466262/101‑1 (MQ=255)
        gCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGTCCCTGCcggggatcggtg                                                                                               >  1:5700/1‑101 (MQ=255)
        gCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatcggtg                                                                                               >  2:1600595/1‑101 (MQ=255)
         cAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatcggtg.a                                                                                             >  1:639688/1‑101 (MQ=255)
          aaTTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCcggtgatcggtg.at                                                                                            <  1:2448610/101‑1 (MQ=255)
          aaTTTACACCGTCTGGGCAGCGGGTTACGAAACGGCAGGCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCCCCGCTGGCCCTGCcggtgatcggtg.at                                                                                            <  1:3591544/101‑1 (MQ=255)
           aTTTACGCCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCggtgatcggtg.atc                                                                                           <  2:3844140/101‑1 (MQ=255)
           aTTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCCCCGCTGGCCCTGCCggtgatcggtg.atc                                                                                           <  1:4950614/101‑1 (MQ=255)
             ttACCCCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAg                                                                                         <  2:270826/101‑1 (MQ=255)
             ttACCCCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAg                                                                                         <  1:3726104/101‑1 (MQ=255)
              tACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGc                                                                                        <  2:1296915/101‑1 (MQ=255)
              tACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGGGATCGGTG.ATCAgg                                                                                        >  2:4837760/1‑100 (MQ=255)
               acacCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCg                                                                                       >  1:182129/1‑101 (MQ=255)
                cacCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCgg                                                                                      >  1:1376203/1‑101 (MQ=255)
                 acCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGt                                                                                     <  1:4832032/101‑1 (MQ=255)
                       tGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTcc                                                                               >  2:4092344/1‑101 (MQ=255)
                        ggTCAGCGGGTTACGAAACGGCAGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTccc                                                                              <  1:2628643/101‑1 (MQ=255)
                           cAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGGGTGGATCGAAACTGCACCCGCTGGTCCTGCCGGTGATTGGTG.ATCAGCGGTCGCTcacccc                                                                           >  2:608452/1‑101 (MQ=38)
                            aGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCCCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCcg                                                                          <  1:3452953/101‑1 (MQ=255)
                            aGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCcg                                                                          >  2:2983726/1‑101 (MQ=255)
                             gCGGGTTACGAAACGGCAGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACTGCTCCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCcgc                                                                         <  1:1944610/101‑1 (MQ=255)
                                ggTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCggg                                                                      >  2:1717411/1‑101 (MQ=255)
                                 gtgaCGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCTCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCgggg                                                                     <  2:826987/98‑1 (MQ=255)
                                   tACGAAGCGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGc                                                                   <  1:2125068/101‑1 (MQ=255)
                                   tACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGc                                                                   <  1:4064809/101‑1 (MQ=255)
                                      gAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTg                                                                <  2:5021319/101‑1 (MQ=255)
                                      gAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTg                                                                >  2:2291211/1‑101 (MQ=255)
                                      gAAACGGCAGTCGCGCCCCCCGCGGGGGCGTGGATCGAACCTGCACCAGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTg                                                                <  2:330801/101‑1 (MQ=255)
                                       tcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.GCTGGCGCGCGCCGCGCCG..GATCTTTGTCTGCGCGGTCGCTCCCCCCGCGGGGGCGTGGAt                                                             >  2:3656624/9‑101 (MQ=2)
                                       aaaCGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTgg                                                               <  1:4452978/101‑1 (MQ=255)
                                        caacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.GCTGGCGCGCGCCGCGCCG..GATCTTTGTCTGCGCGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                            >  1:10638/8‑101 (MQ=2)
                                        caacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.GCTGGCGCGCGCCGCGCCG..GATCTTTGTCTGCGCGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                            <  2:1119595/94‑1 (MQ=2)
                                           ggCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                           >  1:3854686/1‑101 (MQ=255)
                                           ggCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                           <  1:1600727/101‑1 (MQ=255)
                                           cactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.GCTGGCGCGCGCCGCGCCG..GATCTTTGTCTGCGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                         >  2:521330/5‑101 (MQ=2)
                                             cAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGGGGATCGaa                                                         >  1:2983161/1‑101 (MQ=255)
                                              aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                        >  1:247630/1‑101 (MQ=255)
                                              aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.AGCAGCGGTCGCCCCCCCCGCGGGGGCGTGGATCGaaa                                                        <  2:4389555/101‑1 (MQ=255)
                                               gTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.GCTGGCGCGCGCCGCGCCG..GATCTTTGTCTGCGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActc                                                     <  1:4575325/101‑3 (MQ=2)
                                                tCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                                      <  1:1717541/101‑1 (MQ=255)
                                                 cGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                                     >  1:787531/1‑101 (MQ=255)
                                                  gCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATc                                                    >  1:1292418/1‑101 (MQ=255)
                                                   ccccccccGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcg                                                   <  1:518709/99‑1 (MQ=255)
                                                   cTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcg                                                   <  2:2031632/101‑1 (MQ=255)
                                                    tCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGGGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgc                                                  >  2:4060775/1‑101 (MQ=255)
                                                     ccccccGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgcc                                                 >  1:2082614/1‑101 (MQ=255)
                                                     ccccccGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgcc                                                 >  1:2502519/1‑101 (MQ=255)
                                                          cgcgGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCttt                                            <  1:1497987/101‑1 (MQ=255)
                                                           gcgGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTg                                           <  2:3032144/101‑1 (MQ=255)
                                                                ggCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCggg                                      >  2:3548410/1‑100 (MQ=255)
                                                                ggCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGt                                      <  2:3462169/101‑1 (MQ=255)
                                                                 gCGTGGATCGAAACTGCCCCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTc                                     <  1:2286798/101‑1 (MQ=255)
                                                                 gCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTc                                     <  2:3137470/101‑1 (MQ=255)
                                                                  cGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTTTTCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCg                                    >  2:1676202/1‑101 (MQ=255)
                                                                     ggATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACa                                 >  2:4015583/1‑101 (MQ=255)
                                                                     ggATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACa                                 >  2:4778598/1‑101 (MQ=255)
                                                                     ggATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACa                                 >  1:4188883/1‑101 (MQ=255)
                                                                      gATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAg                                >  1:90282/1‑101 (MQ=255)
                                                                       aTCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAgg                               <  1:1157456/101‑1 (MQ=255)
                                                                        tCGAAACTGCCCCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGc                              <  2:4359736/101‑1 (MQ=255)
                                                                          gAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTg                            <  1:2880644/101‑1 (MQ=255)
                                                                          gAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTg                            >  2:1218574/1‑101 (MQ=255)
                                                                           aaaCTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTgg                           >  1:2146742/1‑101 (MQ=255)
                                                                            aaCTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGt                          <  2:3733854/101‑1 (MQ=255)
                                                                             aCTGCCCCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTa                         <  1:3270020/101‑1 (MQ=255)
                                                                             aCTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTa                         >  2:1690843/1‑101 (MQ=255)
                                                                              cTGCCCCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTaa                        <  2:2455350/101‑1 (MQ=255)
                                                                              cTGCCCCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTaa                        <  1:3892476/101‑1 (MQ=255)
                                                                               tGCCCCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAc                       <  2:4243897/101‑1 (MQ=255)
                                                                                gCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTaagg                      >  1:679051/1‑99 (MQ=255)
                                                                                gCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACg                      <  1:2589447/101‑1 (MQ=255)
                                                                                gCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCCACAGGCTGGTAACg                      >  2:1554630/1‑101 (MQ=255)
                                                                                 cACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACgt                     <  2:4811535/101‑1 (MQ=255)
                                                                                  ccccGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACgtg                    <  1:4343591/100‑1 (MQ=255)
                                                                                   cccGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACgtgt                   <  1:3228443/101‑1 (MQ=255)
                                                                                    ccGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTc                  >  1:1097574/1‑101 (MQ=255)
                                                                                     cccTGTCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCg                 <  2:2532249/99‑1 (MQ=255)
                                                                                     cGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCg                 <  2:1097460/101‑1 (MQ=255)
                                                                                      gCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGc                <  2:1376077/101‑1 (MQ=255)
                                                                                       cTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTcccc               >  2:1791523/1‑98 (MQ=255)
                                                                                        tGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTc              >  2:3934627/1‑101 (MQ=255)
                                                                                         ggCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGATGGTAACGTGTCGCTcc             >  2:719928/1‑101 (MQ=255)
                                                                                           cccTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTcccc           >  2:3230512/1‑101 (MQ=255)
                                                                                           cccTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTcccc           <  2:444843/101‑1 (MQ=255)
                                                                                           cccTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTcccc           >  2:4610795/1‑101 (MQ=255)
                                                                                            ccTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTccccc          <  2:90272/101‑1 (MQ=255)
                                                                                            ccTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTccccc          >  1:2344285/1‑101 (MQ=255)
                                                                                            ccTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTccccc          >  1:3297734/1‑101 (MQ=255)
                                                                                             cTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTcccccc         <  1:3042657/101‑1 (MQ=255)
                                                                                              tGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCcg        <  1:4857107/101‑1 (MQ=255)
                                                                                                ccGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCcgcg      >  1:993649/1‑101 (MQ=255)
                                                                                                 cggtgatcggtg.atCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCgg     >  1:202781/1‑101 (MQ=255)
                                                                                                   gtgatcggtg.atcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCgggg   >  1:4410603/1‑101 (MQ=255)
                                                                                                    tgatcggtg.atcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCggggg  >  2:642205/1‑101 (MQ=255)
                                                                                                    |                                                                                                     
CCACCGGAGCAATTTACACCGTCTGGTCAGCGGGTTACGAAACGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTG.ATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGG  >  NC_014034/1395318‑1395518

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG