Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,494,048 G→A 8.3% intergenic (+96/‑146) RCAP_rcc02319 → / → RCAP_rcc02320 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,494,0480G→A 8.3% 20.4 ‑0.3 ‑0.1121intergenic (+96/‑146)RCAP_rcc02319/RCAP_rcc02320hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (6/4):  ref base (57/54):  total (63/58)
Fisher's exact test strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.56e-01

Only 100 of 142 total aligned reads displayed.

AGTAAGGGGGGGCGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AGCAGAAGGGGGGCAAGCG.CCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAC.GGGCGGGGGCGCGGGCATTAACCGAAATTGCA  >  NC_014034/2493948‑2494147
                                                                                                    |                                                                                                         
gtaaggggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgg                                                                                                           <  1:41580/96‑1 (MQ=255)
 taaggggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGAGATTTGCAGCACGATGAAACCGGATGCCTTTCATATTGCCCGCAATATCCCcgg                                                                                                          >  2:2879137/5‑101 (MQ=18)
 aagggggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCggg                                                                                                          <  1:3531791/97‑1 (MQ=255)
  aaggggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggg                                                                                                         <  1:2488527/98‑1 (MQ=255)
   aggggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCggggg                                                                                                        <  1:2500196/99‑1 (MQ=255)
   aggggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCggggg                                                                                                        >  2:1781926/3‑101 (MQ=255)
   aggggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCggggg                                                                                                        <  2:2205881/99‑1 (MQ=255)
    ggggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggggg                                                                                                       >  1:2602730/2‑101 (MQ=255)
     gggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGt                                                                                                      >  1:1407280/1‑101 (MQ=255)
      ggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCcgggggcg                                                                                                     >  1:4000270/1‑99 (MQ=255)
       gggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGa                                                                                                    >  1:2113492/1‑101 (MQ=255)
         gcgcccTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.t                                                                                                 <  2:915472/96‑1 (MQ=255)
          ggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.Tg                                                                                                >  1:1980829/1‑101 (MQ=255)
              cTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.a                                                                                         >  2:3680676/1‑101 (MQ=255)
               tccccccCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aa                                                                                        <  1:1852439/99‑2 (MQ=255)
               tGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..G.aa                                                                                        >  1:3321008/1‑100 (MQ=255)
                gCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aag                                                                                       >  1:301015/1‑99 (MQ=255)
                gCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aag                                                                                       <  2:4031543/101‑3 (MQ=255)
                 cccccG.TCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATCCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGGC.TGGC..G.aagcc                                                                                     >  2:2595949/1‑97 (MQ=18)
                 cccccCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagc                                                                                      <  2:1014895/101‑4 (MQ=255)
                    cccgtCCGCAATCTGACTCCCCCCGTGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccag                                                                                   <  1:713144/97‑7 (MQ=18)
                    cccgtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccag                                                                                   <  1:2631089/97‑7 (MQ=255)
                    cccgtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccag                                                                                   <  1:1478641/97‑7 (MQ=255)
                       gtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccagagg                                                                                >  2:1266907/2‑92 (MQ=255)
                       gtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccagagg                                                                                <  2:4026053/100‑10 (MQ=255)
                             caGCGGTCTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGCa                                                                           <  2:1690685/100‑1 (MQ=18)
                             aaGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGCa                                                                           <  1:1217599/101‑1 (MQ=18)
                             aaGCGCACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGCa                                                                           <  1:2125088/101‑1 (MQ=18)
                                 ggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGCagcgc                                                                       >  2:4233939/1‑97 (MQ=18)
                                  gACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCG.cc                                                                    >  2:3495375/1‑101 (MQ=17)
                                   aCTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCG.ccc                                                                   >  2:181846/1‑101 (MQ=17)
                                    cTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCG.cccc                                                                  >  2:231394/1‑101 (MQ=18)
                                    cTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCG.cccc                                                                  <  1:2553565/101‑1 (MQ=17)
                                     tCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCG.ccccc                                                                 >  2:321368/1‑101 (MQ=18)
                                       cccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCG.cccccct                                                               <  2:640139/101‑3 (MQ=17)
                                       cccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCG.cccccct                                                               <  1:4113550/101‑3 (MQ=17)
                                       cccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCG.cccccct                                                               >  2:3155734/1‑99 (MQ=17)
                                       cccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCG.cccccct                                                               <  2:936538/101‑3 (MQ=17)
                                         cccGGGATAGTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCG.cccccctgc                                                             <  2:2113354/101‑5 (MQ=17)
                                           cGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGA..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGccc                                                            <  1:3748899/101‑1 (MQ=17)
                                            gggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGcccc                                                           <  2:3255155/101‑1 (MQ=17)
                                            gggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGcccc                                                           <  2:3008822/101‑1 (MQ=17)
                                             ggATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AGGCCAAGAGGGGGCAGCG.CCCCCTGCCCcg                                                          >  1:176553/1‑101 (MQ=2)
                                             ggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCcg                                                          <  2:1524554/101‑1 (MQ=17)
                                              gATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AGGCCAAGAGGGGGCAGCG.CCCCCTGCCCcgc                                                         >  1:286327/1‑100 (MQ=2)
                                              gATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AGGCCAAGAGGGGGCAGCG.CCCCCTGCCCcgc                                                         >  1:1383546/1‑100 (MQ=0)
                                              gATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCcgg                                                         >  2:3035939/1‑101 (MQ=18)
                                               atatTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AGGCCAAGAGGGGGCAGCG.CCCCCTGCCCcgcc                                                        >  2:2911175/1‑99 (MQ=2)
                                               atatTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AGGCCAAGAGGGGGCAGCG.CCCCCTGCCCcgcc                                                        >  1:1063102/1‑99 (MQ=2)
                                               atatTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..GAAGCCAGAGGGGGGC.CGCGCCCCCCCGCCCcggc                                                        >  1:1983187/1‑101 (MQ=14)
                                                 atTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCcggcgg                                                      >  1:3984679/1‑101 (MQ=17)
                                                   ttGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AGGCCAAGAGGGGGCAGCG.CCCCCTGCCCcgccggat                                                    <  2:3744366/101‑3 (MQ=2)
                                                   ttGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AGGCCAAGAGGGGGCAGCG.CCCCCTGCCCcgccggat                                                    <  1:582313/101‑3 (MQ=2)
                                                    tGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTcc                                                   >  1:3235319/1‑101 (MQ=17)
                                                     gcagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGCATCCC.CCCCCTGCCCCGGCGGTCCg                                                  <  2:3138149/101‑1 (MQ=18)
                                                     gcagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCg                                                  <  1:3695775/101‑1 (MQ=17)
                                                       agcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTg                                                <  2:1717607/101‑1 (MQ=17)
                                                       agcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATAGCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTg                                                <  1:181884/101‑1 (MQ=17)
                                                         caAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGGCCGTggg                                              >  1:143678/1‑100 (MQ=17)
                                                           aGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTggcgg                                            >  2:625069/1‑101 (MQ=18)
                                                           aGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTggcgg                                            >  1:4224779/1‑101 (MQ=17)
                                                             aTGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGCAGCTC.CCCCCTGCCCCGGCGGTCCGTGgcggcg                                          >  1:2845295/1‑101 (MQ=18)
                                                             aTGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGgcggcg                                          <  1:3195103/101‑1 (MQ=17)
                                                               gAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGGGGCGcc                                        >  2:1265273/1‑101 (MQ=17)
                                                               gAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGCGAA.TAGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGcc                                        >  1:335947/1‑101 (MQ=17)
                                                                aaaCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCt                                       <  1:3680930/101‑1 (MQ=17)
                                                                aaaCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGGCCGGGGCGGCGGCt                                       >  1:2207444/1‑101 (MQ=9)
                                                                  aCCGGATGCCTTCCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCAGGCCGTGGGGGGGCCTTg                                     >  1:932663/1‑101 (MQ=9)
                                                                   ccGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGa                                    >  1:2841566/1‑101 (MQ=17)
                                                                    cGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAc                                   <  1:882183/101‑1 (MQ=17)
                                                                    cGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..GAAGCCAGAGGGGGGCAGGGC.CCCCCTGCCCCGGCGGGCCGCGGCGGCGCCCTGAc                                   >  1:3622924/1‑101 (MQ=17)
                                                                      gATGTCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCGGAGGGGGGCAGGGC.CCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAc.gg                                <  2:1353657/101‑1 (MQ=18)
                                                                      gATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAc.gg                                <  2:2298375/101‑1 (MQ=17)
                                                                       aTGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAc.ggg                               <  1:4122895/101‑1 (MQ=17)
                                                                        tGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAc.gggc                              <  2:1784540/101‑1 (MQ=17)
                                                                        tGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAc.gggc                              >  2:517161/1‑101 (MQ=17)
                                                                         gCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAc.gggcg                             <  2:2046365/101‑1 (MQ=17)
                                                                          ccTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAc.gggcgg                            <  1:1165216/101‑1 (MQ=17)
                                                                           cTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAc.gggcggg                           >  2:3396601/1‑101 (MQ=18)
                                                                            tttCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAC.GGGCgggg                          >  2:3468148/1‑101 (MQ=18)
                                                                              tCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAC.GGGCGGGGgc                        >  2:1319747/1‑101 (MQ=18)
                                                                                 tCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGGC.GGGCGGGGGcgcg                     >  2:1136849/1‑101 (MQ=18)
                                                                                  cTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAC.GGGCGGGGGCGCgg                    >  2:2713873/1‑101 (MQ=18)
                                                                                    tGCCCGAAATATGCGGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAC.GGGCGGGGGCGCGGGc                  <  2:2602539/101‑1 (MQ=18)
                                                                                    tGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAC.GGGCGGGGGCGCGGGc                  <  1:734103/101‑1 (MQ=18)
                                                                                      cccGAAATATCCCGGGGGGTGAA.TGAC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAC.GGGCGGGGGCGCGGGCAt                >  2:65134/1‑101 (MQ=18)
                                                                                        cGAAATTTCCCGGGGGGGTAA.TGGC..G.AAGCCAGGGGGGGGCAGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAC.GGGCGGGGGCGCGGGCACTa              <  1:2414389/101‑1 (MQ=16)
                                                                                        cGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAC.GGGCGGGGGCGCGGGCATTa              <  2:262806/101‑1 (MQ=18)
                                                                                         gAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAC.GGGCGGGGGCGCGGGCATTaa             <  2:3159051/101‑1 (MQ=18)
                                                                                          aagtgtCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAC.GGGCGGGGGCGCGGGCATTAAc            <  1:3190945/96‑1 (MQ=18)
                                                                                             tatCCGGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAC.GGGCGGGGGCGCGGGCATTAACCGa         <  2:1519688/101‑1 (MQ=21)
                                                                                              atCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGG.GAGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTCGAC.GGGGGGGGGCGCGGGCAGTAACCCaa        >  1:640583/1‑101 (MQ=12)
                                                                                               tCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAC.GGGCGGGGGCGCGGGCATTAACCGaaa       >  2:2713553/1‑101 (MQ=21)
                                                                                               tCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAC.GGGCGGGGGCGCGGGCATTAACCGaaa       >  1:1570566/1‑101 (MQ=21)
                                                                                                  cGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAC.GGGCGGGGGCGCGGGCATTAACCGAAATTg    >  2:1127724/1‑101 (MQ=25)
                                                                                                   ggggggTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAC.GGGCGGGGGCGCGGGCATTAACCGAAATTGc   >  1:1912598/1‑101 (MQ=25)
                                                                                                   ggggggTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAC.GGGCGGGGGCGCGGGCATTAACCGAAATTGc   <  1:4170531/101‑1 (MQ=25)
                                                                                                   ggggggTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGAGGTCCGAGGGGGGGCCTTGACTGGGCGGGGG.GCGGGCGTTAACGGaagtggc   >  1:1608576/1‑96 (MQ=2)
                                                                                                   ggggggCGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAC.GGGCGGGGGCGCGGGCCTTAACCGAAATTGc   >  2:1282704/1‑101 (MQ=21)
                                                                                                    gggggTGAA.AGGC..GAAGCCAGAGGGGGGCAGCCC.CCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAC.GGGCGGGGGCGCGGGCATTTACCGAAATTGCa  <  1:2285623/101‑1 (MQ=25)
                                                                                                    |                                                                                                         
AGTAAGGGGGGGCGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AGCAGAAGGGGGGCAAGCG.CCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAC.GGGCGGGGGCGCGGGCATTAACCGAAATTGCA  >  NC_014034/2493948‑2494147

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 23 ≤  ATCG  < 32 ≤  ATCG  < 40 ≤  ATCG