Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,494,003 A→G 7.2% intergenic (+51/‑191) RCAP_rcc02319 → / → RCAP_rcc02320 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,494,0030A→G 7.2% 17.5 ‑1.1 ‑0.0139intergenic (+51/‑191)RCAP_rcc02319/RCAP_rcc02320hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (6/4):  ref base (52/77):  total (58/81)
Fisher's exact test strand distribution p-value = 3.20e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.81e-01

Only 100 of 165 total aligned reads displayed.

GCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAA..GGGG.GGG.CGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AGCAGAAGGGGGGCAAGCG.CCCCCTGCCCCGGCGGTCCGTG  >  NC_014034/2493904‑2494102
                                                                                                       |                                                                                                        
gcgTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTgca                                                                                                          <  1:2111405/101‑1 (MQ=25)
 cgTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTgcag                                                                                                         <  1:3663506/101‑1 (MQ=25)
  gTCTGGCCAAGCAGCGCGGCGGGCGCGCGGGGCGACAGCCTGAGTAA..GGGGTGGG.CGATGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTgcagc                                                                                                        <  2:2670810/101‑1 (MQ=17)
   tCTTGCCCAGCGCCGCGGCGGGCGCGCGGCCCGACAGGCCGAGTAA..GGGG.GGGTCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTgcagca                                                                                                       <  2:4078883/101‑1 (MQ=21)
   tCTGGCCCAGCGGCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GGGG.GGG.GGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGGTATTTgccgca                                                                                                       >  1:2298543/1‑101 (MQ=21)
   tCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTgcagca                                                                                                       >  2:87518/1‑101 (MQ=25)
   tCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTgcagca                                                                                                       <  1:3974864/101‑1 (MQ=25)
     tGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAg                                                                                                     <  1:2475969/101‑1 (MQ=21)
     tGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAg                                                                                                     <  1:1563538/101‑1 (MQ=21)
      ggCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCCAGCGGACTCCCCCCGGGATATTTGCAGCAAGa                                                                                                    >  1:1784659/1‑101 (MQ=21)
      ggCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGGTATTTTCAGCAAGa                                                                                                    >  1:3680077/1‑101 (MQ=21)
        cccAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATg                                                                                                  <  2:573109/101‑1 (MQ=21)
        cccAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATg                                                                                                  <  2:4193359/101‑1 (MQ=21)
          cAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGaa                                                                                                <  1:450733/101‑1 (MQ=21)
          cAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACGCCCCCCGGGATATTTGCACCAAGATGaa                                                                                                >  2:547352/1‑101 (MQ=18)
            gcgcCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAAc                                                                                              <  2:3840071/101‑1 (MQ=21)
              gccgcGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCg                                                                                            <  2:2745107/101‑1 (MQ=21)
              gccgcGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCg                                                                                            >  2:2630897/1‑101 (MQ=21)
              gccgcGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCg                                                                                            <  2:2226720/101‑1 (MQ=21)
               ccgcGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCgg                                                                                           <  1:93742/101‑1 (MQ=21)
                cgcgGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGa                                                                                          <  1:3543301/101‑1 (MQ=21)
                   ggcgggcgCGCGGCACGACAGGCCGAGTAAG.GGGG.GGG.GGCTGCCCCCCGTCCGCCAGCGGGCTCCCCCCGGGATATTTTCAGCAAGATGAAACCGGATGc                                                                                       >  2:2461382/1‑101 (MQ=18)
                   ggcgggcgCGCGGCACGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGc                                                                                       <  1:2661231/101‑1 (MQ=21)
                     cgggcgCGCGGCACGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCt                                                                                     >  2:1314684/1‑101 (MQ=21)
                      gggcgCGCGGCACGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCCGGCTCCCCCCGGGATATTTGCAGCAAGAGGAAACCGGATGCCtt                                                                                    >  2:1852319/1‑101 (MQ=21)
                       ggcgCGCGGCACGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCttt                                                                                   <  2:3798128/101‑1 (MQ=21)
                        gcgcgcGGCACGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTc                                                                                  >  1:2611616/1‑101 (MQ=21)
                        gcgcgcGGCACGACAGGCCGAGGAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTc                                                                                  <  2:1972186/101‑1 (MQ=18)
                         ggcgcgGCACGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCa                                                                                 <  2:1565830/100‑1 (MQ=18)
                         cgcgcgGCACGACAGGCCGAGTAAGGGGGG.GGG.CGCTGCCCCCCGCCCGCAAGCGGACTCCCCCCGGGAATTTTGCAGCAAAATGAAACCGGAGGCCTTTc                                                                                  >  1:1059554/1‑101 (MQ=18)
                         cgcgcgGCACGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCa                                                                                 <  1:3896555/101‑1 (MQ=21)
                         cgcgcgGCACGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCa                                                                                 <  1:337161/101‑1 (MQ=21)
                         cccgcgGCACGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCa                                                                                 <  1:1388558/99‑1 (MQ=18)
                          gcgcgGCACGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGACTTTCAt                                                                                >  2:734027/1‑101 (MQ=18)
                            gcgGCACGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCt                                                                              >  1:3138380/1‑101 (MQ=21)
                            gcgGCACGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCt                                                                              <  2:1549294/101‑1 (MQ=21)
                             cgGCACGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCtt                                                                             <  1:134877/101‑1 (MQ=21)
                               gCACGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGc                                                                           <  1:3274269/101‑1 (MQ=21)
                               gCACGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGc                                                                           >  2:1165122/1‑101 (MQ=21)
                                cACGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGcc                                                                          >  1:1014986/1‑101 (MQ=21)
                                 aCGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGccc                                                                         <  2:2354520/101‑1 (MQ=21)
                                  cGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCg                                                                        >  1:1888872/1‑101 (MQ=21)
                                  cGACAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGAGATTTGCAGCAAGATGAAACCGGACGCCTTTCATCTTGCCCg                                                                        >  1:1524658/1‑101 (MQ=18)
                                     cAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGaaa                                                                     <  2:3183986/101‑1 (MQ=21)
                                     cAGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCCGGGCATTTGCAGAAAGATGAAAACGGATGCCTTTTATCTTGTCCGaca                                                                     >  1:833567/1‑99 (MQ=17)
                                      aGGCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAat                                                                    <  1:2461566/101‑1 (MQ=21)
                                       ggCCGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAata                                                                   <  1:3201927/101‑1 (MQ=21)
                                          cGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATcc                                                                <  1:3496873/101‑1 (MQ=21)
                                          cGAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATcc                                                                >  1:1615732/1‑101 (MQ=18)
                                           gAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATccc                                                               <  1:2594615/101‑1 (MQ=18)
                                           gAGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATccc                                                               <  1:275202/101‑1 (MQ=18)
                                            aGTAAG.GGGG.GGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCg                                                              <  1:1137456/101‑1 (MQ=18)
                                             aggaaggggg.ggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCg                                                              <  2:177793/95‑1 (MQ=18)
                                              gtaaggggg.ggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgg                                                             <  1:41580/96‑1 (MQ=255)
                                               taaggggg.ggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGAGATTTGCAGCACGATGAAACCGGATGCCTTTCATATTGCCCGCAATATCCCcgg                                                            >  2:2879137/5‑101 (MQ=18)
                                                aaggggg.ggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggg                                                           <  1:2488527/98‑1 (MQ=255)
                                                 aggggg.ggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCggggg                                                          <  2:2205881/99‑1 (MQ=255)
                                                 aggggg.ggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCggggg                                                          <  1:2500196/99‑1 (MQ=255)
                                                   gggg.ggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGCAATATCCCggggggg                                                        >  2:1224466/1‑100 (MQ=255)
                                                   gggg.ggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGt                                                        >  1:1407280/1‑101 (MQ=255)
                                                   gggg.ggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGt                                                        <  2:3679823/101‑1 (MQ=255)
                                                     gg.ggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGa                                                      >  1:2113492/1‑101 (MQ=255)
                                                      g.ggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGaa                                                     <  1:3964032/101‑1 (MQ=255)
                                                        ggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.t                                                   >  1:2855095/1‑101 (MQ=255)
                                                              cTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.a                                           >  2:3680676/1‑101 (MQ=255)
                                                               tGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..G.aa                                          >  1:3321008/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.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)
                                                                      cgtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccagag                                   <  1:4031167/99‑9 (MQ=255)
                                                                       gtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccagagg                                  <  2:4026053/100‑10 (MQ=255)
                                                                             aaGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGCa                             <  1:1217599/101‑1 (MQ=18)
                                                                                cGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGCagcg                          >  1:262851/1‑98 (MQ=18)
                                                                                 ggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAGGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGCaccgc                         >  1:1861491/1‑97 (MQ=17)
                                                                                  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:1471921/101‑1 (MQ=17)
                                                                                     tCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCG.ccccc                   >  2:633975/1‑101 (MQ=17)
                                                                                     tCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCG.ccccc                   >  2:321368/1‑101 (MQ=18)
                                                                                       cccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCG.cccccct                 <  2:936538/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.acccccg                 >  2:2285455/1‑99 (MQ=17)
                                                                                         cccGGGATAGTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCG.cccccctgc               <  2:2113354/101‑5 (MQ=17)
                                                                                           cGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGA..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGccc              <  1:3748899/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:76569/101‑2 (MQ=0)
                                                                                              gATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCcgg           >  2:3035939/1‑101 (MQ=18)
                                                                                               atatTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AGGCCAAGAGGGGGCAGCG.CCCCCTGCCCcgcc          >  2:2911175/1‑99 (MQ=2)
                                                                                               atatTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..GAAGCCAGAGGGGGGC.CGCGCCCCCCCGCCCcggc          >  1:1983187/1‑101 (MQ=14)
                                                                                                 atTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCcggcgg        >  1:3984679/1‑101 (MQ=17)
                                                                                                 atTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCcggcgg        >  2:3916301/1‑101 (MQ=18)
                                                                                                   ttGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AGGCCAAGAGGGGGCAGCG.CCCCCTGCCCcgccggat      <  1:582313/101‑3 (MQ=2)
                                                                                                     gcagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGCATCCC.CCCCCTGCCCCGGCGGTCCg    <  2:3138149/101‑1 (MQ=18)
                                                                                                       agcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTg  <  2:1717607/101‑1 (MQ=17)
                                                                                                       agcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATAGCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTg  <  1:181884/101‑1 (MQ=17)
                                                                                                       |                                                                                                        
GCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAA..GGGG.GGG.CGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AGCAGAAGGGGGGCAAGCG.CCCCCTGCCCCGGCGGTCCGTG  >  NC_014034/2493904‑2494102

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