Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,696,185 A→C 100% intergenic (‑113/+112) 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,1850A→C100.0% 39.1 266intergenic (‑113/+112)RCAP_rcc01564/mmsAisobutyryl‑CoA dehydrogenase/methylmalonate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  new base (109/137):  ref base (13/7):  total (122/144)
Fisher's exact test strand distribution p-value = 1.01e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.30e-06

Only 100 of 273 total aligned reads displayed.

AAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAAGGGCAGCTTGCCCCTCATGGCCTTTCGG.CAATTCCACCCGGGAAGGTTGGGCCAAAGGGTGAAACCACCCTGGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCC  >  NC_014034/1696085‑1696286
                                                                                                     |                                                                                                     
aaGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGG.CCAATTCAcc                                                                                                     >  1:1169992/1‑101 (MQ=255)
aaGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGG.CCAATTCAcc                                                                                                     >  2:5272432/1‑101 (MQ=255)
 aGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGG.CCAATTCAccc                                                                                                    >  1:2403507/1‑101 (MQ=255)
  gTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGG.CCAATTCAcccc                                                                                                   >  1:5060482/1‑100 (MQ=255)
  gTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGG.CCAATTCAcccc                                                                                                   <  2:3635451/101‑2 (MQ=37)
  gTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGG.CCAATTCAcccc                                                                                                   <  1:935325/101‑2 (MQ=37)
    gCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGG.CCAATTCACCCcgg                                                                                                 >  2:4121638/1‑101 (MQ=37)
       tGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGat                                                                                               >  1:2555858/1‑100 (MQ=37)
        gAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGatg                                                                                              >  2:4697478/1‑99 (MQ=255)
        gAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGatg                                                                                              >  2:3510539/1‑99 (MQ=255)
          agtcaaCGTCTGTTTAATTTCGAATCCTTGCGCGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGatgtt                                                                                            <  1:3235478/97‑5 (MQ=255)
          aTTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGatgtt                                                                                            >  1:3500385/1‑97 (MQ=25)
           ttGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGatgttt                                                                                           >  1:332376/1‑96 (MQ=255)
            tGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTg                                                                                          >  1:1759181/1‑101 (MQ=255)
            tGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTg                                                                                          >  2:5338113/1‑101 (MQ=255)
             gAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTgg                                                                                         <  1:5628285/101‑1 (MQ=255)
              aaCGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTggg                                                                                        >  1:5725083/1‑101 (MQ=255)
              aaCGTCTGGTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGGCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGCTGTTTggg                                                                                        >  2:5763773/1‑101 (MQ=255)
                cGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGcc                                                                                      >  2:1045481/1‑101 (MQ=25)
                 gtgtGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCa                                                                                     <  1:76336/98‑1 (MQ=19)
                 gTCTTTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCa                                                                                     <  2:4231027/101‑1 (MQ=19)
                 gTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCa                                                                                     >  2:9278/1‑101 (MQ=255)
                  tCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaa                                                                                    <  2:4319786/101‑1 (MQ=255)
                   cTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaag                                                                                   <  1:4130743/101‑2 (MQ=255)
                    tGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaagg                                                                                  <  2:3344486/101‑3 (MQ=19)
                    tGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaagg                                                                                  <  1:1557816/101‑3 (MQ=19)
                     gTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaaggt                                                                                 >  1:4290886/1‑98 (MQ=255)
                     gTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAAGgg                                                                                 >  1:1727236/1‑101 (MQ=255)
                      tttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGCGG..TGCCCCCCCCGGCCCTTCGGCCCATTCACCCCCGGGTGTTTGGGCCAAGggg                                                                                >  2:2063718/1‑101 (MQ=255)
                        tAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaaggtgaa                                                                              >  1:2666089/1‑95 (MQ=255)
                         aaTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaaggtgaaa                                                                             >  2:3585925/1‑94 (MQ=16)
                         aaTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaaggtgaaa                                                                             <  2:5387312/101‑8 (MQ=255)
                         aaTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaaggtgaaa                                                                             >  1:2424577/1‑94 (MQ=16)
                          aTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAac                                                                          <  1:5656327/101‑1 (MQ=16)
                          aTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAac                                                                          >  2:4582163/1‑101 (MQ=255)
                            ttCGAATCCGTGCGCGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacac                                                                        <  1:1574542/101‑3 (MQ=255)
                               gAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacctt                                                                     <  2:2394705/101‑6 (MQ=255)
                               gAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacctt                                                                     >  1:1464273/1‑96 (MQ=12)
                               gAATCCTTGCACGAATTTTGCAAGGAGGGGAG.AGGGCCGC.TGCCCCTCCTCGCCTGCCGGCCAACTCACCCCCGGATGTTTGGGGCAA..GGTGAAacacctt                                                                     >  1:470609/1‑96 (MQ=255)
                                aaTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttg                                                                    >  2:711231/1‑95 (MQ=255)
                                aaTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttg                                                                    <  2:2876140/101‑7 (MQ=12)
                                aaTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTCGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacagcttg                                                                    >  2:792852/1‑95 (MQ=255)
                                 aTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgc                                                                   <  2:545941/101‑8 (MQ=255)
                                  tCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgct                                                                  <  2:4404361/101‑9 (MQ=255)
                                   ccTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgctt                                                                 <  2:4833660/101‑10 (MQ=9)
                                    cTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgcttt                                                                >  1:5080158/1‑91 (MQ=255)
                                    cTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgcttt                                                                <  1:1663861/101‑11 (MQ=255)
                                        cACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCt                                                          <  1:2889633/101‑1 (MQ=255)
                                         tctaaTTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCtt                                                         <  1:4534363/98‑1 (MQ=2)
                                         aCGAGTTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCtt                                                         <  2:2641926/101‑1 (MQ=2)
                                         aCGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCtt                                                         <  2:4620939/101‑1 (MQ=9)
                                          cGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTg                                                        <  1:5175985/101‑1 (MQ=255)
                                          cGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTg                                                        >  1:101160/1‑101 (MQ=255)
                                           gAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTgg                                                       >  1:2460034/1‑101 (MQ=9)
                                            aaTTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCttaga                                                      >  2:4534091/1‑98 (MQ=255)
                                                 tGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaaa                                                 <  2:448122/101‑1 (MQ=255)
                                                 tGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaaa                                                 <  1:3894203/101‑1 (MQ=255)
                                                 tGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaaa                                                 <  1:4084682/101‑1 (MQ=255)
                                                 tGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaaa                                                 >  1:4960620/1‑101 (MQ=9)
                                                  gCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCC.AAGGG.GAAA.CA.CCTTGCTTTTTCTTGGCGCAAAt                                                >  1:376425/1‑101 (MQ=255)
                                                      ggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGc                                            >  2:1906649/1‑101 (MQ=255)
                                                        aggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCct                                          <  2:1464180/101‑1 (MQ=9)
                                                          gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                        <  1:3586142/101‑1 (MQ=255)
                                                          gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                        <  2:1376288/101‑1 (MQ=255)
                                                          gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                        >  2:314333/1‑101 (MQ=9)
                                                               agaggggcgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGc                                     <  1:3006840/91‑1 (MQ=255)
                                                                gaggggcgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGcc                                    >  1:815329/10‑101 (MQ=16)
                                                                  ggggtgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCg                                  <  2:105838/94‑1 (MQ=255)
                                                                  ggggcgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCg                                  <  1:4797495/94‑1 (MQ=255)
                                                                      cgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAgg                              <  2:4064968/98‑1 (MQ=16)
                                                                      cgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAgg                              <  2:2452049/98‑1 (MQ=16)
                                                                      cgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAgg                              <  2:1586700/98‑1 (MQ=16)
                                                                       gctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAggg                             <  2:1773231/99‑1 (MQ=255)
                                                                          gCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCGTCTGGCCCGAAGGGCaa                          >  2:5195249/1‑101 (MQ=255)
                                                                          gCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCaa                          <  1:5786815/101‑1 (MQ=255)
                                                                           ccccTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCaac                         >  1:2721148/1‑100 (MQ=16)
                                                                            cccTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGc                        <  2:2903700/101‑1 (MQ=16)
                                                                            cccTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGc                        <  1:5338403/101‑1 (MQ=16)
                                                                            cccTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAACGGCAAGc                        <  2:1178357/101‑1 (MQ=16)
                                                                             ccTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGcc                       <  1:1512092/101‑1 (MQ=255)
                                                                                cctGGCCTTTCGGCCAATTCCCCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCcgcg                    <  2:4897380/99‑1 (MQ=19)
                                                                                 ctGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCgcgc                   >  2:161873/2‑101 (MQ=255)
                                                                                  tGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAgatcggaa                  >  2:5706634/1‑94 (MQ=16)
                                                                                  tGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGcc                  >  2:2169963/1‑101 (MQ=255)
                                                                                   ggCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGTTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCg                 <  1:4809667/101‑1 (MQ=19)
                                                                                   ggCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCg                 >  1:1394413/1‑101 (MQ=255)
                                                                                    gCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGa                <  2:380442/101‑1 (MQ=19)
                                                                                     ccTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGaa               <  2:2523438/101‑1 (MQ=19)
                                                                                      cTTTCGGCCAATTCCCCCCGGGATGTTTGGGCCAA..GGTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAg              <  1:1336457/101‑1 (MQ=25)
                                                                                      cTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAg              >  2:272152/1‑101 (MQ=255)
                                                                                      cTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAg              <  1:3521573/101‑1 (MQ=19)
                                                                                      cTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTCTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAg              <  2:5010192/101‑1 (MQ=255)
                                                                                       tttCGGCCAATTCCCCCCGGGATGTTTGGGCCAA..GGTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAgg             <  1:314366/101‑1 (MQ=25)
                                                                                       tttCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAgg             <  1:213777/101‑1 (MQ=19)
                                                                                       tttCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAgg             <  1:5422572/101‑1 (MQ=19)
                                                                                       tttCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAgg             >  1:3979311/1‑101 (MQ=255)
                                                                                          cggccAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGgcgc          <  2:3282864/97‑1 (MQ=21)
                                                                                          cggccAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGgcgc          >  1:226557/5‑101 (MQ=25)
                                                                                           ggccAATTCCCCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGAAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGcgcg         <  2:2555716/98‑1 (MQ=255)
                                                                                                 ttCCCCCCGGGATGTTTGGGCCAA..GGTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGccccc   <  1:1045556/101‑1 (MQ=25)
                                                                                                     |                                                                                                     
AAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAAGGGCAGCTTGCCCCTCATGGCCTTTCGG.CAATTCCACCCGGGAAGGTTGGGCCAAAGGGTGAAACCACCCTGGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCC  >  NC_014034/1696085‑1696286

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 27 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG