Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,397,543 G→A 10.1% M475I (ATG→CAA RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,397,5430G→A10.1% 81.9 ‑2.4 ‑0.4326M475I (ATG→ATARCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (11/21):  ref base (143/136):  total (161/165)
Fisher's exact test strand distribution p-value = 9.23e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.71e-01

Only 100 of 362 total aligned reads displayed.

CCCGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCTCCCCCC  >  NC_014034/1397443‑1397644
                                                                                                     |                                                                                                         
 ccGCGGGGGCGTGGATCGAAACACCGTCATTACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATgg                                                                                                          <  1:5250688/101‑1 (MQ=255)
 ccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATgg                                                                                                          >  1:268016/1‑101 (MQ=255)
  cgcgGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGc                                                                                                         <  2:2945543/101‑1 (MQ=255)
   gcgGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGc.c                                                                                                       <  1:1256653/101‑1 (MQ=255)
   cttccgatctGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGc.c                                                                                                       <  2:5761822/92‑1 (MQ=255)
    cgGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGc.c..c                                                                                                    <  2:3362940/101‑1 (MQ=255)
     gggggCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..Cg                                                                                                   >  1:813085/1‑101 (MQ=255)
     gggggCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..Cg                                                                                                   >  2:556348/1‑101 (MQ=255)
            tGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.Agc                                                                                           >  2:3767413/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)
              gATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.Agcgc                                                                                         >  2:25384/1‑101 (MQ=255)
              gATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.Agcgc                                                                                         >  1:3710900/1‑101 (MQ=255)
               aTCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGcgcc                                                                                        >  1:1714358/1‑101 (MQ=255)
               aTCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGcgcc                                                                                        >  1:5405108/1‑101 (MQ=255)
                   aaaCACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCgcgc                                                                                    >  2:3554289/1‑101 (MQ=255)
                    gctcttccgatctCCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCt                                                                                   <  2:410508/92‑1 (MQ=38)
                    aaCACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCt                                                                                   <  1:236997/101‑1 (MQ=255)
                     acacCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTg                                                                                  <  1:1104576/101‑1 (MQ=255)
                         cGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGt                                                                              >  2:5025190/1‑101 (MQ=255)
                          gTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTc                                                                             >  2:2120496/1‑101 (MQ=255)
                           tcatcaCCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCg                                                                            >  2:103393/1‑101 (MQ=255)
                              ccaccacGCAACCCCCCGGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTc                                                                         <  2:5481547/100‑1 (MQ=255)
                                accacGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTccc                                                                       <  1:3321922/101‑1 (MQ=255)
                                    cGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCcg                                                                   <  1:380966/101‑1 (MQ=255)
                                     gCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCcgc                                                                  <  2:2355602/101‑1 (MQ=255)
                                       aaCCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCgg                                                                >  2:5380486/1‑101 (MQ=255)
                                       aaCCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCAGCgg                                                                <  2:4897496/101‑1 (MQ=255)
                                       aaCCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCCCCCCCCGCgg                                                                >  2:4718922/1‑101 (MQ=255)
                                        aCCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGGGCCGCGCTGGCGTCGCTCCCCCCGCggg                                                               >  2:1932552/1‑101 (MQ=255)
                                           gCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGc                                                            <  1:474413/101‑1 (MQ=255)
                                              cTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTg                                                         <  2:4245535/101‑1 (MQ=255)
                                               ggaacaggatgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTAT.GCTGGC.TCGA.C..CGAGTT.CACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGa                                                       <  1:1010202/91‑1 (MQ=2)
                                                gCACCACCGGTCGCTCCCCCCGGGGGGGCGTGGATCGAAACAAC.CCAAAAATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGa                                                       <  1:4810009/101‑1 (MQ=38)
                                                   cgctgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACC.CCGAACCTGTCGG..CGGC.G.AACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                    <  2:682534/95‑1 (MQ=2)
                                                   cgctgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACC.CCGAACCTGTCGG..CGGC.G.AACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                    <  1:766838/95‑1 (MQ=2)
                                                   ccaccGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                    <  1:2222844/101‑1 (MQ=255)
                                                   ccaccGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                    >  2:5314202/1‑101 (MQ=255)
                                                    aggatgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTAT.GCTGGC.TCGA.C..CGAGTT.CACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                  <  2:3980391/96‑1 (MQ=2)
                                                     accGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                  >  1:1779768/1‑101 (MQ=255)
                                                      tgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACC.CCGAACCTGTCGG..CGGC.G.AACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                 <  2:3320313/98‑1 (MQ=2)
                                                       cGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                >  1:516151/1‑101 (MQ=255)
                                                        cgTCGCCCCCCCCGCGGGGGCGTGGATCGAAACACC.AATGCCGAAGC.CAAGGCCGC.AGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                <  2:4198967/100‑1 (MQ=2)
                                                         gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACC.CCGAACCTGTCGG..CGGC.G.AACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                              >  2:1663266/1‑101 (MQ=2)
                                                         gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACC.AATGCCGAAGC.CAAGGCCGC.AGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                               >  1:3595309/1‑101 (MQ=2)
                                                         gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                              <  2:4417526/101‑1 (MQ=255)
                                                         gTCGCACCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATC.C..GCAATC.GGCCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcg                                               <  1:4623812/101‑2 (MQ=2)
                                                          tCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATc                                             >  1:1904913/1‑101 (MQ=255)
                                                           cGCTCCCCCCGCGGGGGCGTGGATCGAAACTAT.GCTGGC.TCGA.C..CGAGTT.CACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGGGGGTCGaaaaggggc                                           >  1:4491237/1‑95 (MQ=2)
                                                           cGCTCCCCCCGCGGGGGCGTGGATCGAAACACC.CCGAACCTGTCGG..CGGC.G.AACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcg                                            <  1:1586032/101‑2 (MQ=2)
                                                           cGCTCCCCCCGCGGGGGCGTGGATCGAAACACC.CCGAACCTGTCGG..CGGC.G.AACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcg                                            >  1:2507597/1‑100 (MQ=2)
                                                           cGCTCCCCCCGCGGGGGCGTGGATCGAAACACC.AATGCCGAAGC.CAAGGCCGC.AGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaa                                             <  1:2082591/101‑3 (MQ=2)
                                                           cGCTCCCCCCGCGGGGGCGTGGATCGAAACACC.AATGCCGAAGC.CAAGGCCGC.AGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaa                                             >  2:1224800/1‑99 (MQ=2)
                                                           cGCCCCCCCCGCGGGGGCGTGGATCGAAACACC.AATGCCGAAGC.CAAGGCCGC.AGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaa                                             <  2:5119525/101‑3 (MQ=0)
                                                            gCTCCCCCCGCGGGGGCGTGGATCGAAACTAT.GCTGGC.TCGA.C..CGAGTT.CACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgggtcg                                          <  1:2812404/101‑7 (MQ=2)
                                                            gCTCCCCCCGCGGGGGCGTGGATCGAAACTAT.GCTGGC.TCGA.C..CGAGTT.CACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgggtcg                                          <  2:2253908/101‑7 (MQ=2)
                                                            gCTCCCCCCGCGGGGGCGTGGATCGAAACACC.AATGCCGAAGC.CAAGGCCGC.AGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaa                                            <  2:4993897/101‑4 (MQ=2)
                                                            gCTCCCCCCGCGGGGGCGTGGATCGAAACACC.AATGCCGAAGC.CAAGGCCGC.AGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaa                                            <  2:3507299/101‑4 (MQ=2)
                                                            gCTCCCCCCGCGGGGGCGTGGATCGAAACACC.AATGCCGAAGC.CAAGGCCGC.AGCAGCG.GCTGCGGTCGCTCCCCCCGCGGAGGCGTGGATCGAAACaaaa                                            <  2:3157138/101‑4 (MQ=2)
                                                            gCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCcg                                           >  1:2574482/1‑101 (MQ=255)
                                                            gCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCcg                                           >  1:1189406/1‑101 (MQ=255)
                                                             cTCCCCCCGCGGGGGCGTGGATCGAAACACC.CCGAACCTGTCGG..CGGC.G.AACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgat                                          <  1:3055915/101‑4 (MQ=2)
                                                             cTCCCCCCGCGGGGGCGTGGATCGAAACACC.CCGAACCTGTCGG..CGGC.G.AACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgat                                          <  1:1499886/101‑4 (MQ=2)
                                                             cTCCCCCCGCGGGGGCGTGGATCGAAACACC.AATGCCGAAGC.CAAGGCCGC.AGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaaa                                           >  1:2312889/1‑97 (MQ=2)
                                                             cTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCcgc                                          <  1:4905232/101‑1 (MQ=255)
                                                              tCCCCCCGCGGGGGCGTGGATCGAAACTAT.GCTGGC.TCGA.C..CGAGTT.CACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgggtcgaa                                        <  2:1792929/101‑9 (MQ=2)
                                                              tCCCCCCGCGGGGGCGTGGATCGAAACACC.CCGAACCT.GT.C..GGCCGCGAACACCGATACCGCGGCGCTCCCCCCGCGGGGGCGGGGATCGAAACATcgttc                                         >  2:2835365/1‑97 (MQ=2)
                                                              cccccccGCGGGGGCGTGGATCGAATCAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCcgcg                                         <  1:4164889/100‑1 (MQ=255)
                                                              cccccccGCGGGGGCGTGGATCGAAACACC.CCGAACCTGTCGG..CGGC.G.AACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatg                                         <  2:4630797/100‑5 (MQ=2)
                                                               ccccccGCGGGGGCGTGGATCGAAACACC.AATGCCGAAGC.CAAGGCCGC.AGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaaagc                                         <  1:5570144/101‑7 (MQ=0)
                                                               ccccccGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCgg                                        >  1:3729148/1‑101 (MQ=255)
                                                                cccccGCGGGGGCGTGGATCGAAACACC.CCGAACCTGTCGG..CGGC.G.AACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatgat                                       >  2:469249/1‑95 (MQ=2)
                                                                cccccGCGGGGGCGTGGATCGAAACACC.CCGAACCTGTCGG..CGGC.G.AACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatgat                                       >  1:3913391/1‑95 (MQ=2)
                                                                cccccGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.CGCGCCGCTCTCACGTCGTTCCTCCCGCGGGTGCGTTTATCTAAActcacgaaga                                       >  2:310365/1‑92 (MQ=21)
                                                                  cccGCGGGGGCGTGGATCGAAACACC.CCGAACCTGTCGG..CGGC.G.AACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatgatca                                     <  1:5667739/101‑9 (MQ=0)
                                                                        ggggCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTg                               <  2:3821029/101‑1 (MQ=255)
                                                                         gggCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGa                              <  2:3645560/101‑1 (MQ=255)
                                                                         gggCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGCCCCCACCCCGCCGGGGGCGGGGATCGAAAGAACCGCGGATCGGCGCGa                              >  2:1353659/1‑101 (MQ=255)
                                                                          ggCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAc                             >  1:2903989/1‑101 (MQ=255)
                                                                            ggTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAccc                           <  1:3552292/100‑1 (MQ=255)
                                                                             gTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccg                          <  2:1398028/101‑1 (MQ=255)
                                                                              tGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgc                         >  1:4111341/1‑101 (MQ=255)
                                                                               ggATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCTCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgcc                        >  2:3849294/1‑101 (MQ=255)
                                                                                 aTCGAATCAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgcccc                      <  1:973383/101‑1 (MQ=255)
                                                                                 aTCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgcccc                      <  1:2158737/101‑1 (MQ=255)
                                                                                     aaaCAAC.CAGAGCATGGC.C..CCCCTC.AGCGCCGCACGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCcgccccgccc                  <  1:644530/101‑1 (MQ=37)
                                                                                       acaac.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCcgcccccccccc                >  1:4264022/1‑99 (MQ=255)
                                                                                       acaac.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCAt                >  2:1541827/1‑101 (MQ=255)
                                                                                        caac.CCGATCATGGC.C..CGCCTG.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATc               >  1:3278616/1‑101 (MQ=255)
                                                                                          ac.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCgg             >  2:2436555/1‑101 (MQ=255)
                                                                                          ac.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCgg             <  2:734497/101‑1 (MQ=255)
                                                                                              cGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGGGACCCGCCCCGCCCATCGGGCg          >  2:2743840/1‑101 (MQ=255)
                                                                                              cGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGAGCGTCGTGGCACGCCACGCCCCTCGGtag          >  2:888457/1‑99 (MQ=25)
                                                                                               gATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGGCGc         >  1:5062284/1‑101 (MQ=255)
                                                                                                aTCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCt        >  2:3248869/1‑101 (MQ=255)
                                                                                                aTCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCt        <  2:2649768/101‑1 (MQ=255)
                                                                                                    tGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCTcccc    >  2:4572133/1‑101 (MQ=255)
                                                                                                    tGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCTcccc    <  1:339128/101‑1 (MQ=255)
                                                                                                     ggC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCTccccc   <  1:3841488/101‑1 (MQ=255)
                                                                                                     ggC.C..CGCCTC.AGCGCCGCGCTGACGTCGCT.CCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCTcccccc  <  1:873876/101‑1 (MQ=255)
                                                                                                     |                                                                                                         
CCCGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTC.AGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCTCCCCCC  >  NC_014034/1397443‑1397644

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