Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,301,060 +C 91.6% 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.→C91.6% 39.2 ‑1.7 ‑1.7167intergenic (+17/‑92)urtA/urtCurea ABC transporter permease UrtA/urea ABC transporter permease UrtC
Reads supporting (aligned to +/- strand):  new base (74/79):  ref base (4/10):  total (78/89)
Fisher's exact test strand distribution p-value = 1.75e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.94e-02

Only 100 of 184 total aligned reads displayed.

TGGATGATCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGGCG.CC.TTTCATCTTGCCCAAAATATCCCGCGGGGGTCCCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGA  >  NC_014034/1300960‑1301152
                                                                                                       |                                                                                             
tGGATGATCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGccg                                                                                                >  2:118492/1‑101 (MQ=255)
tGGATGATCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGccg                                                                                                >  2:4739297/1‑101 (MQ=255)
tGGATGATCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCACCGGAGATTGAGCCGATGC.CC.CCCGccg                                                                                                <  1:750796/101‑1 (MQ=255)
 ggATGATCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGccg.c                                                                                              <  1:1123859/101‑1 (MQ=255)
   atgatCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCG.ccc                                                                                            <  2:4200542/101‑2 (MQ=255)
    tgatCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCG.ccct                                                                                           >  1:4406181/1‑99 (MQ=255)
     gatCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGCGATTGAGCCTATGC.CC.CCCGCCG.CC.ctt                                                                                         >  1:2046735/1‑101 (MQ=255)
     gatCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCG.CC.ctt                                                                                         <  1:4350582/101‑1 (MQ=255)
     gatCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCG.CC.ctt                                                                                         <  2:5619421/101‑1 (MQ=255)
     gatCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCG.CC.ctt                                                                                         >  1:4132835/1‑101 (MQ=255)
      atCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATTATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCG.CC.cttt                                                                                        <  1:779626/101‑2 (MQ=255)
      atCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCG.CC.cttt                                                                                        <  2:4134259/101‑2 (MQ=255)
       tCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCG.CC.ctttc                                                                                       >  1:22312/1‑99 (MQ=255)
       tCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCG.CC.ctttc                                                                                       <  1:5523736/101‑3 (MQ=255)
       tCCTCTTCATCATCATCTTCATCCAGTTCAGCCCGCGCGGCATTATCGCGCTGCGCGGCCGCGCCGCCGGGGATTGAGCCGATGA.CC.CCCGCCG.CC.ctttc                                                                                       >  1:343873/1‑99 (MQ=255)
        ccTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCa                                                                                       <  2:2212568/101‑1 (MQ=255)
         ctctTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCAt                                                                                      >  1:2530497/1‑101 (MQ=255)
           ctTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCt                                                                                    <  1:1390094/101‑1 (MQ=255)
            ttCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCtt                                                                                   <  2:4004894/101‑1 (MQ=255)
            ttCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCtt                                                                                   <  1:3055913/101‑1 (MQ=255)
             tcatcatcaTCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTg                                                                                  <  1:722403/101‑1 (MQ=255)
              catcatcatCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGc                                                                                 >  2:4751159/1‑101 (MQ=255)
               atcatcatcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCGGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGcc                                                                                <  2:5732908/101‑1 (MQ=255)
                tcatcatcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGccc                                                                               >  1:1077185/1‑101 (MQ=255)
                 catcatcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCa                                                                              <  1:3147875/101‑1 (MQ=255)
                  atcatcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCaa                                                                             >  2:2893881/1‑101 (MQ=255)
                   tcatcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCaaa                                                                            <  2:467350/101‑1 (MQ=255)
                   tcatcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCaaa                                                                            >  2:2069551/1‑101 (MQ=255)
                    catcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCaaaa                                                                           <  1:2320717/101‑1 (MQ=255)
                    catcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCaaaa                                                                           <  1:1037226/101‑1 (MQ=255)
                     atcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAat                                                                          <  1:4306864/101‑1 (MQ=255)
                     atcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAat                                                                          <  1:3801675/101‑1 (MQ=255)
                      tcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAata                                                                         <  2:3120349/101‑1 (MQ=255)
                      tcTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAata                                                                         <  1:1378937/101‑1 (MQ=255)
                       cTTCATCCAGTTCCGGCCGCGCCGCATCCTGCCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAatat                                                                        <  2:4970730/101‑1 (MQ=255)
                       cTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAatat                                                                        >  1:3044036/1‑101 (MQ=255)
                       cTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAatat                                                                        >  1:1531248/1‑101 (MQ=255)
                        ttCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATc                                                                       <  2:3384022/101‑1 (MQ=255)
                        ttCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATc                                                                       <  1:948971/101‑1 (MQ=255)
                        ttCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATc                                                                       >  1:1699857/1‑101 (MQ=255)
                         tCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATcc                                                                      <  1:3268519/101‑1 (MQ=255)
                          cATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATccc                                                                     >  1:5685707/1‑101 (MQ=255)
                          cATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATccc                                                                     <  1:2925621/101‑1 (MQ=255)
                           aTCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCcg                                                                    <  2:109736/101‑1 (MQ=255)
                           aTCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCcg                                                                    >  2:119675/1‑101 (MQ=255)
                           aTCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCcg                                                                    >  2:781666/1‑101 (MQ=255)
                           aTCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCcg                                                                    <  1:1024553/101‑1 (MQ=255)
                           aTCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCcg                                                                    <  1:5256638/101‑1 (MQ=255)
                            tCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCcgc                                                                   >  1:5296586/1‑101 (MQ=255)
                             ccAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCcgcg                                                                  <  2:5172664/101‑1 (MQ=255)
                              cAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCgg                                                                 <  2:3228705/101‑1 (MQ=255)
                              cAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCgg                                                                 >  2:2527529/1‑101 (MQ=255)
                              cAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCgg                                                                 <  1:1478356/101‑1 (MQ=255)
                               aGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCggg                                                                <  2:3343098/101‑1 (MQ=255)
                               aGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCggg                                                                >  2:3817800/1‑101 (MQ=255)
                               aGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCggg                                                                >  2:4272101/1‑101 (MQ=255)
                                 ttCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCggggg                                                              >  2:631965/1‑101 (MQ=255)
                                 ttCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCggggg                                                              >  1:3597097/1‑101 (MQ=255)
                                  tCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGt                                                             <  2:22301/101‑1 (MQ=255)
                                   ccgcccgcGCGGCATCAGCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGTTGCGCCTCCCGGC..CC.TTTCATCTTGCCCAAAATACACCGCGGGCGtt                                                            >  2:3578294/1‑100 (MQ=255)
                                    cgcccgcGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGTcc                                                           <  2:2394750/101‑1 (MQ=255)
                                    cgcccgcGCGCCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGTcc                                                           <  1:5767754/101‑1 (MQ=38)
                                      cccgcGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGTccgg                                                         <  2:71250/101‑3 (MQ=255)
                                       ccgcGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGTCCGgg                                                        >  2:2490415/1‑101 (MQ=37)
                                       ccgcGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGGCCGgg                                                        >  1:4611879/1‑101 (MQ=37)
                                        cgcgcgGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGTCCGggg                                                       <  1:5603421/101‑1 (MQ=37)
                                           gcgGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGgcg                                                   <  1:2416847/101‑1 (MQ=25)
                                           gcgGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGgcg                                                   >  1:5098041/1‑101 (MQ=25)
                                            cgGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGgcgc                                                  >  1:5788236/1‑101 (MQ=25)
                                             ggCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGcgcg                                                 <  1:4356054/101‑1 (MQ=25)
                                             ggCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGcgcg                                                 <  1:4272361/101‑1 (MQ=25)
                                              gCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGGCCGggggggcga                                                 >  2:5613970/1‑97 (MQ=25)
                                                atcatcGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAg                                              <  2:2322907/101‑1 (MQ=25)
                                                atcatcGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAg                                              <  2:3626893/101‑1 (MQ=25)
                                                 tcatcGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGc                                             >  1:4161116/1‑101 (MQ=25)
                                                 tcatcGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGc                                             >  2:863247/1‑101 (MQ=25)
                                                 tcatcGCGCGGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGc                                             <  1:5183575/101‑1 (MQ=21)
                                                 acatcGCGCTGCCCGGCCGCGCCGCGCGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGc                                             <  1:1162159/100‑1 (MQ=18)
                                                    tcGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGcccc                                          <  2:316344/101‑1 (MQ=21)
                                                     cgcgCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGccccc                                         <  1:5799659/101‑1 (MQ=21)
                                                     cgcgCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGccccc                                         <  2:1078547/101‑1 (MQ=21)
                                                          ggcgcGGCCGCGCCGCCGGAGTTGGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGCCCCCGgcgc                                    <  2:3718690/100‑1 (MQ=18)
                                                             gcgGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGCCCCCGGcgccgc                                 >  2:5105802/1‑101 (MQ=18)
                                                                gccgcgccgcCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAg                              >  2:3262192/1‑101 (MQ=18)
                                                                  cgcgccgcCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAGAt                            >  1:5331412/1‑101 (MQ=18)
                                                                    cgccgcCGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATcc                          >  1:3674561/1‑101 (MQ=18)
                                                                     gccgccGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATccc                         <  2:2772138/101‑1 (MQ=18)
                                                                     gccgccGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATccc                         >  2:4468154/1‑101 (MQ=18)
                                                                     gccgccGGAGATTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATccc                         >  1:5094515/1‑101 (MQ=18)
                                                                             agaTTGAGCCGATGC.CC.CCCGCCGCCC.TTTCATCTTCCCAAAAAAATCCCGCGGGGG.CCCGGGGGCGCGAAGCCCCCGGCCCCCCCAGACCCCGGGGCccc                 >  1:478918/1‑101 (MQ=12)
                                                                               aTTAAGCCGATGC.CC.CCCGTCG.CCATTTCATCTTTCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATCCCGGAGACCCCg               >  2:181324/1‑101 (MQ=18)
                                                                                     ccGATGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGaca         >  2:354226/1‑101 (MQ=18)
                                                                                      cgaccc.cc.cccGCCC.CCTTTTCATCTTCCCAAAAATACCCCG.GGGGGGCCGGGGGGCGGGAACCCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGacac        <  1:4739578/96‑1 (MQ=255)
                                                                                       gctcc.cc.cccGCCGCCC.TTTCATCTTGCCAAAAATCCCCCGCGGGGGT.CGGGGGCCGCAAAGCCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGAcaca       <  1:2634691/97‑1 (MQ=18)
                                                                                        cccc.cc.cccCCCC.CC.TTTCCTTCTCCCCAAAAATCCCCGGGGGCGGCGGGGGGGCGCGACCCCCCGGGCGCCGCCAGATCCCGGAGACCCCGATGACACAg      <  1:395358/98‑1 (MQ=255)
                                                                                        aTGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAg      >  2:5225650/1‑101 (MQ=18)
                                                                                        aTGC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAg      >  2:4102543/1‑101 (MQ=18)
                                                                                          gC.CC.CCCGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGcc    >  2:2445402/1‑101 (MQ=18)
                                                                                          gC.CC.CCCGCCGCCC.TTTCATCTCGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGcc    >  1:5690994/1‑101 (MQ=18)
                                                                                             cc.cccGCCGCCC.TTTCATCTTGCCCAAAATATCCCGCGGGGGT.CCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGa  >  1:1209208/1‑101 (MQ=18)
                                                                                                       |                                                                                             
TGGATGATCCTCTTCATCATCATCTTCATCCAGTTCCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATGC.CC.CCCGGCG.CC.TTTCATCTTGCCCAAAATATCCCGCGGGGGTCCCGGGGGCGCGAAGCCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGA  >  NC_014034/1300960‑1301152

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