Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 3,439,792 G→A 2.5% intergenic (‑95/+68) RCAP_rcc03240 ← / ← mutB ATPase AAA/methylmalonyl‑CoA mutase large subunit

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140343,439,7920G→A 2.5% 8.6 ‑2.2 0.0375intergenic (‑95/+68)RCAP_rcc03240/mutBATPase AAA/methylmalonyl‑CoA mutase large subunit
Reads supporting (aligned to +/- strand):  new base (3/7):  ref base (205/159):  total (208/167)
Fisher's exact test strand distribution p-value = 1.16e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 471 total aligned reads displayed.

GCTCATGGCGGCGACCCTACACGGCGCGACCAGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCCCGCGCCGCCACGGTTCACGCCGCGACCGTCTCCGGGGTCGCCAGCGCGA  >  NC_014034/3439692‑3439894
                                                                                                     |                                                                                                         
  tCATGGCGGCGACCCTACACGGCGCGACCAGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGGGA.TTTTGGGGAAGAAAAAaggg                                                                                                         >  1:7751015/1‑101 (MQ=255)
     tGGCGGCGACCCTACACGGCGCGACCAGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGgctt                                                                                                      <  1:3174106/101‑1 (MQ=255)
      ggcggcGACCCTACACGGCGCGACCAGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGGTA.TATTGGGTCAGCCAAATGGgcttt                                                                                                     >  2:1679041/1‑101 (MQ=255)
      ggcggcGACCCTACACGGCGCGACCAGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGgcttt                                                                                                     <  2:8435879/101‑1 (MQ=255)
       gcggcgACCCTACACGGCGCGACCAGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGgctttg                                                                                                    <  2:3432869/101‑1 (MQ=255)
       gcggcgACCCTACACGGCGCGACCAGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGgctttg                                                                                                    <  2:6898552/101‑1 (MQ=255)
        cggcgACCCTACACGGCGCGACCAGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGgctttgc                                                                                                   <  1:8775360/101‑1 (MQ=255)
         ggcgACCCTACACGGCGCGACCAGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGGTATTTTGGGGGAGAAAAAAGGGggtttgt                                                                                                   >  1:6343145/1‑100 (MQ=255)
         ggcgACCCTACACGGCGCGACCAGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGgctttgct                                                                                                  <  2:863565/101‑1 (MQ=255)
          gcgACCCTACACGGCGCGACCAGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGgctttgctt                                                                                                 <  2:3768663/101‑1 (MQ=255)
            gACCCTACACGGCGCGACCAGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGGTA.TTTTGGGGAAGATAAAGGGGGTTTGGtttt                                                                                               >  2:7897683/1‑101 (MQ=255)
            gACCCTACACGGCGCGACCAGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGGCTTTGCtttt                                                                                               <  2:2850977/101‑1 (MQ=255)
             aCCCTACCCGGCGCGACCAGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTa                                                                                              <  1:336593/101‑1 (MQ=255)
             aCCCTACACGGCGCGACCAGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGGGA.TTTTGGGGAAGATAAAAGGGGGTTGCttttt                                                                                              >  1:1672717/1‑100 (MQ=255)
             aCCCTACACGGCGCGACCAGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTa                                                                                              <  2:6703929/101‑1 (MQ=255)
              cccTACACGGCGCGACCAGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGGGA.TTTGGGGGAAGATAAAAGGGCTTTGCtttttt                                                                                             >  1:5206304/1‑99 (MQ=255)
               ccTACACGGCGCGACCAGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGGTA.TTTTGGGCAAGATAAAGGGGGTTTGCTTTTAtt                                                                                            >  1:2341124/1‑100 (MQ=255)
                cTACCCGGCGCGACCAGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGAGGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCt                                                                                           <  1:8456197/101‑1 (MQ=255)
                cTACACGGCGCGACCAGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGGAA.TTTTGGGCAAGGTAAAAGGGGTTTGTTTTTTTCt                                                                                           >  1:2205568/1‑101 (MQ=255)
                   cacGGCGCGACCAGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTc                                                                                        <  2:1520589/101‑1 (MQ=255)
                    acGGCGCGACCAGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTcc                                                                                       <  2:560240/101‑1 (MQ=255)
                     cGGCGCGACCAGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTccc                                                                                      <  1:8275740/101‑1 (MQ=255)
                     cGGCGCGACCAGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTccc                                                                                      <  1:809732/101‑1 (MQ=255)
                        cgcgACCAGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCaa                                                                                   <  2:2525773/101‑1 (MQ=255)
                          cgACCAGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAat                                                                                 <  2:3286458/101‑1 (MQ=255)
                          cgACCAGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAat                                                                                 <  1:6823556/101‑1 (MQ=255)
                          cgACCAGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAat                                                                                 <  1:407693/101‑1 (MQ=255)
                            aCCAGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAata.t                                                                              <  1:3724368/101‑1 (MQ=255)
                            aCCAGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAata.t                                                                              <  2:5975025/101‑1 (MQ=255)
                             ccAGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.Tc                                                                             <  1:5897177/101‑1 (MQ=255)
                               aGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.Tccc                                                                           <  2:8034285/101‑1 (MQ=255)
                               aGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.Tccc                                                                           <  1:8646745/101‑1 (MQ=255)
                               aGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGAAA.TTTGGGGGAAGATAAAGGGGCTTTGCTTTTATCTTTTCCCAAATA.Tccc                                                                           >  1:5589487/1‑101 (MQ=255)
                                gCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGGGTTTGCTTTTATCTTTCCCCCAATA.TTCCg                                                                          >  1:7374924/1‑101 (MQ=255)
                                    tCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGGAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TACCggggg                                                                      >  2:4283642/1‑101 (MQ=255)
                                      gAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGGTA.TTTGGGGGAAGATAAAAGGGCTTTGCTTTTTTCTTTTCCCCAATA.TCCCggggggg                                                                    >  2:2330136/1‑99 (MQ=37)
                                      gAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGGTA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCggggggg                                                                    >  1:4937359/1‑99 (MQ=255)
                                      gAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGGTA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTg                                                                    >  2:777416/1‑101 (MQ=255)
                                       aaaGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTTGGGCAAGATAAAGGGGCTTTGCTTTTTTCTTTTCCCCAATATTCCCggggggc                                                                    >  2:8234110/1‑99 (MQ=37)
                                       aaaGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGc                                                                   <  1:4360606/101‑1 (MQ=255)
                                        aaGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..g                                                                <  1:1355211/101‑1 (MQ=255)
                                           aCAGGGATTCGGCCTGCGGCCGATGCCTCCGGCGGGGAAA.TTTGGGGCAAGATAAAGGGGGTTTGCTTTTTTTTTTCCCCAAAAA.TCCCGGGGGGGC..gggg                                                             >  1:5453641/1‑101 (MQ=37)
                                            cAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGGTA.TTTTGGGCAAGATAAAGGGGCTTTGCTTTTTTCTTTCCCCCAAAA.TCCCGGGGGGGC..ggggg                                                            >  2:8848162/1‑101 (MQ=37)
                                             aGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGc                                                           <  1:8542827/101‑1 (MQ=255)
                                               gAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCag                                                         <  1:580298/101‑1 (MQ=255)
                                                aaTTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTTGGGGAAGATAAAAGGGCTTTGCTTTTTTTTTTCCCCAAAAA.TTCCGGGGGGGG..GGGGGCagg                                                        >  1:7906068/1‑100 (MQ=25)
                                                aaTTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGGAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGCGGGGGC..GGGGGCgga                                                        >  2:9019060/1‑101 (MQ=37)
                                                aaTTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTGCCCCAAATA.TCCCGGGGGTGC..GGGGGCaga                                                        <  2:3008037/101‑1 (MQ=255)
                                                   tCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGcc                                                     <  1:9009695/101‑1 (MQ=255)
                                                    cGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGGAAGATAAAAGGGGTTTGCTTTTATCTTTTCCCCAACA.TCCCGGGGGGGC..GGGGGCAGAGccc                                                    >  1:8023470/1‑101 (MQ=37)
                                                     ggCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTTGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTTCCCAAAAA.TCCCGGGGGGGC..GGGGGCAGAGcccc                                                   >  2:2076442/1‑101 (MQ=25)
                                                      gCCTGCGGCCGATGCCTCCGGCGGGGGTA.TTTGGGGCAAGATAAAGGGGGTTTGCTTTTATCTTTCCCCAAATA.TCCCCGGGGGGC..GGGGGGAGAGccccc                                                  >  2:5425203/1‑101 (MQ=37)
                                                      gCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCCAATA.TCCCGGGGGTGC..GGGGGCAGAGccccc                                                  >  2:7110040/1‑101 (MQ=255)
                                                        tgactttcGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATGAA.ACCCATCCCTTTCATCTTGCCCCAAATA.TCCCGGGGGGCCCGGGGGGCAGCGcccccc                                                 <  1:66026/94‑2 (MQ=25)
                                                        tgactttcGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATGAA.ACCCATCCCTTTCATCTTGCCCCAAATA.TCCCGGGGGGCCCGGGGGGCAGCGcccccc                                                 <  1:4760667/94‑2 (MQ=25)
                                                        tgactttcGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATGAA.ACCCATCCCTTTCATCTTGCCCCAAATA.TCCCGGGGGGCCCGGGGGGCAGCGcccccc                                                 <  2:257900/94‑2 (MQ=25)
                                                         tGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCCGGGGGGC..GGGGGCAGAGCCCCcgcc                                               >  2:4489460/1‑100 (MQ=255)
                                                               cGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGGCGAGCCCCcgccgccacg                                         >  2:1898846/1‑94 (MQ=21)
                                                               cGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGGGC..GGGGGGAGAGCCCCcgccgccccg                                         >  2:5026146/1‑94 (MQ=21)
                                                                gATGCCTCCGGCGGGGCTA.CTTGGGGGTAGACAAATGGGCTTTTTTTCTTACTTTCCCCCAATA.TCCCTGGGGGGC..GGGGGC.GGGCCCCC.CGCCGCCAAgg                                      >  2:3196534/1‑101 (MQ=255)
                                                                gATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACgg                                      >  2:4856636/1‑101 (MQ=255)
                                                                 aTGCCTCCGGCGGGGATA.TTTGGGGCAAGATGAA.ACCCCTCCCTTTCATCTTGCCCCAAATA.TCCCGGGGGGCCCGGGGGGCAGCGcccccccccccttcc                                        >  1:8944819/1‑91 (MQ=18)
                                                                 aTGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCCcgccgccccggt                                       >  1:3597661/1‑92 (MQ=255)
                                                                   gCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTc                                   >  1:5037213/1‑101 (MQ=255)
                                                                           cggGGATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTCACGCcgcg                           >  2:4682804/1‑101 (MQ=255)
                                                                            ggggATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTCACGCCGCGa                          >  1:1030055/1‑101 (MQ=255)
                                                                             gggATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCGCC..CGCCGCCACGGTTCACGCCGCGAc                         >  1:615585/1‑101 (MQ=255)
                                                                              ggATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTCACGCCGCGAcc                        >  1:124399/1‑101 (MQ=255)
                                                                               gATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTCACGCCGCGACCg                       >  1:920765/1‑101 (MQ=255)
                                                                               gATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTCACGCCGCGACCg                       >  1:2505487/1‑101 (MQ=255)
                                                                                ata.tTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTCACGCCGCGACCGt                      >  2:6552024/1‑101 (MQ=255)
                                                                                ata.tTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTCACGCCGCGACCGt                      >  2:6227201/1‑101 (MQ=255)
                                                                                 ta.tTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTCACGCCGCGACCGtc                     >  2:3625453/1‑101 (MQ=255)
                                                                                 ta.tTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTCACGCCGCGACCGtc                     >  1:459610/1‑101 (MQ=255)
                                                                                    gttGGGGAAAAAAAAAGGGG.TTTTTTTTTTTTTCTCCCCCAAAA.TCCGGGGGGGGG..GGGGGGCCCGCCCCCCCGCCGCCACGGTTCACGCCGCGACCGtct                    <  1:7132452/100‑1 (MQ=255)
                                                                                      tGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTCACGCCGCGACCGTCTCCg                 >  2:3995226/1‑101 (MQ=255)
                                                                                      tGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTCACGCCGCGACCGTCTCCg                 >  1:2240371/1‑101 (MQ=255)
                                                                                      tGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTCACGCCGCGACCGTCTCCg                 >  2:2704565/1‑101 (MQ=255)
                                                                                      tGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTCACGCCGCGACCGTCTCCg                 >  2:1773490/1‑101 (MQ=255)
                                                                                       ggggCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTCACGCCGCGACCGTCTCCgg                >  2:8288355/1‑101 (MQ=255)
                                                                                        gggCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTCACGCCGCGACCGTCTCCggg               >  1:7057766/1‑101 (MQ=255)
                                                                                           cAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTCACGCCGCGACCGTCTCCGGGGTc            >  1:1865527/1‑101 (MQ=255)
                                                                                           cAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTCACGCCGCGACCGTCTCCGGGGTc            >  2:6940845/1‑101 (MQ=255)
                                                                                            aaGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTCACGCCGCGACCGTCTCCGGGGGCg           >  1:2969677/1‑101 (MQ=255)
                                                                                              gATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTCACGCCGGGACCGTCTCCGGGGTCGcc         >  2:586756/1‑101 (MQ=255)
                                                                                               aTAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTCACGCCGCGACCGTCTCCGGGGTCGCCa        >  2:5503241/1‑101 (MQ=255)
                                                                                               aTAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTCACGCCGCGACCGTCTCCGGGGTCGCCa        >  2:8711388/1‑101 (MQ=255)
                                                                                               aTAAAGGGGCTTTGCTTTTATCTTGCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTCACGCCGCGACCGTCTCCGGGGTCGCCa        >  2:4382164/1‑101 (MQ=25)
                                                                                                tAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTCACGCCGCGACCGTCTCCGGGGTCGCCAg       >  1:2036226/1‑101 (MQ=255)
                                                                                                tAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTCACGCCGCGACCGTCTCCGGGGTCGCCAg       >  1:9122992/1‑101 (MQ=255)
                                                                                                 aaaGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTCACGCCGCGACCGTCTCCGGGGTCGCCAgc      >  2:7315059/1‑101 (MQ=255)
                                                                                                 aaaGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTCACGCCGCGACCGTCTCCGGGGTCGCCAgc      >  1:2846141/1‑101 (MQ=255)
                                                                                                 aaaGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTCACGCCGCGACCGTCTCCGGGGTCGCCAgc      >  2:1272516/1‑101 (MQ=255)
                                                                                                 aaaGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTCACGCCGCGACCGTCTCCGGGGTCGCCAgc      >  2:6439485/1‑101 (MQ=255)
                                                                                                  aaGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTCACGCCGCGACCGTCTCCGGGGTCGCCAgcg     >  2:1125362/1‑101 (MQ=255)
                                                                                                  aaGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTCACGCCGCGACCGTCTCCGGGGTCGCCAgcg     >  2:2240606/1‑101 (MQ=255)
                                                                                                   aGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTCACGCCGCGACCGTCTCCGGGGTCGCCAgcgc    >  2:7537450/1‑101 (MQ=255)
                                                                                                   aGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTCACGCCGCGACCGTCTCCGGGGTCGCCAgcgc    >  1:7718591/1‑101 (MQ=255)
                                                                                                     gggCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTCACGCCGCGACCGTCTCCGGGGTCGCCAGCGcga  >  2:4160766/1‑101 (MQ=255)
                                                                                                     gggCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCC..CGCCGCCACGGTTCACGCCGCGACCGTCTCCGGGGTCGCCAGCGcga  >  2:5806717/1‑101 (MQ=255)
                                                                                                     |                                                                                                         
GCTCATGGCGGCGACCCTACACGGCGCGACCAGCTTTCGAAAGACAGGAATTCGGCCTGCGGCCGATGCCTCCGGCGGGGATA.TTTGGGGCAAGATAAAGGGGCTTTGCTTTTATCTTTCCCCAAATA.TCCCGGGGGTGC..GGGGGCAGAGCCCCCGCGCCGCCACGGTTCACGCCGCGACCGTCTCCGGGGTCGCCAGCGCGA  >  NC_014034/3439692‑3439894

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 27 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG