Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,303,844 Δ1 bp 2.2% intergenic (+31/+126) 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,8440A→. 2.2% 3.7 ‑1.2 ‑1.9182intergenic (+31/+126)urtE/RCAP_rcc01231urea ABC transporter ATP‑binding protein UrtE/UbiH/UbiF/VisC/COQ6 family ubiquinone biosynthesis hydroxylase
Reads supporting (aligned to +/- strand):  new base (0/4):  ref base (65/108):  total (67/115)
Fisher's exact test strand distribution p-value = 2.98e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.21e-02

Only 100 of 193 total aligned reads displayed.

GGTGGCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCTC..CGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCG  >  NC_014034/1303747‑1303919
                                                                                                     |                                                                             
ggTGGCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCAc                                                                              <  1:770578/101‑1 (MQ=255)
 gTGGCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCAcc                                                                             >  2:2705618/1‑101 (MQ=255)
  tGGCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCAccc                                                                            <  2:2082976/101‑1 (MQ=255)
  tGGCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCAccc                                                                            >  2:1170963/1‑101 (MQ=255)
   ggCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCAcccc                                                                           <  1:3724255/101‑1 (MQ=255)
    gCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCAccccc                                                                          <  1:1026270/101‑1 (MQ=39)
     cGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCCcccccc                                                                         >  1:3757657/1‑100 (MQ=255)
     cGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCg                                                                         >  2:1217857/1‑101 (MQ=39)
      gATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTTACCCCCgg                                                                        >  1:1754909/1‑101 (MQ=255)
      gATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCgg                                                                        >  2:1962614/1‑101 (MQ=39)
       aTGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGGTG..CCCCCTCGGCTTCGCCTCATCCCCggg                                                                       >  2:1291626/1‑101 (MQ=255)
       aTGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCggg                                                                       <  2:518430/101‑1 (MQ=39)
         ggAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGat                                                                     <  2:3680532/101‑1 (MQ=38)
         ggAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGat                                                                     <  2:360240/101‑1 (MQ=38)
         ggAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGat                                                                     <  2:622757/101‑1 (MQ=38)
         ggAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCCGGat                                                                     >  1:1069662/1‑101 (MQ=38)
         ggAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCGG..CCCCCTCGGCTTCGCCTCACCCCCGGGat                                                                     >  1:2060014/1‑101 (MQ=255)
          gAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGata                                                                    <  2:2541706/101‑1 (MQ=38)
           aGGGCACGAAGGACATCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCCGGagat                                                                   >  2:2156165/1‑101 (MQ=37)
           aGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGatat                                                                   <  1:238094/101‑1 (MQ=38)
           aGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGatat                                                                   <  1:3949095/101‑1 (MQ=38)
           aGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGatat                                                                   <  2:70487/101‑1 (MQ=38)
           aGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCCGGatat                                                                   >  2:1717776/1‑101 (MQ=37)
            gggCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATAtt                                                                  <  2:3870730/101‑1 (MQ=38)
               cACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTgg                                                               <  2:1809427/101‑1 (MQ=38)
               cACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTgg                                                               <  2:471414/101‑1 (MQ=38)
                aCGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTggg                                                              <  1:3415657/101‑1 (MQ=38)
                  gAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCCCCCCCGGGATATTTGGGGc                                                            >  2:3730921/1‑101 (MQ=38)
                  gAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGc                                                            <  1:3376921/101‑1 (MQ=38)
                  gAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGc                                                            >  1:537155/1‑101 (MQ=38)
                   aaGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCa                                                           <  2:3757404/101‑1 (MQ=38)
                   aaGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCa                                                           <  2:2105852/101‑1 (MQ=38)
                   aaGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCa                                                           >  2:2672922/1‑101 (MQ=38)
                    aGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCaa                                                          <  2:1895797/101‑1 (MQ=37)
                     ggACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAg                                                         <  1:222855/101‑1 (MQ=37)
                     ggACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAg                                                         <  2:720635/101‑1 (MQ=37)
                      gACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGa                                                        >  2:2567570/1‑101 (MQ=37)
                      gACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGa                                                        <  1:1868420/101‑1 (MQ=37)
                        cAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.a                                                     <  1:3025391/101‑1 (MQ=37)
                         aGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGGGGTGCCCCCC..CGGCTTCGCCTCCCCCCCCGGGTATTTGGGGGAAGAT.aa                                                    >  1:1552880/1‑101 (MQ=255)
                          gCCTGTCCCGCGAACCCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGATG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.aaa                                                   <  2:3396055/101‑1 (MQ=25)
                          gCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.aaa                                                   <  1:4429985/101‑1 (MQ=37)
                           ccTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAc                                                  <  1:372764/101‑1 (MQ=37)
                            cTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAc.c                                                <  1:909789/101‑1 (MQ=37)
                            cTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCCGGGTATTTGGGGCAAGAT.AAAc.c                                                >  2:3936609/1‑101 (MQ=37)
                             tGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCCCCCCCGGGATATTTGGGGCAAGAT.AAAc.cc                                               >  2:941612/1‑101 (MQ=37)
                             tGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAc.cc                                               >  1:1415423/1‑101 (MQ=37)
                              gTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAc.ccc                                              >  1:675918/1‑101 (MQ=37)
                               tCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCCGCCCCCC..CGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAA.CCCa                                             >  1:3584021/1‑101 (MQ=37)
                                cccGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCAt                                            <  1:981910/101‑1 (MQ=37)
                                 ccGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATc                                           <  1:2455505/101‑1 (MQ=37)
                                 ccGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATc                                           <  1:2504877/101‑1 (MQ=37)
                                   gcgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCt                                         <  1:1209386/101‑1 (MQ=37)
                                   gcgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCt                                         >  2:2713562/1‑101 (MQ=37)
                                    cgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCtt                                        <  2:2294799/101‑1 (MQ=37)
                                    cgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCtt                                        <  2:3924512/101‑1 (MQ=37)
                                    cgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCGCACCCCCGGGATATTTGGGGCAAGAT.AAAA.CCCATCCtt                                        >  2:2516216/1‑101 (MQ=25)
                                     gAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCttt                                       >  2:549019/1‑101 (MQ=37)
                                      aaaCCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTc                                      >  1:1964968/1‑101 (MQ=37)
                                      aaaCCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTc                                      <  2:1153367/101‑1 (MQ=37)
                                      aaaCCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTc                                      <  1:4183610/101‑1 (MQ=37)
                                      aaaCCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAA.CCCATCATTTc                                      >  2:1634156/1‑101 (MQ=25)
                                       aaCCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCa                                     <  1:560586/101‑1 (MQ=37)
                                       aaCCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCa                                     <  1:2941568/101‑1 (MQ=37)
                                        aCCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCAc                                    <  2:1520131/101‑1 (MQ=37)
                                         cccTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCAcc                                   <  1:2909380/101‑1 (MQ=37)
                                          ccTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCACCt                                  <  1:268955/101‑1 (MQ=37)
                                           cTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCACCtt                                 <  2:958495/101‑1 (MQ=37)
                                            tGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCACCTTg                                <  1:1029661/101‑1 (MQ=37)
                                            tGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCACCTTg                                <  1:2340026/101‑1 (MQ=37)
                                             gcgcgcCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCACCTTGc                               >  1:4069879/1‑101 (MQ=37)
                                             gcgcgcCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCACCTTGc                               <  2:232175/101‑1 (MQ=37)
                                              cgcgcCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCACCTTGcc                              >  2:268817/1‑101 (MQ=37)
                                               gcgcCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCACCTTGCCg                             >  1:381672/1‑101 (MQ=37)
                                               gcgcCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCACCTTGCCg                             >  2:198702/1‑101 (MQ=37)
                                                   cGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCACCTTGCCGCaaa                         >  2:676889/1‑101 (MQ=37)
                                                   cGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCACCTTGCCGCaaa                         <  1:2567733/101‑1 (MQ=37)
                                                   cGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCACCTTGCCGCaaa                         <  1:3252659/101‑1 (MQ=37)
                                                    gTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCACCTTGCCGCAAat                        >  1:945980/1‑101 (MQ=25)
                                                      ggTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCTCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCACCTTGCCGCAAatat                      <  1:4408920/101‑1 (MQ=37)
                                                       gTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCACCTTGCCGCAAATATc                     <  2:2894150/101‑1 (MQ=37)
                                                       gTCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCACCTTGCCGCAAATATc                     <  2:4366074/101‑1 (MQ=37)
                                                        tCAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCACCTTGCCGCAAATATcc                    <  2:3109187/101‑1 (MQ=37)
                                                         cAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCACCTTGCCGCAAATATccc                   <  2:1415316/101‑1 (MQ=37)
                                                         cAGCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCACCTTGCCCCAAATATccc                   >  1:4031346/1‑101 (MQ=25)
                                                           gCGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCCCCCCCGGGGTATTTGGGGCAAGAT.AAAC.CCCCTCCTTTCACCTTGCCGCAAATATCCCgg                 >  2:4200577/1‑101 (MQ=25)
                                                            cGTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCACCTTGCCGCAAATATCCCggg                >  2:3450970/1‑101 (MQ=25)
                                                             gTTTAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCACCTTGCCGCAAATATCCCgggg               <  1:2713752/101‑1 (MQ=37)
                                                                tAGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGt            <  1:3450426/101‑1 (MQ=37)
                                                                 cggacggggggCTCTGC.CCCC..CGGCTTCG.CTC.CCCCCGGGATATTTGCGGCAAGAT.GAAGGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAAt        <  1:3652129/95‑1 (MQ=2)
                                                                 aGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTg           >  1:2803911/1‑101 (MQ=37)
                                                                 aGGCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCCGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTg           >  2:3034117/1‑101 (MQ=25)
                                                                   gCAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGGGaa         >  1:3093844/1‑101 (MQ=25)
                                                                    cAAGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAAt        <  1:2281954/101‑1 (MQ=37)
                                                                     ctgggggCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAAtt       <  1:2093029/99‑1 (MQ=25)
                                                                     aaGGGGGCGCTG..CCCCCTCGGCTTCGCCTCGCCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAAtt       <  1:2705808/101‑1 (MQ=25)
                                                                     aaGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCCGGATATTTGGGGCAAGAT.AAAC.CCCATACCTGTTCCTTGCCACAAATATGCCgggggggcaact       >  1:3837157/1‑95 (MQ=12)
                                                                      aGGGGGCGCTG..CCCCCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTg      <  2:3443337/101‑1 (MQ=37)
                                                                       gggggCTCTGC.CCCC..CGGCTTCG.CTC.CCCCCGGGATATTTGCGGCAAGATGAAGG.GGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGgcca  <  2:4436892/101‑2 (MQ=2)
                                                                       gggggCTCTGC.CCCC..CGGCTTCG.CTC.CCCCCGGGATATTTGCGGCAAGAT.GAAGGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGgcca  <  2:4155742/101‑2 (MQ=2)
                                                                                                     |                                                                             
GGTGGCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCTC..CGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAC.CCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCG  >  NC_014034/1303747‑1303919

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