Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,303,831 +T 97.0% intergenic (+18/+139) 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,8312.→T97.0% 4.6 ‑3.8 0.0134intergenic (+18/+139)urtE/RCAP_rcc01231urea ABC transporter ATP‑binding protein UrtE/UbiH/UbiF/VisC/COQ6 family ubiquinone biosynthesis hydroxylase
Reads supporting (aligned to +/- strand):  new base (49/81):  ref base (4/0):  total (53/81)
Fisher's exact test strand distribution p-value = 2.28e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 167 total aligned reads displayed.

GCTCAAACGCGGCGCGGTGGCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCT.C..CGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGCC  >  NC_014034/1303732‑1303921
                                                                                                       |                                                                                            
gCTCAAACGCGGCGCGGTGGCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.c..c                                                                                            <  1:1235315/101‑1 (MQ=255)
gCTCAAACGCGGCGCGGTGGCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.c..c                                                                                            <  1:1771397/101‑1 (MQ=255)
  tCAAACGCGGCGCGGTGGCGATGGAGTGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.c..ctc                                                                                          >  1:317666/1‑99 (MQ=255)
    aaaCGCGGCGCGGTGGCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.c..ctcgg                                                                                        <  1:992179/101‑5 (MQ=255)
    aaaCGCGGCGCGGTGGCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.c..ctcgg                                                                                        <  2:3919965/101‑5 (MQ=255)
    aaaCGCGGCGCGGTGGCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGAAAGGGGGCGCTG..CCC.c..ctcgg                                                                                        >  1:4499100/1‑97 (MQ=255)
      aCGCGGCGCGGTGGCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.c..ctcggct                                                                                      <  1:1057551/101‑7 (MQ=255)
        gcggcgcggTGGCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.c..ctcggcttc                                                                                    <  2:1994459/101‑9 (MQ=255)
        gcggcgcggTGGCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.c..ctcggcttc                                                                                    <  2:2723933/101‑9 (MQ=255)
        gcggcgcggTGGCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.c..ctcggcttc                                                                                    <  2:3739624/101‑9 (MQ=255)
         cggcgcggTGGCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.c..ctcggcttcg                                                                                   <  1:3325655/101‑10 (MQ=255)
          ggcgcggTGGCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGc                                                                                    >  2:1748276/1‑101 (MQ=255)
          ggcgcggTGGCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCGG..CCC.CCTCGGCTTCGc                                                                                    >  1:3552726/1‑101 (MQ=255)
           gcgcggTGGCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGcc                                                                                   >  1:133698/1‑101 (MQ=255)
             gcggTGGCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTc                                                                                 <  1:1495276/101‑1 (MQ=255)
              cggTGGCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCATGGTCAGCGTTTAGGCAAGGGGTCGGTG..CCC.CCTCGGCTTCGCCTcc                                                                                >  2:3427209/1‑100 (MQ=255)
                 tGGCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCAccc                                                                             <  2:3237410/101‑1 (MQ=255)
                   gCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGATGCCCCC.C..CGGCTTCCCCTCGccccc                                                                           >  2:4407804/1‑101 (MQ=255)
                    cGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCCCCTCACCCCCg                                                                          >  2:2532750/1‑101 (MQ=255)
                     gATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCgg                                                                         >  1:3152089/1‑101 (MQ=39)
                      aTGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCggg                                                                        <  2:1341999/101‑1 (MQ=39)
                      aTGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCggg                                                                        <  1:3225517/101‑1 (MQ=39)
                       tgtaGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGa                                                                       <  2:471683/98‑1 (MQ=38)
                        ggAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGat                                                                      <  2:248994/101‑1 (MQ=38)
                        ggAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGat                                                                      <  1:248800/101‑1 (MQ=38)
                         gAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGata                                                                     <  1:1209026/101‑1 (MQ=38)
                         gAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGata                                                                     >  1:2163331/1‑101 (MQ=38)
                         gAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGata                                                                     <  1:2018953/101‑1 (MQ=38)
                           gggCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATAtt                                                                   <  1:618758/101‑1 (MQ=38)
                            ggCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATAttt                                                                  <  2:1246470/101‑1 (MQ=38)
                             gCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTg                                                                 <  2:1886666/101‑1 (MQ=38)
                              cACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTgg                                                                <  1:2164282/101‑1 (MQ=38)
                              cACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTgg                                                                >  1:350301/1‑101 (MQ=38)
                               aCGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGCGCTTTggg                                                               >  2:4438322/1‑101 (MQ=38)
                                 gAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGc                                                             <  1:1036236/101‑1 (MQ=38)
                                   aGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCaa                                                           <  1:404474/101‑1 (MQ=37)
                                   aGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCaa                                                           <  2:3833843/101‑1 (MQ=37)
                                    ggACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAg                                                          <  2:2915276/101‑1 (MQ=37)
                                    ggACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAg                                                          <  2:4498805/101‑1 (MQ=37)
                                     gACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGa                                                         >  1:2789282/1‑101 (MQ=37)
                                       cAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.a                                                      <  1:2940153/101‑1 (MQ=37)
                                        aGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.aa                                                     <  1:4305326/101‑1 (MQ=37)
                                         gCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.aaa                                                    <  2:3607043/101‑1 (MQ=37)
                                         gCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.aaa                                                    >  2:2575792/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)
                                            tGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAccc                                                 >  1:4448669/1‑101 (MQ=37)
                                             gTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAcccc                                                >  1:3002695/1‑101 (MQ=37)
                                              tCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGCGATATTTGGGGCAAGAT.AAACCCCa                                               >  1:519772/1‑101 (MQ=37)
                                                  gcgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCt                                           <  2:3374429/101‑1 (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                                          <  1:2820522/101‑1 (MQ=37)
                                                   cgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCtt                                          <  2:3151845/101‑1 (MQ=37)
                                                   cgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCtt                                          <  2:3753589/101‑1 (MQ=37)
                                                     aaaCCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTc                                        <  1:4546506/101‑1 (MQ=37)
                                                     aaaCCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCCCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTc                                        >  2:1331851/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                                     >  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)
                                                             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)
                                                                gcCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCa                             <  2:3567447/101‑1 (MQ=37)
                                                                 ccGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCaa                            <  2:3000834/101‑1 (MQ=37)
                                                                   gTGGTCAGCGTTTAGGCAAGGGAGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAat                          <  2:4474152/101‑1 (MQ=37)
                                                                    tGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAata                         <  1:2986212/101‑1 (MQ=37)
                                                                    tGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAata                         <  2:4252110/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:4116336/101‑1 (MQ=37)
                                                                     ggTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAatat                        <  1:1129550/101‑1 (MQ=37)
                                                                      gTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATc                       <  1:2185199/101‑1 (MQ=37)
                                                                      gTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATc                       <  2:568293/101‑1 (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)
                                                                         aGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCg                    <  1:3876810/101‑1 (MQ=37)
                                                                         aGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCg                    >  1:3606527/1‑101 (MQ=25)
                                                                          gCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCgg                   <  1:442255/101‑1 (MQ=37)
                                                                          gCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCgg                   <  1:186243/101‑1 (MQ=37)
                                                                          gCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCgg                   <  2:2272683/101‑1 (MQ=37)
                                                                            ctgacggacggggggCTCTGC.CCC.C..CGGCTTCG.CTC.CCCCCGGGATATTTGCGGCAAGATGAAGGGGCATGATTTCATCTTGCCCGAAATATCCCaggggg               >  1:3167541/11‑101 (MQ=2)
                                                                             tttAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCggggg                >  1:3705785/1‑101 (MQ=37)
                                                                                aGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCgggggggg             >  1:541675/1‑99 (MQ=25)
                                                                                aGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTg             <  1:2882259/101‑1 (MQ=37)
                                                                                aGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTg             <  1:457214/101‑1 (MQ=37)
                                                                                 ggCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGAAAAAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGa            <  1:4653332/101‑1 (MQ=37)
                                                                                    aaGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAAtt         <  2:4271463/101‑1 (MQ=37)
                                                                                    aaGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAAtt         <  1:762851/101‑1 (MQ=37)
                                                                                      gggggCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTgg       >  1:3012954/1‑101 (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.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGgccg    <  2:3069157/101‑1 (MQ=37)
                                                                                                       |                                                                                            
GCTCAAACGCGGCGCGGTGGCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCT.C..CGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGCC  >  NC_014034/1303732‑1303921

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