Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,486,442 +C 3.1% 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 3.1% 39.0 ‑1.2 ‑1.6532intergenic (‑70/+124)RCAP_rcc02312/soxGhypothetical protein/sarcosine oxidase subunit gamma
Reads supporting (aligned to +/- strand):  new base (9/6):  ref base (234/281):  total (244/288)
Fisher's exact test strand distribution p-value = 3.01e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.73e-02

Only 100 of 538 total aligned reads displayed.

GCGCCCCCGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAA.GGG.GCGCTGCCCCTCTCTTG.CCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTCCGCAG  >  NC_014034/2486347‑2486538
                                                                                                                                                                                                   
gcgcCCCCGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCtt                                                                                              >  2:7110840/1‑101 (MQ=255)
gcgcCCCCGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCtt                                                                                              >  1:2658829/1‑101 (MQ=255)
ccccccccGGATTTCGAATCGTCCATGCCGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCtt                                                                                              <  1:8934509/98‑1 (MQ=255)
      ccGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTtg.ccctg                                                                                       >  2:5205381/1‑98 (MQ=255)
      ccGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTtg.ccctg                                                                                       <  2:1651223/101‑4 (MQ=255)
       cGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTtg.ccctgc                                                                                      <  2:6213627/101‑5 (MQ=255)
       cGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTtg.ccctgc                                                                                      >  1:6002175/1‑97 (MQ=255)
         gATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCg                                                                                     >  1:6345304/1‑101 (MQ=255)
          aTTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCgg                                                                                    >  1:7938334/1‑101 (MQ=255)
          aTTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCgg                                                                                    <  2:748458/101‑1 (MQ=255)
          aTTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCgg                                                                                    >  2:8035007/1‑101 (MQ=255)
           tttCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCggg                                                                                   <  2:34089/101‑1 (MQ=255)
           tttCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCggg                                                                                   >  1:1699085/1‑101 (MQ=255)
           tttCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCggg                                                                                   >  1:6893877/1‑101 (MQ=255)
               gAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAg                                                                               >  2:2908323/1‑101 (MQ=38)
               gAATAGTCCATGCTGGATCTCCTTAACCTTCGACTTGCGGCAACCCCGGCGGCAAGGCAAGCGAAAAGGGG.GCCCTG.CCCGCTCTTGCCCTGCCGGGCacg                                                                               >  2:3945080/1‑99 (MQ=255)
                aaTAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGt                                                                              <  2:7903644/101‑1 (MQ=38)
                 aTAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTc                                                                             <  2:3053887/101‑1 (MQ=38)
                  tAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTcc                                                                            >  2:3753327/1‑100 (MQ=38)
                  tAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTcc                                                                            >  1:6720615/1‑100 (MQ=38)
                  tAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCa                                                                            >  2:7651999/1‑101 (MQ=38)
                   aGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCAc                                                                           <  1:3895149/101‑1 (MQ=38)
                   aGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCAc                                                                           <  2:5710498/101‑1 (MQ=38)
                    gTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCAcc                                                                          >  2:7757597/1‑101 (MQ=38)
                    ctCCATGCTGCATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCAcc                                                                          <  1:7359761/100‑1 (MQ=37)
                        aTGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCg                                                                      <  2:6311662/101‑1 (MQ=37)
                          gCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCggg                                                                    >  1:6150329/1‑101 (MQ=37)
                            tGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGat                                                                  >  2:8892138/1‑101 (MQ=37)
                            tGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGat                                                                  <  2:1213702/101‑1 (MQ=37)
                            tGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGat                                                                  <  1:7052498/101‑1 (MQ=37)
                            tGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGat                                                                  >  1:1265085/1‑101 (MQ=37)
                             ggATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGata                                                                 <  2:8824826/101‑1 (MQ=37)
                             ggATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGata                                                                 >  1:5230019/1‑101 (MQ=37)
                               aTCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATAtt                                                               <  2:1828994/101‑1 (MQ=37)
                               aTCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATAtt                                                               <  2:3318163/101‑1 (MQ=37)
                                 ctcCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCa                                                             >  1:8186887/1‑101 (MQ=37)
                                 ctcCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCa                                                             >  2:3905768/1‑101 (MQ=37)
                                 ctcCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCa                                                             >  2:1350857/1‑101 (MQ=37)
                                 ctcCTTGGCCTGCGAGCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCa                                                             <  2:1952518/101‑1 (MQ=25)
                                  tcCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCAt                                                            <  2:8292099/101‑1 (MQ=37)
                                   ccTTTGCCTTCGGTCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGGTATTCATg                                                           >  2:3528153/1‑101 (MQ=25)
                                   ccTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATg                                                           >  2:4428154/1‑101 (MQ=37)
                                    ctgtGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATgg                                                          <  1:2243860/98‑1 (MQ=25)
                                    cTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATgg                                                          >  1:4886948/1‑101 (MQ=37)
                                      ttGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCa                                                        <  2:2658747/101‑1 (MQ=37)
                                      ttGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCa                                                        >  1:1434732/1‑101 (MQ=37)
                                         ccgccGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCTCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGa                                                     <  2:4130120/97‑1 (MQ=37)
                                         ccTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGa                                                     >  2:3124639/1‑101 (MQ=37)
                                           ttCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATg                                                   <  1:882622/101‑1 (MQ=37)
                                           ttCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATg                                                   >  1:7490363/1‑101 (MQ=37)
                                           ttCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATg                                                   >  1:4530554/1‑101 (MQ=37)
                                             cGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGaa                                                 <  1:542733/101‑1 (MQ=37)
                                             cGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGaa                                                 <  1:399978/101‑1 (MQ=37)
                                              gATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAag                                                <  1:5094514/101‑1 (MQ=37)
                                               aTCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAagg                                               <  1:9039153/101‑1 (MQ=37)
                                               aTCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAagg                                               <  1:1898924/101‑1 (MQ=37)
                                                         aGCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGttt                                     <  1:1023667/101‑1 (MQ=37)
                                                          gCCCGGCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTc                                    <  1:1808862/101‑1 (MQ=25)
                                                            ccGGCCGCAGGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCAt                                  >  1:7166334/1‑101 (MQ=25)
                                                              ggCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCt                                <  1:5893924/101‑1 (MQ=37)
                                                              ggCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCt                                >  2:5895066/1‑101 (MQ=25)
                                                               gCCGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCtt                               <  2:6836327/101‑1 (MQ=25)
                                                                ccGCAAGGCAAGCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTg                              >  1:8410909/1‑101 (MQ=25)
                                                                  gcaaggcaagCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGcc                            >  1:862530/1‑101 (MQ=21)
                                                                   caaggcaagCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCa                           >  1:4168135/1‑101 (MQ=21)
                                                                   caaggcaagCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCa                           >  1:3304209/1‑101 (MQ=21)
                                                                    aaggcaagCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCAt                          <  2:2532415/101‑1 (MQ=25)
                                                                    aaggcaagCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCAt                          >  1:2521383/1‑101 (MQ=21)
                                                                     aggcaagCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATa                         <  2:4616182/101‑1 (MQ=25)
                                                                      ggcaagCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATaa                        <  1:2908421/101‑1 (MQ=25)
                                                                      ggcaagCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATaa                        >  2:819701/1‑101 (MQ=21)
                                                                      ggcaagCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATaa                        <  2:6581890/101‑1 (MQ=25)
                                                                       gcaagCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATaaa                       <  2:4323192/101‑1 (MQ=25)
                                                                       gcaagCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATaaa                       >  1:3725866/1‑101 (MQ=21)
                                                                        gaagCGAAAAGGGG.GCTCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAat                      <  2:9153584/100‑1 (MQ=21)
                                                                        caagCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAat                      <  2:500667/101‑1 (MQ=25)
                                                                        caagCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAat                      >  1:5637781/1‑101 (MQ=21)
                                                                          agCGAAAA.GGGTGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAatat                    <  2:1584776/101‑1 (MQ=25)
                                                                           gCGAAAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATc                   <  2:7338656/101‑1 (MQ=21)
                                                                            gaaaan.GGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATccc                 <  1:3064184/99‑1 (MQ=25)
                                                                            gaaaag.GGG.GCGATGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATccc                 >  1:8073372/3‑101 (MQ=21)
                                                                             aaaag.GGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCg                <  1:2774179/100‑1 (MQ=21)
                                                                              aaag.GGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCgg               >  1:4366199/1‑101 (MQ=21)
                                                                              aaag.GGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCgg               >  1:4655274/1‑101 (MQ=21)
                                                                              aaag.GGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCgg               <  2:2356762/101‑1 (MQ=21)
                                                                                aagggg.gCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCggg              >  1:206776/4‑101 (MQ=21)
                                                                                aagggg.gCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCggg              >  2:6095396/4‑101 (MQ=21)
                                                                                    gg.gCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGa          <  1:7652326/101‑1 (MQ=255)
                                                                                    gg.gCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGa          <  2:4742718/101‑1 (MQ=255)
                                                                                     g.gCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAg         <  1:3490998/101‑1 (MQ=255)
                                                                                     g.gCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAg         <  2:8744029/101‑1 (MQ=255)
                                                                                        cgcTGCCCC.CTCTTGCCCTGTCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTc       <  2:8536943/101‑1 (MQ=21)
                                                                                        cgcTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTc       <  2:7167009/101‑1 (MQ=255)
                                                                                         gcTGCCCC.CTCTTGCCCTGCCGGGCAAGTCATCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCTCGGGGGGAGTcc      <  1:6584564/101‑1 (MQ=25)
                                                                                         gcTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTcc      >  2:6432209/1‑101 (MQ=255)
                                                                                          cTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTCCg     <  2:5794803/101‑1 (MQ=255)
                                                                                            tgcccCCTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTCCGc    >  2:7055573/3‑101 (MQ=255)
                                                                                            cgcccCCTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGGccat    <  2:5986200/99‑3 (MQ=21)
                                                                                              cccCCTCTTGCCCTGCCGGGCAAGTCTCCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTCCGCAg  <  2:548206/101‑1 (MQ=21)
                                                                                              cccCCTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTCCGCAg  <  1:3245782/101‑1 (MQ=255)
                                                                                                                                                                                                   
GCGCCCCCGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAAAA.GGG.GCGCTGCCCCTCTCTTG.CCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTCCGCAG  >  NC_014034/2486347‑2486538

Base quality scores:  ATCG  < 3 ≤  ATCG  < 23 ≤  ATCG  < 30 ≤  ATCG  < 35 ≤  ATCG  < 41 ≤  ATCG