Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,494,067 A→G 91.5% intergenic (+115/‑127) RCAP_rcc02319 → / → RCAP_rcc02320 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,494,0670A→G91.5% 22.8 ‑0.6 ‑0.0139intergenic (+115/‑127)RCAP_rcc02319/RCAP_rcc02320hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (60/60):  ref base (7/4):  total (73/66)
Fisher's exact test strand distribution p-value = 5.32e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.99e-01

Only 100 of 149 total aligned reads displayed.

GCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.A..G.CAG..A.AGGGGGGCA.AGCG.C.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCAACAGA  >  NC_014034/2493975‑2494152
                                                                                                     |                                                                                         
gCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AAGG.CAG..A.Gggggg                                                                                    >  1:1601862/1‑101 (MQ=18)
  aaGCGGGCTCCCGCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGCa                                                                                 <  1:2260694/101‑1 (MQ=18)
    gCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGCagc                                                                               >  1:2176995/1‑99 (MQ=18)
    gCGGACTCCCCCCGGGAAATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGCaga                                                                               >  1:4734439/1‑99 (MQ=18)
      ggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAC.TGGC..GAA..G.CCA..G.AGGGGGGCagcgc                                                                             >  1:4638215/1‑97 (MQ=17)
      ggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGCagcgc                                                                             >  1:4676075/1‑97 (MQ=18)
       gACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCG.c.c                                                                        >  1:2725543/1‑101 (MQ=17)
        aCTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCG.c.cc                                                                       >  2:2414556/1‑101 (MQ=18)
        aCTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCG.c.cc                                                                       >  1:283443/1‑101 (MQ=18)
        aCTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCG.c.cc                                                                       >  2:3025943/1‑101 (MQ=18)
         cTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..GCCAG....AGGGGGGC..AGCG.c.ccc                                                                      <  1:402027/101‑1 (MQ=17)
         cTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..GCCAG....AGGGGGGC..AGCG.c.ccc                                                                      <  2:2954777/101‑1 (MQ=17)
         cTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCG.c.ccc                                                                      <  2:196547/101‑1 (MQ=17)
         cTCCCCCCGGGATAGTTGCAGCAAGATGAATCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCG.c.ccc                                                                      <  2:43056/101‑1 (MQ=17)
          tCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCG.c.cccc                                                                     >  1:1565968/1‑101 (MQ=18)
          tCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCG.c.cccc                                                                     >  2:128706/1‑101 (MQ=17)
            cccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCG.c.ccccct                                                                   <  2:848858/101‑3 (MQ=17)
               ccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCG.c.ccccctgcc                                                                <  2:4018385/101‑6 (MQ=17)
                cGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..T.AGGGGGGC..AGCGCC.CCCCTGccc                                                                >  1:2925682/1‑101 (MQ=18)
                 gggATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGC..C.A..A.AGGCCA.AGAGGGGGC.AGCG.C.CCCCTGcccc                                                               <  1:4206111/101‑1 (MQ=0)
                 gggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGcccc                                                               >  1:3588006/1‑101 (MQ=17)
                 gggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGcccc                                                               >  1:3846964/1‑101 (MQ=17)
                  ggATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.A..G.GCC..A.AGAGGGGGC.AGCG.C.CCCCTGCCCcg                                                              >  2:3750535/1‑101 (MQ=2)
                   gATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.A..G.GCC..A.AGAGGGGGC.AGCG.C.CCCCTGCCCcgc                                                             >  1:4404457/1‑100 (MQ=2)
                    atatTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.A..G.GCC..A.AGAGGGGGC.AGCG.C.CCCCTGCCCcgcc                                                            >  1:2418385/1‑99 (MQ=2)
                    atatTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCcggg                                                            >  1:3337824/1‑100 (MQ=17)
                     tatTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.A..G.GCC..A.AGAGGGGGC.AGCG.C.CCCCTGCCCcgccg                                                           >  2:1139753/1‑101 (MQ=2)
                      atTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCG.CACCCCTGCCCcggcgg                                                          >  2:2029738/1‑101 (MQ=18)
                       tttGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCcggcggc                                                         >  2:977702/1‑100 (MQ=17)
                         tGCAGCAATATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTcc                                                       >  1:3303403/1‑101 (MQ=17)
                         tGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTcc                                                       >  2:263673/1‑101 (MQ=17)
                          gcagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCg                                                      >  1:4463377/1‑101 (MQ=17)
                           cagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGt                                                     <  1:4346247/101‑1 (MQ=17)
                            agcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTg                                                    <  1:998655/101‑1 (MQ=17)
                            agcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTg                                                    <  1:2414691/101‑1 (MQ=17)
                              caAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGG.A.AGCGCC.CCCCTGCCCCGGCGGTCCGTggc                                                  >  1:345412/1‑101 (MQ=18)
                                aGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTggcgg                                                >  2:1161814/1‑101 (MQ=18)
                                  aTGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGgcggcg                                              >  2:1882578/1‑101 (MQ=18)
                                    gAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGcc                                            >  1:2626861/1‑101 (MQ=17)
                                     aaaCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGGA.TGGC..GAA..G.CCA..G.AGGGGGGCA.GGGC.C.CCCCTGCCCCCGCGGTCAGAGGCGGCGccg                                           >  1:4463056/1‑100 (MQ=12)
                                     aaaCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGAGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCt                                           >  1:2831091/1‑101 (MQ=18)
                                       aCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTg                                         >  1:812992/1‑101 (MQ=17)
                                         ccgATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAc                                       <  2:1774723/99‑1 (MQ=17)
                                         cGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAc                                       >  1:769470/1‑101 (MQ=18)
                                         cGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAc                                       <  2:4640867/101‑1 (MQ=17)
                                          ggATGCCTTTCATCTTGCCCGAAATGTCCCGGGGGTTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcg                                      <  1:4552032/101‑1 (MQ=17)
                                            aTGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGCAGAGC..C.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGatcgg                                    >  2:857681/1‑97 (MQ=17)
                                            aTGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcggg                                    <  2:1430771/101‑1 (MQ=17)
                                            aTGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcggg                                    <  1:1172984/101‑1 (MQ=17)
                                             tgctttacgtCTTGCCCGACATATTCCGGGGGGCAAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCCTGAcgggg                                   >  2:142639/10‑100 (MQ=9)
                                             tGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcgggc                                   <  1:4037518/101‑1 (MQ=17)
                                               ccTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcgggcgg                                 <  1:1254947/101‑1 (MQ=17)
                                                cTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcgggcggg                                >  1:3228256/1‑101 (MQ=18)
                                                cTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..GAA..G.CCA..G.AGGGGGGCA.GCAC.C.CCCCTGCCCCGGCGGTCCGTGGCGGCGGCGTGAcgggcggg                                >  2:4502618/1‑101 (MQ=17)
                                                 tttCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTTCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCgggg                               >  1:2587620/1‑101 (MQ=17)
                                                 tttCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCgggg                               <  2:944117/101‑1 (MQ=17)
                                                 tttCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGCG..A.A..A.CCA..A.AGGGGGGCC.AGGGCC.CCCCGGCCCCGGCGG.CCGGGGCCGGGCCCTTGCCGGCggggg                              >  1:1221501/1‑101 (MQ=9)
                                                 tttCATCTTCCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCgggg                               <  2:4386003/101‑1 (MQ=17)
                                                  ttCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCggggg                              >  2:1718966/1‑101 (MQ=18)
                                                   tCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCgggggg                             >  2:3581/1‑100 (MQ=18)
                                                   tCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGgc                             <  2:4272800/101‑1 (MQ=18)
                                                   tCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGgc                             <  2:2564012/101‑1 (MQ=18)
                                                   tCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGAGGCGGCGCCTTGACGGGCGGGGgc                             >  2:2012970/1‑101 (MQ=17)
                                                     aTCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCGCGGCGGTCCGTGGCGGCGCCTTGACGGGCgggggggg                           >  1:4009401/1‑98 (MQ=17)
                                                      tCTTTCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGcgcg                          <  2:2176871/101‑1 (MQ=18)
                                                      tCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGcgcg                          >  1:46093/1‑101 (MQ=18)
                                                       cTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCgg                         >  2:1633663/1‑101 (MQ=18)
                                                         tGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGGGCCACGAAGGGCCGGGGCGCgggg                       >  1:835590/1‑100 (MQ=9)
                                                         tGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGc                       <  1:617765/101‑1 (MQ=18)
                                                          gCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCa                      >  2:4636090/1‑101 (MQ=18)
                                                           cccGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCAt                     <  1:172579/101‑1 (MQ=18)
                                                            ccGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGGGCCTTGACGGGCGGGGGCGCGGGCAtt                    >  1:2628535/1‑101 (MQ=18)
                                                            ccGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCAtt                    <  2:3593736/101‑1 (MQ=18)
                                                            ccGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCAtt                    <  2:2304365/101‑1 (MQ=18)
                                                              gAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTaa                  <  1:98413/101‑1 (MQ=18)
                                                              gAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTaa                  <  1:56142/101‑1 (MQ=18)
                                                              gAAATATCCCGGGGGGCGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGGCCCTTCACGGGCGGGGGGGCGGGCATTaa                  >  2:1959993/1‑101 (MQ=12)
                                                               aaaTATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAAc                 <  2:4734183/101‑1 (MQ=18)
                                                               aaaTATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAAc                 <  1:459589/101‑1 (MQ=18)
                                                                aaTATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAAcc                >  1:1127494/1‑101 (MQ=21)
                                                                aaTATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAAcc                <  1:2740830/101‑1 (MQ=21)
                                                                 atatCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCg               >  1:1296638/1‑101 (MQ=21)
                                                                 atatCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCg               <  2:4398025/101‑1 (MQ=21)
                                                                  taccccGGGGGGTGAA.TGGC..G.A..A.GCC..A.GAGGGGGCC.AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGa              <  2:1329142/98‑1 (MQ=21)
                                                                   atCCGGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGaa             <  2:81842/101‑1 (MQ=21)
                                                                   atCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGaa             <  1:2278845/101‑1 (MQ=21)
                                                                   atCCCGGGGGGTGAA.TGGC..GAA..G.CCA..G.AGGGGGG.G.AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGGCGCTTGACGGGCGGGGGCGCGGGCATTAACCGaa             >  1:538066/1‑101 (MQ=16)
                                                                   atCCCGGGGGGGGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGGGCGGGCATTAACCGaa             >  2:2821257/1‑101 (MQ=21)
                                                                        ggggggTGAA.TGGC..GAA..G.CCA..G.AGGGCGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGc        >  1:4254827/1‑101 (MQ=25)
                                                                         gggggTGAA.TGGC..GAA..GCCAG....AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCa       >  1:316332/1‑101 (MQ=25)
                                                                         gggggTGAA.TGGC..G.AG.G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCa       <  2:417052/101‑1 (MQ=25)
                                                                           gggTGAA.TGGC..GAA..GCCAG....AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCAAc     >  1:3130257/1‑101 (MQ=25)
                                                                           gggTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCAAc     <  1:749122/101‑1 (MQ=25)
                                                                           gggTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCAAc     <  1:2590566/101‑1 (MQ=25)
                                                                           gggTGAA.TGGC..G.A..A.GAG..AGAGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCAAc     <  2:71212/101‑1 (MQ=37)
                                                                            ggTGAA.TGGC..GAA..GCCAG....AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGGGCGGGCATTAACCGAAATTGCaaaa    >  1:16817/1‑99 (MQ=21)
                                                                            ggTGAA.TGGC..GAA..GCCAG....AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCAACa    >  2:377256/1‑101 (MQ=25)
                                                                             gTGAA.TGGC..GAA..GCCAG....AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCAACAg   >  1:3642685/1‑101 (MQ=25)
                                                                             gTGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCAACAg   <  1:3047534/101‑1 (MQ=25)
                                                                              tGAA.TGGC..GAA..G.CCA..G.AGGGGGGC..AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCAACAGa  <  1:2132218/101‑1 (MQ=25)
                                                                                                     |                                                                                         
GCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.A..G.CAG..A.AGGGGGGCA.AGCG.C.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCAACAGA  >  NC_014034/2493975‑2494152

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 24 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG