Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,608,042 G→C 13.1% intergenic (+312/‑35) RCAP_rcc01484 → / → murD hypothetical protein/UDP‑N‑acetylmuramoylalanine‑‑D‑glutamate ligase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,608,0420G→C13.1% 19.6 ‑0.8 ‑0.077intergenic (+312/‑35)RCAP_rcc01484/murDhypothetical protein/UDP‑N‑acetylmuramoylalanine‑‑D‑glutamate ligase
Reads supporting (aligned to +/- strand):  new base (4/6):  ref base (16/50):  total (20/57)
Fisher's exact test strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.84e-01

TTCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCCTTGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGGCCGCGCGGCGCGGGCCCAATTCACCCCCGCAGCTTACGGA  >  NC_014034/1607942‑1608068
                                                                                                    |                          
ctCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGCGGCcc                          >  1:5541712/2‑101 (MQ=38)
ctCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGCGGCcc                          <  1:678386/100‑1 (MQ=38)
ctCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGCGGCcc                          <  2:156570/100‑1 (MQ=38)
ctCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGCGGCcc                          >  1:7784027/2‑101 (MQ=38)
ctCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGAGCGCCCCCCTCTTGGTCCGC.TGGCGCGGCcc                          >  1:5863920/2‑101 (MQ=255)
 tCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGCGGCcca                         <  2:4049335/101‑2 (MQ=38)
 tCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGCGGCcca                         >  2:4866936/1‑100 (MQ=38)
 tCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGCGGCcca                         <  2:7953867/101‑2 (MQ=38)
 tCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGCGGCcca                         <  1:5250225/101‑2 (MQ=38)
 tCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGCGGCcca                         <  2:7570935/101‑2 (MQ=38)
 tCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGAGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.CGGCGCGGAcca                         >  2:2502712/1‑100 (MQ=38)
  ccGCCCTTTCCCCGGCGGACCCGGCGCCGGGGT.CGGC.CCCCGCCGCCCCCGCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAAt                      <  1:1921035/101‑1 (MQ=9)
  ccGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGCTC.GCGC.GGCCCaa                       >  2:7057926/1‑101 (MQ=38)
  ccGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGCG.GGCGC.GGCCCaa                       <  1:8129816/101‑1 (MQ=38)
  ccGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGCGGCCCCa                        >  1:2368578/1‑101 (MQ=38)
  ccGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCaa                       <  2:1148519/101‑1 (MQ=38)
  ccGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCaa                       <  1:6809175/101‑1 (MQ=38)
  ccGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCaa                       >  1:2875827/1‑101 (MQ=38)
  ccGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCaa                       >  2:4346557/1‑101 (MQ=38)
   cGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAAt                      <  2:2233097/101‑1 (MQ=37)
   cGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCAACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAAt                      <  2:345376/101‑1 (MQ=37)
    gACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAAtt                     <  1:9068334/101‑1 (MQ=38)
    gACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAAtt                     <  2:6596926/101‑1 (MQ=37)
    gACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAAtt                     <  2:9159796/101‑1 (MQ=38)
    gACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAAtt                     <  2:7405228/101‑1 (MQ=38)
     gccTTTCCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTc                    <  2:4031076/100‑1 (MQ=37)
     aCCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTc                    <  2:4831057/101‑1 (MQ=37)
     aCCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTc                    <  2:8744814/101‑1 (MQ=37)
     aCCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGAGGCGCCGCCCGAAAGGGGGGCTGCGGGCCCCCCACTTGGGCCGC.AGGCG.GGGCCCaagcc                    >  2:8927082/1‑98 (MQ=255)
      ccTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCa                   <  1:5957303/101‑1 (MQ=37)
      ccTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCa                   <  2:4047231/101‑1 (MQ=37)
      ccTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCa                   <  1:273306/101‑1 (MQ=37)
      ccTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCa                   <  1:848760/101‑1 (MQ=37)
      ccTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCa                   <  2:6209851/101‑1 (MQ=37)
       cTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCAc                  >  2:3785440/1‑101 (MQ=38)
        ttttCCCAGCGGATCCCGCGACGGGGTGCGACGGC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCAcc                 >  1:4763555/1‑101 (MQ=37)
        ttttCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCAcc                 <  1:8687340/101‑1 (MQ=37)
         tttCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCAccc                >  1:4117545/1‑101 (MQ=38)
         tttCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCAccc                <  1:761705/101‑1 (MQ=37)
          ttCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCAcccc               <  2:6030505/101‑1 (MQ=37)
          ttCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCAcccc               <  2:7525882/101‑1 (MQ=38)
           tCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGGCCCCCTCTTGGCCCGC.TGGCGCGGCCCCATCTCCcccc               >  2:6887919/1‑101 (MQ=255)
           tCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCAccccc              <  1:7793885/101‑1 (MQ=38)
           tCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCAccccc              <  1:2276010/101‑1 (MQ=38)
           tCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCAccccc              <  2:880224/101‑1 (MQ=38)
           tCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCAccccc              <  1:7006702/101‑1 (MQ=37)
           tCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCAccccc              <  2:3366443/101‑1 (MQ=37)
           tCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCAccccc              >  1:4837367/1‑101 (MQ=38)
           tCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCAccccc              <  1:538/101‑1 (MQ=37)
           tCCCAGCGGATCCCGCGACGGAGTGCGACGCC.TGCGGCGCCGCACGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCAccccc              <  2:5990203/101‑1 (MQ=37)
            cccAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GTCCCAATTCACCCCCg             <  1:7264149/101‑1 (MQ=37)
            cccAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCg             >  2:7165992/1‑101 (MQ=38)
            cccAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCg             >  1:1518052/1‑101 (MQ=38)
            cccAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCg             >  2:7508360/1‑101 (MQ=38)
             ccAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCTCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGc            >  2:7430445/1‑101 (MQ=255)
              cAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGGCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcg           <  1:1278775/101‑2 (MQ=37)
              cAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcg           >  2:2311062/1‑100 (MQ=38)
              cAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcg           <  2:7316183/101‑2 (MQ=38)
              cAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcg           >  1:5736609/1‑100 (MQ=38)
               ggcgtgtCCCGCGACGGTGTGCCAC.CCCTGCGGCGCCCCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgg          <  2:7613014/95‑3 (MQ=19)
               aGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGAGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCGATTCACCCCCGcgg          >  2:5233808/1‑99 (MQ=37)
                gCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcggg         >  1:4554602/1‑98 (MQ=37)
                 cGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcggga        <  2:2744594/101‑5 (MQ=37)
                 cGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcggga        >  1:3003729/1‑97 (MQ=37)
                  ggATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggat       <  2:7670588/101‑6 (MQ=37)
                  ggATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggat       >  2:7411857/1‑96 (MQ=37)
                  ggATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggat       >  1:533914/1‑96 (MQ=37)
                  ggATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggat       <  1:5513982/101‑6 (MQ=37)
                  ggATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggat       <  2:518095/101‑6 (MQ=37)
                   gagcccGCGACGGGGGGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggata      <  1:1921148/98‑7 (MQ=25)
                   gATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggata      <  2:3139401/101‑7 (MQ=37)
                    aTCCCGCGGCGGGGTTCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggatat     <  1:2752091/101‑8 (MQ=25)
                     tCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggatatt    <  1:1527482/101‑9 (MQ=37)
                     tCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggatatt    <  1:3469065/101‑9 (MQ=37)
                     tCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggatatt    <  2:1906782/101‑9 (MQ=37)
                     tCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggatatt    <  2:1858061/101‑9 (MQ=37)
                     tCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCCTCGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggatctt    >  1:5459141/1‑93 (MQ=255)
                     tCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGACCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggatatt    >  2:4479082/1‑93 (MQ=255)
                      cccGCGACGGGGTGCGACGCC.TGCGGCTCCGCCCGAGAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggatattt   <  1:4246234/101‑10 (MQ=25)
                      cccGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggatattt   <  2:2794824/101‑10 (MQ=25)
                      cccGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggatattt   <  2:3087238/101‑10 (MQ=25)
                      cccGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggatattt   >  1:2104161/1‑92 (MQ=25)
                      cccGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggatattt   <  1:8582431/101‑10 (MQ=25)
                       ccGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggatatttc  <  2:5439557/101‑11 (MQ=25)
                       ccGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggatatttc  <  1:9196086/101‑11 (MQ=25)
                       ccGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggatatttc  <  1:8425344/101‑11 (MQ=25)
                                                                                                    |                          
TTCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCCTTGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGGCCGCGCGGCGCGGGCCCAATTCACCCCCGCAGCTTACGGA  >  NC_014034/1607942‑1608068

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 24 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG