Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,494,023 C→A 7.0% intergenic (+71/‑171) RCAP_rcc02319 → / → RCAP_rcc02320 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,494,0230C→A 7.0% 20.4 ‑0.7 ‑0.7143intergenic (+71/‑171)RCAP_rcc02319/RCAP_rcc02320hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (6/4):  ref base (61/72):  total (67/76)
Fisher's exact test strand distribution p-value = 5.16e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.02e-01

Only 100 of 153 total aligned reads displayed.

CGGGCGCGCGGCACGACAGGCCGAGTAA..GGG.GGGGCGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AGCAG..AAGGGGGGCAAGCG.CCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGG  >  NC_014034/2493925‑2494122
                                                                                                     |                                                                                                          
cgggcgCGCGGCACGACAGGCCGAGTAAG.GGG.GGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCt                                                                                                           >  2:1314684/1‑101 (MQ=21)
 gggcgCGCGGCACGACAGGCCGAGTAAG.GGG.GGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCtt                                                                                                          >  2:1217505/1‑101 (MQ=21)
 gggcgCGCGGCACGACAGGCCGAGTAAG.GGG.GGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCtt                                                                                                          <  2:1517101/101‑1 (MQ=21)
 gggcgCGCGGCACGACAGGCCGAGTAAG.GGG.GGGGCGCTGCCCCCCGTCCGCAAGCCGGCTCCCCCCGGGATATTTGCAGCAAGAGGAAACCGGATGCCtt                                                                                                          >  2:1852319/1‑101 (MQ=21)
  ggcgCGCGGCACGACAGGCCGAGTAAG.GGG.GGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCttt                                                                                                         <  2:3798128/101‑1 (MQ=21)
   gcgcgcGGCACGACAGGCCGAGTAAG.GGG.GGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTc                                                                                                        >  1:2611616/1‑101 (MQ=21)
   gcgcgcGGCACGACAGGCCGAGTAAG.GGG.GGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTc                                                                                                        <  2:1834049/101‑1 (MQ=21)
    ggcgcgGCACGACAGGCCGAGTAAG.GGG.GGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCa                                                                                                       <  2:1565830/100‑1 (MQ=18)
    cgcgcgGCACGACAGGCCGAGTAAGGGGG.GGGGCGCTGCCCCCCGCCCGCAAGCGGACTCCCCCCGGGAATTTTGCAGCAAAATGAAACCGGAGGCCTTTc                                                                                                        >  1:1059554/1‑101 (MQ=18)
    cgcgcgGCACGACAGGCCGAGTAAG.GGG.GGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCa                                                                                                       <  1:3896555/101‑1 (MQ=21)
    cgcgcgGCACGACAGGCCGAGTAAG.GGG.GGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCa                                                                                                       <  1:337161/101‑1 (MQ=21)
    cccgcgGCACGACAGGCCGAGTAAG.GGG.GGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCa                                                                                                       <  1:1388558/99‑1 (MQ=18)
     gcgcgGCACGACAGTCCGAGTAA..GGGCGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCAt                                                                                                      <  2:2808832/101‑1 (MQ=18)
     gcgcgGCACGACAGGCCGAGTAAG.GGG.GGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGACTTTCAt                                                                                                      >  2:734027/1‑101 (MQ=18)
       gcgGCACGACAGGCCGAGTAAG.GGG.GGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCt                                                                                                    <  2:1549294/101‑1 (MQ=21)
        cggccggACAGGCCGAGTAAG.GGG.GGGGCCCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCtt                                                                                                   <  1:2413190/95‑1 (MQ=18)
        cgGCACGACAGGCCGAGTAAG.GGG.GGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCtt                                                                                                   <  1:134877/101‑1 (MQ=21)
          gCACGACAGGCCGAGTAAG.GGG.GGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGc                                                                                                 <  1:3274269/101‑1 (MQ=21)
           cACGACAGGCCGAGTAAG.GGG.GGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGcc                                                                                                >  1:1014986/1‑101 (MQ=21)
            aCGACAGGCCGAGTAAG.GGG.GGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGccc                                                                                               <  2:2354520/101‑1 (MQ=21)
             cGACAGGCCGAGTAAG.GGG.GGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGAGATTTGCAGCAAGATGAAACCGGACGCCTTTCATCTTGCCCg                                                                                              >  1:1524658/1‑101 (MQ=18)
                cAGGCCGAGTAAG.GGG.GGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGaaa                                                                                           <  2:3183986/101‑1 (MQ=21)
                cAGGCCGAGTAAG.GGG.GGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGaaa                                                                                           >  2:3814169/1‑101 (MQ=21)
                cAGGCCGAGTAAG.GGG.GGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCCGGGCATTTGCAGAAAGATGAAAACGGATGCCTTTTATCTTGTCCGaca                                                                                           >  1:833567/1‑99 (MQ=17)
                     cGAGTAAG.GGG.GGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATcc                                                                                      <  1:3496873/101‑1 (MQ=21)
                     cGAGTAAG.GGG.GGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATcc                                                                                      <  1:746828/101‑1 (MQ=21)
                      gAGTAAG.GGG.GGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATccc                                                                                     >  2:3928627/1‑101 (MQ=18)
                       aGTAAG.GGG.GGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCg                                                                                    <  1:2718827/101‑1 (MQ=18)
                       aGTAAG.GGG.GGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCg                                                                                    <  1:1137456/101‑1 (MQ=18)
                        aggaagggg.ggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCg                                                                                    <  2:177793/95‑1 (MQ=18)
                         gtaagggg.ggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgg                                                                                   <  1:41580/96‑1 (MQ=255)
                          taagggg.ggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCggg                                                                                  >  1:1814213/5‑101 (MQ=255)
                          taagggg.ggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGAGATTTGCAGCACGATGAAACCGGATGCCTTTCATATTGCCCGCAATATCCCcgg                                                                                  >  2:2879137/5‑101 (MQ=18)
                          aaggggg.ggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCggg                                                                                  <  1:3531791/97‑1 (MQ=255)
                           aagggg.ggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggg                                                                                 <  1:2488527/98‑1 (MQ=255)
                            agggg.ggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCggggg                                                                                <  2:2205881/99‑1 (MQ=255)
                             gggg.ggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggggg                                                                               >  1:2602730/2‑101 (MQ=255)
                              ggg.ggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGCAATATCCCggggggg                                                                              >  2:1224466/1‑100 (MQ=255)
                              ggg.ggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCggggggg                                                                              >  1:2476508/1‑100 (MQ=255)
                              ggg.ggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGt                                                                              >  1:1407280/1‑101 (MQ=255)
                              ggg.ggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGt                                                                              <  2:3679823/101‑1 (MQ=255)
                                g.ggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGa                                                                            >  1:2113492/1‑101 (MQ=255)
                                  ggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGaa                                                                           <  1:3964032/101‑1 (MQ=255)
                                    ggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.Tg                                                                        >  1:1980829/1‑101 (MQ=255)
                                      cgcTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGc                                                                      <  1:4123983/101‑1 (MQ=255)
                                        cTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.a                                                                 >  2:3680676/1‑101 (MQ=255)
                                         tccccccCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aa                                                                <  1:1852439/99‑2 (MQ=255)
                                         tGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGCA.TGGC..G.aa                                                                >  1:123812/1‑100 (MQ=255)
                                          gCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aag                                                               <  2:4031543/101‑3 (MQ=255)
                                          gCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aag                                                               >  1:301015/1‑99 (MQ=255)
                                           cccccG.TCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagcc                                                             <  1:157847/101‑5 (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:1478641/97‑7 (MQ=255)
                                              cccgtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccag                                                           <  1:2631089/97‑7 (MQ=255)
                                                cgtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccagag                                                         <  1:4031167/99‑9 (MQ=255)
                                                 gtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccagagg                                                        >  2:1266907/2‑92 (MQ=255)
                                                 gtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccagagg                                                        <  2:4026053/100‑10 (MQ=255)
                                                        aGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCA..GAGGGGGGCag                                                <  2:1980703/101‑2 (MQ=18)
                                                           ggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCA..GAGGGGGGCagcgc                                             >  2:4233939/1‑97 (MQ=18)
                                                           ggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAGGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCA..GAGGGGGGCaccgc                                             >  1:1861491/1‑97 (MQ=17)
                                                              cTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCA..GAGGGGGGC.AGCG.cccc                                        <  1:2553565/101‑1 (MQ=17)
                                                              cTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCA..GAGGGGGGC.AGCG.cccc                                        >  2:231394/1‑101 (MQ=18)
                                                              cTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCA..GAGGGGGGC.AGCG.cccc                                        <  2:1471921/101‑1 (MQ=17)
                                                               tCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCA..GAGGGGGGC.AGCG.ccccc                                       >  2:633975/1‑101 (MQ=17)
                                                                 cccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCA..GAGGGGGGC.AGCG.cccccct                                     <  2:936538/101‑3 (MQ=17)
                                                                 cccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCA..GAGGGGGGC.AGCG.cccccct                                     >  2:3155734/1‑99 (MQ=17)
                                                                 cccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCA..GAGGGGGGC.AGCG.cccccct                                     <  1:4113550/101‑3 (MQ=17)
                                                                 cccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCA..GAGGGGGGC.AGCG.acccccg                                     >  2:2285455/1‑99 (MQ=17)
                                                                     cGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGA..GAAGCCA..GAGGGGGGC.AGCGCCCCCCTGccc                                  <  1:3748899/101‑1 (MQ=17)
                                                                      gggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCA..GAGGGGGGC.AGCGCCCCCCTGcccc                                 <  2:3255155/101‑1 (MQ=17)
                                                                      gggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCA..GAGGGGGGC.AGCGCCCCCCTGcccc                                 <  2:3008822/101‑1 (MQ=17)
                                                                       ggATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AGGCC..AAGAGGGGGCAGCG.CCCCCTGCCCcg                                >  1:176553/1‑101 (MQ=2)
                                                                        gATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AGGCC..AAGAGGGGGCAGCG.CCCCCTGCCCcgc                               >  1:286327/1‑100 (MQ=2)
                                                                        gATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCA..GAGGGGGGC.AGCGCCCCCCTGCCCcgg                               >  2:3035939/1‑101 (MQ=18)
                                                                         atatTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AGGCC..AAGAGGGGGCAGCG.CCCCCTGCCCcgcc                              <  1:2038261/101‑3 (MQ=2)
                                                                         atatTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCA..GAGGGGGGC.AGCGCCCCCCTGCCCcggc                              >  2:2615727/1‑101 (MQ=18)
                                                                           atTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGC..C.AAACGCCAAGAGGGGGAAGCG.CCCCCTGCCCcgccgg                            >  1:3675887/1‑101 (MQ=2)
                                                                           atTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCA..GAGGGGGGC.AGCGCCCCCCTGCCCcggcgg                            >  1:3984679/1‑101 (MQ=17)
                                                                             ttGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AGGCC..AAGAGGGGGCAGCG.CCCCCTGCCCcgccggat                          <  1:582313/101‑3 (MQ=2)
                                                                              tGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCA..GAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTcc                         >  1:3235319/1‑101 (MQ=17)
                                                                               gcagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCA..GAGGGGGGCATCCC.CCCCCTGCCCCGGCGGTCCg                        <  2:3138149/101‑1 (MQ=18)
                                                                                 agcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCA..GAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTg                      <  2:1717607/101‑1 (MQ=17)
                                                                                  gcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCA..GAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTgg                     <  2:1059460/101‑1 (MQ=17)
                                                                                   caAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..GAAGCCA..GAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGGCCGTggg                    >  1:143678/1‑100 (MQ=17)
                                                                                     aGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCA..GAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTggcgg                  <  1:2538613/101‑1 (MQ=17)
                                                                                     aGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCA..GAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTggcgg                  >  2:625069/1‑101 (MQ=18)
                                                                                       aTGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCA..GAGGGGGGCAGCTC.CCCCCTGCCCCGGCGGTCCGTGgcggcg                >  1:2845295/1‑101 (MQ=18)
                                                                                       aTGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCA..GAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGgcggcg                <  1:3195103/101‑1 (MQ=17)
                                                                                         gAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGCGAA.TAGC..GAAGCCA..GAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGcc              >  1:335947/1‑101 (MQ=17)
                                                                                          aaaCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCA..GAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGGGCCt             >  2:911065/1‑101 (MQ=17)
                                                                                          aaaCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCA..GAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCt             <  1:3680930/101‑1 (MQ=17)
                                                                                          aaaCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..GAAGCCA..GAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGGCCGGGGCGGCGGCt             >  1:2207444/1‑101 (MQ=9)
                                                                                             ccGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCA..GAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGa          >  1:2841566/1‑101 (MQ=17)
                                                                                              cGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCA..GAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAc         <  1:882183/101‑1 (MQ=17)
                                                                                                  tGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCA..GAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcgggc     <  2:1784540/101‑1 (MQ=17)
                                                                                                   gCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCA..GAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcgggcg    <  2:2046365/101‑1 (MQ=17)
                                                                                                    ccTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCA..GAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcgggcgg   <  1:1165216/101‑1 (MQ=17)
                                                                                                     cTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCA..GAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcgggcggg  >  2:3396601/1‑101 (MQ=18)
                                                                                                     |                                                                                                          
CGGGCGCGCGGCACGACAGGCCGAGTAA..GGG.GGGGCGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AGCAG..AAGGGGGGCAAGCG.CCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGG  >  NC_014034/2493925‑2494122

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