Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,696,179 A→C 7.8% 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.8% 34.2 ‑0.8 ‑0.1208intergenic (‑107/+118)RCAP_rcc01564/mmsAisobutyryl‑CoA dehydrogenase/methylmalonate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  new base (9/8):  ref base (78/113):  total (87/121)
Fisher's exact test strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.57e-01

Only 100 of 214 total aligned reads displayed.

CCAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAAGGGCAGCTTGCCCCTCATGGCCTTTC.GG.CAATTCCA.CCCGGGAAGGTTGGGCCAAAGGGTGAAACCACCCTGGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCC  >  NC_014034/1696082‑1696283
                                                                                                   |                                                                                                         
ccAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTc                                                                                                       <  2:2210227/101‑1 (MQ=255)
 cAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTca                                                                                                      >  2:677068/1‑100 (MQ=255)
 cAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTca                                                                                                      >  2:3929199/1‑100 (MQ=255)
  aaaGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTcac                                                                                                     <  1:733058/101‑3 (MQ=255)
   aaGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GG.CCAATTCA.cc                                                                                                  <  2:3590879/101‑1 (MQ=255)
    aGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GG.CCAATTCA.ccc                                                                                                 >  2:926325/1‑101 (MQ=255)
     gTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GG.CCAATTCA.cccc                                                                                                >  2:2670416/1‑100 (MQ=255)
      tGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GG.CCAATTCA.ccccg                                                                                               >  2:4059117/1‑99 (MQ=255)
       gCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GG.CCAATTCA.CCCcgg                                                                                              >  2:2048056/1‑101 (MQ=37)
         cTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCTCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGa                                                                                             >  2:1236542/1‑101 (MQ=255)
         cTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGa                                                                                             >  2:641394/1‑101 (MQ=255)
          tGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGat                                                                                            <  1:1676892/101‑2 (MQ=255)
           gAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGatg                                                                                           <  1:3067310/101‑3 (MQ=25)
            aaTTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGatgt                                                                                          <  2:782878/101‑4 (MQ=255)
            aaTTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGatgt                                                                                          >  1:2305139/1‑98 (MQ=255)
              ttGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGatgttt                                                                                        <  2:2924136/101‑6 (MQ=25)
               tGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGTGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCGC.CCCGGGATGTTTg                                                                                       <  1:838781/101‑1 (MQ=255)
               tGAACGTCTGTTTAATTTCGAATCCTAGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTg                                                                                       <  1:1858763/101‑1 (MQ=255)
                 aaCGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCCATTCAC.CCCGGGATGTTTggg                                                                                     >  1:154355/1‑101 (MQ=255)
                  aCGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGc                                                                                    <  1:3211974/101‑1 (MQ=255)
                   cGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGcc                                                                                   <  1:4059378/101‑1 (MQ=255)
                   cGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGcc                                                                                   >  2:413712/1‑101 (MQ=25)
                     tCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCaa                                                                                 >  2:2330463/1‑101 (MQ=255)
                       tGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCaagg                                                                               >  1:112763/1‑99 (MQ=255)
                        gTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCaaggt                                                                              <  1:4232443/101‑4 (MQ=255)
                        gTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAAGgg                                                                              >  1:3860893/1‑101 (MQ=255)
                          ttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCaaggtga                                                                            <  1:4081391/101‑6 (MQ=16)
                          ttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.CGCCCCTCCTGGCCTTTC.GGCCAATTAAC.CCCGGGATTTTTGGGCCAAGggga                                                                            >  1:356887/1‑100 (MQ=255)
                           tAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCaaggtgaa                                                                           >  1:3470800/1‑95 (MQ=255)
                            aaTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCaaggtgaaa                                                                          <  1:4008949/101‑8 (MQ=255)
                            aaTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCaaggtgaaa                                                                          <  1:4101421/101‑8 (MQ=255)
                               ttCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacac                                                                     >  1:3243274/1‑99 (MQ=12)
                               ttCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacac                                                                     >  2:1742628/1‑99 (MQ=12)
                                 cGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGTGCCAA..GGTGAAacacct                                                                   >  1:1706762/1‑97 (MQ=255)
                                 cGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAAGgggaaacacct                                                                     >  1:65759/1‑93 (MQ=255)
                                  gAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacacctt                                                                  <  1:926410/101‑6 (MQ=255)
                                   aaTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacaccttg                                                                 <  1:457876/101‑7 (MQ=12)
                                   aaTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacaccttg                                                                 <  1:1311932/101‑7 (MQ=12)
                                   aaTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacaccttg                                                                 <  2:1964853/101‑7 (MQ=12)
                                   aaTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacaccttg                                                                 <  2:3495203/101‑7 (MQ=12)
                                    aTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacaccttgc                                                                >  1:2229965/1‑94 (MQ=255)
                                    aTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacaccttgc                                                                <  2:1296326/101‑8 (MQ=255)
                                    aTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacaccttgc                                                                >  2:3054120/1‑94 (MQ=255)
                                     tCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacaccttgct                                                               <  2:2154012/101‑9 (MQ=255)
                                      ccTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacaccttgctt                                                              <  1:3498769/101‑10 (MQ=9)
                                       cTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacaccttgcttt                                                             <  1:3728606/101‑11 (MQ=255)
                                       cTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacaccttgcttt                                                             >  1:2091301/1‑91 (MQ=255)
                                           cACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCt                                                       >  2:3313167/1‑101 (MQ=9)
                                           cACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCt                                                       <  2:2532239/101‑1 (MQ=255)
                                               aaTTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCACTTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTggg                                                   >  1:2325921/1‑100 (MQ=255)
                                               aaTTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGgc                                                   <  1:3182710/101‑1 (MQ=9)
                                               aaTTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCCCCCGGCCTTTCTGC.CAATTCAC.CCCCGGATGTTTGGGCCCA..GGTGGAA.CA.CCTTGCTTTTTCTTGgc                                                   >  2:809034/1‑101 (MQ=255)
                                                 ttttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGgcgc                                                 >  1:2593942/1‑101 (MQ=9)
                                                   ttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTTGCGCaa                                               >  1:3840992/1‑101 (MQ=255)
                                                   ttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaa                                               >  1:2793385/1‑101 (MQ=255)
                                                    tGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaaa                                              <  2:4000554/101‑1 (MQ=255)
                                                    tGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaaa                                              <  2:3720736/101‑1 (MQ=255)
                                                    tGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaaa                                              >  2:51160/1‑101 (MQ=9)
                                                     gCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAAt                                             <  2:2011600/101‑1 (MQ=9)
                                                      cAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATa                                            <  1:1735552/101‑1 (MQ=255)
                                                       aaGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATAc                                           <  2:236406/101‑1 (MQ=255)
                                                        aggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACg                                          <  2:3685589/101‑1 (MQ=9)
                                                        aggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACg                                          >  1:3660745/1‑101 (MQ=255)
                                                         ggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGc                                         >  1:595900/1‑101 (MQ=255)
                                                          gaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGcc                                        >  1:447266/1‑101 (MQ=9)
                                                          gagcGGAGAGGGGC.GC.TGCTCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGCGATGTTTGGGCCAA..GGTTAAAC..CCCTTGCTTTTTCTTGGCGCAAATACGcc                                        <  1:4166880/101‑1 (MQ=9)
                                                           aggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGTCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCct                                       >  2:1387883/1‑101 (MQ=12)
                                                           aggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCct                                       >  2:4073574/1‑101 (MQ=12)
                                                           aggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCct                                       <  2:17495/101‑1 (MQ=9)
                                                             gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                     >  1:150170/1‑101 (MQ=9)
                                                             gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                     <  2:3104716/101‑1 (MQ=12)
                                                             gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                     <  1:214203/101‑1 (MQ=255)
                                                             gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                     <  2:549700/101‑1 (MQ=255)
                                                                   gaggggcgctGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTTGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTTGGGCAAATACGCCTCTGGcc                                 >  1:3758137/10‑101 (MQ=255)
                                                                       ggcgctGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGaa                             >  1:1084837/6‑101 (MQ=16)
                                                                        gcgctGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAg                            <  1:1361192/97‑1 (MQ=255)
                                                                         cgctGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAgg                           <  2:1723791/98‑1 (MQ=16)
                                                                          gctGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAggg                          >  1:2482379/3‑101 (MQ=19)
                                                                          gctGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAggg                          <  1:2434553/99‑1 (MQ=255)
                                                                           ctGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGc                         >  2:2870644/2‑101 (MQ=16)
                                                                            tGCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCa                        <  2:3226963/101‑1 (MQ=16)
                                                                             gCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCaa                       >  2:162784/1‑101 (MQ=19)
                                                                             gCCCCTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCaa                       <  2:3342690/101‑1 (MQ=255)
                                                                              ccccTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGCAGGGCAAg                      >  2:1692870/1‑101 (MQ=16)
                                                                              ccccTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAg                      >  1:1000892/1‑101 (MQ=19)
                                                                              ccccTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAg                      >  2:2977780/1‑101 (MQ=19)
                                                                                 cTCCTGGCCTTTC.GGCCAATTCCC.CCCGGGATGTTTGGGCCAA..GGTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCcg                   <  1:2978002/101‑1 (MQ=21)
                                                                                 cTCCTGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCcg                   <  1:1186680/101‑1 (MQ=255)
                                                                                     tGGCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGcc               <  1:2387071/101‑1 (MQ=19)
                                                                                      ggCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCg              >  2:1057942/1‑101 (MQ=255)
                                                                                      ggCCTTTC.GGCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAAGACGCCTCTGGCCCGAAGGGCAAGCCGCGCCg              <  1:3800518/101‑1 (MQ=255)
                                                                                      ggCCTTTC.GG.CCAATCCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCg              <  2:2414690/101‑1 (MQ=19)
                                                                                              cggccAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGgcgc       >  2:2594750/5‑101 (MQ=255)
                                                                                               ggccAATTCCC.CCCGGGATGTTTGGGCCAA..GGTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGcgcg      <  1:3704033/98‑1 (MQ=25)
                                                                                                gccAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCgg     <  2:3840733/99‑1 (MQ=25)
                                                                                                 ccAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGc    >  1:2693596/2‑101 (MQ=25)
                                                                                                  cAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCTGCGCCGAAGGCGCGGcc   >  1:2508434/1‑101 (MQ=25)
                                                                                                  cAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGcc   <  2:2185880/101‑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)
                                                                                                   |                                                                                                         
CCAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAAGGGCAGCTTGCCCCTCATGGCCTTTC.GG.CAATTCCA.CCCGGGAAGGTTGGGCCAAAGGGTGAAACCACCCTGGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCC  >  NC_014034/1696082‑1696283

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