Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,696,217 Δ1 bp 100% intergenic (‑145/+80) 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,2170C→.100.0% 86.4 178intergenic (‑145/+80)RCAP_rcc01564/mmsAisobutyryl‑CoA dehydrogenase/methylmalonate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  new base (63/89):  ref base (14/12):  total (77/101)
Fisher's exact test strand distribution p-value = 2.86e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.89e-03

Only 100 of 184 total aligned reads displayed.

GCACGAATTTTGCAAGGAGGGGAGAAGGGCAGCTTGCCCCTCATGGCCTTTCGG.C.AATTCCACCCGGGAAGGTTGGGCCAAAGGGTGAAACCACCCTGGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCAGTCGAGCGCCTTGAAGTT  >  NC_014034/1696124‑1696317
                                                                                               |                                                                                                    
 cACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCt                                                                                         <  1:289385/101‑1 (MQ=255)
 cACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCt                                                                                         <  1:2889633/101‑1 (MQ=255)
 cACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCt                                                                                         <  1:2209872/101‑1 (MQ=255)
  aCGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCtt                                                                                        >  2:3291777/1‑101 (MQ=255)
   cGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTg                                                                                       >  1:101160/1‑101 (MQ=255)
   cGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTg                                                                                       <  1:5175985/101‑1 (MQ=255)
   cGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTg                                                                                       <  2:2720343/101‑1 (MQ=255)
      aTTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGgcg                                                                                    <  2:1421258/101‑1 (MQ=255)
       ttttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGgcgc                                                                                   >  1:5073464/1‑101 (MQ=9)
        tttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCa                                                                                  <  1:1568370/101‑1 (MQ=9)
         ttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaa                                                                                 >  1:4806610/1‑101 (MQ=255)
          tGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaaa                                                                                <  1:3894203/101‑1 (MQ=255)
          tGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaaa                                                                                <  2:448122/101‑1 (MQ=255)
          tGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaaa                                                                                >  1:4960620/1‑101 (MQ=9)
           gCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAAt                                                                               >  1:103316/1‑101 (MQ=255)
           gCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAAt                                                                               <  2:667129/101‑1 (MQ=9)
           gCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAAt                                                                               <  2:3500168/101‑1 (MQ=9)
           gCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCC.AAGGG.GAAA.CA.CCTTGCTTTTTCTTGGCGCAAAt                                                                               >  1:376425/1‑101 (MQ=255)
            cAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATa                                                                              <  2:5337219/101‑1 (MQ=255)
            cAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCCTTCGGCC.AATTCACCCCGGGATGTTTGGGCC.AAGGG.GAAA.CA.CCTTGCTTTTTCTTGGCGCAAATa                                                                              >  1:4191496/1‑101 (MQ=255)
             aaGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATAc                                                                             >  1:2315620/1‑101 (MQ=9)
             aaGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATAc                                                                             >  2:4138663/1‑101 (MQ=9)
              aggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACg                                                                            <  1:4121887/101‑1 (MQ=9)
                gaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGcc                                                                          <  1:4582439/101‑1 (MQ=255)
                 gggggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCct                                                                         <  2:4710891/100‑1 (MQ=9)
                 aggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCC.AATTCCCCCCGGGATGTTTTGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATCCTCct                                                                         >  1:792507/1‑101 (MQ=255)
                  ggggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctc                                                                        <  2:484049/101‑1 (MQ=255)
                  ggggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctc                                                                        <  2:4588254/101‑1 (MQ=255)
                   gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                                                       <  1:3586142/101‑1 (MQ=255)
                   gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                                                       >  2:314333/1‑101 (MQ=9)
                   gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                                                       <  2:1376288/101‑1 (MQ=255)
                   gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                                                       <  2:4290626/101‑1 (MQ=255)
                        agaggggcgctGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGc                                                                    >  1:612674/11‑101 (MQ=9)
                        agaggggcgctGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGc                                                                    >  2:389819/11‑101 (MQ=9)
                        agaggggcgctGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGc                                                                    <  1:3006840/91‑1 (MQ=255)
                          aggggcgctGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGccc                                                                  >  1:5508406/9‑101 (MQ=16)
                          aggggcgctGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGccc                                                                  >  2:4356471/9‑101 (MQ=16)
                           ggggcgctGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGTCCCg                                                                 <  2:5414676/94‑1 (MQ=255)
                           ggggcgctGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCg                                                                 <  1:4797495/94‑1 (MQ=255)
                            gggcgctGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGa                                                                <  2:761079/95‑1 (MQ=12)
                             ggcgctGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGaa                                                               <  1:2459546/96‑1 (MQ=255)
                             ggcgctGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGaa                                                               <  2:817268/96‑1 (MQ=255)
                              gcgctGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAg                                                              <  1:411131/97‑1 (MQ=16)
                               cgctGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAgg                                                             >  2:2899745/4‑101 (MQ=16)
                                gctGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAggg                                                            <  2:1773231/99‑1 (MQ=255)
                                gctGCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAggg                                                            >  1:852627/3‑101 (MQ=16)
                                  ttccgatcttGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCa                                                          <  1:5706929/92‑1 (MQ=16)
                                   gCCCCTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCaa                                                         <  1:5002549/101‑1 (MQ=255)
                                    ccccTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAg                                                        <  2:332332/101‑1 (MQ=255)
                                    ccccTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCaac                                                        >  1:2721148/1‑100 (MQ=16)
                                    ccccTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAg                                                        >  1:1362150/1‑101 (MQ=16)
                                    ccccTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAg                                                        <  1:5344937/101‑1 (MQ=255)
                                     cccTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGc                                                       <  1:5435314/101‑1 (MQ=16)
                                     cccTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAACGGCAAGc                                                       <  2:1178357/101‑1 (MQ=16)
                                      ccTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGcc                                                      <  2:2459900/101‑1 (MQ=255)
                                      ccTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGcc                                                      <  1:1512092/101‑1 (MQ=255)
                                       cTCCTGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCcg                                                     >  1:2585213/1‑101 (MQ=16)
                                          ctGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCgcgc                                                  >  2:161873/2‑101 (MQ=255)
                                           tGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGcc                                                 >  2:2169963/1‑101 (MQ=255)
                                           tGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGcc                                                 >  1:1577233/1‑101 (MQ=255)
                                           tGGCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGcc                                                 >  1:2523580/1‑101 (MQ=255)
                                            ggCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCg                                                >  2:5158674/1‑101 (MQ=255)
                                            ggCCTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCg                                                <  2:4529224/101‑1 (MQ=19)
                                              ccTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATCCGCCTCTGGCCCGAAGGGCAAGCCGCGCCGaa                                              <  2:2403373/101‑1 (MQ=255)
                                              ccTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGaa                                              <  2:2523438/101‑1 (MQ=19)
                                               cTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAg                                             <  1:3521573/101‑1 (MQ=19)
                                               cTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAg                                             >  1:2404714/1‑101 (MQ=19)
                                               cTTTCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTCTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAg                                             <  2:5010192/101‑1 (MQ=255)
                                               cTTTCGG.CAAATTCGCCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAg                                             <  2:1668573/101‑1 (MQ=19)
                                                tttCGGCC.AATTCCCCCCGGGATGTTTGGGCCAA..GGTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAgg                                            <  1:314366/101‑1 (MQ=25)
                                                tttCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAgg                                            >  1:1449903/1‑101 (MQ=255)
                                                tttCGGCC.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAgg                                            >  1:3979311/1‑101 (MQ=255)
                                                  tcggcc.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGgcg                                          >  2:3664880/6‑101 (MQ=255)
                                                   cggcc.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGgcgc                                         >  1:226557/5‑101 (MQ=25)
                                                    ggcc.AATTCCCCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGAAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGcgcg                                        <  2:2555716/98‑1 (MQ=255)
                                                       c.AATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGcc                                     <  2:1590470/101‑1 (MQ=25)
                                                          attcaccccGGGATGTTTGGGCCAA..GGTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGcccc                                   <  1:5374652/95‑1 (MQ=25)
                                                          attcaccccGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGcccc                                   >  1:3499306/7‑101 (MQ=255)
                                                          attcaccccGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGcccc                                   >  2:1892533/7‑101 (MQ=255)
                                                          attcaccccGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGcccc                                   >  1:4486051/7‑101 (MQ=255)
                                                          aTTCCCCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGcccc                                   <  1:3609094/101‑1 (MQ=25)
                                                            tCCCCCGGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATCCGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGcccccc                                 <  2:4577258/101‑1 (MQ=255)
                                                            tCCCCCCGGGATGTTTGGGCCAA..GGTGAAAC..CCCTTGTTTTTTCTGGGCGCAAATCCCCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGcccccc                                 <  1:3998092/101‑1 (MQ=255)
                                                             caccccGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCt                                >  2:456789/4‑101 (MQ=25)
                                                              accccGGGATGTTTGGGCCAA..GGTGAAAC..CCCTTGCTTTTTCTTGGCGAAAATACGCCTCTGGCCGGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTg                               <  1:3291977/99‑1 (MQ=255)
                                                               ccccGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGGCCCCCTGa                              >  2:2724426/2‑101 (MQ=255)
                                                                cccGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGAt                             >  1:1948850/1‑101 (MQ=255)
                                                                  cGGGATGTTTGGGCCAA..GGTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCg                           <  1:5760672/101‑1 (MQ=37)
                                                                   gggATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATCCGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGt                          <  2:1116139/101‑1 (MQ=255)
                                                                   gggATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCTGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGt                          >  1:871041/1‑101 (MQ=255)
                                                                   gggATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCgg                          >  2:2004892/1‑100 (MQ=255)
                                                                    ggatgtttGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTg                         >  1:806789/7‑101 (MQ=255)
                                                                     gatgtttGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGTCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGc                        >  1:1214005/6‑101 (MQ=255)
                                                                     gatgtttGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGc                        >  1:4744639/6‑101 (MQ=255)
                                                                      atgtttGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGGGcc                       >  2:5706908/5‑101 (MQ=255)
                                                                       tgtttGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCt                      >  2:544694/4‑101 (MQ=255)
                                                                       tgtttGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCt                      >  2:3572899/4‑101 (MQ=255)
                                                                             tgggccaaggTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATCCGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCAGt                  <  1:1697491/93‑1 (MQ=255)
                                                                                gccaaggTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCAGTCGa               <  2:3062305/96‑1 (MQ=38)
                                                                                    aggTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCAGTCGAgccc           >  2:1499860/2‑99 (MQ=255)
                                                                                               |                                                                                                    
GCACGAATTTTGCAAGGAGGGGAGAAGGGCAGCTTGCCCCTCATGGCCTTTCGG.C.AATTCCACCCGGGAAGGTTGGGCCAAAGGGTGAAACCACCCTGGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCAGTCGAGCGCCTTGAAGTT  >  NC_014034/1696124‑1696317

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