Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,486,442 +C 4.7% intergenic (‑70/+124) RCAP_rcc02312 ← / ← soxG hypothetical protein/sarcosine oxidase subunit gamma

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,486,4420.→C 4.7% 22.4 ‑0.3 ‑0.7192intergenic (‑70/+124)RCAP_rcc02312/soxGhypothetical protein/sarcosine oxidase subunit gamma
Reads supporting (aligned to +/- strand):  new base (5/4):  ref base (85/98):  total (90/102)
Fisher's exact test strand distribution p-value = 7.37e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.82e-01

Only 100 of 201 total aligned reads displayed.

GCGCCCCCGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAA.GGG.G.CGCTGCCCCTCTCTTG.CCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTCCGCAG  >  NC_014034/2486347‑2486538
                                                                                                                                                                                                    
gcgcCCCCGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCtt                                                                                              <  2:3781877/101‑1 (MQ=255)
  gcCCCCGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTtg.c                                                                                           >  2:50842/1‑101 (MQ=255)
    ccccGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTtg.ccc                                                                                         >  2:2169043/1‑100 (MQ=255)
     cccGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTtg.ccct                                                                                        >  1:173204/1‑99 (MQ=255)
      ccGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTtg.ccctg                                                                                       <  2:2227102/101‑4 (MQ=255)
       cGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTtg.ccctgc                                                                                      <  2:2126269/101‑5 (MQ=255)
       cGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTtg.ccctgc                                                                                      >  1:707328/1‑97 (MQ=255)
        ggATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTtg.ccctgcc                                                                                     >  2:232025/1‑96 (MQ=255)
          aTTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCgg                                                                                    >  1:4100010/1‑101 (MQ=255)
           tttCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCggg                                                                                   <  2:4397730/101‑1 (MQ=255)
           tttCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAA.GGGCG.CGCTGCCCC.CTCTTGCCCTGCCggg                                                                                   >  1:2803404/1‑101 (MQ=255)
            ttCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCGGCCCC.CTCTTGCCCTGCCGGGc                                                                                  >  1:787790/1‑101 (MQ=255)
              cGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCaa                                                                                >  1:2190598/1‑101 (MQ=38)
              cGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCaa                                                                                >  1:4439632/1‑101 (MQ=38)
               gAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAg                                                                               >  1:3133686/1‑101 (MQ=38)
               gAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAg                                                                               >  2:1680282/1‑101 (MQ=38)
               gAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAg                                                                               >  1:4242559/1‑101 (MQ=38)
                 aTAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTc                                                                             >  2:1042529/1‑101 (MQ=38)
                 aTAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTc                                                                             <  1:2450844/101‑1 (MQ=38)
                 aTAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTc                                                                             <  1:154904/101‑1 (MQ=38)
                   aGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCGGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCAc                                                                           >  2:2377320/1‑101 (MQ=37)
                   aGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCAc                                                                           <  1:3453116/101‑1 (MQ=38)
                     tCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCAccc                                                                         <  1:3732450/101‑1 (MQ=37)
                        aTGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCg                                                                      <  2:2545654/101‑1 (MQ=37)
                        aTGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCg                                                                      <  2:263945/101‑1 (MQ=37)
                         tGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCgg                                                                     >  2:3765163/1‑101 (MQ=37)
                          gCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCggg                                                                    >  2:362349/1‑101 (MQ=37)
                          gCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCggg                                                                    <  2:1235027/101‑1 (MQ=37)
                          gCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCggg                                                                    <  1:1321418/101‑1 (MQ=37)
                           cTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGa                                                                   >  2:3455030/1‑101 (MQ=37)
                            tGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGat                                                                  <  2:2204311/101‑1 (MQ=37)
                            tGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGat                                                                  >  1:3810140/1‑101 (MQ=37)
                            tGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGat                                                                  <  1:2985294/101‑1 (MQ=37)
                             ggATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGata                                                                 >  1:3043526/1‑101 (MQ=37)
                               aTCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATAtt                                                               <  1:2916392/101‑1 (MQ=37)
                                tctcCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTc                                                              <  1:189189/101‑1 (MQ=37)
                                tctcCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAA.GGG.GACGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCAACACCGGGATAtta                                                              >  2:1168411/1‑100 (MQ=37)
                                 ctcCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCa                                                             <  1:434432/101‑1 (MQ=37)
                                     tttGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGCGATTCATGGc                                                         >  2:1863889/1‑101 (MQ=25)
                                     tttGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTTATGGc                                                         >  1:4640250/1‑101 (MQ=25)
                                       tGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCaa                                                       <  2:3150516/101‑1 (MQ=37)
                                        gCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAg                                                      <  1:2171322/101‑1 (MQ=37)
                                          cTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGAt                                                    <  1:4524187/101‑1 (MQ=37)
                                           ttCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGAta                                                   >  2:3174661/1‑100 (MQ=37)
                                             cGATCTGCGGCAGGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGaa                                                 <  2:4099881/101‑1 (MQ=25)
                                             cGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGaa                                                 <  2:4518351/101‑1 (MQ=37)
                                              gATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAag                                                >  1:2074535/1‑101 (MQ=37)
                                               aTCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAagg                                               <  2:1248628/101‑1 (MQ=37)
                                               aTCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAagg                                               <  2:1713300/101‑1 (MQ=37)
                                                tCTGCGGCAAGCCCGGCAGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAagga                                              >  2:1352992/1‑101 (MQ=25)
                                                gcTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAagga                                              <  2:1379741/100‑1 (MQ=25)
                                                 cTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAaggag                                             <  1:3100335/101‑1 (MQ=37)
                                                   gCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGc                                           <  1:863636/101‑1 (MQ=37)
                                                   gCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGACGAAGGAGGc                                           >  1:4406185/1‑101 (MQ=25)
                                                    cGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCg                                          >  2:2431851/1‑101 (MQ=25)
                                                     ggCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCgt                                         <  2:2517622/101‑1 (MQ=37)
                                                      gCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCgtg                                        >  2:3260924/1‑101 (MQ=25)
                                                       cAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCgtgt                                       >  1:3898138/1‑101 (MQ=25)
                                                        aaGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGtt                                      <  2:111754/101‑1 (MQ=37)
                                                        aaGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGtt                                      >  1:254223/1‑101 (MQ=25)
                                                        aaGCCCGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGtt                                      >  1:3654398/1‑101 (MQ=25)
                                                           cccGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCa                                   <  1:2670098/101‑1 (MQ=37)
                                                           cccGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCa                                   <  2:329890/101‑1 (MQ=37)
                                                            ccGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCAt                                  <  2:1743066/101‑1 (MQ=25)
                                                            ccGGCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCAt                                  <  2:2803188/101‑1 (MQ=25)
                                                              ggCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCt                                >  2:1409368/1‑101 (MQ=25)
                                                              ggCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCt                                >  1:1671972/1‑101 (MQ=25)
                                                               gCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCtt                               <  1:1980023/101‑1 (MQ=21)
                                                               gCCGCAAGGCAAGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATAtt                               >  1:4457612/1‑101 (MQ=21)
                                                                ccGCAAGGCACGCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTg                              <  2:1130919/101‑1 (MQ=21)
                                                                  gcaaggcaagCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGcc                            >  1:2260654/1‑101 (MQ=21)
                                                                    aaggcaagCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCAt                          >  1:3150087/1‑101 (MQ=21)
                                                                    aaggcaagCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCAt                          <  2:2295345/101‑1 (MQ=21)
                                                                    aaggcaagCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCAt                          <  1:3722181/101‑1 (MQ=25)
                                                                    aaggcaagCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCAt                          >  2:1191268/1‑101 (MQ=21)
                                                                    aaggcaagCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGGGTTTCATCTTGCCAt                          >  1:1167005/1‑101 (MQ=21)
                                                                     aggcaagCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATa                         <  1:1409197/101‑1 (MQ=25)
                                                                     aggcaagCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATa                         <  1:3231394/101‑1 (MQ=25)
                                                                     aggcaagCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATa                         >  1:3399976/1‑101 (MQ=21)
                                                                     aggcaagCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATAGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATa                         >  1:1144593/1‑101 (MQ=21)
                                                                        caagCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAat                      >  1:3387021/1‑101 (MQ=21)
                                                                        caagCGAAAAGGGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAat                      <  1:4155558/101‑1 (MQ=21)
                                                                            gaaaag.GGG.G.CGCTGCTCC.CTCTTGCCCTGCCGGGCAAGTCGCCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATccc                 <  2:3703747/99‑1 (MQ=21)
                                                                            gaaaag.GGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATccc                 <  1:50848/99‑1 (MQ=21)
                                                                             aaaag.GGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGCACATTCATTGCAAGATGAAGAGAGCGTGTTTCATCTGGTCAAAAttataacg                >  2:3219497/2‑97 (MQ=12)
                                                                              aaag.GGG.G.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCgg               >  1:1002947/1‑101 (MQ=21)
                                                                               aaaggtg.g.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCgg               <  2:444468/95‑1 (MQ=21)
                                                                                 agggg.g.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCgggg             <  2:4242271/99‑1 (MQ=21)
                                                                                  gggg.g.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCggggg            >  2:2910855/2‑101 (MQ=21)
                                                                                   ggg.g.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCgggggg           >  2:1587474/1‑101 (MQ=255)
                                                                                   ggg.g.CGCTGCCCC.CTCGTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCgggggg           <  2:4194864/101‑1 (MQ=21)
                                                                                    gg.g.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGa          <  1:4043239/101‑1 (MQ=255)
                                                                                    gg.g.CGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGa          <  2:775884/101‑1 (MQ=255)
                                                                                       g.cgcTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGt        <  2:4018386/101‑1 (MQ=255)
                                                                                         cgcTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTc       <  1:1617886/101‑1 (MQ=255)
                                                                                          gcTGCCCC.CTCTTGCCCTGCCGGGCAAGTCCCCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTcc      <  1:632316/101‑1 (MQ=21)
                                                                                             tgcccCCTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTCCGc    <  2:2081338/99‑1 (MQ=255)
                                                                                             tgcccCCTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTCCGc    <  1:2113035/99‑1 (MQ=255)
                                                                                              gcccCCTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTCCGCa   <  2:1813326/100‑1 (MQ=255)
                                                                                              gcccCCTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTCCGCa   >  2:1532199/2‑101 (MQ=255)
                                                                                                                                                                                                    
GCGCCCCCGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAA.GGG.G.CGCTGCCCCTCTCTTG.CCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTCCGCAG  >  NC_014034/2486347‑2486538

Base quality scores:  ATCG  < 3 ≤  ATCG  < 21 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG