Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,494,014 A→G 6.5% intergenic (+62/‑180) RCAP_rcc02319 → / → RCAP_rcc02320 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,494,0140A→G 6.5% 19.4 0.0 ‑0.1168intergenic (+62/‑180)RCAP_rcc02319/RCAP_rcc02320hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (4/7):  ref base (57/100):  total (61/107)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.63e-01

Only 100 of 184 total aligned reads displayed.

AGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.A.A.GGGGGGGCGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AGCAGA.AGGGGGGCAAGCG.CC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTG  >  NC_014034/2493915‑2494113
                                                                                                      |                                                                                                          
aGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGaaa                                                                                                            <  2:3113942/101‑1 (MQ=21)
  cgccgcGGCGGGCGCGCGGCACGACAGGCCGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGGTATTTGCAGCCAAATTAAAcc                                                                                                          >  1:4142835/1‑101 (MQ=18)
  cgccgcGGCGGGCGCGCGGCACGACAGGCCGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAAcc                                                                                                          <  2:1418470/101‑1 (MQ=21)
  cgccgcGCCGGGCGCGCGGCACGACAGGCCGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAAcc                                                                                                          <  2:2107580/101‑1 (MQ=18)
   gccgcGGCGGGCGCGCGGCACGACAGGCCGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCg                                                                                                         <  1:2286702/101‑1 (MQ=21)
      gcggcgGGCGCGCGGGACGACAGGCCGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAt                                                                                                      <  2:2440153/101‑1 (MQ=18)
      gcggcgGGCGCGCGGCACGACAGGCCGAGTAA.A.GGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGag                                                                                                      >  2:1337476/1‑100 (MQ=18)
      gcggcgGGCGCGCGGCACGACAGGCCGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAt                                                                                                      <  1:2455584/101‑1 (MQ=21)
         gcgggcgCGCGGCACGACAGGCCGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGcc                                                                                                   <  1:1811761/101‑1 (MQ=21)
          cgggcgCGCGGCACGACAGGCCGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCt                                                                                                  <  2:1168021/101‑1 (MQ=21)
           gggcgCGCGGCACGACAGGCCGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCtt                                                                                                 >  2:628191/1‑101 (MQ=21)
            ggcgCGCGGCACGACAGGCCGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGGCATTTGCAGCAAGATGAAACCGGTTGCCttt                                                                                                >  2:1014575/1‑101 (MQ=18)
              ggcgcgGCGCGACAGGCCGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCa                                                                                              <  2:3910641/100‑1 (MQ=18)
              cgcgcgGCACGACAGGCCGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCa                                                                                              <  2:171363/101‑1 (MQ=21)
              cgcgcgGCACGACAGGCCGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCa                                                                                              <  1:3510733/101‑1 (MQ=21)
              cgcgcgGCACGACAGGCCGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCa                                                                                              <  1:2563115/101‑1 (MQ=21)
                cgcgGCACGACAGGCCGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATc                                                                                            <  2:420997/101‑1 (MQ=21)
                 gcgGCACGACAGGCCGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCt                                                                                           <  1:736814/101‑1 (MQ=21)
                  cgGCACGACAGGCCGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCtt                                                                                          <  1:1235454/101‑1 (MQ=21)
                   ggCACGACAGGCCGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTg                                                                                         <  2:2326949/101‑1 (MQ=21)
                   ggCACGACAGGCCGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTg                                                                                         <  2:1014675/101‑1 (MQ=21)
                    gCACGACAGGCCGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGc                                                                                        >  2:2713331/1‑101 (MQ=21)
                     cACGACAGGCCGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGcc                                                                                       >  2:3568212/1‑101 (MQ=21)
                      aCGACAGGCCGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGccc                                                                                      >  2:1249258/1‑101 (MQ=21)
                      aCGACAGGCCGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGccc                                                                                      <  1:1075202/101‑1 (MQ=21)
                       cGACGGGCCGAGT.AGG.GGGGGGGCGCTCCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATAAAACCGGATGCCTTTCATCTTGCCCg                                                                                     <  2:2023637/101‑1 (MQ=18)
                       cGACAGGCCGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCg                                                                                     <  2:2359697/101‑1 (MQ=21)
                       cGACAGGCCGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTATCATCTTGCCCg                                                                                     <  1:2350473/101‑1 (MQ=18)
                         aCAGGCCGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGaa                                                                                   <  2:3266273/101‑1 (MQ=21)
                          cAGGCCGATT.AGG.GGGGGGGCGCTGCCCCCCGTCCGCAATCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGaaa                                                                                  <  2:1865859/101‑1 (MQ=18)
                            ggCCGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAata                                                                                <  1:1698919/101‑1 (MQ=21)
                            ggCCGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAata                                                                                <  2:1477498/101‑1 (MQ=21)
                            ggCCGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAata                                                                                >  2:4043230/1‑101 (MQ=21)
                            ggCCGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAata                                                                                <  1:269917/101‑1 (MQ=21)
                             gCCGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAatat                                                                               <  1:2503497/101‑1 (MQ=21)
                             gCCGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAatat                                                                               <  1:1991105/101‑1 (MQ=21)
                             atctaGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAatat                                                                               <  1:4384874/97‑1 (MQ=18)
                              ccGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATc                                                                              <  1:1463652/101‑1 (MQ=21)
                               cGAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATcc                                                                             <  1:3699225/101‑1 (MQ=18)
                                gAGT.A.AGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATccc                                                                            <  1:3869050/101‑1 (MQ=18)
                                       aggggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCggggg                                                                       >  1:281857/3‑101 (MQ=255)
                                         gggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGt                                                                     <  1:701182/101‑1 (MQ=255)
                                          ggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTg                                                                    <  2:1735060/101‑1 (MQ=255)
                                              ggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.Tg                                                               <  2:257003/101‑1 (MQ=255)
                                               gcgcTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGACATTTGCAGCAAGATGAAAACGGATGCCTTTCATCTTGCCCGAACTATACCGGGGGGGGAG.Tgg                                                              >  1:883252/1‑101 (MQ=18)
                                                cgcTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGc                                                             <  1:2374334/101‑1 (MQ=255)
                                                 gcTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..g                                                          <  1:3097965/101‑1 (MQ=255)
                                                   tGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aa                                                       <  1:2144125/101‑2 (MQ=255)
                                                    gCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTTCCCGAAATATCTCTGTGGGGTAA.TGGT..G.atc                                                      >  1:1313118/1‑99 (MQ=18)
                                                    gCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aag                                                      <  1:1081080/101‑3 (MQ=255)
                                                     cccccG.TCCGCAAGCGGACTCCCCCCGGGATATTTGCCGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGAGAA.TGGC..G.aagcc                                                    >  1:2946572/1‑97 (MQ=18)
                                                      ccccG.TCCTCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATTCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagcca                                                   <  1:1331605/101‑6 (MQ=18)
                                                      ccccG.TCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagcca                                                   <  1:1897539/101‑6 (MQ=255)
                                                        cccgtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccag                                                  <  1:769318/97‑7 (MQ=255)
                                                        cccgtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccag                                                  >  1:3570867/5‑95 (MQ=255)
                                                          cgtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccagag                                                >  1:3401873/3‑93 (MQ=255)
                                                           gtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..G.aagccagagg                                               >  1:1827773/2‑92 (MQ=21)
                                                            tCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccagaggg                                              >  2:2511595/1‑91 (MQ=255)
                                                            tCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..G.aagccagaggg                                              >  2:3007661/1‑91 (MQ=18)
                                                                 aaGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGCa                                         >  2:2302328/1‑101 (MQ=18)
                                                                 aaGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..GAAGCCAG.AGGGGGGCa                                         >  1:4364755/1‑101 (MQ=18)
                                                                  aGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGCag                                        <  2:413502/101‑2 (MQ=18)
                                                                  aGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGCag                                        <  1:2028967/101‑2 (MQ=18)
                                                                   gCGGACTCCCCCCGGGATTTTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGCagc                                       <  2:4145191/101‑3 (MQ=18)
                                                                   gCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGCagc                                       >  2:450441/1‑99 (MQ=18)
                                                                        cTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGC.AGCG.cc.cc                               >  2:688702/1‑101 (MQ=17)
                                                                            ccccGGGATATGTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGC.AGCG.cc.cccctg                           <  2:1322038/101‑4 (MQ=17)
                                                                              ccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGC.AGCG.cc.cccctgcc                         <  2:2951645/101‑6 (MQ=17)
                                                                               cGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGC.AGCGACC.CCCTGccc                         >  2:1214419/1‑101 (MQ=17)
                                                                                gggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGC.AGCGCCC.CCCTGcccc                        >  1:576951/1‑101 (MQ=17)
                                                                                 ggATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AGGCCA.AGAGGGGGCAGCG.CC.CCCTGCCCcg                       >  2:1447132/1‑101 (MQ=0)
                                                                                 ggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGAGAA.GGGC..GAAGCCAG.AGGGGGGC.AGCGCCC.CCCTGCCCcg                       >  1:4009471/1‑101 (MQ=17)
                                                                                 ggATANNTGCAGAAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGC.AGCGCCC.CCCTGCCCcg                       <  1:1594892/101‑1 (MQ=17)
                                                                                  gATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AGGCCA.AGAGGGGGCAGCG.CC.CCCTGCCCcgc                      >  2:1683206/1‑100 (MQ=2)
                                                                                  gATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AGGCCA.AGAGGGGGCAGCG.CC.CCCTGCCCcgc                      <  1:2793530/101‑2 (MQ=0)
                                                                                  gATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGC.AGCGCCC.CCCTGCCCcgg                      >  1:2613171/1‑101 (MQ=17)
                                                                                   atatTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AGGCCA.AGAGGGGGCAGCG.CC.CCCTGCCCcgcc                     <  1:1155047/101‑3 (MQ=2)
                                                                                   atatTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AGGCCA.AGAGGGGGCAGCG.CC.CCCTGCCCcgcc                     <  2:2819530/101‑3 (MQ=2)
                                                                                    tatTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGC.AGCGCCC.CCCTGCCCcggcg                    >  2:3118786/1‑101 (MQ=18)
                                                                                    tatTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AAGCCA.GAGGGGGC.AGCGCCC.CCCTGCCCcggcgg                   >  2:225289/1‑101 (MQ=17)
                                                                                    tatTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGCGAA.TGGC..GAAGCCAG.AGGGGGGC.AGCGCCC.CCCTGCACcggcg                    >  2:341801/1‑101 (MQ=17)
                                                                                     atTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AGGCCA.AGAGGGGGCAGCG.CC.CCCTGCCCcgccgg                   <  2:337533/101‑1 (MQ=2)
                                                                                     atTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AAGACAGAGGGGGGC.AGCGCCC.CCCTGCCCcggcgg                   >  1:477110/1‑101 (MQ=18)
                                                                                     atTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..GAAGCCAG.AGGGGGGCAGCAC.CC.CCCTGCCCcggcgg                   >  1:501559/1‑101 (MQ=17)
                                                                                       ttGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AGGCCA.AGAGGGGGCAGCG.CC.CCCTGCCCcgccggat                 >  2:91130/1‑99 (MQ=2)
                                                                                       ttGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTc                 <  2:3706895/101‑1 (MQ=17)
                                                                                        tGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGC.AGCG.CCGCCCTGCCCCGGCGGTcc                >  2:1376656/1‑101 (MQ=18)
                                                                                          cagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGTGC.AGCGCCC.CCCTGCCCCGGCGGTCCGt              >  1:1490453/1‑101 (MQ=18)
                                                                                          cagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCCGCGGGCCgg              >  2:3847046/1‑100 (MQ=14)
                                                                                           agcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTg             >  2:1075039/1‑101 (MQ=18)
                                                                                           agcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTg             >  1:3046297/1‑101 (MQ=17)
                                                                                            gcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTgg            >  1:1966347/1‑101 (MQ=17)
                                                                                             caAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTggc           <  1:3185265/101‑1 (MQ=17)
                                                                                              aaGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTggcg          <  2:1110206/101‑1 (MQ=17)
                                                                                               aGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTggcgg         <  2:1229934/101‑1 (MQ=17)
                                                                                                gATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAA.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTggcggc        <  1:2171140/101‑1 (MQ=18)
                                                                                                  tGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGgcgc      >  2:4269411/1‑101 (MQ=17)
                                                                                                    aaaCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCt    <  1:4261496/101‑1 (MQ=17)
                                                                                                      aCCGGATGCCTTTCATCTTGCCCGAAATATCCTGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTg  >  2:3138641/1‑101 (MQ=18)
                                                                                                      aCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTg  <  1:2335769/101‑1 (MQ=17)
                                                                                                      |                                                                                                          
AGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.A.A.GGGGGGGCGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AGCAGA.AGGGGGGCAAGCG.CC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTG  >  NC_014034/2493915‑2494113

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