Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,396,967 T→C 15.5% R283R (CGT→CGC RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,396,9670T→C15.5% 133.4 ‑4.6 ‑0.9530R283R (CGT→CGCRCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (38/20):  ref base (173/198):  total (245/285)
Fisher's exact test strand distribution p-value = 1.06e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.15e-01

Only 100 of 552 total aligned reads displayed.

CCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTG  >  NC_014034/1396867‑1397067
                                                                                                              |                                                                                                     
 cGGTACAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTg                                                                                                      >  2:1644415/1‑101 (MQ=255)
  ggTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTg.g                                                                                                    >  1:1651645/1‑101 (MQ=255)
    tCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCGGC.CTG..CGTG.GTc                                                                                                  >  2:402632/1‑101 (MQ=255)
    tCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTc                                                                                                  >  1:5672399/1‑101 (MQ=255)
    tCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTc                                                                                                  >  1:2444277/1‑101 (MQ=255)
    tCCAGAGATAGACCGCGCGGACCTCCTCAGACCTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTc                                                                                                  <  2:5248607/101‑1 (MQ=255)
       agagATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCt                                                                                               >  1:148271/1‑101 (MQ=255)
         agaTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTcc                                                                                             <  2:541266/101‑1 (MQ=255)
            tAGACCGCGCGGCCCTCCTCAGACGTCGGTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGTG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTccccc                                                                                          >  1:5068089/1‑101 (MQ=255)
            tAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTccccc                                                                                          <  1:4821969/101‑1 (MQ=255)
             aGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTcccccc                                                                                         >  1:3363613/1‑101 (MQ=255)
                    gcgGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCggggg                                                                                  <  2:2060205/101‑1 (MQ=255)
                     cgGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGc                                                                                 <  1:3417197/101‑1 (MQ=255)
                       gACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCCC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCgg                                                                               >  2:2152926/1‑100 (MQ=255)
                        gcctcctCAGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTg                                                                              <  1:2795884/100‑1 (MQ=255)
                        aCCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTg                                                                              <  1:6110576/101‑1 (MQ=255)
                         ctcgcggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG..TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTg                                                                              <  1:3973994/91‑1 (MQ=0)
                         ctcgcggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG..TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTg                                                                              <  1:95678/91‑1 (MQ=0)
                         cctcctCAGACGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTgg                                                                             <  2:5715011/101‑1 (MQ=255)
                          ggtcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTgg                                                                             >  1:6745831/11‑101 (MQ=9)
                          attggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                                            <  2:2468786/91‑1 (MQ=0)
                            tggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                          <  2:6433906/93‑1 (MQ=2)
                            tggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                          >  1:5774372/9‑101 (MQ=2)
                            tggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                          >  2:3516712/9‑101 (MQ=2)
                            tcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                                           >  1:1279405/9‑101 (MQ=9)
                             ggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGTGCGTGGATCg                                                                         <  2:6705523/94‑1 (MQ=0)
                             ggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                         <  2:3216650/94‑1 (MQ=2)
                             cggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG..TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                          >  2:4607307/7‑101 (MQ=2)
                             caaccctgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                          <  2:140106/94‑1 (MQ=9)
                              tcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGTGGCGTGGATCGa                                                                        <  2:103894/101‑1 (MQ=255)
                              ggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG..TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                         <  2:3392758/96‑1 (MQ=2)
                              ccAGACGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                        <  2:231277/100‑1 (MQ=255)
                              aacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                         >  1:364820/7‑101 (MQ=9)
                              aacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                         >  1:1954229/7‑101 (MQ=9)
                               cAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                       <  1:4741295/101‑1 (MQ=255)
                                ttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG..TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                       <  1:5204680/98‑1 (MQ=2)
                                ttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG..TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                       <  2:6580384/98‑1 (MQ=2)
                                aGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                      <  1:702851/101‑1 (MQ=255)
                                 ccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                     <  1:491709/98‑1 (MQ=2)
                                  gcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG..TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                     <  2:5003180/100‑1 (MQ=2)
                                  cggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                                                    <  1:6963133/99‑2 (MQ=2)
                                  cggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                                                    >  1:6210111/3‑100 (MQ=2)
                                  ccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                       <  2:604060/100‑1 (MQ=2)
                                   tgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                                                    <  1:5313027/100‑2 (MQ=9)
                                   cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                      >  1:6413968/1‑101 (MQ=2)
                                   cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACac                                                                   >  2:2687303/1‑100 (MQ=255)
                                   cGTCCCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                                                   <  1:592315/101‑1 (MQ=255)
                                    gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                     >  2:3620972/1‑101 (MQ=2)
                                    gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActac                                                                  <  2:1457818/101‑4 (MQ=2)
                                    gTCGCCCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                     <  1:1169098/101‑1 (MQ=2)
                                     tCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCa                                                                 <  1:222824/101‑1 (MQ=255)
                                     tCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcg                                                                  <  1:789161/101‑4 (MQ=9)
                                      cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACac                                                                   <  1:2715877/101‑2 (MQ=2)
                                      cGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAg                                                                <  2:1130326/101‑1 (MQ=255)
                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACacc                                                                  <  1:557487/101‑3 (MQ=2)
                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG..TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcccgcc                                                                <  2:2199894/101‑6 (MQ=2)
                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGc                                                               <  1:3894899/101‑1 (MQ=255)
                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacacc                                                               >  1:3067454/1‑95 (MQ=2)
                                       gCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggc                                                                >  2:3843686/1‑96 (MQ=9)
                                        cTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGcc                                                              >  2:1852078/1‑100 (MQ=255)
                                        cTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCt                                                              <  2:2853730/101‑1 (MQ=255)
                                        atctccccGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCa                                                                 <  2:7083470/97‑1 (MQ=2)
                                         tCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG..TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGGGGATCGAAAcccgccac                                                              >  1:2496554/1‑94 (MQ=2)
                                         tCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcct                                                              >  1:5472416/1‑94 (MQ=9)
                                          ccccccGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaat                                                               <  2:4688640/101‑3 (MQ=2)
                                          ccccccGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CC.AGGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaat                                                               <  2:3173614/101‑3 (MQ=2)
                                          ccccccGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacaccggg                                                            >  2:893597/1‑92 (MQ=2)
                                           cccccGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaatg                                                              <  1:3151241/101‑4 (MQ=2)
                                             cccGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaatgcc                                                            <  1:2098123/101‑6 (MQ=0)
                                               cgcgGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGtt                                                       <  1:4362662/101‑1 (MQ=255)
                                                  gggggCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTggc                                                    >  1:3208237/1‑100 (MQ=255)
                                                   ggggCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGGGGATCGAAACATCAGCTCGGCGTTgggc                                                   >  1:6359208/1‑99 (MQ=255)
                                                     ggCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAg                                                 >  1:2869020/1‑101 (MQ=255)
                                                      gCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGc                                                <  2:6971360/101‑1 (MQ=255)
                                                       cGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGcc                                               <  2:1305590/101‑1 (MQ=255)
                                                        gTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCg                                              >  2:59166/1‑101 (MQ=255)
                                                        gTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCg                                              >  1:4452202/1‑101 (MQ=255)
                                                        gTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCg                                              >  2:4231118/1‑101 (MQ=255)
                                                          ggATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGa                                            >  2:2983762/1‑101 (MQ=255)
                                                            aTCGAACCA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACg                                          <  1:1852150/101‑1 (MQ=255)
                                                            aTCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACg                                          <  1:2167717/101‑1 (MQ=255)
                                                              cGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGAt                                        >  2:4089102/1‑101 (MQ=255)
                                                              cGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGAt                                        >  1:3796274/1‑101 (MQ=255)
                                                               gAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATg                                       <  2:61725/101‑1 (MQ=255)
                                                                aaaCA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCCC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAACCATCAGCTCGGCGTTGGTCAGCCCGACGATGt                                      >  2:5800637/1‑101 (MQ=38)
                                                                   cA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCg                                   <  2:3056437/101‑1 (MQ=255)
                                                                         t.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCt                              <  2:4905796/101‑1 (MQ=255)
                                                                             gcg.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTccc                           >  2:5600468/1‑101 (MQ=255)
                                                                                  c.cGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCcg                       <  2:1983307/101‑1 (MQ=255)
                                                                                    cGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCcgc                      <  1:4186446/101‑1 (MQ=255)
                                                                                     gCA..TCGGCG.TCCGG.CTG..CGTG.GTCCGCCCCCCCGCGGGGGCGTGGAGCGAAACCTCCGCCCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCcgcg                     <  1:2295269/101‑1 (MQ=37)
                                                                                           cGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCggggg                 >  2:6508168/1‑101 (MQ=255)
                                                                                           cGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCggggg                 >  1:462086/1‑101 (MQ=255)
                                                                                            ggCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGc                <  2:2514821/101‑1 (MQ=255)
                                                                                              cT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGt              >  2:1521011/1‑101 (MQ=255)
                                                                                                    gC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATc         <  2:3933235/101‑1 (MQ=255)
                                                                                                    gC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATc         <  1:6518074/101‑1 (MQ=255)
                                                                                                     c.cTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATCg        <  1:1600350/101‑1 (MQ=255)
                                                                                                             gTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt   >  2:4421296/1‑101 (MQ=255)
                                                                                                             gTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt   >  1:1003564/1‑101 (MQ=255)
                                                                                                              |                                                                                                     
CCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTG  >  NC_014034/1396867‑1397067

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