Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,696,177 +C 95.0% intergenic (‑105/+120) RCAP_rcc01564 ← / ← mmsA isobutyryl‑CoA dehydrogenase/methylmalonate‑semialdehyde dehydrogenase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,696,1771.→C95.0% 58.9 0.0 ‑1.1456intergenic (‑105/+120)RCAP_rcc01564/mmsAisobutyryl‑CoA dehydrogenase/methylmalonate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  new base (217/216):  ref base (12/11):  total (229/227)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.30e-02

Only 100 of 491 total aligned reads displayed.

GCTCCTCCCAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAAGGGCAGCTTGCCCCTCATGGCCTTTC.GG.CAATTCCA.CCCGGGAAGGTTGGGCCAAAGGGTGAAACCACCCTGGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGC  >  NC_014034/1696075‑1696278
                                                                                                        |                                                                                                      
gcTCCTCCCAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.gg                                                                                                         >  1:3754836/1‑101 (MQ=255)
  tcctccCAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GG.cc                                                                                                      <  2:2735116/101‑2 (MQ=255)
     tccCAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GG.ccaat                                                                                                   >  1:5748622/1‑97 (MQ=255)
      cccAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GG.ccaatt                                                                                                  >  2:5671478/1‑96 (MQ=255)
      cccAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GG.ccaatt                                                                                                  <  2:3092476/101‑6 (MQ=255)
       ccAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTc                                                                                                  <  2:5150027/101‑1 (MQ=255)
         aaaGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTcac                                                                                                >  2:8330812/1‑99 (MQ=255)
          aaGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GG.CCAATTCA.cc                                                                                             >  1:1831856/1‑101 (MQ=255)
            gTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GG.CCAATTCA.cccc                                                                                           <  1:4470687/101‑2 (MQ=37)
            gTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GG.CCAATTCA.cccc                                                                                           >  1:4037620/1‑100 (MQ=255)
              gCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GG.CCAATTCA.CCCcgg                                                                                         >  2:8146330/1‑101 (MQ=37)
               ccTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCggg                                                                                         <  2:1735253/101‑1 (MQ=37)
                cTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGa                                                                                        <  2:2348854/101‑1 (MQ=255)
                 tGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGat                                                                                       >  2:7722963/1‑100 (MQ=37)
                  gAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGatg                                                                                      <  1:4350597/101‑3 (MQ=255)
                  gAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGatg                                                                                      >  2:6194318/1‑99 (MQ=255)
                   aaTTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGatgt                                                                                     >  2:4274420/1‑98 (MQ=255)
                   aaTTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGatgt                                                                                     >  2:4088540/1‑98 (MQ=255)
                      tGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTg                                                                                  <  1:2194374/101‑1 (MQ=255)
                      tGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTg                                                                                  >  2:724668/1‑101 (MQ=255)
                      tGAACGTCTGTTTAATGTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTg                                                                                  <  2:1987588/101‑1 (MQ=255)
                       gAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTgg                                                                                 <  2:4443020/101‑1 (MQ=255)
                        aaCGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTggg                                                                                <  1:7498158/101‑1 (MQ=25)
                        aaCGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTggg                                                                                >  1:1979111/1‑101 (MQ=255)
                         aCGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGc                                                                               >  2:6916556/1‑101 (MQ=255)
                          cGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGcc                                                                              >  1:846411/1‑101 (MQ=25)
                          cGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGcc                                                                              <  1:5125584/101‑1 (MQ=255)
                           gTCTGTTTAATTTCGACTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCa                                                                             <  2:3789565/101‑1 (MQ=21)
                           gTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCa                                                                             <  1:5979220/101‑1 (MQ=25)
                             cTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCaag                                                                           <  1:4232545/101‑2 (MQ=255)
                             cTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCaag                                                                           >  2:54869/1‑100 (MQ=21)
                              tGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCaagg                                                                          <  2:7609857/101‑3 (MQ=19)
                              tGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCaagg                                                                          <  1:4869644/101‑3 (MQ=19)
                              tGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGGTTGGGCCaagg                                                                          >  1:5535577/1‑99 (MQ=255)
                               gTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCaaggt                                                                         <  2:811296/101‑4 (MQ=255)
                                tttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCaaggtg                                                                        >  1:937686/1‑97 (MQ=19)
                                tttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCaaggtg                                                                        <  1:3171981/101‑5 (MQ=255)
                                tttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCaaggtg                                                                        <  1:5671685/101‑5 (MQ=255)
                                 ttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCaaggtga                                                                       <  2:4869992/101‑6 (MQ=16)
                                  tAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCaaggtgaa                                                                      <  1:3432965/101‑7 (MQ=255)
                                  tAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCaaggtgaa                                                                      <  2:1178228/101‑7 (MQ=255)
                                   aaTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCaaggtgaaa                                                                     <  2:5052087/101‑8 (MQ=255)
                                    aTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAac                                                                  >  2:4623506/1‑101 (MQ=255)
                                     tttCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAaca                                                                 <  1:1253615/101‑2 (MQ=255)
                                      ttCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacac                                                                >  2:3585672/1‑99 (MQ=12)
                                         gAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacacctt                                                             <  1:6711735/101‑6 (MQ=255)
                                         gAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacacctt                                                             >  2:6108718/1‑96 (MQ=12)
                                         gAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacacctt                                                             >  1:9120559/1‑96 (MQ=12)
                                          aaTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacaccttg                                                            >  2:705191/1‑95 (MQ=255)
                                             ccTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacaccttgctt                                                         <  2:6684122/101‑10 (MQ=9)
                                              cTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacaccttgcttt                                                        >  1:74018/1‑91 (MQ=255)
                                              cTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacaccttgcttt                                                        <  2:801511/101‑11 (MQ=255)
                                                  cACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCt                                                  <  1:1117682/101‑1 (MQ=255)
                                                  cACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCt                                                  <  1:7956894/101‑1 (MQ=255)
                                                   aCGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCtt                                                 <  1:2400185/101‑1 (MQ=9)
                                                   aCGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCtt                                                 <  1:7628558/101‑1 (MQ=9)
                                                    cGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTg                                                >  2:1406863/1‑101 (MQ=255)
                                                     gAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTgg                                               <  1:7024941/101‑1 (MQ=255)
                                                     gAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTgg                                               <  1:4395264/101‑1 (MQ=255)
                                                      aaTTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGgc                                              <  1:6888984/101‑1 (MQ=9)
                                                        ttttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGgcgc                                            >  1:9000840/1‑101 (MQ=9)
                                                         tttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCa                                           >  2:8637723/1‑101 (MQ=255)
                                                          ttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGACAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaa                                          >  1:3331828/1‑101 (MQ=255)
                                                           tGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaaa                                         <  2:92221/101‑1 (MQ=255)
                                                           tGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaaa                                         >  2:1655709/1‑101 (MQ=9)
                                                             cAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATa                                       <  2:4537870/101‑1 (MQ=255)
                                                              aaGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATAc                                      <  1:1327004/101‑1 (MQ=255)
                                                               aggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACg                                     >  1:3911271/1‑101 (MQ=255)
                                                               aggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACg                                     >  1:2963891/1‑101 (MQ=255)
                                                               aggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACg                                     >  1:8133973/1‑101 (MQ=255)
                                                                ggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGc                                    <  1:5364500/101‑1 (MQ=255)
                                                                 gaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAACTACGcc                                   >  2:250904/1‑101 (MQ=2)
                                                                    gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                <  2:8804757/101‑1 (MQ=255)
                                                                         agaggggcgctGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGc                             <  1:4623641/91‑1 (MQ=255)
                                                                           aggggcgctGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGccc                           >  2:4742247/9‑101 (MQ=12)
                                                                           aggggcgctGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGccc                           <  1:1478243/93‑1 (MQ=255)
                                                                            ggggcgctGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCg                          <  1:8699939/94‑1 (MQ=16)
                                                                                cgctGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAgg                      <  1:2787232/98‑1 (MQ=16)
                                                                                 gctGCCCCGCCTGGCCTTTCGGG.CAATTCCC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAggg                     <  1:1640930/99‑1 (MQ=19)
                                                                                   tGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCa                   <  2:642616/101‑1 (MQ=16)
                                                                                   tGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCa                   <  2:5926367/101‑1 (MQ=16)
                                                                                   tGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCa                   <  1:7824082/101‑1 (MQ=16)
                                                                                    gCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCaa                  <  1:616045/101‑1 (MQ=19)
                                                                                    gCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCaa                  >  2:221261/1‑101 (MQ=255)
                                                                                    gCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCaa                  >  2:3496701/1‑101 (MQ=19)
                                                                                      cccTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGc                <  2:3256364/101‑1 (MQ=16)
                                                                                      cccTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGc                <  1:4742396/101‑1 (MQ=16)
                                                                                      cccTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGc                <  1:5827454/101‑1 (MQ=16)
                                                                                        cTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCcg              >  2:2682931/1‑101 (MQ=255)
                                                                                         tcctGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCcgc             >  1:2575379/4‑101 (MQ=255)
                                                                                         tcctGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGCCCCGAAGGGCAAGAcgc             >  1:360605/4‑101 (MQ=16)
                                                                                          cctGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCcgcg            >  1:992224/3‑101 (MQ=255)
                                                                                            tGGCCTTTC.GGCCAATTCCC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGcc          <  1:6194557/101‑1 (MQ=25)
                                                                                            tGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGcc          >  2:5294485/1‑101 (MQ=255)
                                                                                             ggCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCg         >  1:8336171/1‑101 (MQ=255)
                                                                                             ggCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCg         <  1:6583477/101‑1 (MQ=19)
                                                                                              gCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGa        >  1:4658545/1‑101 (MQ=255)
                                                                                                 tttC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAgg     >  2:536109/1‑101 (MQ=255)
                                                                                                  ttC.GGCCAATTCCC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGgc    <  2:1831803/101‑1 (MQ=25)
                                                                                                    c.GG.CCAATCCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGgcgc  <  2:9014228/101‑1 (MQ=19)
                                                                                                        |                                                                                                      
GCTCCTCCCAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAAGGGCAGCTTGCCCCTCATGGCCTTTC.GG.CAATTCCA.CCCGGGAAGGTTGGGCCAAAGGGTGAAACCACCCTGGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGC  >  NC_014034/1696075‑1696278

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