Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,486,442 +C 1.6% 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 1.6% 5.4 ‑1.1 ‑1.0244intergenic (‑70/+124)RCAP_rcc02312/soxGhypothetical protein/sarcosine oxidase subunit gamma
Reads supporting (aligned to +/- strand):  new base (3/1):  ref base (110/130):  total (113/131)
Fisher's exact test strand distribution p-value = 3.39e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.92e-02

Only 100 of 250 total aligned reads displayed.

GCGCCCCCGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAA.GGG.GCGCTGCCCCTCTCTTG.CCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTCCGCAG  >  NC_014034/2486347‑2486538
                                                                                                                                                                                                    
ccgcCCCCGGATTTCGAATAGTCCATTCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCtt                                                                                              <  1:2856697/100‑1 (MQ=255)
 cgcCCCCGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTtg                                                                                             >  2:2504734/1‑101 (MQ=255)
   cccccGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCAGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTtg.cc                                                                                          >  1:4245961/1‑101 (MQ=255)
     cccGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTtg.ccct                                                                                        >  1:1302685/1‑99 (MQ=255)
     cccGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTtg.ccct                                                                                        <  2:4385163/101‑3 (MQ=255)
      ccGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTtg.ccctg                                                                                       >  2:24641/1‑98 (MQ=255)
       cGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTtg.ccctgc                                                                                      >  2:4541832/1‑97 (MQ=255)
        ggATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTtg.ccctgcc                                                                                     <  1:3126640/101‑6 (MQ=255)
        ggATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTtg.ccctgcc                                                                                     >  2:3868257/1‑96 (MQ=255)
         gATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCg                                                                                     >  2:4091050/1‑101 (MQ=255)
         gATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCg                                                                                     >  1:1199452/1‑101 (MQ=255)
           tttCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCggg                                                                                   >  1:919132/1‑101 (MQ=255)
             tCGAATAGTCCATGCTGGATCTCCTTTGGCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCa                                                                                 >  2:359321/1‑101 (MQ=38)
             tCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCa                                                                                 <  2:4337394/101‑1 (MQ=38)
              cGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAATCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCaa                                                                                <  2:4634372/101‑1 (MQ=38)
              cGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCaa                                                                                <  2:2227753/101‑1 (MQ=38)
               gAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAg                                                                               >  1:4023474/1‑101 (MQ=38)
               gAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAg                                                                               >  2:1024221/1‑101 (MQ=38)
                 aTAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTc                                                                             <  2:1728938/101‑1 (MQ=38)
                 aTAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTc                                                                             >  1:1214823/1‑101 (MQ=38)
                 aTAGTCCATGCTGGATCTCCTTTGCCTTCGATCTCCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTc                                                                             <  2:4050520/101‑1 (MQ=38)
                  tAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCa                                                                            >  2:2318431/1‑101 (MQ=38)
                   aGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GGGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCAc                                                                           >  2:1121563/1‑101 (MQ=38)
                   aGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCAc                                                                           <  2:2084081/101‑1 (MQ=38)
                     tCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCAccc                                                                         >  2:1093151/1‑101 (MQ=37)
                     tCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCAccc                                                                         >  2:3780465/1‑101 (MQ=37)
                        aTGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCg                                                                      <  2:1784138/101‑1 (MQ=37)
                          gCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCggg                                                                    >  2:1583364/1‑101 (MQ=37)
                            tGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGat                                                                  <  1:1869106/101‑1 (MQ=37)
                            tGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAA.GGGTGCGCTGCCCC.CTGTTGCCCTGCCGGGCAAGTCACCCCCGGGat                                                                  >  1:300800/1‑101 (MQ=255)
                             ggATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGata                                                                 >  1:4351233/1‑101 (MQ=37)
                             ggATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGata                                                                 >  1:1238550/1‑101 (MQ=37)
                             ggATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGata                                                                 <  1:336308/101‑1 (MQ=37)
                             ggATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGata                                                                 >  1:2627044/1‑101 (MQ=37)
                                tctcCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTc                                                              >  1:1691976/1‑101 (MQ=37)
                                tctcCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTc                                                              <  2:4753808/101‑1 (MQ=37)
                                  tcCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCAt                                                            <  1:2503639/101‑1 (MQ=37)
                                  tcCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCAt                                                            >  1:2760288/1‑101 (MQ=37)
                                  tcCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCCGGAAATTCAt                                                            >  1:2949022/1‑101 (MQ=25)
                                  tcCTTTGCCTACTATCTACGTCTAGCCCGACCGCAAGGCAAGCGACAAC.GGG.GCGCTGCCCC.CTCTTGCCCTGCCGGTCAATTCTCCCCCGGGATATTCAt                                                            <  2:780071/101‑1 (MQ=12)
                                   ccTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATg                                                           >  1:3750078/1‑101 (MQ=37)
                                     tttGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGc                                                         <  1:3579985/101‑1 (MQ=37)
                                     tttGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGATGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGc                                                         >  1:4695732/1‑101 (MQ=37)
                                      ttGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCa                                                        >  2:3168230/1‑101 (MQ=37)
                                       tGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCaa                                                       <  2:337088/101‑1 (MQ=37)
                                        gCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAg                                                      <  1:4628645/101‑1 (MQ=37)
                                        gCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAg                                                      >  2:3125813/1‑101 (MQ=37)
                                        gCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAg                                                      <  2:2626892/101‑1 (MQ=37)
                                         ccTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGa                                                     <  1:1001390/101‑1 (MQ=37)
                                         ccTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGa                                                     <  1:3083533/101‑1 (MQ=37)
                                          cTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGAt                                                    <  1:1625586/101‑1 (MQ=37)
                                          cTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGAt                                                    >  2:1066680/1‑101 (MQ=37)
                                           ttCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATg                                                   <  1:945695/101‑1 (MQ=37)
                                             cGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGaa                                                 <  2:2139092/101‑1 (MQ=37)
                                             cGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGaa                                                 >  1:4192402/1‑101 (MQ=37)
                                              gATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAag                                                <  1:488244/101‑1 (MQ=37)
                                              gATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAag                                                >  2:3935534/1‑101 (MQ=37)
                                               aTCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAagg                                               <  1:2086193/101‑1 (MQ=37)
                                                  tGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAaggagg                                            <  1:1558740/101‑1 (MQ=37)
                                                    cGGCAAGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCg                                          <  1:4551444/101‑1 (MQ=37)
                                                         aGCCCGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGttt                                     <  1:4134911/101‑1 (MQ=37)
                                                           cccGGCCGCAAGGCAAGCGA.AAAGGGG.GCCCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCa                                   <  2:1334566/101‑1 (MQ=37)
                                                            ccGGCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCAt                                  >  2:1682973/1‑101 (MQ=25)
                                                              ggCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCt                                >  1:329219/1‑101 (MQ=25)
                                                               gCCGCAAGGCCACCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCCCCCCCCGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCtt                               <  1:2349077/101‑1 (MQ=21)
                                                               gCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCtt                               <  2:1826899/101‑1 (MQ=25)
                                                               gCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCtt                               <  1:4091292/101‑1 (MQ=25)
                                                               gCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCtt                               <  2:2545801/101‑1 (MQ=25)
                                                               gCCGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCtt                               <  2:1219468/101‑1 (MQ=25)
                                                                ccGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTg                              <  1:2178431/101‑1 (MQ=25)
                                                                ccGCAAGGCAAGCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATAGCAAGATGAAGGAGGCGTGTTTCATCTTg                              >  2:690125/1‑101 (MQ=25)
                                                                  gcaaggcaagCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGcc                            <  1:468432/101‑1 (MQ=25)
                                                                  gcaaggcaagCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGcc                            >  1:2560109/1‑101 (MQ=21)
                                                                   caaggcaagCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCa                           >  1:1956822/1‑101 (MQ=21)
                                                                     aggcaagCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATa                         <  1:1460852/101‑1 (MQ=21)
                                                                     aggcaagCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATa                         <  1:4495301/101‑1 (MQ=25)
                                                                     aggcaagCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATa                         <  2:3024018/101‑1 (MQ=25)
                                                                     aggcaagCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATa                         <  1:3677498/101‑1 (MQ=21)
                                                                      ggcaagCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATaa                        <  2:1985094/101‑1 (MQ=25)
                                                                       gcaagCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATaaa                       <  2:1199362/101‑1 (MQ=25)
                                                                          agCGA.ACAGGGG.GCGGTGCCCC.CCCTGGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAatat                    <  2:2978101/101‑1 (MQ=18)
                                                                          agCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGGAGGAGGCGTGTTTCATCTTGCCATAAatat                    <  2:3438804/101‑1 (MQ=21)
                                                                          agCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAatat                    <  1:1038650/101‑1 (MQ=21)
                                                                          agCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAatat                    >  1:1521012/1‑101 (MQ=21)
                                                                           gCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATc                   <  1:223038/101‑1 (MQ=21)
                                                                           gCGA.AAAGGGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATc                   <  2:4449041/101‑1 (MQ=21)
                                                                           cgaa.aag.GGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATcc                  <  1:13913/98‑1 (MQ=21)
                                                                           cgaa.aag.GGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATcc                  >  2:3955092/4‑101 (MQ=21)
                                                                            gaa.aag.GGG.GCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATACccc                 <  2:4023232/99‑1 (MQ=21)
                                                                                   gggg.gCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCggggg            <  2:1433966/100‑1 (MQ=21)
                                                                                    ggg.gCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCgggggg           <  1:1765198/101‑1 (MQ=255)
                                                                                    ggg.gCGCTACCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCgggggg           <  2:2843299/101‑1 (MQ=21)
                                                                                      g.gCGCTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAg         <  2:425341/101‑1 (MQ=255)
                                                                                        gcgcTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGt        <  1:3868499/101‑1 (MQ=255)
                                                                                         cgcTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTc       >  1:3732785/1‑101 (MQ=255)
                                                                                          gcTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTcc      >  2:2910247/1‑101 (MQ=255)
                                                                                          gcTGCCCC.CTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGGcc      >  2:326782/1‑101 (MQ=21)
                                                                                             tgcccCCTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTCCGc    >  2:3452214/3‑101 (MQ=255)
                                                                                               cccCCTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTCCGCAg  >  1:708816/1‑101 (MQ=255)
                                                                                               cccCCTCTTGCCCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTCCGCAg  <  2:3555344/101‑1 (MQ=255)
                                                                                                                                                                                                    
GCGCCCCCGGATTTCGAATAGTCCATGCTGGATCTCCTTTGCCTTCGATCTGCGGCAAGCCCGGCCGCAAGGCAAGCGA.AAA.GGG.GCGCTGCCCCTCTCTTG.CCTGCCGGGCAAGTCACCCCCGGGATATTCATGGCAAGATGAAGGAGGCGTGTTTCATCTTGCCATAAATATCCCGGGGGGAGTCCGCAG  >  NC_014034/2486347‑2486538

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