Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,397,690 +G 2.2% coding (1572/1578 nt) RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,397,6901.→G 2.2% 7.3 ‑1.7 ‑0.0223coding (1572/1578 nt)RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (1/4):  ref base (121/97):  total (122/101)
Fisher's exact test strand distribution p-value = 1.79e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.28e-01

Only 100 of 261 total aligned reads displayed.

GATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGCGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGCGTGGA  >  NC_014034/1397590‑1397790
                                                                                                                                                                                                         
gATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTtgttctgttc                                                                                                      >  1:3425769/1‑101 (MQ=255)
 aTCGAAACATCCGCGGATCGGCGTGCCCCGCCCCGCCCATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTcc                                                                                                     <  2:3613971/101‑1 (MQ=255)
     ttccgatctcgGATCGGCGTGCCCCGCCCCGCCCATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTa                                                                                                 <  2:202051/92‑1 (MQ=38)
     aaaCATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTa                                                                                                 >  1:2721873/1‑101 (MQ=255)
      aaCATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTNTTGTTCTGTNNCgnnnn                                                                                                >  2:1850527/1‑97 (MQ=255)
      aaCACCCGCGGATCGGGGTGCCCCGCCCCGCCCATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTaa                                                                                                <  1:4470090/101‑1 (MQ=255)
       aCATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAAc                                                                                               <  2:539793/101‑1 (MQ=255)
        cATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACa                                                                                              >  2:4446311/1‑101 (MQ=255)
        cATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACa                                                                                              >  2:2554434/1‑101 (MQ=255)
         aTCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAg                                                                                             <  2:2524106/101‑1 (MQ=255)
            cgcgGATCGGCGTGACCCGCCCCGCCCATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTc                                                                                          <  1:1229291/101‑1 (MQ=255)
             gcgGAGCGGCGTGACCCGCCCCGCCCATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCg                                                                                         <  2:3380374/101‑1 (MQ=255)
               ggATCGGCGTGACCCGCCCCGCCCATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGCCTTGTTCTGTTCCGGCAACAGTTCGCt                                                                                       >  2:2429101/1‑101 (MQ=255)
                 aTCGGCGTGACCCGCCCCGCCCATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTcc                                                                                     >  2:3455468/1‑101 (MQ=255)
                    ggCGTGACCCGCCCCGCCCATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTccccc                                                                                  >  2:1729168/1‑101 (MQ=255)
                    ggCGTGACCCGCCCCGCCCATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTccccc                                                                                  >  2:2597590/1‑101 (MQ=255)
                    ggCGTGACCCGCCCCGCCCATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGCTCTTGTTCTGTTCCGGTAATAGGTCGCTccccc                                                                                  >  1:4339196/1‑101 (MQ=255)
                       gTGACCCGCCCCGCCCATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCcgc                                                                               >  2:4609492/1‑101 (MQ=255)
                        ggccccgccccgcCCATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCcgcg                                                                              <  2:2458510/98‑1 (MQ=255)
                         gACCCGCCCCGCCCATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCgg                                                                             >  1:1774567/1‑101 (MQ=255)
                            ccgccccgccCATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCggggg                                                                          <  1:2050508/101‑3 (MQ=255)
                              gccccgcccATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGggcg                                                                        >  2:408589/1‑101 (MQ=255)
                              gccccgcccATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAATACAGAGTTCCGGCCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGggcg                                                                        >  2:1775060/1‑101 (MQ=255)
                               ccccgcccATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGt                                                                       >  2:2221204/1‑101 (MQ=255)
                                cccgcccATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTTTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTg                                                                      >  1:2290058/1‑101 (MQ=255)
                                cccgcccATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGNGCGTg                                                                      >  1:1000238/1‑101 (MQ=255)
                                cccgcccATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTg                                                                      >  2:651626/1‑101 (MQ=255)
                                 ccgcccATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTgg                                                                     <  1:4050089/101‑1 (MQ=255)
                                   gcccATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                   <  1:1800029/101‑1 (MQ=255)
                                    cccATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                  <  1:2173160/101‑1 (MQ=255)
                                     ccATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                 <  1:3916087/101‑1 (MQ=255)
                                     ccATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                 <  2:3160805/101‑1 (MQ=255)
                                      cATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                <  2:2648792/101‑1 (MQ=255)
                                         cGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                             <  1:2630642/101‑1 (MQ=255)
                                          ggTCGCTCCCCCCGTGGGGGCGTGGATAGAAACTACATAGTTCCGCTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                                            <  2:2490752/101‑1 (MQ=37)
                                          ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                                            >  2:3664076/1‑101 (MQ=255)
                                          ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                                            >  2:4346889/1‑101 (MQ=255)
                                           gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACtt                                                           <  1:1941827/101‑1 (MQ=255)
                                             cGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGt                                                         <  1:335001/101‑1 (MQ=255)
                                                tCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTg                                                      >  2:3116734/1‑101 (MQ=255)
                                                 ccccccGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGa                                                     <  1:2841096/101‑1 (MQ=255)
                                                     ccGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGAcgcg                                                 <  1:505148/101‑1 (MQ=255)
                                                      cgcgGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACgcgc                                                >  1:1421842/1‑101 (MQ=255)
                                                       gcgGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCt                                               >  2:954915/1‑101 (MQ=255)
                                                        cgGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCtc                                              >  1:3199903/1‑101 (MQ=255)
                                                        cgGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCtc                                              <  2:1948239/101‑1 (MQ=255)
                                                        cgGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCtc                                              >  2:1248618/1‑101 (MQ=255)
                                                        cgGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCtc                                              >  1:2436638/1‑101 (MQ=255)
                                                        cgGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCtc                                              >  2:168613/1‑101 (MQ=255)
                                                         gggggCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCtcc                                             >  1:22345/1‑100 (MQ=255)
                                                         gggggCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCtca                                             <  2:4686915/101‑1 (MQ=255)
                                                         gggggCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCACCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCtca                                             >  1:4106228/1‑101 (MQ=255)
                                                          ggggCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCtcat                                            <  2:4967708/101‑1 (MQ=255)
                                                          ggggCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCtcat                                            <  2:2002437/101‑1 (MQ=255)
                                                           gggCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCtcatc                                           <  2:2381454/101‑1 (MQ=255)
                                                           gggCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCtcatc                                           >  2:4752388/1‑101 (MQ=255)
                                                             gCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAg                                         <  1:3096963/101‑1 (MQ=255)
                                                               gTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCa                                       >  2:1029228/1‑101 (MQ=255)
                                                               gTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCa                                       <  1:623022/101‑1 (MQ=255)
                                                                 ggATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACg                                     >  1:4740030/1‑101 (MQ=255)
                                                                 ggATCGAAACGACAGAGTTCCGGTCTTGTTCTGTGCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACg                                     <  2:4264287/101‑1 (MQ=255)
                                                                  gATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGt                                    <  2:2892965/101‑1 (MQ=255)
                                                                  gATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGt                                    >  1:3624922/1‑101 (MQ=255)
                                                                   aTCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTc                                   <  2:1054127/101‑1 (MQ=255)
                                                                   aTCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTc                                   <  2:729019/101‑1 (MQ=255)
                                                                      gAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAg                                <  1:587567/101‑1 (MQ=255)
                                                                       aaaCTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGa                               >  1:1169007/1‑101 (MQ=255)
                                                                       aaaCTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGa                               >  1:3565070/1‑101 (MQ=255)
                                                                       aaaCTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGa                               <  1:997417/101‑1 (MQ=255)
                                                                       aaaCTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGa                               >  1:4893761/1‑101 (MQ=255)
                                                                           tACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCa                           >  2:149178/1‑101 (MQ=255)
                                                                             cAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGt                         >  2:1398491/1‑101 (MQ=255)
                                                                              agagTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTc                        <  1:1099327/101‑1 (MQ=255)
                                                                               gagTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCg                       <  2:1613055/101‑1 (MQ=255)
                                                                                aggtCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGc                      <  1:2131118/98‑1 (MQ=255)
                                                                                 gTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCt                     >  2:4183129/1‑101 (MQ=255)
                                                                                 gTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCt                     <  2:3995275/101‑1 (MQ=255)
                                                                                   tCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTcc                   >  1:3801792/1‑101 (MQ=255)
                                                                                    ccGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTccc                  >  1:77836/1‑101 (MQ=255)
                                                                                     nggNCNTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTcccc                 <  1:1308944/100‑1 (MQ=255)
                                                                                     cGGTCTTGTTCTGTTCCGGTAACAGGTCGCTTCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTcccc                 >  1:4506962/1‑101 (MQ=255)
                                                                                      gttCTTGTTCTGTTCGGGTAACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTccccc                <  2:4397673/99‑1 (MQ=255)
                                                                                       gTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTcccccc               <  2:282637/101‑1 (MQ=255)
                                                                                        tCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCcg              <  1:2554615/101‑1 (MQ=255)
                                                                                         cTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCcgc             >  2:4803647/1‑101 (MQ=255)
                                                                                         cTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCcgc             >  1:3656744/1‑101 (MQ=255)
                                                                                          ttGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCcgcg            >  1:1570582/1‑101 (MQ=255)
                                                                                           tgttctgttcCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCgg           >  1:1583044/1‑101 (MQ=255)
                                                                                           tgttctgttcCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCgg           >  2:1312722/1‑101 (MQ=255)
                                                                                           tgttctgttcCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCgg           >  1:4156354/1‑101 (MQ=255)
                                                                                           tgttctgttcCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCgg           >  2:4882663/1‑101 (MQ=255)
                                                                                            gttctgttcCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCggg          >  1:1340627/1‑101 (MQ=255)
                                                                                             ttctgttcCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGACGCTCCCCCCGCgggg         >  1:838111/1‑101 (MQ=255)
                                                                                              tctgttcCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCggggg        >  2:1167559/1‑101 (MQ=255)
                                                                                              tctgttcCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCggggg        >  2:4993117/1‑101 (MQ=255)
                                                                                               ctgttcCGGTAACAGGTCGGTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGc       <  2:212962/101‑1 (MQ=255)
                                                                                                 gttcCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGCGt     >  1:4073987/1‑101 (MQ=255)
                                                                                                 gttcCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGCGt     <  2:3888883/101‑1 (MQ=255)
                                                                                                 gttcCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGCGt     <  1:4000871/101‑1 (MQ=255)
                                                                                                   tcCGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGCGTgg   <  1:3455733/101‑1 (MQ=255)
                                                                                                                                                                                                         
GATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTCCGGTCTTGTTCTGTTCCGGTAACAGGTCGCTCCCCCCGCGGGCGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGCGTGGA  >  NC_014034/1397590‑1397790

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