Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,608,040 Δ1 bp 87.0% intergenic (+310/‑37) 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,0400G→.87.0% 29.6 ‑1.3 ‑1.077intergenic (+310/‑37)RCAP_rcc01484/murDhypothetical protein/UDP‑N‑acetylmuramoylalanine‑‑D‑glutamate ligase
Reads supporting (aligned to +/- strand):  new base (16/51):  ref base (4/6):  total (20/57)
Fisher's exact test strand distribution p-value = 2.75e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.07e-01

TCTTCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCCTTGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGGCCGCGCGGCGCGGGCCCAATTCACCCCCGCAGCTTACGGA  >  NC_014034/1607940‑1608068
                                                                                                    |                            
gtctCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCgc.tggcgcggc                            <  1:7089966/98‑2 (MQ=38)
gtctCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCgc.tggcgcggc                            >  1:7741775/4‑100 (MQ=38)
gtctCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCgc.tggcgcggc                            <  2:5599256/98‑2 (MQ=38)
gtctCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCgc.tggcgcggc                            <  2:1101812/98‑2 (MQ=38)
 tctCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCgc.tggcgcggcc                           <  1:7720653/99‑3 (MQ=38)
 tctCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCgc.tggcgcggcc                           <  2:6590769/99‑3 (MQ=38)
  ctCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGCGGCcc                          <  1:678386/100‑1 (MQ=38)
  ctCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGCGGCcc                          >  1:7784027/2‑101 (MQ=38)
  ctCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGCGGCcc                          >  1:5541712/2‑101 (MQ=38)
  ctCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGCGGCcc                          <  2:156570/100‑1 (MQ=38)
  ctCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGAGCGCCCCCCTCTTGGTCCGC.TGGCGCGGCcc                          >  1:5863920/2‑101 (MQ=255)
   tCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGCGGCcca                         <  2:7953867/101‑2 (MQ=38)
   tCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGCGGCcca                         <  2:7570935/101‑2 (MQ=38)
   tCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGCGGCcca                         >  2:4866936/1‑100 (MQ=38)
   tCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGCGGCcca                         <  2:4049335/101‑2 (MQ=38)
   tCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGCGGCcca                         <  1:5250225/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                       >  2:4346557/1‑101 (MQ=38)
    ccGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCaa                       >  1:2875827/1‑101 (MQ=38)
    ccGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCaa                       <  1:6809175/101‑1 (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:7405228/101‑1 (MQ=38)
      gACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAAtt                     <  2:9159796/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                   <  2:4047231/101‑1 (MQ=37)
        ccTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCa                   <  2:6209851/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                   <  1:5957303/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:761705/101‑1 (MQ=37)
           tttCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCAccc                >  1:4117545/1‑101 (MQ=38)
            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:4837367/1‑101 (MQ=38)
             tCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCAccccc              <  2:880224/101‑1 (MQ=38)
             tCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCAccccc              <  1:538/101‑1 (MQ=37)
             tCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCAccccc              <  1:7793885/101‑1 (MQ=38)
             tCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCAccccc              <  1:2276010/101‑1 (MQ=38)
             tCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCAccccc              <  1:7006702/101‑1 (MQ=37)
             tCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCAccccc              <  2:3366443/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             >  2:7508360/1‑101 (MQ=38)
              cccAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCg             >  1:1518052/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           >  1:5736609/1‑100 (MQ=38)
                cAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcg           <  2:7316183/101‑2 (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:518095/101‑6 (MQ=37)
                    ggATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggat       <  1:5513982/101‑6 (MQ=37)
                    ggATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggat       >  2:7411857/1‑96 (MQ=37)
                    ggATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggat       >  1:533914/1‑96 (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    <  2:1858061/101‑9 (MQ=37)
                       tCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggatatt    <  2:1906782/101‑9 (MQ=37)
                       tCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggatatt    <  1:1527482/101‑9 (MQ=37)
                       tCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggatatt    <  1:3469065/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   >  1:2104161/1‑92 (MQ=25)
                        cccGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggatattt   <  1:8582431/101‑10 (MQ=25)
                        cccGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggatattt   <  2:2794824/101‑10 (MQ=25)
                        cccGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggatattt   <  2:3087238/101‑10 (MQ=25)
                         ccGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggatatttc  <  1:9196086/101‑11 (MQ=25)
                         ccGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggatatttc  <  1:8425344/101‑11 (MQ=25)
                         ccGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggatatttc  <  2:5439557/101‑11 (MQ=25)
                                                                                                    |                            
TCTTCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCCTTGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGGCCGCGCGGCGCGGGCCCAATTCACCCCCGCAGCTTACGGA  >  NC_014034/1607940‑1608068

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