Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,486,442 +C 3.3% 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.3% 28.5 ‑0.3 ‑0.2368intergenic (‑70/+124)RCAP_rcc02312/soxGhypothetical protein/sarcosine oxidase subunit gamma
Reads supporting (aligned to +/- strand):  new base (5/7):  ref base (179/177):  total (184/184)
Fisher's exact test strand distribution p-value = 7.71e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.60e-01

Only 100 of 381 total aligned reads displayed.

GCGCCCCCGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AA.GGGGCGCTGCCCCTCTCTTG.CCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTCCGCAG  >  NC_014034/2486347‑2486538
                                                                                                                                                                                                   
 cgcCCCCGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTtg                                                                                             >  1:6468123/1‑101 (MQ=255)
 cgcCCCCGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTtg                                                                                             <  2:5928495/101‑1 (MQ=255)
 cgcCCCCGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTtg                                                                                             <  1:2187718/101‑1 (MQ=255)
   cccccGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTtg.cc                                                                                          <  2:3301547/101‑1 (MQ=255)
    ccccGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTtg.ccc                                                                                         >  2:3010803/1‑100 (MQ=255)
     cccGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTtg.ccct                                                                                        <  1:6079832/101‑3 (MQ=255)
      ccGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTtg.ccctg                                                                                       <  2:1948833/101‑4 (MQ=255)
      ccGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTtg.ccctg                                                                                       <  1:2608391/101‑4 (MQ=255)
       cGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTtg.ccctgc                                                                                      >  1:811702/1‑97 (MQ=255)
         gATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCg                                                                                     >  2:3605849/1‑101 (MQ=255)
         gATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCg                                                                                     >  2:6456211/1‑101 (MQ=255)
          aTTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCgg                                                                                    >  2:4466909/1‑101 (MQ=255)
           tttCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCGC.CTCTTGCCCTGCCggg                                                                                   >  1:1492864/1‑101 (MQ=255)
           tttCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTtg.cccttccggg                                                                                  <  1:4683909/101‑9 (MQ=255)
           tttCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCggg                                                                                   <  2:1450917/101‑1 (MQ=255)
              cGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCaa                                                                                >  2:5429839/1‑101 (MQ=38)
              cGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCaa                                                                                <  2:5843785/101‑1 (MQ=38)
                aaTAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGt                                                                              <  1:6541419/101‑1 (MQ=38)
                aaTAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGt                                                                              >  2:4252737/1‑101 (MQ=38)
                aaTAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGt                                                                              >  2:2831375/1‑101 (MQ=38)
                aaTAGTCCATGCTGGATCTCCTTTGCCGTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGt                                                                              <  2:3715975/101‑1 (MQ=38)
                  tAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCa                                                                            >  1:2470304/1‑101 (MQ=38)
                   aGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCACGTCAc                                                                           >  2:1080056/1‑101 (MQ=38)
                    gTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCAcc                                                                          >  1:5377818/1‑101 (MQ=38)
                       cATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCAccccc                                                                       >  1:519376/1‑101 (MQ=37)
                        aTGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCTCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCg                                                                      <  2:4858444/101‑1 (MQ=37)
                        aTGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCg                                                                      >  1:4217952/1‑101 (MQ=37)
                          gCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCggg                                                                    >  1:6414388/1‑101 (MQ=37)
                           cTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGa                                                                   >  1:2037799/1‑101 (MQ=37)
                            tGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CCCTTGCCCTGCCGGGCAAGTCACCCCCgggcc                                                                  >  1:1527870/1‑99 (MQ=37)
                              gATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGatat                                                                >  1:4349351/1‑101 (MQ=37)
                               aTCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCCTGCCCTGCCGGGCAAGTCACCCCCGGGATAtt                                                               >  2:3897100/1‑101 (MQ=37)
                                   ccTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATg                                                           <  2:3093439/101‑1 (MQ=37)
                                    ttttGCCTGCGAACTGCGGCAAGGCCGACAGCAAGGAAAGCGAA.AAGGGGGCGCTTCCCC.CTCTTGATCTGCCGGGCAAGTCACCCCCGGGATATTCATgg                                                          <  1:3588665/100‑1 (MQ=255)
                                    cTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATgg                                                          <  1:5644841/101‑1 (MQ=37)
                                     tttGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGc                                                         >  1:5742442/1‑101 (MQ=37)
                                      ttGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCa                                                        >  1:5006962/1‑101 (MQ=37)
                                      ttGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCa                                                        <  1:58040/101‑1 (MQ=37)
                                      ttGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTAATGGCa                                                        >  1:3296793/1‑101 (MQ=25)
                                      ttGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGCATATTCATGGCa                                                        >  1:7081129/1‑101 (MQ=37)
                                        gCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAg                                                      <  1:2580503/101‑1 (MQ=37)
                                        gCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAg                                                      >  2:6911239/1‑101 (MQ=37)
                                          cTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGAt                                                    <  2:1632448/101‑1 (MQ=37)
                                           ttCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTTCCGGGCAAGTCAACCCAGGGATATTCATGGCAAGATg                                                   >  2:5691591/1‑101 (MQ=37)
                                            tCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAACAG.GGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGa                                                  <  1:6273220/101‑1 (MQ=25)
                                              gATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAag                                                <  2:7057428/101‑1 (MQ=37)
                                               aTCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAagg                                               >  1:3096064/1‑101 (MQ=37)
                                                tCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAagga                                              <  2:3501473/101‑1 (MQ=37)
                                                   gCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGc                                           >  1:5403900/1‑101 (MQ=25)
                                                   gCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGc                                           >  1:6898320/1‑101 (MQ=25)
                                                    cGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCg                                          <  1:3386880/101‑1 (MQ=37)
                                                    cGGCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCg                                          >  1:5988286/1‑101 (MQ=25)
                                                     ggCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCgt                                         <  2:2256865/101‑1 (MQ=37)
                                                     ggCAAGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCgt                                         <  1:3173363/101‑1 (MQ=37)
                                                        aaGCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCgtggt                                      >  1:2973361/1‑99 (MQ=25)
                                                          gCCCGGCCGCCAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAAGAGGCGTGTTTc                                    <  1:5079449/101‑1 (MQ=21)
                                                          gCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTc                                    <  2:6232961/101‑1 (MQ=25)
                                                          gCCCGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTc                                    <  2:7103357/101‑1 (MQ=25)
                                                           cccGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCa                                   <  1:4130014/101‑1 (MQ=37)
                                                            ccGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCAt                                  <  1:6633836/101‑1 (MQ=25)
                                                            ccGGCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCAt                                  <  2:259742/101‑1 (MQ=25)
                                                              ggCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCt                                <  2:2470210/101‑1 (MQ=37)
                                                               gCCGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCtt                               <  2:5004680/101‑1 (MQ=25)
                                                                ccGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTg                              <  1:4462696/101‑1 (MQ=25)
                                                                ccGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTg                              >  2:1767770/1‑101 (MQ=25)
                                                                 cGCAAGGCAAGCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGc                             >  1:2138855/1‑101 (MQ=21)
                                                                  gcaaggcaagCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTgtc                            <  1:1759242/101‑3 (MQ=25)
                                                                  gcaaggcaagCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGcc                            <  2:1558588/101‑1 (MQ=25)
                                                                  gcaaggcaagCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGcc                            <  1:6799818/101‑1 (MQ=25)
                                                                    aaggcaagCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCAt                          <  2:519364/101‑1 (MQ=25)
                                                                    aaggcaagCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCAt                          >  1:3715191/1‑101 (MQ=21)
                                                                      ggcaagCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATaa                        <  1:4854527/101‑1 (MQ=25)
                                                                      ggcaagCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATaa                        <  1:1912339/101‑1 (MQ=25)
                                                                      ggcaagCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATaa                        <  2:3957730/101‑1 (MQ=25)
                                                                       gcaagCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATaaa                       <  2:4412920/101‑1 (MQ=25)
                                                                        caagCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAat                      >  2:5603124/1‑101 (MQ=21)
                                                                        caagCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAat                      <  1:4499800/101‑1 (MQ=21)
                                                                         aagCGAA.AAGGGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGGGCTTCATCTTGCCATAAata                     >  1:7025617/1‑101 (MQ=25)
                                                                           cgaaa.ag.GGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATcc                  <  1:3946100/98‑1 (MQ=21)
                                                                           cgaaa.ag.GGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATcc                  <  2:4700452/98‑1 (MQ=18)
                                                                            gaaa.ag.GGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATccc                 <  1:1192927/99‑1 (MQ=21)
                                                                             aaa.ag.GGGGCTCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCg                <  2:911783/100‑1 (MQ=18)
                                                                             aaa.ag.GGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCg                <  1:1320310/100‑1 (MQ=21)
                                                                             aaa.ag.GGGGCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCg                >  2:1848156/2‑101 (MQ=21)
                                                                                 aagggggCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCggg              >  1:6575902/4‑101 (MQ=21)
                                                                                     gggCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGa          <  1:1844547/101‑1 (MQ=255)
                                                                                      ggCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAg         <  2:160238/101‑1 (MQ=255)
                                                                                      ggCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAg         <  2:6771845/101‑1 (MQ=255)
                                                                                      ggCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAg         >  1:2415674/1‑101 (MQ=255)
                                                                                       gcgcTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGt        <  2:2223329/101‑1 (MQ=255)
                                                                                         gcTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGGGTcc      >  1:1622877/1‑101 (MQ=21)
                                                                                         gcTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTcc      >  1:1048342/1‑101 (MQ=255)
                                                                                         gcTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTcc      >  2:812807/1‑101 (MQ=255)
                                                                                         gcTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTcc      >  1:3511388/1‑101 (MQ=255)
                                                                                         gcTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTcc      <  2:4325267/101‑1 (MQ=255)
                                                                                            tgcccCCTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTCCGc    <  1:6182547/99‑1 (MQ=255)
                                                                                            tgcccCCTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTCCGc    <  2:5803125/99‑1 (MQ=255)
                                                                                            tgcccCCTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTCCGc    >  2:2779128/3‑101 (MQ=255)
                                                                                             ccccCCTCTTGCCCTGCCGGGCAAGTCCCCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTCCGCa   <  1:2765877/101‑1 (MQ=21)
                                                                                              cccCCTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTCCGCAg  >  2:3351565/1‑101 (MQ=255)
                                                                                                                                                                                                   
GCGCCCCCGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGAA.AA.GGGGCGCTGCCCCTCTCTTG.CCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTCCGCAG  >  NC_014034/2486347‑2486538

Base quality scores:  ATCG  < 3 ≤  ATCG  < 24 ≤  ATCG  < 31 ≤  ATCG  < 35 ≤  ATCG  < 41 ≤  ATCG