Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,494,003 A→G 6.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 6.2% 16.8 ‑0.3 ‑0.1163intergenic (+51/‑191)RCAP_rcc02319/RCAP_rcc02320hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (4/6):  ref base (47/105):  total (52/111)
Fisher's exact test strand distribution p-value = 7.26e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.53e-01

Only 100 of 187 total aligned reads displayed.

CGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AA.GGGGGGGCGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AG.CAGA.AGGGGGGCAAGCG.CC.CCCTGCCCCGGCGGTCCGTG  >  NC_014034/2493905‑2494102
                                                                                                    |                                                                                                           
 gTCTGGCCCCGCTCCGCGGCGGGCGGGCGGCACGGCAGGCGGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTgcagc                                                                                                           <  1:1836389/101‑1 (MQ=18)
 gTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTgcagc                                                                                                           >  1:1110294/1‑101 (MQ=25)
  tCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTgcagca                                                                                                          <  1:3238582/101‑1 (MQ=25)
   cTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCaa                                                                                                         <  2:638159/101‑1 (MQ=25)
    tGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAg                                                                                                        <  1:1683569/101‑1 (MQ=21)
    tGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGGGCGCCGCCCCCCCGCCGCCAGCGGGCTCCCCCCCGGCTATTTgcacccag                                                                                                        >  2:1093764/1‑96 (MQ=255)
        ccAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGGTATTTGCCGCAAGATGa                                                                                                    >  1:3156763/1‑101 (MQ=21)
        ccAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGa                                                                                                    <  2:596979/101‑1 (MQ=21)
         cAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGGNAGTTGCCGCAATATgga                                                                                                   >  2:1594786/1‑99 (MQ=18)
         cAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGaa                                                                                                   <  1:1466599/101‑1 (MQ=21)
         cAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGaa                                                                                                   <  1:287907/101‑1 (MQ=21)
         cAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGaa                                                                                                   <  2:2714774/101‑1 (MQ=21)
          aGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGaaa                                                                                                  <  2:3113942/101‑1 (MQ=21)
            cgccgcGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGGTATTTGCAGCCAAATTAAAcc                                                                                                >  1:4142835/1‑101 (MQ=18)
            cgccgcGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAAcc                                                                                                <  1:1177442/101‑1 (MQ=21)
            cgccgcGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAAcc                                                                                                <  2:1418470/101‑1 (MQ=21)
            cgccgcGCCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAAcc                                                                                                <  2:2107580/101‑1 (MQ=18)
                gcggcgGGCGCGCGGGACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAt                                                                                            <  2:2440153/101‑1 (MQ=18)
                gcggcgGGCGCGCGGCACGACAGGCCGAGTAAA.GGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGag                                                                                            >  2:1337476/1‑100 (MQ=18)
                gcggcgGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAt                                                                                            <  2:112598/101‑1 (MQ=21)
                  ggcgggcgCGCGGAACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGc                                                                                          >  2:1723956/1‑101 (MQ=18)
                   gcgggcgCGCGGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGGCCCCCGTCCGCAAGCGGACTCCCCCCCGGATATTTGCAGCAAGATGaaaaacgaggtc                                                                                         >  1:3065595/1‑92 (MQ=17)
                   gcgggcgCGCGGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGcc                                                                                         <  1:1811761/101‑1 (MQ=21)
                    cgggcgCGCGGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCt                                                                                        <  2:1168021/101‑1 (MQ=21)
                     gggcgCGCGGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCtt                                                                                       >  2:4365982/1‑101 (MQ=21)
                      ggcgCGCGGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGGCATTTGCAGCAAGATGAAACCGGTTGCCttt                                                                                      >  2:1014575/1‑101 (MQ=18)
                        ggcgcgGCGCGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCa                                                                                    <  2:3910641/100‑1 (MQ=18)
                        cgcgcgGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCa                                                                                    <  1:2563115/101‑1 (MQ=21)
                          cgcgGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATc                                                                                  <  2:420997/101‑1 (MQ=21)
                          cgcgGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATc                                                                                  >  1:1322135/1‑101 (MQ=21)
                           gcgGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCt                                                                                 <  1:602055/101‑1 (MQ=21)
                           gcgGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCt                                                                                 <  1:736814/101‑1 (MQ=21)
                            cgGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCtt                                                                                <  2:3551263/101‑1 (MQ=21)
                            cgGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCtt                                                                                <  1:1235454/101‑1 (MQ=21)
                             ggCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTg                                                                               <  2:4069210/101‑1 (MQ=21)
                             ggCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTg                                                                               <  2:2326949/101‑1 (MQ=21)
                              gCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGc                                                                              >  2:2713331/1‑101 (MQ=21)
                                aCGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGccc                                                                            >  2:1249258/1‑101 (MQ=21)
                                aCGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGccc                                                                            <  2:3086890/101‑1 (MQ=21)
                                aCGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGccc                                                                            <  1:1075202/101‑1 (MQ=21)
                                 cGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCg                                                                           <  2:2359697/101‑1 (MQ=21)
                                 cGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCg                                                                           <  2:2992867/101‑1 (MQ=21)
                                  gACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGa                                                                          >  1:2348065/1‑101 (MQ=21)
                                   aCAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGaa                                                                         <  2:3266273/101‑1 (MQ=21)
                                     aGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCGAAATGAAACCGGGTGCCTTTCATCTTGCCCGAAat                                                                       >  1:3626024/1‑101 (MQ=18)
                                      ggCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAata                                                                      <  2:1477498/101‑1 (MQ=21)
                                       gCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAatat                                                                     <  1:2503497/101‑1 (MQ=21)
                                       gCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAatat                                                                     <  1:1991105/101‑1 (MQ=21)
                                       atctaGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAatat                                                                     <  1:4384874/97‑1 (MQ=18)
                                        ccGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATc                                                                    <  2:183852/101‑1 (MQ=18)
                                         cgggtaagggggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATc                                                                    <  1:3026218/92‑1 (MQ=18)
                                          gAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATccc                                                                  <  1:3869050/101‑1 (MQ=18)
                                             gtaaggggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgg                                                                <  1:2456452/96‑1 (MQ=255)
                                              taaggggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCggg                                                               <  1:72196/97‑1 (MQ=255)
                                                      gggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.t                                                      >  1:1378452/1‑101 (MQ=255)
                                                       ggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.Tg                                                     <  2:257003/101‑1 (MQ=255)
                                                        gcgcTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGACATTTGCAGCAAGATGAAAACGGATGCCTTTCATCTTGCCCGAACTATACCGGGGGGGGAG.Tgg                                                    >  1:883252/1‑101 (MQ=18)
                                                         cgcTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGc                                                   <  2:1030138/101‑1 (MQ=255)
                                                            tGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aa                                             <  1:2144125/101‑2 (MQ=255)
                                                            tGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aa                                             <  1:3568460/101‑2 (MQ=255)
                                                             gCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aag                                            <  1:1081080/101‑3 (MQ=255)
                                                              cccccG.TCCGCAAGCGGACTCCCCCCGGGATATTTGCCGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGAGAA.TGGC..G.aagcc                                          >  1:2946572/1‑97 (MQ=18)
                                                              cccccCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagc                                           <  1:2808866/101‑4 (MQ=255)
                                                                 cccgtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccag                                        <  1:3438156/97‑7 (MQ=255)
                                                                 cccgtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccag                                        <  1:769318/97‑7 (MQ=255)
                                                                   cgtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccagag                                      <  1:573662/99‑9 (MQ=255)
                                                                    gtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..G.aagccagagg                                     >  1:1827773/2‑92 (MQ=21)
                                                                     tCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccagaggg                                    >  2:2511595/1‑91 (MQ=255)
                                                                     tCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccagaggg                                    <  1:1724067/101‑11 (MQ=255)
                                                                          aaGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGCa                              >  2:2302328/1‑101 (MQ=18)
                                                                          aaGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..GAAG.CCAG.AGGGGGGCa                              >  1:4364755/1‑101 (MQ=18)
                                                                           aGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGCag                             <  2:4028881/101‑2 (MQ=18)
                                                                           aGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGCag                             <  2:413502/101‑2 (MQ=18)
                                                                           aGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGCag                             <  1:2028967/101‑2 (MQ=18)
                                                                                 cTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCG.cc.cc                    >  2:688702/1‑101 (MQ=17)
                                                                                  tCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGCAATATCCCGGGGGGGGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCG.cc.ccc                   >  1:2395499/1‑101 (MQ=17)
                                                                                  tCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG..AGGGGGGC.AGCG.cc.ccc                   <  2:4142574/101‑1 (MQ=17)
                                                                                   ccccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCG.cc.cccc                  >  1:2365263/1‑100 (MQ=17)
                                                                                   ccccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCG.cc.cccc                  <  2:1952332/101‑2 (MQ=17)
                                                                                    cccccGGGATATTTGCAGCAAGATGAAACCGGATCCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCG.cc.cccct                 <  1:19700/101‑3 (MQ=17)
                                                                                     ccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAT.TGGC..GAAG.CCGG.GGGGGGGCGGGCC.cc.ccccc                 >  2:3136668/1‑99 (MQ=17)
                                                                                     ccccGGGATATGTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCG.cc.cccctg                <  2:1322038/101‑4 (MQ=17)
                                                                                        cGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCGACC.CCCTGccc              >  2:1214419/1‑101 (MQ=17)
                                                                                         gggATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCA.AGAGGGGGCAGCG.CC.CCCTGcccc             >  1:2784255/1‑101 (MQ=0)
                                                                                         gggATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCA.AGAGGGGGCAGCG.CC.CCCTGcccc             <  1:2186554/101‑1 (MQ=0)
                                                                                          ggATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCA.AGAGGGGGCAGCG.CC.CCCTGCCCcg            >  2:1447132/1‑101 (MQ=0)
                                                                                          ggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGAGAA.GGGC..GAAG.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCcg            >  1:4009471/1‑101 (MQ=17)
                                                                                          ggATANNTGCAGAAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCcg            <  1:1594892/101‑1 (MQ=17)
                                                                                           gATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCcgg           >  1:2613171/1‑101 (MQ=17)
                                                                                            atatTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCA.AGAGGGGGCAGCG.CC.CCCTGCCCcgcc          <  2:2819530/101‑3 (MQ=2)
                                                                                             tatTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AA.GCCA.GAGGGGGC.AGCGCCC.CCCTGCCCcggcgg        >  2:225289/1‑101 (MQ=17)
                                                                                              atTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCA.AGAGGGGGCAGCG.CC.CCCTGCCCcgccgg        <  2:337533/101‑1 (MQ=2)
                                                                                              atTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AA.GACAGAGGGGGGC.AGCGCCC.CCCTGCCCcggcgg        >  1:477110/1‑101 (MQ=18)
                                                                                              anNTGCNGCAAGATGAANGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCA.AGAGGGGGCAGCG.CC.CCCTGCCCcgccgg        <  2:1707762/101‑1 (MQ=2)
                                                                                                ttGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCA.AGAGGGGGCAGCG.CC.CCCTGCCCcgccggat      >  2:91130/1‑99 (MQ=2)
                                                                                                ttGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTc      <  2:3706895/101‑1 (MQ=17)
                                                                                                 tGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCG.CCGCCCTGCCCCGGCGGTcc     >  2:1376656/1‑101 (MQ=18)
                                                                                                  gcagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCg    >  2:3410585/1‑101 (MQ=18)
                                                                                                   cagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGTGC.AGCGCCC.CCCTGCCCCGGCGGTCCGt   >  1:1490453/1‑101 (MQ=18)
                                                                                                   cagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCCGCGGGCCg  >  2:3847046/1‑100 (MQ=14)
                                                                                                    |                                                                                                           
CGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AA.GGGGGGGCGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AG.CAGA.AGGGGGGCAAGCG.CC.CCCTGCCCCGGCGGTCCGTG  >  NC_014034/2493905‑2494102

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