Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,696,183 C→T 7.4% intergenic (‑111/+114) 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,1830C→T 7.4% 32.8 ‑0.9 ‑1.1204intergenic (‑111/+114)RCAP_rcc01564/mmsAisobutyryl‑CoA dehydrogenase/methylmalonate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  new base (8/7):  ref base (78/111):  total (86/118)
Fisher's exact test strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.84e-02

Only 100 of 216 total aligned reads displayed.

CCAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAAGGGCAGCTTGCCCCTCATGGCCTTTC.G.G.CAATTCCA.CCCGGGAAGGTTGGGCCAAAGGGTGAAACCACCCTGGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCC  >  NC_014034/1696082‑1696285
                                                                                                        |                                                                                                       
 cAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTca                                                                                                        >  2:677068/1‑100 (MQ=255)
  aaaGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTcac                                                                                                       <  1:733058/101‑3 (MQ=255)
   aaGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.G.CCAATTCA.cc                                                                                                    <  2:3590879/101‑1 (MQ=255)
    aGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.G.CCAATTCA.ccc                                                                                                   >  2:926325/1‑101 (MQ=255)
    aGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.G.CCAATTCA.ccc                                                                                                   >  1:249543/1‑101 (MQ=255)
      tGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.G.CCAATTCA.ccccg                                                                                                 >  2:4059117/1‑99 (MQ=255)
       gCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.G.CCAATTCA.CCCcgg                                                                                                <  2:2642149/101‑1 (MQ=255)
       gCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.G.CCAATTCA.CCCcgg                                                                                                >  2:2048056/1‑101 (MQ=37)
         cTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCTCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGa                                                                                               >  2:1236542/1‑101 (MQ=255)
         cTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGa                                                                                               >  2:641394/1‑101 (MQ=255)
          tGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGat                                                                                              <  1:1676892/101‑2 (MQ=255)
           gAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGatg                                                                                             <  1:3067310/101‑3 (MQ=25)
            aaTTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGatgt                                                                                            <  2:782878/101‑4 (MQ=255)
            aaTTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGatgt                                                                                            >  1:2305139/1‑98 (MQ=255)
              ttGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGatgttt                                                                                          <  2:2924136/101‑6 (MQ=25)
               tGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGTGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCGC.CCCGGGATGTTTg                                                                                         <  1:838781/101‑1 (MQ=255)
               tGAACGTCTGTTTAATTTCGAATCCTAGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTg                                                                                         <  1:1858763/101‑1 (MQ=255)
                  aCGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGc                                                                                      <  1:3211974/101‑1 (MQ=255)
                   cGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGcc                                                                                     >  2:413712/1‑101 (MQ=25)
                   cGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGcc                                                                                     >  2:2594102/1‑101 (MQ=25)
                   cGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGcc                                                                                     <  1:2670622/101‑1 (MQ=255)
                   cGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGcc                                                                                     <  1:4059378/101‑1 (MQ=255)
                       tGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCaagg                                                                                 >  1:112763/1‑99 (MQ=255)
                        gTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCaaggt                                                                                <  1:4232443/101‑4 (MQ=255)
                        gTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAAGgg                                                                                >  1:3860893/1‑101 (MQ=255)
                          ttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCaaggtga                                                                              <  1:4081391/101‑6 (MQ=16)
                           tAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCaaggtgaa                                                                             <  1:214508/101‑7 (MQ=255)
                           tAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCaaggtgaa                                                                             >  2:2143540/1‑95 (MQ=255)
                            aaTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCaaggtgaaa                                                                            <  1:4008949/101‑8 (MQ=255)
                            aaTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCaaggtgaaa                                                                            <  1:4101421/101‑8 (MQ=255)
                            aaTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCaaggtgaaa                                                                            <  1:3794677/101‑8 (MQ=255)
                             aTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAG.GGGGGGGC.TGCCCCTCCTGGCCTTTC.GCT.CAATTTAC.CCCGGGATGATTGGGCCaaggtggaac                                                                           >  1:3463996/1‑93 (MQ=255)
                               ttCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacac                                                                       >  1:3243274/1‑99 (MQ=12)
                               ttCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacac                                                                       >  2:1742628/1‑99 (MQ=12)
                                 cGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGTGCCAA..GGTGAAacacct                                                                     >  1:1706762/1‑97 (MQ=255)
                                 cGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAAGgggaaacacct                                                                       >  1:65759/1‑93 (MQ=255)
                                 cGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacacct                                                                     >  2:2170299/1‑97 (MQ=255)
                                 cGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacacct                                                                     <  2:3330481/101‑5 (MQ=255)
                                  gAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacacctt                                                                    <  1:926410/101‑6 (MQ=255)
                                   aaTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacaccttg                                                                   <  1:1311932/101‑7 (MQ=12)
                                   aaTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacaccttg                                                                   <  2:1964853/101‑7 (MQ=12)
                                   aaTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacaccttg                                                                   <  2:3495203/101‑7 (MQ=12)
                                   aaTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacaccttg                                                                   <  1:457876/101‑7 (MQ=12)
                                    aTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacaccttgc                                                                  <  2:1296326/101‑8 (MQ=255)
                                    aTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacaccttgc                                                                  >  2:3054120/1‑94 (MQ=255)
                                       cTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacaccttgcttt                                                               >  1:2091301/1‑91 (MQ=255)
                                       cTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacaccttgcttt                                                               <  2:3216506/101‑11 (MQ=255)
                                       cTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacaccttgcttt                                                               <  1:3728606/101‑11 (MQ=255)
                                          gCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTc                                                          >  2:3743203/1‑101 (MQ=255)
                                           cACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCCC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCt                                                         <  1:2616245/101‑1 (MQ=255)
                                           cACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCt                                                         <  2:2532239/101‑1 (MQ=255)
                                               aaTTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCACTTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTggg                                                     >  1:2325921/1‑100 (MQ=255)
                                               aaTTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGgc                                                     <  1:3182710/101‑1 (MQ=9)
                                               aaTTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCCCCCGGCCTTTCTG.C.CAATTCAC.CCCCGGATGTTTGGGCCCA..GGTGGAA.CA.CCTTGCTTTTTCTTGgc                                                     >  2:809034/1‑101 (MQ=255)
                                                    tGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaaa                                                <  2:3720736/101‑1 (MQ=255)
                                                    tGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaaa                                                >  2:51160/1‑101 (MQ=9)
                                                    tGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaaa                                                <  2:4000554/101‑1 (MQ=255)
                                                     gCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAAt                                               <  1:3811798/101‑1 (MQ=9)
                                                      cAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATa                                              <  1:3495228/101‑1 (MQ=255)
                                                      cAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATa                                              <  1:1735552/101‑1 (MQ=255)
                                                       aaGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATAc                                             <  2:236406/101‑1 (MQ=255)
                                                        aggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCCA..GGCGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATTCg                                            >  2:3860695/1‑101 (MQ=255)
                                                        aggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACg                                            <  2:3685589/101‑1 (MQ=9)
                                                        aggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACg                                            >  1:3660745/1‑101 (MQ=255)
                                                         ggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGc                                           >  1:595900/1‑101 (MQ=255)
                                                          gaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGcc                                          <  2:2325753/101‑1 (MQ=255)
                                                          gaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGcc                                          <  1:413763/101‑1 (MQ=255)
                                                          gaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGcc                                          >  1:447266/1‑101 (MQ=9)
                                                           aggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GTCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCct                                         >  2:1387883/1‑101 (MQ=12)
                                                           aggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCct                                         <  2:17495/101‑1 (MQ=9)
                                                           aggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCct                                         >  2:4073574/1‑101 (MQ=12)
                                                             gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                       <  2:3104716/101‑1 (MQ=12)
                                                             gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                       >  1:150170/1‑101 (MQ=9)
                                                             gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                       <  1:1560730/101‑1 (MQ=255)
                                                                   gaggggcgctGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTTGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTTGGGCAAATACGCCTCTGGcc                                   >  1:3758137/10‑101 (MQ=255)
                                                                    aggggcgctGCCCCTCCTGGCCTTTC.G.GCCAATTCTC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTTCTTTTTCTTGGCGCAAATACGCCTCTGGccc                                  <  2:1706647/93‑1 (MQ=255)
                                                                    aggggcgctGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGccc                                  <  1:2714615/93‑1 (MQ=255)
                                                                     ggggcgctGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCg                                 <  2:3243033/94‑1 (MQ=255)
                                                                       ggcgctGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGaa                               >  1:1084837/6‑101 (MQ=16)
                                                                        gcgctGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAg                              <  1:1361192/97‑1 (MQ=255)
                                                                        gcgctGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAg                              <  2:536233/97‑1 (MQ=16)
                                                                        gcgctGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTAGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAg                              <  1:3559757/97‑1 (MQ=255)
                                                                         cgctGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAgg                             <  1:2858079/98‑1 (MQ=16)
                                                                         cgctGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAgg                             <  2:1723791/98‑1 (MQ=16)
                                                                          gctGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAggg                            <  1:2434553/99‑1 (MQ=255)
                                                                             gCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCaa                         >  2:2191513/1‑101 (MQ=19)
                                                                             gCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCaa                         >  2:162784/1‑101 (MQ=19)
                                                                              ccccTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAg                        >  1:1000892/1‑101 (MQ=19)
                                                                              ccccTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAg                        >  2:2977780/1‑101 (MQ=19)
                                                                              ccccTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTATTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAg                        >  2:2824071/1‑101 (MQ=16)
                                                                                 cTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCcg                     <  1:1186680/101‑1 (MQ=255)
                                                                                     tGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGcc                 <  1:2387071/101‑1 (MQ=19)
                                                                                      ggCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCg                >  2:1057942/1‑101 (MQ=255)
                                                                                      ggCCTTTC.G.G.CCAATCCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCg                <  2:2414690/101‑1 (MQ=19)
                                                                                       gCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGa               >  1:2248752/1‑101 (MQ=255)
                                                                                         cTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCC.AAGGG.GAAAC..CCCTTGCTTTTTCTTGGCGCAAATCCGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAg             <  1:2564489/101‑1 (MQ=255)
                                                                                                ggccAATTCCC.CCCGGGATGTTTGGGCCAA..GGTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGcgcg        <  1:3704033/98‑1 (MQ=25)
                                                                                                   cAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGcc     >  2:1878215/1‑101 (MQ=255)
                                                                                                    aaTTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGccc    >  2:855119/1‑101 (MQ=255)
                                                                                                     aTTCCC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGcccc   <  2:2025703/101‑1 (MQ=25)
                                                                                                        |                                                                                                       
CCAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAAGGGCAGCTTGCCCCTCATGGCCTTTC.G.G.CAATTCCA.CCCGGGAAGGTTGGGCCAAAGGGTGAAACCACCCTGGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCC  >  NC_014034/1696082‑1696285

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