Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,303,892 G→C 3.8% intergenic (+79/+78) 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,8920G→C 3.8% 2.8 ‑3.2 ‑0.1106intergenic (+79/+78)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 (45/57):  total (49/57)
Fisher's exact test strand distribution p-value = 4.27e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.47e-01

Only 100 of 120 total aligned reads displayed.

GCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG.CCT.C..CGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCC....GCCAAGATGGGG..G..G.C.CCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCTGTCAGCGCATCCCCAAG  >  NC_014034/1303794‑1303985
                                                                                                       |                                                                                                       
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:2139799/101‑1 (MQ=37)
  gcCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG.CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCa                                                                                                       <  2:3567447/101‑1 (MQ=37)
   ccGTGGTCAGCGTTTAGGCAAGGGGGCGCTG.CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCaa                                                                                                      <  1:2874406/101‑1 (MQ=37)
   ccGTGGTCAGCGTTTAGGCAAGGGGGCGCTG.CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCaa                                                                                                      <  2:3000834/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)
     gTGGTCAGCGTTTAGGCAAGGGAGCGCTG.CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAat                                                                                                    <  2:4474152/101‑1 (MQ=37)
      tGGTCAGCGTTTAGGCAAGGGGGCGCTG.CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAata                                                                                                   <  2:4252110/101‑1 (MQ=37)
      tGGTCAGCGTTTAGGCAAGGGGGCGCTG.CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAata                                                                                                   <  1:3162787/101‑1 (MQ=37)
      tGGTCAGCGTTTAGGCAAGGGGGCGCTG.CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAata                                                                                                   <  1:2986212/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                                                                                                 >  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: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                                                                                             <  2:2272683/101‑1 (MQ=37)
            gCGTTTAGGCAAGGGGGCGCTG.CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCgg                                                                                             <  1:186243/101‑1 (MQ=37)
            gCGTTTAGGCAAGGGGGCGCTG.CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCgg                                                                                             <  1:442255/101‑1 (MQ=37)
              gTTTAGGCAAGGGGGCGCTG.CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCgggg                                                                                           <  1:447572/101‑1 (MQ=37)
              ctgacggacggggggCTCTGCCCC.C..CGGCTTCG.CTC.CCCCCGGGATATTTGCGGCAAGATGAAGGGGCATGATTTCATCTTGCCCGAAATATCCCaggggg                                                                                         >  2:2585515/11‑101 (MQ=2)
              ctgacggacggggggCTCTGCCCC.C..CGGCTTCG.CTC.CCCCCGGGATATTTGCGGCAAGATGAAGGGGCATGATTTCATCTTGCCCGAAATATCCCaggggg                                                                                         >  1:3167541/11‑101 (MQ=2)
               tttAGGCAAGGGGGCGCTG.CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCggggg                                                                                          >  1:3705785/1‑101 (MQ=37)
                 tAGGCAAGGGGGCGCTG.CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCCAATATCCCggggggc                                                                                        >  2:567525/1‑100 (MQ=25)
                  aGGCAAGGGGGCGCTG.CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCgggggggg                                                                                       >  1:541675/1‑99 (MQ=25)
                   ggCAAGGGGGCGCTG.CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGAAAAAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGa                                                                                      <  1:4653332/101‑1 (MQ=37)
                    gCAAGGGGGCGCTG.CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGaa                                                                                     <  1:1192725/101‑1 (MQ=37)
                      aaGGGGGCGCTG.CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAAtt                                                                                   >  2:3190409/1‑101 (MQ=25)
                      aaGGGGGCGCTG.CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAAtt                                                                                   <  2:4271463/101‑1 (MQ=37)
                      aaGGGGGCGCTG.CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAAtt                                                                                   <  1:762851/101‑1 (MQ=37)
                      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)
                           ggCGCTG.CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGgcc....g                                                                          <  2:3069157/101‑1 (MQ=37)
                                                                    tGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGgggcggg                                    >  2:169142/1‑101 (MQ=14)
                                                                         cAAGAT.AACCCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCgggccgca                               <  1:2988281/101‑1 (MQ=9)
                                                                          aaGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCgggccgcag                              >  1:3519800/1‑101 (MQ=14)
                                                                          aaGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCgggccgcag                              <  2:3513518/101‑1 (MQ=14)
                                                                          aaGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGGGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCgggccgcag                              >  1:3361931/1‑101 (MQ=9)
                                                                           aGAT.AAACCCCATCCTTTCACCTTGCCGCAAATACCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCgggccgcagg                             <  2:2228843/101‑1 (MQ=9)
                                                                            gAT.AAACCCCATCCTTTCCCCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCgggccgcaggg                            <  1:696988/101‑1 (MQ=9)
                                                                            gAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCCG.CGCCCCCTTTTCGGGGCgggccgccagg                            >  2:2901542/1‑97 (MQ=9)
                                                                             aT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGGGAATTGGCCGAAGGCCAGGAGGGGG..GCAG.CGCCCCCTTTTAGGGGCgggccgcagggc                           >  1:1387025/1‑101 (MQ=9)
                                                                              t.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCgggccgcagggcc                          <  2:3334533/101‑1 (MQ=14)
                                                                              t.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCgggccgcagggcc                          >  1:2445950/1‑101 (MQ=14)
                                                                                aagccctgtCCTTTCGCCTTGCCGCAAAGGTTCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCgggccgcagggccg                         <  1:2545422/93‑1 (MQ=2)
                                                                                aaaCCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCgggccgcagggccg                         >  2:235268/1‑101 (MQ=14)
                                                                                 aaCCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCgggccgcagggccgc                        >  1:951809/1‑101 (MQ=14)
                                                                                  aCCCCATCCTTTCTCCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCgggccgcagggccgca                       <  1:2532924/101‑1 (MQ=9)
                                                                                   ccccATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCAc                      <  2:211985/101‑1 (MQ=14)
                                                                                   ccccATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCAc                      >  1:625181/1‑101 (MQ=14)
                                                                                   ccccATCCTTTCACCTTGCCGCAAATATCCCGGGGGGGGAATTGGCCGAAGGCCAAGAGGGGGTCA..G.CGCGCCCTTTTCGGGGCGGGCCGCAGGGCCGCAc                      >  2:2830202/1‑101 (MQ=9)
                                                                                    cccATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCAcc                     >  2:4394458/1‑101 (MQ=14)
                                                                                     ccATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCAccc                    >  2:468793/1‑101 (MQ=14)
                                                                                      cATCCTTTCACCTTGCCGCAAATCTCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCt                   <  1:1934072/101‑1 (MQ=9)
                                                                                      cATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCt                   <  1:681799/101‑1 (MQ=14)
                                                                                       ctCCTTTCTCGTTGCGGTATATATCCTGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGT..GCAG.CGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCCCCCTg                  <  2:4553359/100‑1 (MQ=2)
                                                                                         ccTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCTGTc                >  1:244201/1‑101 (MQ=14)
                                                                                         ccTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCTGTc                >  2:2763581/1‑101 (MQ=14)
                                                                                          cTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCTGTCa               >  2:3139949/1‑101 (MQ=14)
                                                                                          cTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCTGTCa               <  2:3417652/101‑1 (MQ=14)
                                                                                           tttCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCTGTCAg              <  2:3606259/101‑1 (MQ=14)
                                                                                           tttCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCTGTCAg              >  1:1922202/1‑101 (MQ=14)
                                                                                            ttCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCTGTCAgc             >  1:3545773/1‑101 (MQ=14)
                                                                                             tccccTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCTGTCAgcg            <  1:2659521/98‑1 (MQ=9)
                                                                                             tCACCTAGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCGGGGCGCAGGGCCGCACCCGGTCAgcg            >  2:1984083/1‑101 (MQ=9)
                                                                                              cACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCTGTCAgcgc           >  1:4546583/1‑101 (MQ=17)
                                                                                              cACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCTGTCAgcgc           >  2:2595452/1‑101 (MQ=17)
                                                                                               aCCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCGGGGCGCAGGGCCGGACCCTGTAAgggcc          >  2:1834988/1‑100 (MQ=2)
                                                                                                ccTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCTGTCAGCGCAt         <  1:2833865/101‑1 (MQ=17)
                                                                                                ccTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCTGTCAGCGCAt         >  2:138254/1‑101 (MQ=17)
                                                                                                ccTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCTGTCAGCGCAt         <  2:2129937/101‑1 (MQ=17)
                                                                                                ccTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCTGTCAGCGCAt         <  1:1950939/101‑1 (MQ=17)
                                                                                                 cTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGACAAGAGGGGG..GCAGAG.CCCCCTTTTCGGGGCGGGGCGCCGGGCCGCCCCCCGTCAGCGCATc        >  1:1346634/1‑101 (MQ=2)
                                                                                                  ttGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCACCCTTTTCGGGGCGGGCCACAGGGCCGCACCCTGTCAGCGCATcc       >  1:3921635/1‑101 (MQ=9)
                                                                                                  ttGCCGCAAATATCCCGGGGGGCGAATTGGCCGAAGGCCAAGAGGGGG..GGAG.CGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCTGTCAGCGCATcc       >  2:3620349/1‑101 (MQ=12)
                                                                                                  ttGCCGAAAATAGCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCTGTCAGCGCATcc       <  1:2308745/101‑1 (MQ=9)
                                                                                                    gCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCTGTCAGCGCATcccc     >  1:4168821/1‑101 (MQ=17)
                                                                                                     ccGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCTGTCAGCGCATCCCCa    >  2:2137204/1‑101 (MQ=17)
                                                                                                      cGCAAATATCCGGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCTGTCAGCGCATCCCCaa   <  1:240545/101‑1 (MQ=12)
                                                                                                      cGCAAATATCCCGGGGGTTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCTGTCAGCGCATCCCCaa   <  2:3034055/101‑1 (MQ=12)
                                                                                                      cGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCTGTCAGCGCATCCCCaa   <  2:2367764/101‑1 (MQ=17)
                                                                                                      cGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCTGTCAGCGCATCCCCaa   >  2:1536373/1‑101 (MQ=17)
                                                                                                      cGCAAATATCCCGGGGGGGGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCGGCACGCAGGGCCGCACCCCGTAAGGGCATCCCCaa   >  2:2675182/1‑101 (MQ=2)
                                                                                                       gCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCTGTCAGCGCATCCCCAAg  >  2:4312132/1‑101 (MQ=17)
                                                                                                       gCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAG.CGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCTGTCAGCGCATCCCCAAg  <  1:3847456/101‑1 (MQ=17)
                                                                                                       |                                                                                                       
GCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG.CCT.C..CGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCC....GCCAAGATGGGG..G..G.C.CCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCTGTCAGCGCATCCCCAAG  >  NC_014034/1303794‑1303985

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 24 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG