Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,303,872 C→G 3.8% intergenic (+59/+98) 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,8720C→G 3.8% 2.2 ‑3.8 ‑0.0105intergenic (+59/+98)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/63):  total (42/63)
Fisher's exact test strand distribution p-value = 2.34e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.41e-01

Only 100 of 120 total aligned reads displayed.

CTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCT.C..CGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCC....GCCAAGATGGGG..G..GC.CCCCCTTTTCGGGGCGGGCCGCAGGGCCGCAC  >  NC_014034/1303775‑1303965
                                                                                                       |                                                                                                      
 tGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAccc                                                                                                      >  1:4448669/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:1345128/101‑1 (MQ=37)
       gcgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCt                                                                                                >  1:324403/1‑101 (MQ=37)
       gcgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCt                                                                                                <  1:3385630/101‑1 (MQ=37)
        cgACACCCTGCGCGCAGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCtt                                                                                               <  2:2820085/101‑1 (MQ=25)
        cgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCtt                                                                                               <  2:3151845/101‑1 (MQ=37)
        cgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCtt                                                                                               >  2:3809545/1‑101 (MQ=37)
        cgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCtt                                                                                               <  1:2985676/101‑1 (MQ=37)
        cgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCtt                                                                                               <  1:2820522/101‑1 (MQ=37)
        cgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCtt                                                                                               <  2:40047/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                                                                                             >  2:2915590/1‑101 (MQ=37)
          aaaCCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTc                                                                                             <  2:133684/101‑1 (MQ=37)
          aaaCCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTc                                                                                             <  1:4546506/101‑1 (MQ=37)
          aaaCCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTc                                                                                             >  2:672659/1‑101 (MQ=37)
           aaCCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTcc                                                                                            >  2:2933032/1‑100 (MQ=25)
            aCCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCAc                                                                                           <  2:350269/101‑1 (MQ=37)
             cccTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCAcc                                                                                          >  1:3069401/1‑101 (MQ=37)
             cccTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCAcc                                                                                          >  2:762776/1‑101 (MQ=37)
              ccTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCt                                                                                         >  2:1356017/1‑101 (MQ=37)
               tgcagcgcCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCtt                                                                                        >  2:2187067/5‑101 (MQ=25)
               cTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCtt                                                                                        <  1:2513536/101‑1 (MQ=37)
               cTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCtt                                                                                        <  1:3809830/101‑1 (MQ=37)
                tGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTg                                                                                       >  2:2185049/1‑101 (MQ=37)
                  cgcggCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGcc                                                                                     <  2:2964949/101‑1 (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)
                      ccGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCaa                                                                                 <  2:3000834/101‑1 (MQ=37)
                      ccGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCaa                                                                                 <  1:2874406/101‑1 (MQ=37)
                      ccGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCaa                                                                                 >  1:1585085/1‑101 (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                                                                              <  2:4252110/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:3435157/1‑101 (MQ=25)
                           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)
                              aGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCg                                                                         >  1:3606527/1‑101 (MQ=25)
                              aGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCg                                                                         <  1:3876810/101‑1 (MQ=37)
                               gCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCgg                                                                        <  1:186243/101‑1 (MQ=37)
                               gCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCgg                                                                        <  2:2272683/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)
                                  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)
                                      ggCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGAAAAAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGa                                                                 <  1:4653332/101‑1 (MQ=37)
                                       gCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGaa                                                                <  1:1192725/101‑1 (MQ=37)
                                         aaGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAAtt                                                              <  1:762851/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)
                                          aGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTg                                                             <  1:3012741/101‑1 (MQ=37)
                                           gggggCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTgg                                                            >  1:3012954/1‑101 (MQ=37)
                                            ggggCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGgc                                                           <  2:3256091/101‑1 (MQ=37)
                                             gggCTCTG..CCCCC..CGGCTTCG.CTC.CCCCCGGGATATTTGCGGCAAGATGAAGGGGCATGATTTCATCTTGCCCGAAATATCCCAGGGGG.GAATTGgccaaa                                                       >  1:3204317/1‑98 (MQ=2)
                                             gggCTCTG..CCCCC..CGGCTTCG.CTC.CCCCCGGGATATTTGCGGCAAGATGAAGGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGgccaaa                                                       >  1:2983703/1‑98 (MQ=2)
                                             gggCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGgcc                                                          >  2:4053354/1‑101 (MQ=37)
                                                                                        tGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGgggcggg                >  2:169142/1‑101 (MQ=14)
                                                                                         ggggCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGgggcgggc               <  2:1302747/101‑1 (MQ=14)
                                                                                          gggCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCgggcc              >  2:4068038/1‑101 (MQ=14)
                                                                                             cAAGAT.AACCCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCgggccgca           <  1:2988281/101‑1 (MQ=9)
                                                                                              aaGAT.AACCCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCgggccgcag          <  2:324371/101‑1 (MQ=9)
                                                                                              aaGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCgggccgcag          <  2:3513518/101‑1 (MQ=14)
                                                                                              aaGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCgggccgcag          >  1:3519800/1‑101 (MQ=14)
                                                                                              aaGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGGGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCgggccgcag          >  1:3361931/1‑101 (MQ=9)
                                                                                               aGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCgggccgcagg         >  1:4163937/1‑101 (MQ=14)
                                                                                               aGAT.AAACCCCATCCTTTCACCTTGCCGCAAATACCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCgggccgcagg         <  2:2228843/101‑1 (MQ=9)
                                                                                                gAT.AAACCCCATCCTTTCCCCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCgggccgcaggg        <  1:696988/101‑1 (MQ=9)
                                                                                                gAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCCGCGCCCCCTTTTCGGGGCgggccgccagg        >  2:2901542/1‑97 (MQ=9)
                                                                                                 aT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGGGAATTGGCCGAAGGCCAGGAGGGGG..GCAGCGCCCCCTTTTAGGGGCgggccgcagggc       >  1:1387025/1‑101 (MQ=9)
                                                                                                  t.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCgggccgcagggcc      >  1:2445950/1‑101 (MQ=14)
                                                                                                  t.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCgggccgcagggcc      <  2:3334533/101‑1 (MQ=14)
                                                                                                    aaaCCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCgggccgcagggccg     >  2:1990618/1‑101 (MQ=14)
                                                                                                     aaCCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCgggccgcagggccgc    >  1:951809/1‑101 (MQ=14)
                                                                                                     aaCCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCgggccgcagggccgc    >  2:2594452/1‑101 (MQ=14)
                                                                                                      aCCCCATCCTTTCTCCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCgggccgcagggccgca   <  1:2532924/101‑1 (MQ=9)
                                                                                                       ccccATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCAc  <  2:211985/101‑1 (MQ=14)
                                                                                                       ccccATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCAc  >  1:625181/1‑101 (MQ=14)
                                                                                                       ccccATCCTTTCACCTTGCCGCAAATATCCCGGGGGGGGAATTGGCCGAAGGCCAAGAGGGGGTCA..GCGCGCCCTTTTCGGGGCGGGCCGCAGGGCCGCAc  >  2:2830202/1‑101 (MQ=9)
                                                                                                       |                                                                                                      
CTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCT.C..CGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCC....GCCAAGATGGGG..G..GC.CCCCCTTTTCGGGGCGGGCCGCAGGGCCGCAC  >  NC_014034/1303775‑1303965

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