Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 3,504,126 +G 1.9% intergenic (‑19/+64) recO ← / ← RCAP_rcc03309 DNA repair protein RecO/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140343,504,1261.→G 1.9% 3.0 ‑3.4 ‑0.0217intergenic (‑19/+64)recO/RCAP_rcc03309DNA repair protein RecO/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (3/0):  ref base (67/145):  total (72/145)
Fisher's exact test strand distribution p-value = 3.35e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.15e-01

Only 100 of 267 total aligned reads displayed.

TGGCCGTGAACACCTCGACAATCGCCGAGCTTTCGCCGTGCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGTCTTTCAGTCCGTCTTCATCCAGAAGCGTGCGG  >  NC_014034/3504026‑3504226
                                                                                                       |                                                                                                     
tGGCCGTGAACACCTCGACAATCGCCGAGCTTTCGCCGTGCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.cc                                                                                                        <  1:514464/101‑1 (MQ=255)
 ggCCGTGAACACCTCGACAATCGCCGAGCTTTCGCCGTGCAGGCGGGTCGAAAGGATGGGGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.GC..g                                                                                                     >  1:4663670/1‑101 (MQ=255)
  gCCGTGAACACCTCGACAATCGCCGAGCTTTCGCCGTGCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..gg                                                                                                    <  1:415312/101‑1 (MQ=255)
   ccGTGAACACCTCGACAATCGCCGAGCTTTCGCCGTGCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.GC..ggg                                                                                                   >  2:1595877/1‑101 (MQ=255)
   ccGTGAACACCTCGACAATCGCCGAGCTTTCGCCGTGCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..ggg                                                                                                   >  2:210852/1‑101 (MQ=255)
    cGTGAACACCTCGACAATCGCCGAGCTTTCGCCGTGCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..gggg                                                                                                  <  2:2478994/101‑1 (MQ=255)
    cGTGAACACCTCGACAATCGCCGAGCTTTCGCCGTGCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..cggg                                                                                                  >  2:2994520/1‑101 (MQ=255)
      tGAACACCTCGACAATCGCCGAGCTTTCGCCGTGCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAAGGG.CC..ggggg                                                                                                 >  1:3918237/1‑101 (MQ=255)
      tGAACACCTCGACAATCGCCGAGCTTTCGCCGTGCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..gggggg                                                                                                >  2:685121/1‑100 (MQ=255)
        tacacCTCGACAAGCGCCGAGCTTTCGCCGTGCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGttt                                                                                              <  1:3049413/100‑1 (MQ=255)
         acacCTCGACAATCGCCGAGCTTTCGCCGTGCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGtttt                                                                                             <  2:456815/101‑1 (MQ=255)
              tCGACAATCGCCGAGCTTTCGCCGTGCAGGCGGGACGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGGTTTGGAcc                                                                                        >  2:3137932/1‑101 (MQ=255)
                 aCAATCGCCGAGCTTTCGCCGTGCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCAccccc                                                                                     <  2:874022/101‑1 (MQ=255)
                 aCAATCGCCGAGCTTTCGCCGTGCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCAccccc                                                                                     <  1:231738/101‑1 (MQ=255)
                  cAATCGCCGAGCTTTCGCCGTGCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.TG.CC..GGGGGTTTTGCACCCCCg                                                                                    <  1:2701829/101‑1 (MQ=255)
                   aaTCGCCGAGCTTTCGCCGTGCAGGCGGGGCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.GC..GGGGGGTTTGCACCCCCgg                                                                                   >  2:2157895/1‑101 (MQ=255)
                    aTCGCCGAGCTTTCGCCGTGCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCggg                                                                                  <  2:4387370/101‑1 (MQ=255)
                        ccGAGCTTTCGCCGTGCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.GC..GGGGGGTTTGCACCCCCGGGcccc                                                                              >  1:2087389/1‑101 (MQ=255)
                        ccGAGCTTTCGCCGTGCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGcccc                                                                              <  1:1678590/101‑1 (MQ=255)
                        ccGAGCTTTCGCCGTGCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGGTTTGCACCCCCGGGcccc                                                                              >  1:2961169/1‑101 (MQ=255)
                         cGAGCTTTCGCCGTGCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGTccccc                                                                             <  1:3138132/101‑1 (MQ=255)
                         cGAGCTTTCGCCGTGCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGccccc                                                                             <  2:2721726/101‑1 (MQ=255)
                          gAGCTTTCGCCGTGCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGTAACCCCCGGCCCCCg                                                                            >  2:1282900/1‑101 (MQ=255)
                          gAGCTTTCGCCGTGCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCg                                                                            <  2:1679548/101‑1 (MQ=255)
                           aGCTTTCGCCGTGCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGc                                                                           <  2:1756378/101‑1 (MQ=255)
                           aGCTTTCGCCGTGCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGc                                                                           <  1:483550/101‑1 (MQ=255)
                            gCTTTCGCCGTGCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAAGGG.CC..GGGGGTTTTGACCCCCCGGGccccccc                                                                           >  1:3384337/1‑99 (MQ=255)
                                 cGCCGTGCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGata                                                                     <  1:4177338/101‑1 (MQ=255)
                                  gCCGTGCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGatat                                                                    <  1:2904487/101‑1 (MQ=255)
                                   ccGTGCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATAtt                                                                   <  2:734209/101‑1 (MQ=255)
                                     gTGCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTg                                                                 <  1:567712/101‑1 (MQ=255)
                                     gTGCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTg                                                                 <  1:3241129/101‑1 (MQ=255)
                                      tGCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGa                                                                <  2:2647978/101‑1 (MQ=255)
                                      tGCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGa                                                                <  2:3908486/101‑1 (MQ=255)
                                       gCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAg                                                               <  1:3722673/101‑1 (MQ=255)
                                       gCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAg                                                               <  2:3912992/101‑1 (MQ=255)
                                         aGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGGTTTGCACCCCCGGGCCCCCGCCGGATATTTGAgga                                                             >  1:1476514/1‑100 (MQ=255)
                                          ggCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCa                                                            <  1:221672/101‑1 (MQ=255)
                                           gCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCaa                                                           <  2:1559822/101‑1 (MQ=255)
                                            cGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAg                                                          <  1:1149085/101‑1 (MQ=255)
                                             gggTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGTTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGa                                                         <  2:4318179/101‑1 (MQ=255)
                                              ggTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGAt                                                        <  2:61472/101‑1 (MQ=255)
                                               gTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATa                                                       <  1:810717/101‑1 (MQ=255)
                                               gTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATa                                                       <  1:652438/101‑1 (MQ=255)
                                               gTCGAAAGGATGGTGCCGCTGCCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CT..GTGGGTTTTGCCCCCCCCGGCCCCCCCACGGCGTTTGAGGGAAGATa                                                       >  2:4475234/1‑101 (MQ=255)
                                                tCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGGTTTTCACCCCCCGGCCCCCGCAGGGTATTTGAGGCAAGATaa                                                      >  1:4285020/1‑101 (MQ=255)
                                                 cGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGGTTTGCAACCCCCGGCCCCCCCAGGATATTTTAGGCAAGATaaa                                                     >  1:1036740/1‑101 (MQ=255)
                                                   aaaGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGGTTTGCACCCCCCGGCCCCCGCCGGATATTTGAGGCAAGATaaaaa                                                   >  2:2623981/1‑100 (MQ=255)
                                                     aGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGGTTTGCACCCCCGGGCCCCCCCAGGGTATTTTAGGCAAGATAAAAggg                                                 >  2:4133393/1‑99 (MQ=255)
                                                      ggATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GGCCC..GGGGGGTTTGCCCCCCCGGGCCCCC.CCCGATCTTTTACGCAAGAAaaaaagga                                                >  1:4364705/1‑97 (MQ=255)
                                                       gATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGaa                                               >  2:2369378/1‑101 (MQ=255)
                                                       gATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGaa                                               <  1:403702/101‑1 (MQ=255)
                                                       gATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGaa                                               <  1:1544992/101‑1 (MQ=255)
                                                       gATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGGTTTGCACCCCCGGGCCCCCCCAGGATATTTGAGGCAAGATAAAAGAGaa                                               >  1:763101/1‑101 (MQ=255)
                                                        aTGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGaaa                                              <  2:1446934/101‑1 (MQ=255)
                                                        aTGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGaaa                                              <  2:2796997/101‑1 (MQ=255)
                                                         tGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAg                                             <  1:1282991/101‑1 (MQ=255)
                                                          ggTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..CGGGGGTTTGCACCCCCGGGCCCCCGCAGGATATTTTAGGCAAGATAAAAGCGAAAGc                                            >  2:3735052/1‑101 (MQ=255)
                                                           gTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CCG.GGGGGGTTTGACCCCCCGGGCCCCCCCGGGGACTTTGGGGCAAAAAaaaaagaaaaac                                            >  1:4400896/1‑93 (MQ=255)
                                                            tGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCg                                          <  2:3475617/101‑1 (MQ=255)
                                                            tGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCg                                          <  1:3551517/101‑1 (MQ=255)
                                                                 cTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTc                                     <  2:521707/101‑1 (MQ=255)
                                                                 cTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTc                                     >  2:3407945/1‑101 (MQ=255)
                                                                   gTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CCG.GGGGGTTTTGACCCCCCGGGCCCCC.CCGGGTATTTTGGGCCAGAAAAAAAAGAAAACCcggggccg                                   >  2:198177/1‑97 (MQ=255)
                                                                    tCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CCCGGGGGGTTTTTCCCCCCCGCCCCCCCCCGGGGTTTTTGGGGAAAAAAAAAAGAAAAACCGggggcgg                                    >  2:902500/1‑98 (MQ=255)
                                                                      ttGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGTCt                                >  1:2136140/1‑101 (MQ=255)
                                                                      ttGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGTCt                                <  1:1442402/101‑1 (MQ=255)
                                                                       tGCCAGTCCATGGTGCCCGTGACGAA.GGCCC..GGGGGTTTTTGCCCCCCGGGCCCCCCCCGGAAATTTTGGGGAAAAAAAAAAGAAAACCCgggggtggcct                                >  1:3796413/1‑94 (MQ=255)
                                                                       tGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGTCtt                               <  2:3249247/101‑1 (MQ=255)
                                                                          cAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGTCTTTCa                            <  1:4382080/101‑1 (MQ=255)
                                                                          cAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGACTTTCa                            >  1:3849301/1‑101 (MQ=255)
                                                                           aGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGTCTTTCAg                           <  2:289635/101‑1 (MQ=255)
                                                                           aGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGTCTTTCAg                           <  1:4547637/101‑1 (MQ=255)
                                                                            gTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGTCTTTCAGt                          <  2:1816455/101‑1 (MQ=255)
                                                                             tCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGTCTTTCAGTc                         <  2:530357/101‑1 (MQ=255)
                                                                              ccATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGTCTTTCAGTcc                        <  2:3599392/101‑1 (MQ=255)
                                                                                  ggTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGTCTTTCAGTCCGTCt                    <  2:2036639/101‑1 (MQ=255)
                                                                                  ggTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGTCTTTCAGTCCGTCt                    <  2:3831259/101‑1 (MQ=255)
                                                                                  ggTGCCCGTGACGAA.GG.CC..GGGGGGTTTGCACCCCCGGGCCCCCGCAGGGTATTTGAGGCAAGATAAAAGCGAAAGCCGGGGGCAGTCTTTCAGTCCGTCt                    >  1:2836383/1‑101 (MQ=255)
                                                                                     gCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGTCTTTCAGTCCGTCTTCa                 <  1:2417179/101‑1 (MQ=255)
                                                                                     gCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGTCTTTCAGTCCGTCTTCa                 <  2:3706454/101‑1 (MQ=255)
                                                                                     gCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGTCTTTCAGTCCGTCTTCa                 <  1:3079950/101‑1 (MQ=255)
                                                                                     gCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGTCTTTCAGTCCGTCTTCa                 <  2:1018551/101‑1 (MQ=255)
                                                                                     gCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGTCTTTCAGTCCGTCTTCa                 <  1:3447572/101‑1 (MQ=255)
                                                                                      cccGTGACGAA.GG.CC..GGGGGTTGTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGTCTTTCAGTCCGTCTTCAt                <  2:2176332/101‑1 (MQ=255)
                                                                                         gTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGTCTTTCAGTCCGTCTTCATCCa             <  2:4337689/101‑1 (MQ=255)
                                                                                           gACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGTCTTTCAGTCCGTCTTCATCCAGa           <  1:4580845/101‑1 (MQ=255)
                                                                                             cGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGTCTTTCAGTCCGTCTTCATCCAGAAg         >  2:2284821/1‑101 (MQ=255)
                                                                                              gAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGTCTTTCAGTCCGTCTTCATCCAGAAGc        <  2:1474435/101‑1 (MQ=255)
                                                                                              gAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGTCTTTCAGTCCGTCTTCATCCAGAAGc        <  2:1218500/101‑1 (MQ=255)
                                                                                              gAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGTCTTTCAGTCCGTCTTCATCCAGAAGc        <  2:1887063/101‑1 (MQ=255)
                                                                                               aa.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGTCTTTCAGTCCGTCTTCATCCAGAAGCg       >  2:2851780/1‑101 (MQ=255)
                                                                                                a.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGTCTTTCAGTCCGTCTTCATCCAGAAGCgg      >  1:1540575/1‑100 (MQ=255)
                                                                                                a.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGTCTTTCAGTCCGTCTTCATCCAGAAGCGt      <  2:1403777/101‑1 (MQ=255)
                                                                                                a.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGTCTTTCAGTCCGTCTTCATCCAGAAGCGt      <  2:4635738/101‑1 (MQ=255)
                                                                                                  gg.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGTCTTTCAGTCCGTCTTCATCCAGAAGCGTg     <  1:2646177/101‑1 (MQ=255)
                                                                                                   g.CC..GGGGGTTTTGCCCCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGTCTTTCAGTCCGTCTTCATCCAGAAGCGTGc    <  1:2369511/101‑1 (MQ=255)
                                                                                                   g.CC..GGGGGTTTTGCCCCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGTCTTTCAGTCCGTCTTCATCCAGAAGCGTGc    <  2:3090670/101‑1 (MQ=255)
                                                                                                     cc..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGTCTTTCAGTCCGTCTTCATCCAGAAGCGTGCg   <  1:1842266/101‑1 (MQ=255)
                                                                                                      c..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGTCTTTCAGTCCGTCTTCATCCAGAAGCGTGCgg  <  1:4213757/101‑1 (MQ=255)
                                                                                                       |                                                                                                     
TGGCCGTGAACACCTCGACAATCGCCGAGCTTTCGCCGTGCAGGCGGGTCGAAAGGATGGTGCCGCTGTCTTGCCAGTCCATGGTGCCCGTGACGAA.GG.CC..GGGGGTTTTGCACCCCCGGGCCCCCGCAGGATATTTGAGGCAAGATAAAAGCGAAAGCCGGGGTCAGTCTTTCAGTCCGTCTTCATCCAGAAGCGTGCGG  >  NC_014034/3504026‑3504226

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 27 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG