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% 43.0 368intergenic (‑113/+112)RCAP_rcc01564/mmsAisobutyryl‑CoA dehydrogenase/methylmalonate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  new base (168/176):  ref base (13/11):  total (181/187)
Fisher's exact test strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.92e-09

Only 100 of 375 total aligned reads displayed.

AAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAAGGGCAGCTTGCCCCTCATGGCCTTTCGG.CAATTCCACCCGGGAAGGTTGGGCCAAAGGGTGAAACCACCCTGGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCC  >  NC_014034/1696085‑1696286
                                                                                                     |                                                                                                     
aaGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGG.CCAATTCAcc                                                                                                     <  1:6262574/101‑1 (MQ=255)
  gTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGG.CCAATTCAcccc                                                                                                   <  2:5294909/101‑2 (MQ=37)
  gTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGG.CCAATTCAcccc                                                                                                   <  1:1324883/101‑2 (MQ=37)
  gTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGG.CCAATTCAcccc                                                                                                   >  1:2256036/1‑100 (MQ=255)
   tGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGG.CCAATTCAccccg                                                                                                  <  2:4168244/101‑3 (MQ=255)
    gCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGG.CCAATTCACCCcgg                                                                                                 >  2:5549146/1‑101 (MQ=37)
    gCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGG.CCAATTCACCCcgg                                                                                                 >  1:5647293/1‑101 (MQ=37)
      cTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGa                                                                                                <  2:5670029/101‑1 (MQ=255)
       tGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGat                                                                                               <  1:5336769/101‑2 (MQ=255)
       tGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGat                                                                                               <  2:3693553/101‑2 (MQ=255)
       tGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGat                                                                                               <  1:906203/101‑2 (MQ=255)
       tGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGan                                                                                               >  1:6392415/1‑100 (MQ=37)
         aaTTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGatgt                                                                                             >  1:4777744/1‑98 (MQ=255)
          aTTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGatgtt                                                                                            >  2:1533805/1‑97 (MQ=25)
          aTTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGatgtt                                                                                            >  2:1278140/1‑97 (MQ=25)
          aTTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGatgtt                                                                                            <  2:1037517/101‑5 (MQ=255)
           ttGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGatgttt                                                                                           <  1:6195232/101‑6 (MQ=25)
           ttGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGatgttt                                                                                           >  1:5738109/1‑96 (MQ=255)
             gAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTgg                                                                                         <  1:6634658/101‑1 (MQ=255)
             gAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTgg                                                                                         >  1:1965379/1‑101 (MQ=25)
               aCGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGatttttggcc                                                                                       >  2:296137/1‑92 (MQ=255)
               aCGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGc                                                                                       <  2:3149199/101‑1 (MQ=255)
                 gTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCa                                                                                     >  2:1113242/1‑101 (MQ=255)
                 gTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCa                                                                                     <  2:1847489/101‑1 (MQ=25)
                  tCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaa                                                                                    >  1:5846125/1‑101 (MQ=255)
                   cTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaag                                                                                   >  2:3789027/1‑100 (MQ=21)
                   cTGTTTAATTTCGAATCCTTGCACGAATTATGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaag                                                                                   <  2:4297687/101‑2 (MQ=255)
                     gTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaaggt                                                                                 >  1:4130634/1‑98 (MQ=255)
                     gTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaaggt                                                                                 >  1:1014122/1‑98 (MQ=255)
                      tttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaaggtg                                                                                >  2:1660511/1‑97 (MQ=19)
                      tttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaaggtg                                                                                >  2:786928/1‑97 (MQ=19)
                      tttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaaggtg                                                                                <  1:6135653/101‑5 (MQ=255)
                       ttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaaggtga                                                                               <  2:2676532/101‑6 (MQ=16)
                       ttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaaggtga                                                                               >  1:696536/1‑96 (MQ=255)
                       ttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaaggtga                                                                               >  1:3865330/1‑96 (MQ=255)
                        tAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaaggtgaa                                                                              >  1:2859564/1‑95 (MQ=255)
                         aaTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaaggtgaaa                                                                             >  2:1627516/1‑94 (MQ=16)
                         aaTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaaggtgaaa                                                                             >  1:1810263/1‑94 (MQ=16)
                            ttCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacac                                                                        <  2:7119290/101‑3 (MQ=255)
                            ttCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacac                                                                        <  1:4011842/101‑3 (MQ=255)
                               gAATCCTTGCACGAATTTTGCAAGGAGGTGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacctt                                                                     <  1:4972298/101‑6 (MQ=255)
                               gAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacctt                                                                     <  2:1789753/101‑6 (MQ=255)
                               gAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacctt                                                                     >  1:1055208/1‑96 (MQ=12)
                                aaTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttg                                                                    <  1:4889842/101‑7 (MQ=12)
                                    cTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgcttt                                                                <  2:5200087/101‑11 (MQ=255)
                                       gCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTc                                                           >  1:4932168/1‑101 (MQ=255)
                                         aCGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCtt                                                         <  2:4608307/101‑1 (MQ=9)
                                         aCGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCtt                                                         <  2:4812976/101‑1 (MQ=9)
                                          cGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTg                                                        >  2:5143478/1‑101 (MQ=255)
                                          cGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTg                                                        <  2:2430718/101‑1 (MQ=255)
                                          cGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTg                                                        <  2:5005772/101‑1 (MQ=255)
                                           gAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTgg                                                       >  1:5537181/1‑101 (MQ=9)
                                           gAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTgg                                                       >  1:3606297/1‑101 (MQ=9)
                                             aTTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGgcg                                                     >  1:1229799/1‑101 (MQ=255)
                                                ttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaa                                                  >  1:2166478/1‑101 (MQ=255)
                                                ttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaa                                                  >  2:221816/1‑101 (MQ=255)
                                                 tGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaaa                                                 >  1:515276/1‑101 (MQ=9)
                                                  gCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAAt                                                >  1:1775946/1‑101 (MQ=255)
                                                   cAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATa                                               <  1:1524836/101‑1 (MQ=255)
                                                     aggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACg                                             <  1:1660572/101‑1 (MQ=9)
                                                     aggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACg                                             <  2:4130439/101‑1 (MQ=9)
                                                     aggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACg                                             <  2:6718758/101‑1 (MQ=9)
                                                     aggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACg                                             >  1:2715983/1‑101 (MQ=255)
                                                      ggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGc                                            <  2:6378909/101‑1 (MQ=255)
                                                      ggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGc                                            <  1:1113277/101‑1 (MQ=255)
                                                       gaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGcc                                           <  1:2941516/101‑1 (MQ=255)
                                                       gaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGcc                                           >  1:4465210/1‑101 (MQ=9)
                                                        aggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCct                                          >  2:1395345/1‑101 (MQ=12)
                                                         ggggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctc                                         <  2:2863963/101‑1 (MQ=255)
                                                         ggggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctc                                         >  2:615866/1‑101 (MQ=12)
                                                         ggggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctc                                         <  1:1556075/101‑1 (MQ=255)
                                                          gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                        >  1:7133078/1‑101 (MQ=9)
                                                               agaggggcgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGc                                     >  2:6032570/11‑101 (MQ=9)
                                                                gaggggcgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGcc                                    >  2:3219292/10‑101 (MQ=12)
                                                                gaggggcgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGcc                                    >  2:4754991/10‑101 (MQ=16)
                                                                  ggggcgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCg                                  <  2:146182/94‑1 (MQ=255)
                                                                  ggggcgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCg                                  >  1:3298151/8‑101 (MQ=12)
                                                                    ggcgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGaa                                <  2:2912799/96‑1 (MQ=255)
                                                                     gcgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAg                               <  2:2354298/97‑1 (MQ=255)
                                                                      cgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAgg                              >  2:4445869/4‑101 (MQ=16)
                                                                       gctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAggg                             >  2:274122/3‑101 (MQ=16)
                                                                         cgatcttcctGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCa                           <  2:6847739/92‑1 (MQ=16)
                                                                          gCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCaa                          >  2:3564591/1‑101 (MQ=16)
                                                                          gCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCaa                          <  1:2574605/101‑1 (MQ=255)
                                                                          gCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCaa                          <  2:5752235/101‑1 (MQ=19)
                                                                           ccccTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAg                         <  2:2587903/101‑1 (MQ=255)
                                                                           ccccTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAg                         <  1:2371461/101‑1 (MQ=255)
                                                                              cTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCcg                      <  2:656112/101‑1 (MQ=19)
                                                                              cTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCcg                      >  1:2810338/1‑101 (MQ=19)
                                                                               tccgggCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTGTGGCCCGAAGGGCAAgatcg                     >  1:6848093/5‑97 (MQ=255)
                                                                                cctGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCcgcg                    <  2:334117/99‑1 (MQ=19)
                                                                                   ggCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCg                 >  2:1941426/1‑101 (MQ=255)
                                                                                     ccTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGaa               >  1:3095794/1‑101 (MQ=255)
                                                                                     ccTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGaa               <  1:3670572/101‑1 (MQ=19)
                                                                                      ctttcgtccAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAg              <  2:3998544/93‑1 (MQ=19)
                                                                                      cTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAg              >  2:958249/1‑101 (MQ=255)
                                                                                      cTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAg              <  1:4252521/101‑1 (MQ=19)
                                                                                           ggccAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGcgcg         <  1:3218286/98‑1 (MQ=25)
                                                                                                aTTCCCCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGcccc    <  2:5860035/101‑1 (MQ=25)
                                                                                                  tCCCCCCGGGATGTTTGGGCCAA..GGTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGcccccc  <  1:4673506/101‑1 (MQ=25)
                                                                                                     |                                                                                                     
AAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAAGGGCAGCTTGCCCCTCATGGCCTTTCGG.CAATTCCACCCGGGAAGGTTGGGCCAAAGGGTGAAACCACCCTGGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCC  >  NC_014034/1696085‑1696286

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