Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,397,539 T→C 15.3% I474T (ATC→GCC)  RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,397,5390T→C15.3% 96.2 ‑1.8 ‑1.2340I474T (ATC→ACC) RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (14/24):  ref base (127/127):  total (161/179)
Fisher's exact test strand distribution p-value = 1.64e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.00e-02

Only 100 of 365 total aligned reads displayed.

TCCCCCCGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCT  >  NC_014034/1397439‑1397638
                                                                                                     |                                                                                                        
  cccccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCa                                                                                                        >  1:5469949/1‑101 (MQ=255)
  cccccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCa                                                                                                        <  1:5096106/101‑1 (MQ=255)
   ccccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCAt                                                                                                       <  2:141552/101‑1 (MQ=255)
   ccccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCAt                                                                                                       <  2:3172504/101‑1 (MQ=255)
    cccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATg                                                                                                      <  1:667877/101‑1 (MQ=255)
    cccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATg                                                                                                      <  1:4831626/101‑1 (MQ=255)
     ccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATgg                                                                                                     <  1:5040128/101‑1 (MQ=255)
       gcgGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGc.c                                                                                                  <  1:1256653/101‑1 (MQ=255)
       gcgGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGc.c                                                                                                  <  2:3512484/101‑1 (MQ=255)
         gggggCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..Cg                                                                                              >  1:813085/1‑101 (MQ=255)
         gggggCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..Cg                                                                                              >  1:4137189/1‑101 (MQ=255)
         gggggCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATAATGGC.C..Cg                                                                                              >  1:892680/1‑101 (MQ=255)
          ggggCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGc                                                                                             <  2:3999705/101‑1 (MQ=255)
           gggCGTGGATGGAAACACCGTCATCACCACGCAACCGCCCTGCACCGGCGGTAGGTCACCCCGCGGGGGCGTGGATCGAAACAAA.CCGATCATGGC.C..CGcc                                                                                            >  2:4809728/1‑101 (MQ=37)
           gggCGTGGATCGAAACCCCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGcc                                                                                            <  2:5794635/101‑1 (MQ=255)
            ggCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCt                                                                                           <  1:5395527/101‑1 (MQ=255)
              cGTGGATCCAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.a                                                                                        <  1:3395240/101‑1 (MQ=255)
                tGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.Agc                                                                                      >  1:4335459/1‑101 (MQ=255)
                 ggATCGAAACACCGTCATCCCCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.Agcg                                                                                     <  2:1696517/101‑1 (MQ=255)
                 ggATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.Agcg                                                                                     <  1:1385918/101‑1 (MQ=255)
                    tCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGcgccg                                                                                  >  2:1021237/1‑101 (MQ=255)
                     cGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGcgccg.c                                                                                <  2:5013727/101‑1 (MQ=255)
                      gAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCcg.cg                                                                               <  1:2379563/101‑1 (MQ=255)
                       aaaCACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCg.cgc                                                                              >  2:339084/1‑101 (MQ=255)
                         acacCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTg                                                                            <  1:1104576/101‑1 (MQ=255)
                          cacCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGa                                                                           <  2:263783/101‑1 (MQ=255)
                          cacCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGa                                                                           >  1:2859937/1‑101 (MQ=255)
                           cccGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGAc                                                                          <  1:3801799/100‑1 (MQ=255)
                           acCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGAc                                                                          >  2:2523140/1‑101 (MQ=255)
                            ccGTCATCACCACGCACCCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACg                                                                         <  2:3767670/101‑1 (MQ=255)
                                catcaCCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGc                                                                     <  1:769372/101‑1 (MQ=255)
                                 atcaCCACGCAACCGCCCTGCACCCCCGGTCCCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCt                                                                    <  1:1481508/101‑1 (MQ=255)
                                 atcaCCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCt                                                                    <  1:3822898/101‑1 (MQ=255)
                                  tcaCCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTc                                                                   <  1:2169665/101‑1 (MQ=255)
                                  tcaCCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTc                                                                   <  2:2056374/101‑1 (MQ=255)
                                   cactacGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTcc                                                                  >  2:4322439/1‑101 (MQ=255)
                                    accacGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTccc                                                                 >  2:290964/1‑101 (MQ=255)
                                    accacGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTccc                                                                 >  1:523486/1‑101 (MQ=255)
                                    accacGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTccc                                                                 <  1:3321922/101‑1 (MQ=255)
                                         gCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCcgc                                                            <  2:2355602/101‑1 (MQ=255)
                                               gCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGc                                                      <  1:474413/101‑1 (MQ=255)
                                               gCCCGGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGc                                                      <  1:3957730/101‑1 (MQ=255)
                                                 ccTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGt                                                    >  2:3666177/1‑101 (MQ=255)
                                                  cTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTg                                                   <  2:127250/101‑1 (MQ=255)
                                                   ggaacaggatgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTAT.GCTGGC.TCGA.C..CGAGTT.CACGGCG.GGCTATCGTCGCTCCCCCCGCGGGGGCGTGGa                                                 <  1:1010202/91‑1 (MQ=2)
                                                    gCACCACCGGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGa                                                 <  2:5194685/101‑1 (MQ=255)
                                                     caccacCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGAt                                                <  1:1679812/101‑1 (MQ=255)
                                                       tcgagcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACC.AATGCCGAAGC.CAAGGCCGC.AGCAGCG..GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATc                                               <  2:5621194/95‑1 (MQ=0)
                                                       ccaccGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCg                                              <  1:2222844/101‑1 (MQ=255)
                                                       ccaccGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCg                                              >  2:5314202/1‑101 (MQ=255)
                                                        accgcgTCGCTCCCCCCGCGGGGGCGGGGATCGAAACATCGATGATCAAATC.C..GCAATC.GGCCGGG.CTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                              <  1:1553335/96‑1 (MQ=0)
                                                         cccGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                            <  2:4689789/100‑1 (MQ=255)
                                                         accGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                            >  1:1779768/1‑101 (MQ=255)
                                                          agcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACC.AATGCCGAAGC.CAAGGCCGC.AGCAGCG..GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                            <  1:2670299/98‑1 (MQ=2)
                                                            cgTCGCCCCCCCCGCGGGGGCGTGGATCGAAACACC.AATGCCGAAGC.CAAGGCCGC.AGCAGCG..GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                          <  2:4198967/100‑1 (MQ=2)
                                                             gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACC.CCGAACCTGTCGG..CGGC.G.AACACCG.ATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                        >  2:1663266/1‑101 (MQ=2)
                                                             gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                        <  2:4417526/101‑1 (MQ=255)
                                                              tCGCTCCCCCCGCGGGGGCGTGGATCGAAACACC.AATGCCGAAGC.CAAGGCCGC.AGCAGCG..GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaa                                        <  1:4196028/101‑2 (MQ=2)
                                                              tCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATc                                       >  1:1904913/1‑101 (MQ=255)
                                                               cGCTCCCCCCGCGGGGGCGTGGATCGAAACTAT.GCTGGC.TCGA.C..CGAGTT.CACGGCG.GGCTATCGTCGCTCCCCCCGCGGGGGCGGGGGTCGaaaaggggc                                     >  1:4491237/1‑95 (MQ=2)
                                                               cGCTCCCCCCGCGGGGGCGTGGATCGAAACACC.CCGAACCTGTCGG..CGGC.G.AACACCG.ATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcg                                      <  2:137027/101‑2 (MQ=2)
                                                               cGCTCCCCCCGCGGGGGCGTGGATCGAAACACC.AATGCCGAAGC.CAAGGCCGC.AGCAGCG..GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaa                                       <  1:2082591/101‑3 (MQ=2)
                                                                gCTCCCCCCGCGGGGGCGTGGATCGAAACTAT.GCTGGC.TCGA.C..CGAGTT.CACGGCG.GGCTATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgggtcg                                    <  1:5590525/101‑7 (MQ=2)
                                                                gCTCCCCCCGCGGGGGCGTGGATCGAAACACC.AATGCCGAAGC.CAAGGCCGC.AGCAGCG..GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaa                                      <  1:2802386/101‑4 (MQ=2)
                                                                gCTCCCCCCGCGGGGGCGTGGATCGAAACACC.AATGCCGAAGC.CAAGGCCGC.AGCAGCG..GCTGCGGTCGCTCCCCCCGCGGAGGCGTGGATCGAAACaaaa                                      <  2:3157138/101‑4 (MQ=2)
                                                                gCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCcg                                     >  1:1189406/1‑101 (MQ=255)
                                                                 cTCCCCCCGCGGGGGCGTGGATCGAAACACC.CCGAACCTGTCGG..CGGC.G.AACACCGATGCCG.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgat                                    <  2:3322691/101‑4 (MQ=2)
                                                                 cTCCCCCCGCGGGGGCGTGGATCGAAACACC.AATGCCGAAGC.CAAGGCCGC.AGCAGCG..GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaaa                                     <  2:5069185/101‑5 (MQ=2)
                                                                 cTCCCCCCGCGGGGGCGTGGATCGAAACACC.AATGCCGAAGC.CAAGGCCGC.AGCAGCG..GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaaa                                     >  1:2312889/1‑97 (MQ=2)
                                                                  tCCCCCCGCGGGGGCGTGGATCGAAACTAT.GCTGGC.TCGA.C..CGAGTT.CACGGCG.GGCTATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgggtcgaa                                  <  2:1792929/101‑9 (MQ=2)
                                                                  tCCCCCCGCGGGGGCGTGGATCGAAACACC.CCGAACCT.GT.C..GGCCGCGAACACCG.ATACCGCGGCGCTCCCCCCGCGGGGGCGGGGATCGAAACATcgttc                                   >  2:2835365/1‑97 (MQ=2)
                                                                   ccccccGCGGGGGCGTGGATCGAAACTAT.GCTGGC.TCGA.C..CGAGTT.CACGGCG.GGCTATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgggtcgaag                                 <  2:3833086/101‑10 (MQ=2)
                                                                    cccccGCGGGGGCGTGGATCGAAACACC.CCGAACCTGTCGG..CGGC.G.AACACCG.ATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatgat                                 <  1:3565673/101‑7 (MQ=2)
                                                                     ccccGCGGGGGCGTTGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGAt                                <  2:3852165/101‑1 (MQ=255)
                                                                      cccGCGGGGGCGTGGATCGAAACACC.CCGAACCTGTCGG..CGGC.G.AACACCG.ATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatgatca                               <  1:5667739/101‑9 (MQ=0)
                                                                      cccGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATc                               <  2:550667/101‑1 (MQ=255)
                                                                       ccGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCg                              >  2:4692683/1‑101 (MQ=255)
                                                                       ccGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCg                              <  2:2187822/101‑1 (MQ=255)
                                                                           gggggCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGt                          >  2:5447002/1‑101 (MQ=255)
                                                                            ggggCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTg                         <  2:1941575/101‑1 (MQ=255)
                                                                               gCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcc                      <  2:3090147/101‑1 (MQ=255)
                                                                                cGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAccc                     <  1:3666405/101‑1 (MQ=255)
                                                                                 gTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccg                    <  2:4953732/101‑1 (MQ=255)
                                                                                  tGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgc                   >  1:4624084/1‑101 (MQ=255)
                                                                                  tGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgc                   >  2:4110810/1‑101 (MQ=255)
                                                                                   ggATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccccc                  >  2:723943/1‑101 (MQ=255)
                                                                                   agATCGAAACAAC.CGGATCATGCC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCTGATCGGCGTGAcccgcc                  <  2:3275296/100‑1 (MQ=255)
                                                                                     aTCAAAACACC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgcccc                <  2:1167460/101‑1 (MQ=255)
                                                                                      tCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgccccg               >  1:303692/1‑101 (MQ=255)
                                                                                         aaaCAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCcgccccgccc            >  1:371742/1‑101 (MQ=255)
                                                                                         aaaCAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCcgccccgccc            >  1:3850196/1‑101 (MQ=255)
                                                                                           acacc.CCGATCAGGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCAt          <  1:2523282/101‑1 (MQ=255)
                                                                                           acaac.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCAt          >  2:2408911/1‑101 (MQ=255)
                                                                                           acaac.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCAt          >  2:1541827/1‑101 (MQ=255)
                                                                                            caac.CCGATCATGGC.C..CGCCTG.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATc         >  1:3278616/1‑101 (MQ=255)
                                                                                             aac.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCg        >  2:4622675/1‑101 (MQ=255)
                                                                                              ac.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCgg       <  2:734497/101‑1 (MQ=255)
                                                                                                  cGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCg    >  2:2663564/1‑101 (MQ=255)
                                                                                                  cGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGAGCGTCGTGGCACGCCACGCCCCTCGGtag    >  2:888457/1‑99 (MQ=25)
                                                                                                    aTCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCt  <  2:3598319/101‑1 (MQ=255)
                                                                                                     |                                                                                                        
TCCCCCCGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCG.CGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCT  >  NC_014034/1397439‑1397638

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