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% 16.6 202intergenic (‑113/+112)RCAP_rcc01564/mmsAisobutyryl‑CoA dehydrogenase/methylmalonate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  new base (76/115):  ref base (7/4):  total (83/119)
Fisher's exact test strand distribution p-value = 2.06e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.59e-04

Only 100 of 207 total aligned reads displayed.

AAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAAGGGCAGCTTGCCCCTCATGGCCTTTCG.G.CAATTCCACCCGGGAAGGTTGGGCCAAAGGGTGAAACCACCCTGGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCC  >  NC_014034/1696085‑1696285
                                                                                                      |                                                                                                    
aaGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.G.CCAATTCAcc                                                                                                    >  1:1591287/1‑101 (MQ=255)
 aGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.G.CCAATTCAccc                                                                                                   >  1:204179/1‑101 (MQ=255)
 aGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.G.CCAATTCAccc                                                                                                   <  1:1958533/101‑1 (MQ=255)
  gTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.G.CCAATTCAcccc                                                                                                  >  1:490578/1‑100 (MQ=255)
  gTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.G.CCAATTCAcccc                                                                                                  >  2:4273049/1‑100 (MQ=255)
  gTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.G.CCAATTCAcccc                                                                                                  <  1:1261833/101‑2 (MQ=37)
   tGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.G.CCAATTCAccccg                                                                                                 >  2:2570260/1‑99 (MQ=255)
     ccTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCggg                                                                                                <  2:575909/101‑1 (MQ=37)
      cTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGa                                                                                               >  1:2136846/1‑101 (MQ=255)
       tGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GACAATTCACCCCGGGat                                                                                              >  2:2531876/1‑100 (MQ=255)
        tctttGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGatg                                                                                             <  2:2056095/98‑3 (MQ=25)
         aaTTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGatgt                                                                                            <  2:81073/101‑4 (MQ=255)
           ttGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGatgttt                                                                                          >  1:2409147/1‑96 (MQ=255)
            tGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTg                                                                                         <  1:2823804/101‑1 (MQ=255)
            cgAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTg                                                                                         >  1:2432460/2‑101 (MQ=255)
             gAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTgg                                                                                        >  2:4292133/1‑101 (MQ=25)
             gAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTgg                                                                                        <  2:2238252/101‑1 (MQ=255)
             gAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTgg                                                                                        <  1:4360063/101‑1 (MQ=255)
               aCGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGc                                                                                      >  1:2653741/1‑101 (MQ=255)
                cGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGcc                                                                                     >  1:4292204/1‑101 (MQ=25)
                cGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGcc                                                                                     <  2:3632103/101‑1 (MQ=255)
                 gTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCa                                                                                    <  2:671681/101‑1 (MQ=25)
                 gTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCa                                                                                    >  2:4168684/1‑101 (MQ=255)
                 gTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCa                                                                                    >  2:1175606/1‑101 (MQ=255)
                   cTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGCTGTTTGGGGCaag                                                                                  >  2:1649210/1‑100 (MQ=19)
                   cTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCaag                                                                                  <  1:4190558/101‑2 (MQ=255)
                   cTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCaag                                                                                  <  1:309206/101‑2 (MQ=255)
                    tGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAG.AGGGGGGC.TGCCCCTCCTGGCCCTTCGCC.CAACTCACCCCGGGATGTTTGGGTCaggg                                                                                 >  1:216287/1‑98 (MQ=255)
                    cgTTGAAGTTCGAAGCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGGCCCTCCTGGCCTTTCG.GCCAATTCCCCCCGGGATGTTTGGGACaagg                                                                                 <  2:194412/100‑3 (MQ=16)
                     gTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCaaggt                                                                                >  2:3569322/1‑98 (MQ=255)
                      tttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGcccaggtg                                                                               >  1:1016181/1‑95 (MQ=19)
                      tttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCaaggtg                                                                               <  1:808522/101‑5 (MQ=255)
                      tttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCaaggtg                                                                               <  1:2871673/101‑5 (MQ=255)
                      tttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCaaggtg                                                                               <  1:3938231/101‑5 (MQ=255)
                       ttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCaaggtga                                                                              <  1:2240324/101‑6 (MQ=16)
                       ttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCaaggtga                                                                              <  1:4215441/101‑6 (MQ=16)
                        tAATTTCGAATCTTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCaaggtgaa                                                                             <  2:1708772/101‑7 (MQ=255)
                        tAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCaaggtgaa                                                                             <  2:1618056/101‑7 (MQ=255)
                        tAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCaaggtgaa                                                                             <  2:2120018/101‑7 (MQ=255)
                         aaTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCaaggtgaaa                                                                            <  1:3431316/101‑8 (MQ=255)
                          aTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAac                                                                         <  2:2837212/101‑1 (MQ=16)
                           tttCGAATCCTTGCACGAATTTTGCAAGGAGGGGAG.AGGGGGGC.TGCCCCCCCTGGCCTTTCG.GCCAATTCATCCCCGGGTGTTTGGGCCAAGgggaaaaa                                                                          >  1:117431/1‑96 (MQ=255)
                             tCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacc                                                                      >  2:3529503/1‑98 (MQ=255)
                             tCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacc                                                                      >  1:2768825/1‑98 (MQ=255)
                              cGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacct                                                                     >  2:3693851/1‑97 (MQ=255)
                               gAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacctt                                                                    <  1:260799/101‑6 (MQ=255)
                                aaTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttg                                                                   <  2:3597375/101‑7 (MQ=12)
                                aaTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttg                                                                   <  1:4141894/101‑7 (MQ=12)
                                  tCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgct                                                                 <  2:20199/101‑9 (MQ=255)
                                  tCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgct                                                                 >  2:1856447/1‑93 (MQ=12)
                                   ccTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgctt                                                                <  1:4033649/101‑10 (MQ=9)
                                   ccTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgctt                                                                <  2:2202451/101‑10 (MQ=9)
                                    cTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgcttt                                                               <  1:571807/101‑11 (MQ=255)
                                        cACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCt                                                         <  1:2635474/101‑1 (MQ=255)
                                          cGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTg                                                       <  1:923827/101‑1 (MQ=255)
                                           gAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTgg                                                      <  2:1074620/101‑1 (MQ=255)
                                            aaTTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGgc                                                     <  1:4426896/101‑1 (MQ=9)
                                              ttttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGGTGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGgcgc                                                   >  1:1685177/1‑101 (MQ=9)
                                               tttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCa                                                  <  1:3147993/101‑1 (MQ=9)
                                                  gCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAAt                                               <  2:3218898/101‑1 (MQ=9)
                                                  gCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAAt                                               >  1:4349251/1‑101 (MQ=255)
                                                    aaGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATAc                                             <  1:776373/101‑1 (MQ=255)
                                                    aaGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATAc                                             >  2:1462852/1‑101 (MQ=9)
                                                     aggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACg                                            >  2:438421/1‑101 (MQ=255)
                                                      ggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGc                                           <  1:3694104/101‑1 (MQ=255)
                                                       gaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGcc                                          <  1:1354249/101‑1 (MQ=255)
                                                       gaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGcc                                          >  2:344212/1‑101 (MQ=9)
                                                       gaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGcc                                          <  1:997578/101‑1 (MQ=255)
                                                          gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                       >  1:3065993/1‑101 (MQ=9)
                                                          gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                       <  2:3065772/101‑1 (MQ=255)
                                                                gaggggcgctGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGcc                                   <  1:4168945/92‑1 (MQ=9)
                                                                gaggggcgctGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGcc                                   <  2:3136530/92‑1 (MQ=9)
                                                                 aggggcgcttccccTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGccc                                  <  2:1759983/91‑1 (MQ=255)
                                                                 aggggcgctGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGccc                                  <  1:4241725/93‑1 (MQ=255)
                                                                 aggggcgctGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGccc                                  >  2:46779/9‑101 (MQ=16)
                                                                     gcgctGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAg                              <  1:2348519/97‑1 (MQ=255)
                                                                     gcgctGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAg                              <  1:1649316/97‑1 (MQ=255)
                                                                      cgctGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAgg                             <  2:3560932/98‑1 (MQ=16)
                                                                       gctGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAggg                            >  1:2198103/3‑101 (MQ=16)
                                                                        ctGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGc                           >  1:1511444/2‑101 (MQ=16)
                                                                        ctGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGc                           >  2:1163031/2‑101 (MQ=16)
                                                                         tGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCa                          >  2:3732371/1‑101 (MQ=19)
                                                                         tGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCa                          <  1:1792210/101‑1 (MQ=16)
                                                                         tGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCa                          >  2:132586/1‑101 (MQ=255)
                                                                           ccccTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAg                        <  2:909483/101‑1 (MQ=255)
                                                                           ccccTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAg                        <  2:3371320/101‑1 (MQ=255)
                                                                            cccTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGc                       <  2:628199/101‑1 (MQ=16)
                                                                            cccTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGc                       >  1:2066798/1‑101 (MQ=255)
                                                                               tcctGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCcgc                    >  2:3470728/4‑101 (MQ=255)
                                                                                  tGGCCTTTCG.GCCAATTCCCCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGcc                 <  1:3535242/101‑1 (MQ=25)
                                                                                   ggCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCg                >  2:3266173/1‑101 (MQ=255)
                                                                                   ggCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCg                >  1:4136000/1‑101 (MQ=255)
                                                                                     ccTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGaa              <  1:255591/101‑1 (MQ=19)
                                                                                       tttCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAgg            >  1:612002/1‑101 (MQ=19)
                                                                                          tcggccAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGgcg          <  1:761875/96‑1 (MQ=19)
                                                                                           cggccAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGgcgc         >  1:3464967/5‑101 (MQ=255)
                                                                                            ggccAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGcgcg        >  2:3793398/4‑101 (MQ=255)
                                                                                            ggccAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGcgcg        >  2:3989351/4‑101 (MQ=255)
                                                                                               cAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGcc     >  1:1977741/1‑101 (MQ=255)
                                                                                                aaTTCCCCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGccc    <  1:4248842/101‑1 (MQ=25)
                                                                                                      |                                                                                                    
AAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAAGGGCAGCTTGCCCCTCATGGCCTTTCG.G.CAATTCCACCCGGGAAGGTTGGGCCAAAGGGTGAAACCACCCTGGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCC  >  NC_014034/1696085‑1696285

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