Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,397,542 T→A 9.9% M475K (ATG→CAA)  RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,397,5420T→A 9.9% 89.5 ‑2.4 ‑0.8331M475K (ATG→AAG) RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (11/21):  ref base (148/144):  total (162/169)
Fisher's exact test strand distribution p-value = 9.45e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.42e-01

Only 100 of 362 total aligned reads displayed.

CCCCGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCTCCCC  >  NC_014034/1397442‑1397642
                                                                                                    |                                                                                                       
 cccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATg                                                                                                        <  2:1769706/101‑1 (MQ=255)
  ccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATgg                                                                                                       >  1:268016/1‑101 (MQ=255)
    gcgGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGc.c                                                                                                    <  1:1256653/101‑1 (MQ=255)
    gcgGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGc.c                                                                                                    <  1:776627/101‑1 (MQ=255)
    cttccgatctGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGc.c                                                                                                    <  2:5761822/92‑1 (MQ=255)
     cgGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGc.c..c                                                                                                 <  1:5223908/101‑1 (MQ=255)
     cgGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGc.c..c                                                                                                 <  1:4080749/101‑1 (MQ=255)
     cgGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGc.c..c                                                                                                 <  2:3362940/101‑1 (MQ=255)
      gggggCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..Cg                                                                                                >  2:556348/1‑101 (MQ=255)
         ggCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCt                                                                                             <  2:326086/101‑1 (MQ=255)
          gCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTc                                                                                            >  2:2178349/1‑101 (MQ=255)
            gTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAg                                                                                          >  1:4151734/1‑101 (MQ=255)
             tGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAgc                                                                                         >  2:3767413/1‑101 (MQ=255)
             tGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAgc                                                                                         <  1:3711068/101‑1 (MQ=255)
             tGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAgc                                                                                         <  1:5565396/101‑1 (MQ=255)
              ggATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAgcg                                                                                        <  1:1385918/101‑1 (MQ=255)
              ggATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAgcg                                                                                        >  2:2391789/1‑101 (MQ=255)
                aTCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGcgcc                                                                                      >  1:1714358/1‑101 (MQ=255)
                 tCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGcgccg                                                                                     >  2:1021237/1‑101 (MQ=255)
                    aaaCACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCgcgc                                                                                  >  2:3554289/1‑101 (MQ=255)
                     gctcttccgatctCCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCt                                                                                 <  2:410508/92‑1 (MQ=38)
                     aaCACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCt                                                                                 <  1:236997/101‑1 (MQ=255)
                      acacCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTg                                                                                <  1:1104576/101‑1 (MQ=255)
                       cacCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGa                                                                               <  2:263783/101‑1 (MQ=255)
                          cGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGt                                                                            >  2:5025190/1‑101 (MQ=255)
                           gTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACgcc                                                                           >  2:2835101/1‑99 (MQ=255)
                            tcatcaCCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCg                                                                          <  2:2475402/101‑1 (MQ=255)
                            tcatcaCCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCg                                                                          >  1:836117/1‑101 (MQ=255)
                              atcaCCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCt                                                                        <  1:3822898/101‑1 (MQ=255)
                               ccaccacGCAACCCCCCGGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTc                                                                       <  2:5481547/100‑1 (MQ=255)
                                 accacGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTccc                                                                     <  1:3321922/101‑1 (MQ=255)
                                    acGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTcccccc                                                                  <  1:5364708/101‑1 (MQ=255)
                                    acGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTcccccc                                                                  >  2:2313764/1‑101 (MQ=255)
                                     cGCAACCGCCCTGCACCCCCGGTCCCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCcg                                                                 <  2:1346350/101‑1 (MQ=255)
                                        aaCCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCgg                                                              >  2:5380486/1‑101 (MQ=255)
                                        aaCCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCAGCgg                                                              <  2:4897496/101‑1 (MQ=255)
                                        aaCCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCCCCCCCCGCgg                                                              >  2:4718922/1‑101 (MQ=255)
                                               cTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTg                                                       <  2:4245535/101‑1 (MQ=255)
                                                ggaacaggatgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGa                                                     <  1:1010202/91‑1 (MQ=2)
                                                    cgctgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                  <  2:682534/95‑1 (MQ=2)
                                                    ccaccGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                  >  2:5314202/1‑101 (MQ=255)
                                                    ccaccGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                  <  1:2222844/101‑1 (MQ=255)
                                                     ccccGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                 <  1:596845/99‑1 (MQ=255)
                                                     aggatgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                <  2:3980391/96‑1 (MQ=2)
                                                      accGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                >  1:1779768/1‑101 (MQ=255)
                                                       tgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                               <  2:3320313/98‑1 (MQ=2)
                                                         ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCGCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                             <  1:929925/101‑1 (MQ=255)
                                                         cgTCGCCCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                              <  2:4198967/100‑1 (MQ=2)
                                                          gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                            <  2:4417526/101‑1 (MQ=255)
                                                           tCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATc                                           <  1:5110389/101‑1 (MQ=2)
                                                           tCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaa                                            <  1:4196028/101‑2 (MQ=2)
                                                           tCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATc                                           >  1:1904913/1‑101 (MQ=255)
                                                            cGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcg                                          >  1:2507597/1‑100 (MQ=2)
                                                            cGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcg                                          <  1:1586032/101‑2 (MQ=2)
                                                            cGCTCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaa                                           >  2:1495369/1‑99 (MQ=2)
                                                            cGCTCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaa                                           <  1:2082591/101‑3 (MQ=2)
                                                            cGCTCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaa                                           >  2:1224800/1‑99 (MQ=2)
                                                            cGCCCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaa                                           <  2:5119525/101‑3 (MQ=0)
                                                             gCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgggtcg                                        <  1:4781901/101‑7 (MQ=2)
                                                             gCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgggtcg                                        <  1:2812404/101‑7 (MQ=2)
                                                             gCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgggtcg                                        <  1:5590525/101‑7 (MQ=2)
                                                             gCTCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaa                                          <  2:3507299/101‑4 (MQ=2)
                                                             gCTCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaa                                          <  2:4993897/101‑4 (MQ=2)
                                                             gCTCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaa                                          >  2:1553241/1‑98 (MQ=2)
                                                             gCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCcg                                         >  1:2574482/1‑101 (MQ=255)
                                                             gCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCcg                                         >  2:3074379/1‑101 (MQ=255)
                                                             gCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCcg                                         >  1:1189406/1‑101 (MQ=255)
                                                              cTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgat                                        <  1:1499886/101‑4 (MQ=2)
                                                              cTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgat                                        <  1:3055915/101‑4 (MQ=2)
                                                              cTCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaaa                                         >  1:2312889/1‑97 (MQ=2)
                                                               cccccccGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatg                                       <  2:4630797/100‑5 (MQ=2)
                                                                ccccccGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCgg                                      >  1:3733821/1‑101 (MQ=255)
                                                                ccccccGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGCCTTGGACCGAAACATCCGCgg                                      >  2:313427/1‑101 (MQ=255)
                                                                 cccccGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatgat                                     >  1:702922/1‑95 (MQ=2)
                                                                 cccccGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatgat                                     >  1:433869/1‑95 (MQ=2)
                                                                 cccccGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatgat                                     >  2:469249/1‑95 (MQ=2)
                                                                  ccccGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGAt                                    <  2:1969322/101‑1 (MQ=255)
                                                                       cgGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCg                               >  1:4616252/1‑101 (MQ=255)
                                                                         ggggCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTg                             <  2:3821029/101‑1 (MQ=255)
                                                                          gggCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGa                            <  2:3645560/101‑1 (MQ=255)
                                                                           ggCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAc                           >  1:2903989/1‑101 (MQ=255)
                                                                             ggTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAccc                         <  1:3552292/100‑1 (MQ=255)
                                                                             cGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAccc                         >  2:1898290/1‑101 (MQ=255)
                                                                             cGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAccc                         >  2:2660237/1‑101 (MQ=255)
                                                                              gTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccg                        <  2:1398028/101‑1 (MQ=255)
                                                                                ggATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCTCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgcc                      >  2:3849294/1‑101 (MQ=255)
                                                                                 gATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAAATCCGCGGATCGGCGTGAcccgccc                     >  1:4658313/1‑101 (MQ=255)
                                                                                  ctCGAAACAACCCGGTCATGGC.C..CGCCTCAGCGCCGCGCTGACGGCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgcccc                    <  2:2874002/100‑1 (MQ=255)
                                                                                  aTCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgcccc                    <  1:2158737/101‑1 (MQ=255)
                                                                                      aaaCAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCcgccccgccc                >  1:3850196/1‑101 (MQ=255)
                                                                                      aaaCAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCGCCGACCGGCGTGACCcgccccgccc                >  2:3159004/1‑101 (MQ=255)
                                                                                       aacaacCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCa               >  2:1680940/1‑101 (MQ=255)
                                                                                         caacCCGATCATGGC.C..CGCCTGAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATc             >  1:3278616/1‑101 (MQ=255)
                                                                                           acCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCgg           <  2:734497/101‑1 (MQ=255)
                                                                                              cGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGAGCGTCGTGGCACGCCACGCCCCTCGGtag        >  2:888457/1‑99 (MQ=25)
                                                                                                   aTGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCTccc   >  1:4835748/1‑101 (MQ=255)
                                                                                                   aTGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCTccc   <  1:5025471/101‑1 (MQ=255)
                                                                                                   aTGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCTccc   >  1:3611044/1‑101 (MQ=255)
                                                                                                    tGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCTcccc  >  2:4572133/1‑101 (MQ=255)
                                                                                                    tGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCTcccc  <  1:339128/101‑1 (MQ=255)
                                                                                                    |                                                                                                       
CCCCGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCTCCCC  >  NC_014034/1397442‑1397642

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