Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,696,179 A→C 4.4% 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 4.4% 40.8 0.0 ‑0.1469intergenic (‑107/+118)RCAP_rcc01564/mmsAisobutyryl‑CoA dehydrogenase/methylmalonate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  new base (11/10):  ref base (227/221):  total (238/231)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.55e-01

Only 100 of 480 total aligned reads displayed.

CCAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAAGGGCAGCTTGCCCCTCATGGCCTTTC.G.G.CAATTCCA.CCCGGGAAGGTTGGGCCAAAGGGTGAAACCACCCTGGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCC  >  NC_014034/1696082‑1696283
                                                                                                    |                                                                                                         
 cAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCTG.C.CAATTca                                                                                                      >  1:5553873/1‑100 (MQ=255)
  aaaGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTcac                                                                                                     <  1:5286059/101‑3 (MQ=255)
  aaaGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTcac                                                                                                     <  1:4653035/101‑3 (MQ=255)
  aaaGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTcac                                                                                                     <  1:7921472/101‑3 (MQ=255)
  aaaGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTcac                                                                                                     >  2:8330812/1‑99 (MQ=255)
   aaGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.G.CCAATTCA.cc                                                                                                  >  1:8107541/1‑101 (MQ=255)
   aaGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.G.CCAATTCA.cc                                                                                                  >  1:6684390/1‑101 (MQ=255)
    aGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.G.CCAATTCA.ccc                                                                                                 <  2:2536666/101‑1 (MQ=255)
     gTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.G.CCAATTCA.cccc                                                                                                <  1:645836/101‑2 (MQ=37)
       gCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.G.CCAATTCA.CCCcgg                                                                                              >  2:759173/1‑101 (MQ=37)
        ccTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCggg                                                                                              <  2:1735253/101‑1 (MQ=37)
         cTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGa                                                                                             <  2:2348854/101‑1 (MQ=255)
         cTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGa                                                                                             <  2:5858495/101‑1 (MQ=255)
           gAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTACC.CCCGGGatg                                                                                           >  2:4232428/1‑99 (MQ=255)
            aaTTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGatgt                                                                                          <  1:1333073/101‑4 (MQ=255)
             aTTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGatgtt                                                                                         >  1:7000019/1‑97 (MQ=25)
               tGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTg                                                                                       >  1:5926602/1‑101 (MQ=255)
               tGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTg                                                                                       <  2:357226/101‑1 (MQ=255)
               tGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTg                                                                                       <  1:2194374/101‑1 (MQ=255)
               tGAACGTCTGTTTAATGTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTg                                                                                       <  2:1987588/101‑1 (MQ=255)
                 aaCGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTggg                                                                                     <  2:3042882/101‑1 (MQ=255)
                    gTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCa                                                                                  >  1:8629527/1‑101 (MQ=255)
                      cTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCaag                                                                                <  1:6806929/101‑2 (MQ=255)
                        gTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCaaggt                                                                              <  1:4046330/101‑4 (MQ=255)
                         tttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCaaggtg                                                                             <  1:3171981/101‑5 (MQ=255)
                         tttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCaaggtg                                                                             >  1:937686/1‑97 (MQ=19)
                          ttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCaaggtga                                                                            >  1:7179833/1‑96 (MQ=255)
                           tAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCaaggtgaa                                                                           <  1:3432965/101‑7 (MQ=255)
                           tAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCaaggtgaa                                                                           <  2:4820245/101‑7 (MQ=255)
                            aaTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCaaggtgaaa                                                                          <  1:5807632/101‑8 (MQ=255)
                             aTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAac                                                                       <  1:6180845/101‑1 (MQ=16)
                              tttCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAaca                                                                      <  1:3918602/101‑2 (MQ=255)
                                 cGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacacct                                                                   >  1:474095/1‑97 (MQ=255)
                                  gAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacacctt                                                                  >  1:9120559/1‑96 (MQ=12)
                                  gAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacacctt                                                                  <  2:4512623/101‑6 (MQ=255)
                                   aaTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacaccttg                                                                 >  1:3768502/1‑95 (MQ=255)
                                    aTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacaccttgc                                                                >  2:7690978/1‑94 (MQ=255)
                                     tCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacaccttgct                                                               <  1:4390118/101‑9 (MQ=255)
                                       cTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTTGGCCAA..GGTGAAacaccttgcttt                                                             >  2:4543699/1‑91 (MQ=255)
                                       cTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAacaccttgcttt                                                             <  2:6073741/101‑11 (MQ=255)
                                           cACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCt                                                       <  1:1117682/101‑1 (MQ=255)
                                           cACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCt                                                       <  1:1745840/101‑1 (MQ=255)
                                            aCGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCtt                                                      <  1:2400185/101‑1 (MQ=9)
                                              gAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTgg                                                    <  2:1375216/101‑1 (MQ=255)
                                              gAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTgg                                                    >  1:4839558/1‑101 (MQ=9)
                                                aTTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGgcg                                                  >  2:3445312/1‑101 (MQ=255)
                                                aTTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGgcg                                                  >  2:4674996/1‑101 (MQ=255)
                                                 ttttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGgcgc                                                 >  1:9000840/1‑101 (MQ=9)
                                                   ttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaa                                               >  2:5131253/1‑101 (MQ=255)
                                                   ttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGACAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaa                                               >  1:3331828/1‑101 (MQ=255)
                                                    tGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaaa                                              <  1:186579/101‑1 (MQ=255)
                                                    tGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaaa                                              >  2:1655709/1‑101 (MQ=9)
                                                    tGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaaa                                              <  2:92221/101‑1 (MQ=255)
                                                     gCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAAt                                             >  1:5488121/1‑101 (MQ=255)
                                                     gCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAAt                                             <  1:4277184/101‑1 (MQ=9)
                                                     gCAAGGAGGGGAGAGGGGC..CCTGCCCCGCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAAt                                             <  2:7048738/101‑1 (MQ=9)
                                                      cAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATa                                            <  2:474087/101‑1 (MQ=255)
                                                      cAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATa                                            >  2:6671757/1‑101 (MQ=255)
                                                      cAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATa                                            >  2:8145966/1‑101 (MQ=255)
                                                       aaGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATAc                                           >  2:4399481/1‑101 (MQ=9)
                                                       aaGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATAc                                           <  2:3255654/101‑1 (MQ=255)
                                                       aaGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATAc                                           <  1:8331175/101‑1 (MQ=255)
                                                       aaGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATAc                                           >  1:6150031/1‑101 (MQ=9)
                                                        aggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACg                                          >  1:97208/1‑101 (MQ=255)
                                                         ggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGc                                         >  1:4434701/1‑101 (MQ=255)
                                                           aggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCct                                       <  2:6914665/101‑1 (MQ=9)
                                                           aggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCct                                       <  1:5721986/101‑1 (MQ=9)
                                                           aggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCct                                       <  1:1576521/101‑1 (MQ=9)
                                                           aggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCct                                       >  2:7878682/1‑101 (MQ=12)
                                                             gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                     >  1:4540461/1‑101 (MQ=9)
                                                             gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                     <  2:1142843/101‑1 (MQ=255)
                                                             gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                     <  2:8804757/101‑1 (MQ=255)
                                                             gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                     <  2:8629141/101‑1 (MQ=255)
                                                                  agaggggcgctGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGc                                  <  2:5330525/91‑1 (MQ=255)
                                                                        gcgctGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAg                            <  1:8457017/97‑1 (MQ=255)
                                                                         cgctGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAgg                           <  1:2787232/98‑1 (MQ=16)
                                                                         cgctGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCCAAgg                           >  1:2942640/4‑101 (MQ=16)
                                                                           ctGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGc                         <  1:1264021/100‑1 (MQ=255)
                                                                            tGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCa                        <  2:642616/101‑1 (MQ=16)
                                                                            tGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCa                        >  2:4078980/1‑101 (MQ=19)
                                                                            tGCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCa                        <  2:5614868/101‑1 (MQ=16)
                                                                             gCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCaa                       >  1:7503979/1‑101 (MQ=19)
                                                                             gCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCaa                       <  2:3768393/101‑1 (MQ=255)
                                                                             gCCCCTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCaa                       >  2:221261/1‑101 (MQ=255)
                                                                               cccTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGc                     >  1:7420814/1‑101 (MQ=255)
                                                                               cccTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGc                     <  2:6181029/101‑1 (MQ=16)
                                                                               cccTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGc                     >  2:5220212/1‑101 (MQ=19)
                                                                               cccTCCTGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGc                     <  1:7723289/101‑1 (MQ=16)
                                                                                  tcctGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCcgc                  >  2:5395481/4‑101 (MQ=255)
                                                                                  tcctGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCcgc                  >  1:2575379/4‑101 (MQ=255)
                                                                                  tcctGGCCTTTC.G.GCCAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGCCCCGAAGGGCAAGAcgc                  >  1:360605/4‑101 (MQ=16)
                                                                                       gCCTTTC.GCC.CAATTCCC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGa             <  2:8089587/101‑1 (MQ=19)
                                                                                           ttC.G.GCCAATTCCC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGgc         <  2:1831803/101‑1 (MQ=25)
                                                                                             c.G.G.CCAATCCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGgcgc       <  2:9014228/101‑1 (MQ=19)
                                                                                               cggccAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGgcgc       <  2:4984176/97‑1 (MQ=25)
                                                                                                 gccAATCCCC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCgg     <  2:2445342/99‑1 (MQ=25)
                                                                                                  ccAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGc    >  1:4976106/2‑101 (MQ=255)
                                                                                                   cAATTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGcc   >  2:2695705/1‑101 (MQ=255)
                                                                                                    aaTTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAAC..CCCTTGCTTTTTCTGGGCGCAAATCCGCCTCTGGCCCGAAGGCCAAGCCGCGCCGAAGGCGCGGccc  <  1:2019150/101‑1 (MQ=255)
                                                                                                    aaTTCAC.CCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGAGCAAGCCGCGCCGAAGGCCAGGccc  >  2:7234286/1‑101 (MQ=19)
                                                                                                    |                                                                                                         
CCAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAAGGGCAGCTTGCCCCTCATGGCCTTTC.G.G.CAATTCCA.CCCGGGAAGGTTGGGCCAAAGGGTGAAACCACCCTGGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCC  >  NC_014034/1696082‑1696283

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