Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,303,868 +G 3.9% intergenic (+55/+102) urtE → / ← RCAP_rcc01231 urea ABC transporter ATP‑binding protein UrtE/UbiH/UbiF/VisC/COQ6 family ubiquinone biosynthesis hydroxylase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,303,8681.→G 3.9% 6.6 ‑3.6 ‑0.0102intergenic (+55/+102)urtE/RCAP_rcc01231urea ABC transporter ATP‑binding protein UrtE/UbiH/UbiF/VisC/COQ6 family ubiquinone biosynthesis hydroxylase
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (38/60):  total (42/60)
Fisher's exact test strand distribution p-value = 2.63e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.87e-01

Only 100 of 119 total aligned reads displayed.

CAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCT.C..CGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCC....GCCAAGATGGGGG..GC.CCCCCTTTTCGGGGCGGGCCGCAGGGCC  >  NC_014034/1303771‑1303961
                                                                                                       |                                                                                                    
cAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.a                                                                                                     <  2:434002/101‑1 (MQ=37)
 aGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.aa                                                                                                    <  1:4305326/101‑1 (MQ=37)
 aGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.aa                                                                                                    <  2:1034436/101‑1 (MQ=37)
  gCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.aaa                                                                                                   >  2:2575792/1‑101 (MQ=37)
  gCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.aaa                                                                                                   <  2:3607043/101‑1 (MQ=37)
  gCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.aaa                                                                                                   <  1:3667202/101‑1 (MQ=37)
    cTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAcc                                                                                                 >  2:2833649/1‑101 (MQ=37)
    cTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGCGAAATTTGGGGCAAGAT.aaaac                                                                                                 >  1:2367926/1‑99 (MQ=25)
     tGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAccc                                                                                                >  1:3005976/1‑101 (MQ=37)
     tGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAccc                                                                                                >  2:3876525/1‑101 (MQ=37)
      gTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAcccc                                                                                               >  2:447535/1‑101 (MQ=37)
      gTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAcccc                                                                                               >  1:3002695/1‑101 (MQ=37)
       tCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCa                                                                                              <  1:4521256/101‑1 (MQ=37)
       tCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGCGATATTTGGGGCAAGAT.AAACCCCa                                                                                              >  1:519772/1‑101 (MQ=37)
        cccGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCAt                                                                                             >  1:2116037/1‑101 (MQ=37)
           gcgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCt                                                                                          <  1:3385630/101‑1 (MQ=37)
           gcgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCt                                                                                          >  1:324403/1‑101 (MQ=37)
           gcgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCt                                                                                          <  1:1345128/101‑1 (MQ=37)
           gcgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCt                                                                                          <  2:3374429/101‑1 (MQ=37)
            cgACACCCTGCGCGCAGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCtt                                                                                         <  2:2820085/101‑1 (MQ=25)
            cgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCtt                                                                                         <  2:3753589/101‑1 (MQ=37)
            cgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCtt                                                                                         <  1:2985676/101‑1 (MQ=37)
            cgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCtt                                                                                         <  2:40047/101‑1 (MQ=37)
            cgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCtt                                                                                         <  1:2820522/101‑1 (MQ=37)
            cgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCtt                                                                                         <  2:4375600/101‑1 (MQ=37)
              aaaCCCTGCGCGGCGTGGTCAGCGTTTAGGCAAGGGGGCGCTGCCCCC.C..CGGCTTCGCCTCACCCCCGGGCGATTTGGGGCAAGAT.AAACCCCATCCTTTc                                                                                       >  1:2520620/1‑101 (MQ=37)
              aaaCCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTc                                                                                       <  1:4546506/101‑1 (MQ=37)
              aaaCCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTc                                                                                       >  2:2915590/1‑101 (MQ=37)
              aaaCCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTc                                                                                       <  2:133684/101‑1 (MQ=37)
              aaaCCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTc                                                                                       >  2:672659/1‑101 (MQ=37)
              aaaCCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCCCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTc                                                                                       >  2:1331851/1‑101 (MQ=37)
               aaCCCTGCGCGCCGTTGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCGCCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCa                                                                                      <  1:4612286/101‑1 (MQ=25)
               aaCCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTcc                                                                                      >  2:2933032/1‑100 (MQ=25)
                aCCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCAc                                                                                     <  2:350269/101‑1 (MQ=37)
                 cccTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCAcc                                                                                    >  2:762776/1‑101 (MQ=37)
                 cccTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCAcc                                                                                    >  1:3069401/1‑101 (MQ=37)
                   tgcagcgcCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCtt                                                                                  >  2:2187067/5‑101 (MQ=25)
                   cTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCtt                                                                                  <  1:3809830/101‑1 (MQ=37)
                   cTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCtt                                                                                  <  1:2513536/101‑1 (MQ=37)
                    tGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTg                                                                                 >  2:2185049/1‑101 (MQ=37)
                     gcgcgcCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGc                                                                                <  2:4118282/101‑1 (MQ=37)
                      cgcggCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGcc                                                                               <  2:2964949/101‑1 (MQ=25)
                       gcgcCGTGGTCAGCGTTTAGGCAAGGGGGCGCTGCCCCC.C..CGGCTTCGCCTCACCCCCGGGGTATTTTTGGCAAGAT.AAACCCCATCCTTTCTCCTTTCCg                                                                              >  1:1575536/1‑101 (MQ=25)
                       gcgcCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCg                                                                              >  1:3830234/1‑101 (MQ=37)
                        cgcCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGc                                                                             <  1:489274/101‑1 (MQ=37)
                         gcCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCa                                                                            <  2:3567447/101‑1 (MQ=37)
                         gcCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCa                                                                            <  2:2139799/101‑1 (MQ=37)
                          ccGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCaa                                                                           <  1:2874406/101‑1 (MQ=37)
                          ccGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCaa                                                                           <  2:3000834/101‑1 (MQ=37)
                           cGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCCCCCCCCGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCaaa                                                                          >  1:777301/1‑101 (MQ=37)
                            gTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGGTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAat                                                                         >  2:2545248/1‑101 (MQ=37)
                            gTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAat                                                                         <  2:2241231/101‑1 (MQ=37)
                             tGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAata                                                                        <  1:2986212/101‑1 (MQ=37)
                             tGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAata                                                                        <  1:3162787/101‑1 (MQ=37)
                             gggTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCGCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAata                                                                        <  2:4362535/100‑1 (MQ=25)
                              gggcAGGGTTTAGGGCAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAatat                                                                       <  1:1531611/98‑1 (MQ=21)
                              ggTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAatat                                                                       <  1:1129550/101‑1 (MQ=37)
                              ggTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAatat                                                                       <  1:4116336/101‑1 (MQ=37)
                               gTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATc                                                                      <  2:568293/101‑1 (MQ=37)
                               gTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATc                                                                      <  1:2185199/101‑1 (MQ=37)
                                tCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATcc                                                                     >  1:4187963/1‑101 (MQ=37)
                                 cAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATccc                                                                    <  2:3555021/101‑1 (MQ=37)
                                 cAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATccc                                                                    >  1:3110031/1‑101 (MQ=37)
                                 cAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATccc                                                                    <  1:664412/101‑1 (MQ=37)
                                   gCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCgg                                                                  <  1:442255/101‑1 (MQ=37)
                                   gCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCgg                                                                  <  1:186243/101‑1 (MQ=37)
                                     gTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCgggg                                                                <  1:447572/101‑1 (MQ=37)
                                     ctgacggacggggggCTCTGC.CCC.C..CGGCTTCG.CTC.CCCCCGGGATATTTGCGGCAAGATGAAGGGGCATGATTTCATCTTGCCCGAAATATCCCaggggg                                                              >  2:2585515/11‑101 (MQ=2)
                                     ctgacggacggggggCTCTGC.CCC.C..CGGCTTCG.CTC.CCCCCGGGATATTTGCGGCAAGATGAAGGGGCATGATTTCATCTTGCCCGAAATATCCCaggggg                                                              >  1:3167541/11‑101 (MQ=2)
                                      tttAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCggggg                                                               >  1:3705785/1‑101 (MQ=37)
                                       ttAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCgggggg                                                              <  2:4448374/101‑1 (MQ=37)
                                        tAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCCAATATCCCggggggc                                                             >  2:567525/1‑100 (MQ=25)
                                         aGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCgggggggg                                                            >  1:541675/1‑99 (MQ=25)
                                         aGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTg                                                            <  1:457214/101‑1 (MQ=37)
                                         aGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTg                                                            <  1:2882259/101‑1 (MQ=37)
                                           gCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGaa                                                          <  1:1192725/101‑1 (MQ=37)
                                             aaGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAAtt                                                        <  2:4271463/101‑1 (MQ=37)
                                             aaGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAAtt                                                        >  2:3190409/1‑101 (MQ=25)
                                             aaGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCGATCACCTTGCCGCAAATATCCCGGGGGGTGAAtt                                                        >  1:4005449/1‑101 (MQ=25)
                                             aaGGGGGCGCTG..CCC.CCTCGGCTTCCCCTCCCCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAAtt                                                        <  1:3413847/101‑1 (MQ=25)
                                               gggggCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTgg                                                      >  1:3012954/1‑101 (MQ=37)
                                                ggggCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGgc                                                     <  2:3256091/101‑1 (MQ=37)
                                                 gggCTCTG..CCCCC..CGGCTTCG.CTC.CCCCCGGGATATTTGCGGCAAGATGAAGGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGgccaaa                                                 >  1:2983703/1‑98 (MQ=2)
                                                 gggCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGgcc                                                    >  2:4053354/1‑101 (MQ=37)
                                                  ggCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGgcc....g                                               <  2:3069157/101‑1 (MQ=37)
                                                                                            tGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGGGCAGCGCCCCCTTTTCGgggcggg            >  2:169142/1‑101 (MQ=14)
                                                                                             ggggCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGGGCAGCGCCCCCTTTTCGgggcgggc           <  2:1302747/101‑1 (MQ=14)
                                                                                              gggCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGGGCAGCGCCCCCTTTTCGGGGCgggcc          >  2:4068038/1‑101 (MQ=14)
                                                                                                 cAAGAT.AACCCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGGGCAGCGCCCCCTTTTCGGGGCgggccgca       <  1:2988281/101‑1 (MQ=9)
                                                                                                  aaGAT.AACCCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGGGCAGCGCCCCCTTTTCGGGGCgggccgcag      <  2:324371/101‑1 (MQ=9)
                                                                                                  aaGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGGGCAGCGCCCCCTTTTCGGGGCgggccgcag      >  1:3519800/1‑101 (MQ=14)
                                                                                                  aaGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGGGCAGCGCCCCCTTTTCGGGGCgggccgcag      <  2:3513518/101‑1 (MQ=14)
                                                                                                  aaGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGGGAATTGGCCGAAGGCCAAGAGGGGGGCAGCGCCCCCTTTTCGGGGCgggccgcag      >  1:3361931/1‑101 (MQ=9)
                                                                                                   aGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGGGCAGCGCCCCCTTTTCGGGGCgggccgcagg     >  1:4163937/1‑101 (MQ=14)
                                                                                                   aGAT.AAACCCCATCCTTTCACCTTGCCGCAAATACCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGGGCAGCGCCCCCTTTTCGGGGCgggccgcagg     <  2:2228843/101‑1 (MQ=9)
                                                                                                    gAT.AAACCCCATCCTTTCCCCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGGGCAGCGCCCCCTTTTCGGGGCgggccgcaggg    <  1:696988/101‑1 (MQ=9)
                                                                                                     aT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGGGAATTGGCCGAAGGCCAGGAGGGGGGCAGCGCCCCCTTTTAGGGGCgggccgcagggc   >  1:1387025/1‑101 (MQ=9)
                                                                                                      t.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGGGCAGCGCCCCCTTTTCGGGGCgggccgcagggcc  <  2:3334533/101‑1 (MQ=14)
                                                                                                      t.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGGGCAGCGCCCCCTTTTCGGGGCgggccgcagggcc  >  1:2445950/1‑101 (MQ=14)
                                                                                                      t.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGGGCAGCGCCCCCTTTTCGGGGCgggccgcagggcc  >  1:3835737/1‑101 (MQ=14)
                                                                                                       |                                                                                                    
CAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCT.C..CGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCC....GCCAAGATGGGGG..GC.CCCCCTTTTCGGGGCGGGCCGCAGGGCC  >  NC_014034/1303771‑1303961

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