Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,396,539 T→A 15.1% C141S (TGC→AGC)  RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,396,5390T→A15.1% 103.5 ‑1.6 ‑0.3238C141S (TGC→AGC) RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (14/22):  ref base (102/100):  total (116/122)
Fisher's exact test strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.26e-01

Only 100 of 254 total aligned reads displayed.

TCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGG  >  NC_014034/1396439‑1396639
                                                                                                    |                                                                                                     
tCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                                                                                       <  1:961432/101‑1 (MQ=255)
tCGCTCCCCCCGCGCGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                                                                                       <  2:2817391/101‑1 (MQ=255)
  gCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGc                                                                                                     <  2:1546214/101‑1 (MQ=255)
    tCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCtt                                                                                                   <  2:1201937/101‑1 (MQ=255)
    cccccccGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCtt                                                                                                   <  1:1483258/100‑1 (MQ=255)
     ccccccGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCttt                                                                                                  <  1:109823/101‑1 (MQ=255)
     ccccccGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCttt                                                                                                  >  1:1418578/1‑101 (MQ=255)
     ccccccGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCttt                                                                                                  >  1:4020715/1‑101 (MQ=255)
        cccGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.Gt                                                                                              <  2:2176291/101‑1 (MQ=255)
        cccGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.Gt                                                                                              <  1:2011050/101‑1 (MQ=255)
            cgGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGa                                                                                          <  2:3959987/101‑1 (MQ=255)
             gggggCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGAc                                                                                         >  1:3481857/1‑101 (MQ=255)
                 gCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGc                                                                                     <  1:3104114/101‑1 (MQ=255)
                   gTGGATCTAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCa                                                                                   <  1:709101/101‑1 (MQ=255)
                   gTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCa                                                                                   >  2:3824226/1‑101 (MQ=255)
                   gTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCa                                                                                   <  2:893189/101‑1 (MQ=255)
                   gTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCa                                                                                   >  2:1095347/1‑101 (MQ=255)
                    tGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAg                                                                                  <  2:3748515/101‑1 (MQ=255)
                     ggATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAgc                                                                                 <  2:2016854/101‑1 (MQ=255)
                     ggATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAgc                                                                                 <  2:3727252/101‑1 (MQ=255)
                      gATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAgcc                                                                                >  2:127670/1‑101 (MQ=255)
                       aTCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAgccg                                                                               >  1:781544/1‑101 (MQ=255)
                         cGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGGCAgccgcc                                                                             >  1:587087/1‑101 (MQ=255)
                            aaCGGACCTGGAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCAAAACTGCTTTC.GTGGGACGGGCCAGccgccccc                                                                          >  1:3654518/1‑99 (MQ=255)
                            aaCGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAACCTGCTTTC.GTGGGACTGGCCAGccgcccgc                                                                          >  1:280697/1‑101 (MQ=255)
                            aaCGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGAACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGccgcccgc                                                                          <  1:2950112/101‑1 (MQ=255)
                             aCGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCcgcg                                                                         >  2:475075/1‑101 (MQ=255)
                              cGGCCCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCgcgc                                                                        <  2:1745843/101‑1 (MQ=255)
                              cGGACCTGAAATTCCGCGGGAAACTCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCgcgc                                                                        <  1:2574129/101‑1 (MQ=255)
                                 aCCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCaa                                                                     <  2:3131673/101‑1 (MQ=255)
                                 aCCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCaa                                                                     >  2:3500584/1‑101 (MQ=255)
                                   cTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGt                                                                   >  1:408501/1‑101 (MQ=255)
                                   cTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGt                                                                   >  1:2112474/1‑101 (MQ=255)
                                    tGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTc                                                                  <  1:839496/101‑1 (MQ=255)
                                        aTTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTc                                                              <  1:3911975/101‑1 (MQ=255)
                                        aTTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTc                                                              >  1:1764399/1‑101 (MQ=255)
                                        aTTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGCCGCTc                                                              <  1:270410/101‑1 (MQ=255)
                                          tCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTccc                                                            >  1:2366114/1‑101 (MQ=255)
                                              cgGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCcg                                                        <  2:2439207/101‑1 (MQ=255)
                                              cgGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCCCTccccccc                                                        >  2:3151363/1‑100 (MQ=255)
                                                ggAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTGCCCCcgcg                                                      >  2:3206959/1‑101 (MQ=255)
                                                 gagaCGCGGCGGCGAGCGAGTCGCTCCCCCCGGGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCgg                                                     <  2:1886236/98‑1 (MQ=38)
                                                 gAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCgg                                                     <  2:2924059/101‑1 (MQ=255)
                                                 gAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCgg                                                     <  1:2861485/101‑1 (MQ=255)
                                                  aaaCGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCggg                                                    >  1:1518230/1‑101 (MQ=255)
                                                   aaCGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCgggg                                                   <  1:1376638/101‑1 (MQ=255)
                                                     cgcgGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGc                                                 >  1:758223/1‑101 (MQ=255)
                                                        ggcggcgAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTg                                              >  2:147744/1‑101 (MQ=255)
                                                         gcggcgAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGGgg                                             >  2:199098/1‑101 (MQ=255)
                                                          cggcgAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGa                                            <  2:1408839/101‑1 (MQ=255)
                                                          cggcgAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGa                                            >  1:1921959/1‑101 (MQ=255)
                                                          cggcgAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGa                                            >  2:3887248/1‑101 (MQ=255)
                                                           ggcgAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGAt                                           <  1:3236311/101‑1 (MQ=255)
                                                           gcctcgCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAACCTGCTTTC.GTGGGACTCGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGAt                                           <  2:433268/96‑1 (MQ=255)
                                                            gggacggggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGAt                                           <  1:3070347/93‑1 (MQ=2)
                                                            gggacggagTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGAt                                           <  2:3074343/93‑1 (MQ=2)
                                                            gcgAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATc                                          <  1:3373078/101‑1 (MQ=255)
                                                              gacggagTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGATCg                                         <  1:63788/95‑1 (MQ=2)
                                                              gacggagTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGATCg                                         <  2:777493/95‑1 (MQ=2)
                                                              gAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                        >  1:3479219/1‑101 (MQ=255)
                                                               gcggagTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                        <  2:2557425/96‑1 (MQ=2)
                                                               aGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                       >  2:1671555/1‑101 (MQ=255)
                                                                gCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                      <  2:4044722/101‑1 (MQ=255)
                                                                  gcGTCGCGCCCCCCGCGTGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACg                                    <  1:1887697/100‑1 (MQ=255)
                                                                    gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGAt                                  <  2:4192719/101‑1 (MQ=255)
                                                                    gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGagatcggaag                                   >  1:46186/1‑92 (MQ=2)
                                                                      cGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGAttt                                <  2:936574/101‑1 (MQ=255)
                                                                      cGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGAAActgtg                                 >  2:566769/1‑97 (MQ=2)
                                                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTtc                               >  2:2374115/1‑101 (MQ=255)
                                                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTtc                               >  1:2542809/1‑101 (MQ=255)
                                                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTtc                               <  2:2632148/101‑1 (MQ=255)
                                                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGAAActgtgc                                >  1:4226149/1‑96 (MQ=2)
                                                                        cTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGCACTGGCCAGCCCCCCGCGCCAAGTCGCTCCCCCCGCCGGGGCGCCGATCGAAACGATTtcg                              >  1:973087/1‑101 (MQ=37)
                                                                        cTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGAAActgtgcc                               <  1:361102/101‑7 (MQ=2)
                                                                         cccccccGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTtcgt                             <  2:3257949/100‑1 (MQ=255)
                                                                          ccccccGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTtcgtc                            >  2:4121428/1‑101 (MQ=255)
                                                                          ccccccGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTtcgtc                            >  2:1128107/1‑101 (MQ=255)
                                                                          ccccccGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTtcgtc                            <  1:93156/101‑1 (MQ=255)
                                                                          ccccccGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGAAActgtgccgg                             <  1:1143593/101‑9 (MQ=2)
                                                                          ccccccGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGAAActgtgccgg                             >  1:1183267/1‑93 (MQ=2)
                                                                            ccccGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTAGc                          >  2:3921989/1‑101 (MQ=255)
                                                                            ccccGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGAAActgtgccggtt                           <  1:3733122/101‑11 (MQ=2)
                                                                              ccGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCtt                        >  2:1607721/1‑101 (MQ=255)
                                                                              ccGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCtt                        <  1:509974/101‑1 (MQ=255)
                                                                              ccGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCtt                        >  2:509971/1‑101 (MQ=255)
                                                                               cgcgGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTa                       <  2:1021764/101‑1 (MQ=255)
                                                                                 cgGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTacg                     <  2:3991475/101‑1 (MQ=255)
                                                                                 cgGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTacg                     <  2:3039800/101‑1 (MQ=255)
                                                                                  gggggCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTacga                    >  2:3410091/1‑101 (MQ=255)
                                                                                   ggggCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTacgac                   >  1:1037435/1‑101 (MQ=255)
                                                                                      gCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGcc                <  1:1607830/101‑1 (MQ=255)
                                                                                        gTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTa              >  2:3663005/1‑101 (MQ=255)
                                                                                         tGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTAt             <  1:2250318/101‑1 (MQ=255)
                                                                                         tGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTAt             >  1:3142320/1‑101 (MQ=255)
                                                                                          ggATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATc            <  1:1664517/101‑1 (MQ=255)
                                                                                          ggATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATc            >  2:2229406/1‑101 (MQ=255)
                                                                                          ggATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATAGAAACGATTTAGCCGCTTACGACGCCTATc            >  1:3810044/1‑101 (MQ=255)
                                                                                           gATTGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCt           <  1:2123366/101‑1 (MQ=255)
                                                                                           gATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACGATTTCGTCCCCTACGACGCCTCTCt           >  2:3605997/1‑101 (MQ=38)
                                                                                                aaaCTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGAttt      <  1:1966419/101‑1 (MQ=255)
                                                                                                    |                                                                                                     
TCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGG  >  NC_014034/1396439‑1396639

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