Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,301,060 +C 89.0% intergenic (+17/‑92) urtA → / → urtC urea ABC transporter permease UrtA/urea ABC transporter permease UrtC

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,301,0601.→C89.0% 41.0 0.0 ‑1.7127intergenic (+17/‑92)urtA/urtCurea ABC transporter permease UrtA/urea ABC transporter permease UrtC
Reads supporting (aligned to +/- strand):  new base (47/66):  ref base (6/8):  total (53/74)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.92e-02

Only 100 of 136 total aligned reads displayed.

TGGATGATCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGGCG.CCTTTCATCTTGCCCAAAATATCCCGC.GGGGGTCCCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGAT  >  NC_014034/1300960‑1301153
                                                                                                     |                                                                                              
tGGATGATCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGccg                                                                                                 >  1:3733849/1‑101 (MQ=255)
 ggATGATCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGccg.c                                                                                               >  2:1371421/1‑101 (MQ=255)
    tgatCCTCGTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCG.ccct                                                                                            >  2:1180402/1‑99 (MQ=255)
     gatCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCG.CCctt                                                                                           >  1:2657144/1‑101 (MQ=255)
     gatCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCG.CCctt                                                                                           <  1:918549/101‑1 (MQ=255)
     gatCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCG.CCctt                                                                                           <  1:1598520/101‑1 (MQ=255)
     gatCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGATGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCG.CCctt                                                                                           <  1:1825182/101‑1 (MQ=255)
      atCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCG.CCcttt                                                                                          <  1:2588682/101‑2 (MQ=255)
      atCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCG.CCcttt                                                                                          <  1:3833341/101‑2 (MQ=255)
      atCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCG.CCcttt                                                                                          <  1:3949867/101‑2 (MQ=255)
      atCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCG.CCcttt                                                                                          >  2:2259680/1‑100 (MQ=255)
       tCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCG.CCctttc                                                                                         >  1:2897893/1‑99 (MQ=255)
        ccTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCa                                                                                         >  1:3122072/1‑101 (MQ=255)
        ccTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCa                                                                                         <  1:3291526/101‑1 (MQ=255)
        ccTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCa                                                                                         >  1:1028340/1‑101 (MQ=255)
         ctctTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCAt                                                                                        <  1:283190/101‑1 (MQ=255)
         ctctTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCAt                                                                                        <  2:2944401/101‑1 (MQ=255)
         ctctTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCAt                                                                                        <  1:850752/101‑1 (MQ=255)
          tctTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATc                                                                                       >  1:3676435/1‑101 (MQ=255)
            ttCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCtt                                                                                     >  2:3163762/1‑101 (MQ=255)
            ttCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCtt                                                                                     <  2:313800/101‑1 (MQ=255)
            ttCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCtt                                                                                     >  1:1555928/1‑101 (MQ=255)
            ttCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCtt                                                                                     <  2:1867885/101‑1 (MQ=255)
             tcatcatcaTCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTg                                                                                    <  1:344185/101‑1 (MQ=255)
              catcatcatCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGc                                                                                   >  2:3184726/1‑101 (MQ=255)
              catcatcatCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGc                                                                                   >  2:3845183/1‑101 (MQ=255)
               atcatcatcTTCATCCAGTACCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGcc                                                                                  <  1:2484008/101‑1 (MQ=255)
                 catcatcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCa                                                                                <  1:533694/101‑1 (MQ=255)
                 catcatcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCa                                                                                <  1:2390781/101‑1 (MQ=255)
                 catcatcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCa                                                                                <  1:1611642/101‑1 (MQ=255)
                 catcatcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCa                                                                                >  2:678799/1‑101 (MQ=255)
                 catcatcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCa                                                                                <  1:1797334/101‑1 (MQ=255)
                  atcatcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCaa                                                                               >  2:895554/1‑101 (MQ=255)
                   tcatcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCaaa                                                                              >  2:1903052/1‑101 (MQ=255)
                    catcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGcagatcg                                                                             >  2:782439/1‑95 (MQ=255)
                    catcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCaaaa                                                                             >  1:1657408/1‑101 (MQ=255)
                     atcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAat                                                                            >  1:1108722/1‑101 (MQ=255)
                     atcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAat                                                                            >  2:1995151/1‑101 (MQ=255)
                      tcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAata                                                                           <  1:1577447/101‑1 (MQ=255)
                      tcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAata                                                                           <  1:722201/101‑1 (MQ=255)
                      ccTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAata                                                                           <  2:2416526/100‑1 (MQ=255)
                       cTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAatat                                                                          <  1:1029147/101‑1 (MQ=255)
                        ttCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATc                                                                         <  1:2909071/101‑1 (MQ=255)
                         tCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATcc                                                                        <  2:2649082/101‑1 (MQ=255)
                          cATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATccc                                                                       >  2:1960567/1‑101 (MQ=255)
                          cATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATccc                                                                       <  1:2445193/101‑1 (MQ=255)
                          cATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATccc                                                                       <  2:4151602/101‑1 (MQ=255)
                          cATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATccc                                                                       >  1:2977149/1‑101 (MQ=255)
                          cATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATccc                                                                       <  1:2033852/101‑1 (MQ=255)
                          cATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATccc                                                                       <  2:2604095/101‑1 (MQ=255)
                            ccgaggtCCGTCCGCGCGGCATCGTCGCGCTGCGCGGCGGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCcgc                                                                     <  2:1544687/95‑1 (MQ=37)
                             ccAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCcgc.g                                                                   >  2:1898265/1‑101 (MQ=255)
                             ccAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCcgc.g                                                                   >  1:2797108/1‑101 (MQ=255)
                             ccAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCcgc.g                                                                   <  2:2517357/101‑1 (MQ=255)
                              cAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.gg                                                                  >  1:621705/1‑101 (MQ=255)
                                gTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.gggg                                                                >  1:807560/1‑101 (MQ=255)
                                gTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.gggg                                                                >  2:2669399/1‑101 (MQ=255)
                                 ttCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.ggggg                                                               >  1:1256137/1‑101 (MQ=255)
                                   ccgcccgcGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCACAATATCCCGC.GGGGGTc                                                             >  1:2505091/1‑101 (MQ=38)
                                   ccgcccgcGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGTc                                                             <  2:3121841/101‑1 (MQ=255)
                                    cgcccgcGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGTcc                                                            <  1:1180496/101‑1 (MQ=255)
                                    cgcccgcGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGTcc                                                            <  1:1216131/101‑1 (MQ=255)
                                    cgcccgcGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGTcc                                                            <  1:3164001/101‑1 (MQ=255)
                                      cccgcGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGTccgg                                                          <  2:1820261/101‑3 (MQ=255)
                                      cccgcGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGTccgg                                                          >  2:1966121/1‑99 (MQ=255)
                                      cccgcGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGTccgg                                                          <  2:3403200/101‑3 (MQ=255)
                                       ccgcGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGTCCGgg                                                         <  1:2171823/101‑1 (MQ=37)
                                        cgcgcgGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGTCCGggg                                                        <  2:2944062/101‑1 (MQ=37)
                                        cgcgcgGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGTCCGggg                                                        <  1:3215391/101‑1 (MQ=37)
                                          cgcgGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGTCCGggggc                                                      <  1:3224575/101‑2 (MQ=37)
                                           ccgGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGgcg                                                    <  2:2069296/100‑1 (MQ=25)
                                            cgGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGgcgc                                                   <  2:4145544/101‑1 (MQ=25)
                                            cgGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGgcgc                                                   <  2:2662265/101‑1 (MQ=25)
                                               catcatCGCGCTGCGCGGCCGTTCCTCCGGATATGGAGCCTATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGGCGCGaa                                                <  2:1009576/101‑1 (MQ=18)
                                                atcatcGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGGCGCGAAg                                               <  1:3951820/101‑1 (MQ=25)
                                                atcatcGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGGCGCGAAg                                               <  2:1269987/101‑1 (MQ=25)
                                                  catcGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGGCGCGAAGcc                                             <  2:2473398/101‑1 (MQ=25)
                                                    tcGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGGCGCGAAGcccc                                           <  1:51144/101‑1 (MQ=21)
                                                       cgcTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGGCGCGAAGCCCCCgg                                        >  1:2143630/1‑101 (MQ=21)
                                                         cTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGGCGCGAAGCCCCCGgcg                                      >  2:3291572/1‑101 (MQ=21)
                                                           gcgcGGCCGCGCCGCGGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGGCGCGAAGCCCCCGGcgcc                                    <  2:4159770/101‑1 (MQ=18)
                                                           gcgcGGCCGCGCCGCCGGAGATTGAGCCGATCCCCCCCGCCGCCCTTTCATCTTGCCCAAAATACCCCGG.GGGGGT.CCGGGGGCGCGAAGCCCCCGGcgcc                                    <  1:1898385/101‑1 (MQ=18)
                                                             gcgGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGGCGCGAAGCCCCCGGcgccgc                                  >  2:236584/1‑101 (MQ=18)
                                                             gcgGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGGCGCGAAGCCCCCGGcgccgc                                  >  2:1807558/1‑101 (MQ=18)
                                                               ggCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCa                                >  2:3818404/1‑101 (MQ=18)
                                                               ggCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCa                                >  1:456038/1‑101 (MQ=18)
                                                                    cgccgcCGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATcc                           >  2:2618648/1‑101 (MQ=18)
                                                                     gccgccGGAGATTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTCCCCAAAATACGCCGC.GGCGGT.CCGGGGGCGCGAAGCCCCCCGCGCCGCCAGACccc                          >  1:3555097/1‑101 (MQ=18)
                                                                              gaTTGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATCCCGGAGAcccc                 >  2:3538515/1‑101 (MQ=18)
                                                                              gaTTGAGCCGACCCCCCCCCCCC.CCTTTTCCCCTGGCCAAAA.ATCCCCCGGGGGGGGCGGGGGGCGCGACGCCCCCGGCGCCGCCAGATCCCGGAGAcccc                 <  1:1371521/101‑1 (MQ=255)
                                                                               attgccccGCTGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGTGT.CCGGGGGCGTGAAGCCCCCGGCGCCGCCAGATCCCGGAGACCCCg                <  2:2106927/95‑1 (MQ=17)
                                                                                ttGAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGGCGCGAAGCCCCCGGTGCGGCCAGATCACGGAGACCCCgg               >  2:2389124/1‑100 (MQ=18)
                                                                                  gAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATCCCGGAGACCCCGATg             <  1:1903172/101‑1 (MQ=18)
                                                                                  gAGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATCCCGGAGACCCCGATg             >  2:1197086/1‑101 (MQ=18)
                                                                                  gAGCCGAGGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATCCCGGAGACCCCGATg             <  2:2796889/101‑1 (MQ=18)
                                                                                   aGCCGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGa            >  1:3948614/1‑101 (MQ=18)
                                                                                      cGATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGacac         >  2:3297222/1‑101 (MQ=18)
                                                                                       gATGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGAcaca        >  1:3058778/1‑101 (MQ=18)
                                                                                        aTGCCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAg       <  2:3486536/101‑1 (MQ=18)
                                                                                             ccccGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGAt  <  1:2542402/101‑1 (MQ=18)
                                                                                                     |                                                                                              
TGGATGATCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCGGCG.CCTTTCATCTTGCCCAAAATATCCCGC.GGGGGTCCCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGAT  >  NC_014034/1300960‑1301153

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