Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,696,179 A→C 7.4% intergenic (‑107/+118) 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,1790A→C 7.4% 43.3 ‑2.3 ‑0.0268intergenic (‑107/+118)RCAP_rcc01564/mmsAisobutyryl‑CoA dehydrogenase/methylmalonate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  new base (13/7):  ref base (111/137):  total (124/144)
Fisher's exact test strand distribution p-value = 1.03e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.92e-01

Only 100 of 277 total aligned reads displayed.

CCAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAAGGGCAGCTTGCCCCTCATGGCCTTT.CG.G.CAATTCCACCCGGGAAGGTTGGGCCAAAGGGTGAAACCACCCTGGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCC  >  NC_014034/1696082‑1696283
                                                                                                    |                                                                                                        
ccAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTc                                                                                                      >  1:1727383/1‑101 (MQ=255)
   aaGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.G.CCAATTCAcc                                                                                                  >  2:5272432/1‑101 (MQ=255)
     gTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.G.CCAATTCAcccc                                                                                                <  2:1056625/101‑2 (MQ=37)
     gTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.G.CCAATTCAcccc                                                                                                >  1:5060482/1‑100 (MQ=255)
      tGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.G.CCAATTCAccccg                                                                                               <  1:1151131/101‑3 (MQ=255)
          tGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGat                                                                                            >  1:2555858/1‑100 (MQ=37)
           gAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGatg                                                                                           >  2:3510539/1‑99 (MQ=255)
           gAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGatg                                                                                           >  1:1668672/1‑99 (MQ=255)
           gAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGatg                                                                                           >  2:4697478/1‑99 (MQ=255)
            aaTTTATCGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAG.AGGGTCGC.TCCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGatgt                                                                                          <  2:4616024/101‑4 (MQ=255)
            aaTTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGatgt                                                                                          >  1:649354/1‑98 (MQ=255)
             agtcaaCGTCTGTTTAATTTCGAATCCTTGCGCGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGatgtt                                                                                         <  1:3235478/97‑5 (MQ=255)
             aTTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGatgtt                                                                                         >  1:425845/1‑97 (MQ=25)
             aTTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGatgtt                                                                                         >  1:3500385/1‑97 (MQ=25)
              ttGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGatgttt                                                                                        <  2:3651499/101‑6 (MQ=255)
              ttGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGatgttt                                                                                        >  1:332376/1‑96 (MQ=255)
               tGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTg                                                                                       >  1:1759181/1‑101 (MQ=255)
               tGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTg                                                                                       >  2:1485453/1‑101 (MQ=255)
                gAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCgggctgttagg                                                                                      >  2:4877616/1‑93 (MQ=25)
                 aaCGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTggg                                                                                     >  1:5725083/1‑101 (MQ=255)
                  aCGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGc                                                                                    <  2:2068822/101‑1 (MQ=255)
                   cGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGcc                                                                                   <  2:1414954/101‑1 (MQ=255)
                    gTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCa                                                                                  >  2:9278/1‑101 (MQ=255)
                    gTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCa                                                                                  >  2:1848175/1‑101 (MQ=255)
                     tCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCaa                                                                                 >  2:1934612/1‑101 (MQ=255)
                     tCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCaa                                                                                 >  2:4469657/1‑101 (MQ=255)
                      cTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.G.ATGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCaag                                                                                >  2:410019/1‑100 (MQ=19)
                       tGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCaagg                                                                               <  1:251494/101‑3 (MQ=19)
                         tttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCaaggtg                                                                             <  1:3778367/101‑5 (MQ=255)
                          ttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCaaggtga                                                                            <  1:1214201/101‑6 (MQ=16)
                          ttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCaaggtga                                                                            <  2:2849156/101‑6 (MQ=16)
                           tAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCaaggtgaa                                                                           >  1:2666089/1‑95 (MQ=255)
                            aaTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCaaggtgaaa                                                                          >  2:3585925/1‑94 (MQ=16)
                            aaTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCaaggtgaaa                                                                          >  2:1336369/1‑94 (MQ=16)
                            aaTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCaaggtgaaa                                                                          <  2:5387312/101‑8 (MQ=255)
                             aTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAac                                                                       <  1:5656327/101‑1 (MQ=16)
                                tCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacc                                                                    >  2:5760377/1‑98 (MQ=255)
                                tCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacc                                                                    <  1:974318/101‑4 (MQ=12)
                                  gAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacctt                                                                  <  2:2394705/101‑6 (MQ=255)
                                   aaTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGCC.CAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttg                                                                 >  2:382282/1‑95 (MQ=255)
                                   aaTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttg                                                                 <  2:425801/101‑7 (MQ=12)
                                   aaTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttg                                                                 <  1:4048124/101‑7 (MQ=12)
                                   aaTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttg                                                                 >  2:711231/1‑95 (MQ=255)
                                   aaTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTCGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacagcttg                                                                 >  2:792852/1‑95 (MQ=255)
                                    aTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgc                                                                <  2:545941/101‑8 (MQ=255)
                                     tCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgct                                                               <  1:277363/101‑9 (MQ=255)
                                      ccTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgctt                                                              >  2:2540259/1‑92 (MQ=255)
                                          gCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTc                                                        <  1:4745554/101‑1 (MQ=255)
                                           cACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCt                                                       <  1:289385/101‑1 (MQ=255)
                                            tctaaTTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCtt                                                      <  1:4534363/98‑1 (MQ=2)
                                            cgatcttttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCtt                                                      <  1:1522735/96‑1 (MQ=2)
                                            aCGAGTTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCtt                                                      <  2:2641926/101‑1 (MQ=2)
                                            aCGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCtt                                                      >  2:3291777/1‑101 (MQ=255)
                                             cGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTg                                                     >  1:101160/1‑101 (MQ=255)
                                             cGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTg                                                     <  1:5175985/101‑1 (MQ=255)
                                                 ttttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGgcgc                                                 >  1:5073464/1‑101 (MQ=9)
                                                  tttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCa                                                <  1:1568370/101‑1 (MQ=9)
                                                   ttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaa                                               >  2:4568820/1‑101 (MQ=255)
                                                   ttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaa                                               >  1:4806610/1‑101 (MQ=255)
                                                    tGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaaa                                              >  1:4960620/1‑101 (MQ=9)
                                                     gCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCCG.C.CAATTTACCCCGGGAGGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAAt                                             >  2:1155929/1‑101 (MQ=255)
                                                      cAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATa                                            <  2:5337219/101‑1 (MQ=255)
                                                      cAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATa                                            <  1:2215698/101‑1 (MQ=255)
                                                      cAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATa                                            <  2:347716/101‑1 (MQ=255)
                                                       aaGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATAc                                           >  2:4138663/1‑101 (MQ=9)
                                                        aggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACg                                          <  1:4121887/101‑1 (MQ=9)
                                                         ggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGc                                         >  1:2305146/1‑101 (MQ=255)
                                                          gaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGcc                                        >  1:390964/1‑101 (MQ=9)
                                                          gaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGcc                                        <  1:4582439/101‑1 (MQ=255)
                                                             gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                     <  1:3586142/101‑1 (MQ=255)
                                                             gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                     >  2:314333/1‑101 (MQ=9)
                                                                  agaggggcgctGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGc                                  <  1:3006840/91‑1 (MQ=255)
                                                                   gaggggcgctGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGcc                                 >  1:1590564/10‑101 (MQ=16)
                                                                   gaggggcgctGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGcc                                 <  2:1624173/92‑1 (MQ=9)
                                                                   gaggggcgctGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGcc                                 >  1:815329/10‑101 (MQ=16)
                                                                    aggggcgctGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGccc                                >  1:5508406/9‑101 (MQ=16)
                                                                       ggcgctGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGaa                             <  1:2459546/96‑1 (MQ=255)
                                                                         cgctGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAgg                           >  2:2899745/4‑101 (MQ=16)
                                                                          gctGCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAggg                          >  2:433399/3‑101 (MQ=16)
                                                                             gCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCGTCTGGCCCGAAGGGCaa                       >  2:5195249/1‑101 (MQ=255)
                                                                             gCCCCTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCaa                       <  1:5786815/101‑1 (MQ=255)
                                                                              ccccTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAg                      <  1:1485546/101‑1 (MQ=255)
                                                                              ccccTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAg                      >  1:1362150/1‑101 (MQ=16)
                                                                               cccTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGc                     <  1:869452/101‑1 (MQ=16)
                                                                               cccTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGc                     <  1:5338403/101‑1 (MQ=16)
                                                                               cccTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGc                     <  2:2903700/101‑1 (MQ=16)
                                                                                ccTCCTGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGcc                    <  2:2459900/101‑1 (MQ=255)
                                                                                     tGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAgatcggaa               >  2:5706634/1‑94 (MQ=16)
                                                                                     tGGCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGcc               >  2:2169963/1‑101 (MQ=255)
                                                                                       gCCTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGa             <  2:4489525/101‑1 (MQ=19)
                                                                                        ccTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGaa            <  2:4792125/101‑1 (MQ=19)
                                                                                         cTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAg           >  1:2404714/1‑101 (MQ=19)
                                                                                         cTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAg           <  1:3521573/101‑1 (MQ=19)
                                                                                         cTTT.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTCTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAg           <  2:5010192/101‑1 (MQ=255)
                                                                                          ttt.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAgg          <  1:213777/101‑1 (MQ=19)
                                                                                          ttt.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAgg          >  1:312558/1‑101 (MQ=255)
                                                                                          ttt.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAgg          >  1:1449903/1‑101 (MQ=255)
                                                                                          ttt.CG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAgg          <  1:5422572/101‑1 (MQ=19)
                                                                                                ggccAATTCCCCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGAAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGcgcg      <  2:2555716/98‑1 (MQ=255)
                                                                                                   cAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGcc   <  2:1590470/101‑1 (MQ=25)
                                                                                                    |                                                                                                        
CCAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAAGGGCAGCTTGCCCCTCATGGCCTTT.CG.G.CAATTCCACCCGGGAAGGTTGGGCCAAAGGGTGAAACCACCCTGGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCC  >  NC_014034/1696082‑1696283

Base quality scores:  ATCG  < 3 ≤  ATCG  < 19 ≤  ATCG  < 28 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG