Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,397,228 C→A 4.0% S370S (TCC→TCA RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,397,2280C→A 4.0% 10.9 ‑3.8 ‑0.1198S370S (TCC→TCARCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (1/7):  ref base (112/78):  total (113/85)
Fisher's exact test strand distribution p-value = 2.20e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.95e-01

Only 100 of 216 total aligned reads displayed.

AAACGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTCCCCCCGCGGGAGCGTGGATCG  >  NC_014034/1397128‑1397328
                                                                                                      |                                                                                                     
aaaCGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATcc                                                                                                       <  1:3146884/101‑1 (MQ=255)
aaaCGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATcc                                                                                                       >  2:3037254/1‑101 (MQ=255)
aaaCGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGAGTCTGCCA.CGAC.GCCCTGCACCAGATcc                                                                                                       >  2:2650/1‑101 (MQ=255)
 aaCGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCcg                                                                                                      <  2:1973916/101‑1 (MQ=255)
 aaCGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCcg                                                                                                      >  1:1916182/1‑101 (MQ=255)
  aCGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCcgc                                                                                                     >  1:2242446/1‑101 (MQ=255)
    gTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCgcgc                                                                                                   >  2:243482/1‑101 (MQ=255)
    gTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCgcgc                                                                                                   >  2:4624959/1‑101 (MQ=255)
     tCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGcc                                                                                                  <  1:1718540/101‑1 (MQ=255)
      ccAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.g                                                                                                <  2:61105/101‑1 (MQ=255)
       cAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.Gt                                                                                               <  1:4625212/101‑1 (MQ=255)
        aGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTc                                                                                              <  2:3531063/101‑1 (MQ=255)
         gCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCg                                                                                             >  1:4313532/1‑101 (MQ=255)
            cgcgCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTc                                                                                          <  1:832878/101‑1 (MQ=255)
              cgcCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTccc                                                                                        <  1:3240935/101‑1 (MQ=255)
              cgcCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTccc                                                                                        >  1:1900238/1‑101 (MQ=255)
                ccACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTccccc                                                                                      <  2:2529536/101‑1 (MQ=255)
                 cacaTGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTcccccc                                                                                     <  2:408907/101‑1 (MQ=255)
                   caTGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCcgc                                                                                   <  1:3562592/101‑1 (MQ=255)
                   caTGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCcgc                                                                                   <  1:2016102/101‑1 (MQ=255)
                    aTGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCcgcg                                                                                  >  2:1265797/1‑101 (MQ=255)
                     tGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCgg                                                                                 <  2:2621619/101‑1 (MQ=255)
                     tGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCgg                                                                                 <  1:1877147/101‑1 (MQ=255)
                      gTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCggg                                                                                >  1:3634401/1‑101 (MQ=255)
                        cTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGTCTGCCA.CGAC.GGCCTCCACCATATCCGCGGC.GCCGCTCCCCCCGCggggg                                                                              >  2:3717516/1‑101 (MQ=37)
                           cGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCgg                                                                           >  2:75198/1‑100 (MQ=255)
                           cGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCgg                                                                           >  1:4487708/1‑100 (MQ=255)
                            gACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAA.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTg                                                                          >  2:1876260/1‑101 (MQ=255)
                             gccGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTgg                                                                         <  1:3037442/100‑1 (MQ=255)
                               cGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                                       <  1:1444410/101‑1 (MQ=255)
                                gTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                      >  2:3598571/1‑101 (MQ=255)
                                 tGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGTGGGCGTGGATCg                                                                     <  1:2577142/101‑1 (MQ=255)
                                 ggcgcCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                     <  1:3362447/100‑1 (MQ=255)
                                  gcggCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGCTGCTGACACCGCC.AACAAATATGTGGTCAACACT.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                     <  1:4562464/101‑1 (MQ=2)
                                  gcggCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGCTGCTGACACCGCC.AACAAATATGTGGTCAACACT.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                     <  1:1591119/101‑1 (MQ=2)
                                    gcCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTGTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCGCCCCCCGCGGGGGCGTGGATCGaaa                                                                  <  1:3196984/101‑1 (MQ=255)
                                     gcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGCTGCTGACACCGCC.AACAAATATGTGGTCAACACT.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                  <  2:337615/100‑1 (MQ=2)
                                     ccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAac                                                                 >  1:3444015/1‑101 (MQ=255)
                                     ccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAac                                                                 <  1:1115594/101‑1 (MQ=255)
                                      cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAaca                                                                <  2:3914230/101‑1 (MQ=255)
                                      cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAaca                                                                <  1:4213506/101‑1 (MQ=255)
                                      cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAaca                                                                <  2:4474970/101‑1 (MQ=255)
                                       gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAacac                                                               >  1:3324558/1‑101 (MQ=255)
                                       gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAacac                                                               <  2:2242322/101‑1 (MQ=255)
                                       gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAacac                                                               >  1:1636018/1‑101 (MQ=255)
                                       gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAacac                                                               >  2:4216297/1‑101 (MQ=255)
                                       gTCGCTCCCCCCGCGGGGGCGTGGACCGAAACCTATGAGGGG.TGGCGGCCCGGTCGTGGATCC.CGCCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaca                                                                <  1:1963394/101‑1 (MQ=0)
                                        tCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAcc                                                              >  2:4663004/1‑101 (MQ=255)
                                        tCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GACGCTCCCCCCGCGGGGGCGTGGATCGAAACAcc                                                              >  1:4071832/1‑101 (MQ=255)
                                          gCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACcgc                                                            >  2:4382313/1‑101 (MQ=255)
                                          gCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACcgc                                                            <  2:793657/101‑1 (MQ=255)
                                            tCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCgcgc                                                          >  1:3968426/1‑101 (MQ=255)
                                             ccccccGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCa                                                         >  2:115182/1‑101 (MQ=255)
                                              cccccGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGTCGTGGATCGAAACACCGCGCAc                                                        <  2:688099/101‑1 (MQ=255)
                                              cccccGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCAc                                                        <  2:1451917/101‑1 (MQ=255)
                                              cccccGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCAc                                                        >  1:4008530/1‑101 (MQ=255)
                                                cccGGGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCAcgc                                                      <  2:586130/101‑1 (MQ=255)
                                                    cgGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGa                                                  >  2:3706405/1‑101 (MQ=255)
                                                    cgGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGa                                                  <  2:2715843/101‑1 (MQ=255)
                                                     gggggCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGaa                                                 <  2:4573648/101‑1 (MQ=255)
                                                     gggggCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGaa                                                 <  1:2333804/101‑1 (MQ=255)
                                                      ggggCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGaag                                                >  1:1819223/1‑100 (MQ=255)
                                                      ggggCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAAc                                                >  2:3137701/1‑101 (MQ=255)
                                                       gggCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACg                                               >  2:3753936/1‑101 (MQ=255)
                                                        ggCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACgg                                              <  2:3232208/101‑1 (MQ=255)
                                                          cGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCg                                            >  2:4170146/1‑101 (MQ=255)
                                                          cGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCg                                            >  2:2052277/1‑101 (MQ=255)
                                                             ggATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCa                                         >  1:771139/1‑101 (MQ=255)
                                                               aTCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGc                                       >  1:1330050/1‑101 (MQ=255)
                                                                tCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGcc                                      <  1:2037716/101‑1 (MQ=255)
                                                                tCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGcc                                      <  2:1506648/101‑1 (MQ=255)
                                                                 cGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGccc                                     >  1:2729149/1‑101 (MQ=255)
                                                                 cGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGccc                                     >  1:598870/1‑101 (MQ=255)
                                                                  gAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGcccc                                    <  2:4313285/101‑1 (MQ=255)
                                                                  gAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGcccc                                    <  1:297238/101‑1 (MQ=255)
                                                                  gAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGcccc                                    <  1:917770/101‑1 (MQ=255)
                                                                   aaaCGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCcg                                   >  2:3430653/1‑101 (MQ=255)
                                                                   aaaCGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGCCAGccccc                                   >  1:2404391/1‑100 (MQ=255)
                                                                     aCGTGTCTGCCA.CGAC.GCCCTGCACTAGATACGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGCTCGAAACACCGCGCACGCGCGCACCGCCTCAGCCCcgcg                                 >  2:3257089/1‑101 (MQ=37)
                                                                       gtgtCTGCCC.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGaa                               <  1:2640068/101‑1 (MQ=255)
                                                                       gtgtCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGaa                               >  1:1002314/1‑101 (MQ=255)
                                                                        tgtCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAAc                              >  1:3748118/1‑101 (MQ=255)
                                                                        tgtCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAAc                              >  1:4527640/1‑101 (MQ=255)
                                                                          tCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGt                            >  2:1606953/1‑101 (MQ=255)
                                                                          tCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGt                            <  1:636198/101‑1 (MQ=255)
                                                                           cTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTc                           >  2:1767154/1‑101 (MQ=255)
                                                                           cTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTc                           <  2:1635918/101‑1 (MQ=255)
                                                                            taacA.CGAA.GCCCTGCACCAGATCCGTGCC.GTCGCCTCCGCCGCGGAGGCGCTGCTCCCACGACCGCGCACGCGCGTAGCGCGGCAGTCCCCCGAACGCCg                          >  2:3901043/4‑101 (MQ=255)
                                                                            tGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCg                          >  1:1557235/1‑101 (MQ=255)
                                                                            tGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCg                          >  2:4369188/1‑101 (MQ=255)
                                                                              ccA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCt                        >  2:905689/1‑101 (MQ=255)
                                                                                    ac.gCCCTGCCCCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTccccc                   <  2:120558/101‑1 (MQ=255)
                                                                                     c.gCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTcccccc                  <  2:1916070/101‑1 (MQ=255)
                                                                                       gCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTCCCCCcg                 <  1:4369435/101‑1 (MQ=255)
                                                                                           tGCACCAGATCCGCGCT.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTCCCCCCGCggg             >  1:723438/1‑101 (MQ=255)
                                                                                           tGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTCCCCCCGCggg             <  1:217218/101‑1 (MQ=255)
                                                                                           tGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTCCCCCCGCggg             <  2:2014887/101‑1 (MQ=255)
                                                                                              aCCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTCCCCCCGCGGGAGc          >  1:2517041/1‑101 (MQ=255)
                                                                                               ccAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTCCCCCCGCGGGAGCg         <  2:3050192/101‑1 (MQ=255)
                                                                                                   aTCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTCCCCCCGCGGGAGCGTGGa     <  2:2728995/101‑1 (MQ=255)
                                                                                                      |                                                                                                     
AAACGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGAC.GCCCTGCACCAGATCCGCGCC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTCCCCCCGCGGGAGCGTGGATCG  >  NC_014034/1397128‑1397328

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