Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,494,015 C→G 4.1% intergenic (+63/‑179) RCAP_rcc02319 → / → RCAP_rcc02320 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,494,0150C→G 4.1% 18.2 ‑0.3 ‑0.0290intergenic (+63/‑179)RCAP_rcc02319/RCAP_rcc02320hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (6/6):  ref base (120/158):  total (126/164)
Fisher's exact test strand distribution p-value = 7.68e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.21e-01

Only 100 of 317 total aligned reads displayed.

GCGCCGCGGCGGGCGCGCGGCACGACAGGCC.GAGTAA..GGGGGGG.CGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AG.CAGAAGGGGGGCAAGCG.CCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGA  >  NC_014034/2493916‑2494114
                                                                                                       |                                                                                                         
gcgcCGCGGCGGGCGCGCGGCACGACAGGCC.GAGTAAG.GGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAAc                                                                                                           <  2:2947535/101‑1 (MQ=21)
 cgatcgggcgggcgGGGGGCACGAGGGGCG.GATTGG..GGGGGGGCCGCCACCCCCCGTCCGAAAGCGGCCTCCCCCCGGGATTTTTGCAGCAAGATGAAAcc                                                                                                          <  2:3360875/95‑1 (MQ=12)
  gccgcGGCGGGCGCGCGGCACGACAGGCC.GAGTAAG.GGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCg                                                                                                         <  1:6574276/101‑1 (MQ=21)
   ccgcGGCGGGCGCGCGGCACGACAGGCC.GAGTAAG.GGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCgg                                                                                                        <  1:1632771/101‑1 (MQ=21)
    cgggcggggcgCGCGGC.CGGCCGGCC.GGGTAAGGGGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGa                                                                                                       <  2:4891201/95‑1 (MQ=14)
     gcggcgGGCGCGCGGCACGACAGGCCGGGGAAA..GGGGGGG.GGCTGCCCCCCGTCCGGCAGCGGGCTCCCCCCCGGCAATTTGCCGCACAAAGAAACCGGaa                                                                                                      >  1:1170684/1‑100 (MQ=255)
       ggcgggcgCGCGGCACGACAGGCC.GAGTAAG.GGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGc                                                                                                    <  1:3737208/101‑1 (MQ=21)
          gggcgCGCGGCACGACAGGCC.GAGTAA..GGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCttt                                                                                                <  1:5570151/101‑1 (MQ=255)
           ggcgCGCGGCACGACAGGCC.GAGTAAG.GGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCttt                                                                                                >  1:5868744/1‑101 (MQ=21)
            gcgcgcGGCACGACAGGCC.GAGTAAG.GGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTc                                                                                               <  2:166973/101‑1 (MQ=21)
             cgcgcgGCACGACAGGCC.GAGTAAG.GGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCa                                                                                              <  2:1249204/101‑1 (MQ=21)
              gcgggGCCCGACAGGCC.GAGTAAG.GGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCAt                                                                                             <  1:6050482/101‑1 (MQ=18)
              gcgcgGCACGACAGGCC.GAGTAAG.GGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCAt                                                                                             >  2:1136266/1‑101 (MQ=21)
              gcgcgGCACGACAGGCC.GAGTAAG.GGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCAt                                                                                             <  1:5963658/101‑1 (MQ=21)
              gcgcgGCACGACAGGCC.GAGTAAG.GGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCAt                                                                                             <  1:4183156/101‑1 (MQ=21)
               cgcgGCACGACAGGCC.GAGTAAG.GGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATc                                                                                            <  1:4354939/101‑1 (MQ=21)
                gcgGCACGACAGGCC.GAGTAAG.GGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCt                                                                                           <  1:458653/101‑1 (MQ=21)
                gcgGCACGACAGGCC.GAGTAAG.GGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCt                                                                                           <  2:3165688/101‑1 (MQ=21)
                 cgGCAGGACAGGCC.GAGTAAG.GGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCtt                                                                                          <  1:1844552/101‑1 (MQ=18)
                 cgGCACGACAGGCC.GAGTAAG.GGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCtt                                                                                          <  1:3289486/101‑1 (MQ=21)
                   gCACGACAGGCC.GAGTAAG.GGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGc                                                                                        >  2:2053341/1‑101 (MQ=21)
                    cACGACAGGCC.GAGTAAG.GGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGcc                                                                                       <  2:1102071/101‑1 (MQ=21)
                     ccGACAGGCC.GAGTAAG.GGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGccc                                                                                      <  2:2847368/100‑1 (MQ=21)
                      cGACAGGCC.GAGTAAG.GGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCg                                                                                     <  2:3126064/101‑1 (MQ=21)
                      cGACAGGCC.GAGTAAG.GGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCg                                                                                     <  2:4467361/101‑1 (MQ=21)
                       gACAGGCC.GAGTAAG.GGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGGCCGa                                                                                    >  1:3212653/1‑101 (MQ=18)
                         cAGGCC.GAGTAAG.GGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGaaa                                                                                  <  1:6101113/101‑1 (MQ=21)
                         cAGGCC.GAGTAAG.GGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGaaa                                                                                  <  1:3644490/101‑1 (MQ=21)
                         cAGGCC.GAGTAAG.GGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGaaa                                                                                  <  2:6934405/101‑1 (MQ=21)
                          aGGCC.GAGTAAG.GGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAat                                                                                 <  1:2055620/101‑1 (MQ=21)
                           ggCC.GAGTAAG.GGGGGGG.GGATGCCCCCCGTTCGCAAGCGGACTCCCCCCGGGACATTTGCAGCAAGCTGAAACCGGACGCGTTTCATCTTGCCCGAAatc                                                                                >  2:6427716/1‑100 (MQ=18)
                           ggCC.GAGTAAG.GGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAata                                                                                <  2:2762090/101‑1 (MQ=21)
                            gCC.GAGTAAG.GGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAatat                                                                               <  1:4924078/101‑1 (MQ=21)
                             cc.GATTAAG.GGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATc                                                                              <  2:40144/101‑1 (MQ=18)
                             cc.GAGTAAG.GGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATc                                                                              <  2:5613774/101‑1 (MQ=21)
                             cc.GAGTAAG.GGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATc                                                                              <  2:535315/101‑1 (MQ=21)
                                gAGTAAG.GGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATccc                                                                            >  1:5242784/1‑101 (MQ=18)
                                cgactaagggggggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAatatac                                                                             >  1:120244/9‑99 (MQ=18)
                                    taagggggggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCggg                                                                         >  2:2995194/5‑101 (MQ=255)
                                     aagggggggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggg                                                                        >  2:707962/4‑101 (MQ=255)
                                     aagggggggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggg                                                                        <  1:2951550/98‑1 (MQ=255)
                                      agggggggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCggggg                                                                       >  1:1594922/3‑101 (MQ=255)
                                      agggggggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATACCggggg                                                                       >  1:1715591/3‑101 (MQ=255)
                                       gggggggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggggg                                                                      >  1:3122284/2‑101 (MQ=255)
                                         gggggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggggggg                                                                    >  1:3456816/1‑99 (MQ=255)
                                         gggggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggggggg                                                                    >  2:6196447/1‑99 (MQ=255)
                                         gggggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTg                                                                    <  1:1978697/101‑1 (MQ=255)
                                         gggggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTg                                                                    <  1:6943153/101‑1 (MQ=255)
                                          ggggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGa                                                                   >  1:6687036/1‑101 (MQ=255)
                                          ggggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGa                                                                   >  2:4706240/1‑101 (MQ=255)
                                          ggggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGa                                                                   <  1:914235/101‑1 (MQ=255)
                                           gggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGaa                                                                  >  2:1473154/1‑101 (MQ=255)
                                           gggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGaa                                                                  >  2:4137507/1‑101 (MQ=255)
                                             gg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.Tg                                                               >  2:3557005/1‑101 (MQ=255)
                                             gg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.Tg                                                               <  2:530791/101‑1 (MQ=255)
                                              g.cgcTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.Tgg                                                              <  1:977624/101‑1 (MQ=255)
                                              g.cgcTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTCAA.Tgg                                                              <  2:3512697/101‑1 (MQ=255)
                                                cgcTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGc                                                             <  2:1639586/101‑1 (MQ=255)
                                                  cTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..G.a                                                        >  1:5416968/1‑101 (MQ=255)
                                                    gCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aag                                                      <  2:3645725/101‑3 (MQ=255)
                                                      ccccG.TCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGACCGAAATATCCCGGGGGGTGAA.TGGC..G.aagcca                                                   >  2:6712796/1‑96 (MQ=21)
                                                        cccgtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccag                                                  <  1:2613364/97‑7 (MQ=255)
                                                         ccgtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccaga                                                 <  1:5334996/98‑8 (MQ=255)
                                                                 aaGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGCa                                         <  1:4841224/101‑1 (MQ=18)
                                                                   gCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGCagc                                       <  1:3524300/101‑3 (MQ=18)
                                                                    cGGCCTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGCagcg                                      <  2:4391520/101‑4 (MQ=18)
                                                                    cGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGCagcg                                      >  2:5775738/1‑98 (MQ=18)
                                                                    cGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGCagcg                                      >  1:381771/1‑98 (MQ=18)
                                                                    cGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGCagcg                                      >  2:5033816/1‑98 (MQ=18)
                                                                        cTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGC.AGCG.cccc                                <  2:252539/101‑1 (MQ=17)
                                                                        cTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGCaccgcccc                                  >  2:5968552/1‑94 (MQ=18)
                                                                        cTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCG.cccc                                <  2:718542/101‑1 (MQ=17)
                                                                          ccccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCG.cccccc                              >  1:2407036/1‑100 (MQ=17)
                                                                          ccccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCG.cccccc                              >  1:884095/1‑100 (MQ=17)
                                                                           cccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCG.cccccct                             <  2:5687321/101‑3 (MQ=17)
                                                                                gggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGcccc                         <  2:4096890/101‑1 (MQ=17)
                                                                                gggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGcccc                         >  2:1766453/1‑101 (MQ=17)
                                                                                gggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGcccc                         >  1:6612003/1‑101 (MQ=17)
                                                                                  gATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCcgg                       >  2:4306352/1‑101 (MQ=18)
                                                                                   atatTTGCAGCAAGATGAACCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCcggc                      <  1:6198708/101‑1 (MQ=17)
                                                                                    tatTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCcggcg                     >  1:6142998/1‑101 (MQ=17)
                                                                                     atTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCAAGAGGGGGCAGCG.CCCCCTGCCCcgccgg                    >  1:3954242/1‑101 (MQ=2)
                                                                                     atTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGCGGGCAGCCC.CCCCCTGCCCcggcgg                    >  2:6129407/1‑101 (MQ=18)
                                                                                       ttGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCAAGAGGGGGCAGCG.CCCCCTGCCCcgccggat                  <  2:6476340/101‑3 (MQ=2)
                                                                                        tGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTcc                 <  1:707982/101‑1 (MQ=17)
                                                                                        tGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTcc                 <  2:2441416/101‑1 (MQ=17)
                                                                                          cagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AA.GCCAGGGGGGGGCAGCGCCCCCCTGCCCCGGCGGTCCGt               <  1:1021714/101‑1 (MQ=18)
                                                                                              aaGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGCCCGTggcg           >  1:3383889/1‑101 (MQ=17)
                                                                                               aGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGGGGTCCGGgccgg          >  1:2858290/1‑101 (MQ=9)
                                                                                                gATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTggcggc         <  1:6883905/101‑1 (MQ=17)
                                                                                                  tGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGA..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGgcgc       >  2:4546162/1‑101 (MQ=18)
                                                                                                   gAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGcc      >  1:334723/1‑101 (MQ=18)
                                                                                                   gAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGcc      <  2:6581054/101‑1 (MQ=17)
                                                                                                    aaaCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCt     >  2:2019594/1‑101 (MQ=18)
                                                                                                    aaaCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGCCCGTCGCGGCGCCt     >  2:1856317/1‑101 (MQ=14)
                                                                                                     aaCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCtt    >  1:2830502/1‑101 (MQ=18)
                                                                                                     aaCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGACCCCCTGCCCCGGCGGTCCGTGGCGGCGCCtt    >  2:4633550/1‑101 (MQ=18)
                                                                                                      aCCGGATGCCTTTCATCTTGCCCGAAATATGCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTg   <  1:229743/101‑1 (MQ=17)
                                                                                                      aCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTg   <  1:4103792/101‑1 (MQ=17)
                                                                                                      aCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTg   <  2:375861/101‑1 (MQ=17)
                                                                                                       |                                                                                                         
GCGCCGCGGCGGGCGCGCGGCACGACAGGCC.GAGTAA..GGGGGGG.CGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AG.CAGAAGGGGGGCAAGCG.CCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGA  >  NC_014034/2493916‑2494114

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG