Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 243,610 C→T 1.5% intergenic (+39/+81) mogA → / ← ubiG molybdopterin biosynthesis mog protein/3‑demethylubiquinone‑9 3‑O‑methyltransferase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034243,6100C→T 1.5% 2.6 ‑2.9 ‑0.1266intergenic (+39/+81)mogA/ubiGmolybdopterin biosynthesis mog protein/3‑demethylubiquinone‑9 3‑O‑methyltransferase
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (127/135):  total (131/135)
Fisher's exact test strand distribution p-value = 5.75e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.76e-01

Only 100 of 289 total aligned reads displayed.

GGGCCGGGCATATCGACACCGATCCGGCGAAAATCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCG  >  NC_014034/243510‑243709
                                                                                                        |                                                                                                      
gggCCGGGCATATCGACACCGATCCGGCGAAAATCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCAccccc                                                                                                        <  2:1797519/101‑1 (MQ=255)
  gCCGGGCATATCGACACCGATCCGGCGAAAATCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCgg                                                                                                      <  2:3527955/101‑1 (MQ=255)
    cGGGCATATCGACACCGATCCGGCGAAAATCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGATGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGa                                                                                                    >  1:4377664/1‑101 (MQ=255)
      ggCATATCGACACCGATCCGGCGAAAATCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGata                                                                                                  >  1:3135112/1‑101 (MQ=255)
          tatCGACATCGATCCGGCGAAAATCGCGGCGTTTCGACCGAAGACGAAATAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGGTATTTg                                                                                              <  1:2839372/101‑1 (MQ=38)
          tatCGACACCGATCCGGCGAAAATCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTg                                                                                              <  1:3978514/101‑1 (MQ=255)
          tatCGACACCGATCCGGCGAAAATCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTg                                                                                              <  1:1394120/101‑1 (MQ=255)
          tatCGACACCGATCCGGCGAAAATCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTg                                                                                              <  1:454005/101‑1 (MQ=255)
             cGACACCGATCCGGCGAAAATCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAgg                                                                                           <  1:1497070/101‑1 (MQ=255)
             cGACACAGATCCGGCGAAAATCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAgg                                                                                           <  1:865377/101‑1 (MQ=255)
                cacCGATCCGGCGAAAATCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCaa                                                                                        >  1:2774167/1‑101 (MQ=255)
                 acCGATCCGGCGAAAATCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAg                                                                                       <  1:1819610/101‑1 (MQ=255)
                    gATCCGGCGAAAATCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.a                                                                                   >  2:548318/1‑101 (MQ=255)
                    gATCCGGCGAAAATCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.a                                                                                   <  2:2206245/101‑1 (MQ=255)
                     aTCCGGCGAAAATCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.aa                                                                                  <  1:3577018/101‑1 (MQ=255)
                      tCCGGCGAAAATCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.aaa                                                                                 <  1:21805/101‑1 (MQ=255)
                       ccGGCGAAAATCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.g                                                                               >  1:3242163/1‑101 (MQ=255)
                         ggCGAAAATCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGa                                                                             >  2:1479240/1‑101 (MQ=255)
                          gCGAAAATCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAg                                                                            >  1:3895310/1‑101 (MQ=255)
                              aaaTCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTg                                                                        <  1:2528061/101‑1 (MQ=255)
                               aaTCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGt                                                                       <  1:5053053/101‑1 (MQ=255)
                               aaTCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGt                                                                       >  2:4662980/1‑101 (MQ=255)
                               aaTCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCAACCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGt                                                                       >  2:3627691/1‑101 (MQ=255)
                                aTCGCGGTGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGtt                                                                      <  2:268250/101‑1 (MQ=255)
                                aTCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCGCCCCCGGGATAGTTGAGGCAAGAT.AAA.GGAGCCTGtt                                                                      <  2:158720/101‑1 (MQ=255)
                                aTCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGtt                                                                      <  1:1221650/101‑1 (MQ=255)
                                aTCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGtt                                                                      <  2:4973271/101‑1 (MQ=255)
                                 tCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGttt                                                                     <  1:1039544/101‑1 (MQ=255)
                                  cgcgGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCGCTC.TGGC..TTCGCC.ATTCACACCCCGGATATTTGAGGCAAGAT.AAA.GGAGCCTGtttt                                                                    >  1:4322195/1‑101 (MQ=255)
                                  cgcgGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGtttt                                                                    <  1:1057115/101‑1 (MQ=255)
                                    cggcgTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTAt                                                                  <  2:4919958/101‑1 (MQ=255)
                                    cggcgTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTAt                                                                  <  2:2374963/101‑1 (MQ=255)
                                    cggcgTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTAt                                                                  <  2:1132647/101‑1 (MQ=255)
                                    cggcgTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTAt                                                                  <  1:4935869/101‑1 (MQ=255)
                                     ggcgTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATc                                                                 <  2:634534/101‑1 (MQ=255)
                                     ggcgTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATc                                                                 <  2:4270120/101‑1 (MQ=255)
                                      gcgTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCt                                                                <  1:4144229/101‑1 (MQ=255)
                                        gTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTg                                                              >  2:3226546/1‑101 (MQ=255)
                                        gTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTg                                                              >  2:2891589/1‑101 (MQ=255)
                                         tttCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGc                                                             <  1:3417905/101‑1 (MQ=255)
                                            cGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.AGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTc                                                          <  1:3511076/101‑1 (MQ=255)
                                             gCCCGAAGACGAAGTAAAAGGGGGCGCTCCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCa                                                         <  1:3847225/101‑1 (MQ=255)
                                              cccGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCaa                                                        <  2:3206193/101‑1 (MQ=255)
                                               ccGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCaaa                                                       <  1:240511/101‑1 (MQ=255)
                                                cGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.t                                                     <  2:3426508/101‑1 (MQ=255)
                                                 gAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.Tg                                                    >  2:3575222/1‑101 (MQ=255)
                                                  aaGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGt                                                   >  1:2284258/1‑101 (MQ=255)
                                                  aaGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGt                                                   >  1:4893149/1‑101 (MQ=255)
                                                  aaGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGt                                                   >  1:4068812/1‑101 (MQ=255)
                                                    gACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTcc                                                 <  1:1914412/101‑1 (MQ=255)
                                                    gACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTcc                                                 <  2:487527/101‑1 (MQ=255)
                                                      cGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCg                                               <  2:3297187/101‑1 (MQ=255)
                                                       gAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCgg                                              <  1:4630242/101‑1 (MQ=255)
                                                        aaGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCggg                                             >  1:3810552/1‑101 (MQ=255)
                                                         aGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCgggg                                            <  1:306068/101‑1 (MQ=255)
                                                           tAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGc                                          >  2:4353821/1‑101 (MQ=255)
                                                           gggcaGGGGGCGCTGCCCCCTCTTGGCCTTTGGCCAATTCACCCCTGGGATATTTCGGGCAAGATGAAA.TCATGCCCCTTCATCTTGCCGCAAA.TATCCCg                                               >  1:3647499/5‑101 (MQ=0)
                                                           gggcaGGGGGCGCTGCCCCCTCTTGGCCTTTGGCCAATTCACCCCTGGGATATTTCGGGCAAGATGAAA.TCATGCCCCTTCATCTTGCCGCAAA.TATCCCg                                               >  2:2059290/5‑101 (MQ=0)
                                                            aaaaGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGcc                                         <  2:4377347/101‑1 (MQ=255)
                                                            aaaaGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGcc                                         >  2:1407155/1‑101 (MQ=255)
                                                                gggggCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTggg                                     >  1:1192276/1‑101 (MQ=255)
                                                                 ggggCAATGCCCCGTC.TGGC..TTCGCT.ATTCTTCCCGGGGGTATTTGAGGCAGGAT.AAA.GGAGGCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTgggg                                    <  1:419654/101‑1 (MQ=21)
                                                                  gggCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTggggg                                   <  2:3374492/101‑1 (MQ=255)
                                                                  gggCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTggggg                                   >  1:1861488/1‑101 (MQ=255)
                                                                   ggCGCTGCCCCCTCTTGGCCTTTGGCCAATTCACCCCTGGGATATTTCGGGCAAGATGAAA.TCATGCCCCTTCATCTTGCCGCAAA.TATCCCggggggagc                                       >  2:1429211/1‑97 (MQ=2)
                                                                      gcTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTGGGGGCAgc                               >  2:3705955/1‑101 (MQ=255)
                                                                      gcTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTGGGGGCAgc                               >  2:2844657/1‑101 (MQ=255)
                                                                         gCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTGGGGGCAGCGcc                            >  2:2974169/1‑101 (MQ=255)
                                                                          cccccTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTGGGGGCAGCGccc                           >  2:3842010/1‑101 (MQ=255)
                                                                           ccccTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTGGGGGCAGCGcccc                          >  1:907412/1‑101 (MQ=255)
                                                                             ccTC.TGGC..TCCGCC.ATTCCCCCCGGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAT.TGCCCCGGGGGCCTGGGGGCAGCGCCCCCa                        <  2:891065/101‑1 (MQ=38)
                                                                               tc.tGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGc                      <  2:398116/101‑1 (MQ=255)
                                                                                c.tGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGcc                     >  2:2288448/1‑101 (MQ=255)
                                                                                  tGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGGAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCt                    >  2:4181296/1‑101 (MQ=255)
                                                                                   ggC..TTCGCC.ATTCCCCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCtt                   <  1:1570754/101‑1 (MQ=255)
                                                                                   ggC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCtt                   <  2:2435210/101‑1 (MQ=255)
                                                                                   ggC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCtt                   >  2:737650/1‑101 (MQ=255)
                                                                                   ggC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCtt                   >  1:2083339/1‑101 (MQ=255)
                                                                                    gC..TTCGCC.ATTCCCCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTa                  <  2:1261246/101‑1 (MQ=255)
                                                                                    gC..TTCCCC.CTTCCCCCCCGGGGTTTTTGGGGGAAAAA.AAAGGGAGCTTGTTTT.TTTTGCCCCAAAGTGCCCGGGGGGCCTGGGGGCAGCGCCCCCAGCCtt                   <  1:1751980/101‑1 (MQ=255)
                                                                                          cGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGt              >  1:257721/1‑101 (MQ=255)
                                                                                           gCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGtt             >  1:1672559/1‑101 (MQ=255)
                                                                                           gCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGtt             >  1:2449926/1‑101 (MQ=255)
                                                                                             c.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTc           >  2:3786972/1‑101 (MQ=255)
                                                                                             c.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTc           >  2:1651171/1‑101 (MQ=255)
                                                                                             c.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTc           >  2:4588565/1‑101 (MQ=255)
                                                                                               aTTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCa          >  2:4791420/1‑101 (MQ=255)
                                                                                               aTTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCa          >  1:3333716/1‑101 (MQ=255)
                                                                                               aTTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCa          <  1:607476/101‑1 (MQ=255)
                                                                                               aTTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCa          <  1:291215/101‑1 (MQ=255)
                                                                                                ttCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCAc         >  2:2744844/1‑101 (MQ=255)
                                                                                                 tCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGACTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACg        >  1:3372880/1‑101 (MQ=255)
                                                                                                  ccccccgGGGATATTTAGGGAAAGAA.AAG.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGc       <  1:4410877/99‑1 (MQ=25)
                                                                                                  cACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGc       >  1:4262442/1‑101 (MQ=255)
                                                                                                   aCCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGAGTTTCACGCa      >  2:2552207/1‑101 (MQ=255)
                                                                                                     ccccGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAgg    <  2:4068509/101‑1 (MQ=255)
                                                                                                     ccccGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAgg    <  2:3895018/101‑1 (MQ=255)
                                                                                                      cccGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAggc   <  1:5007767/101‑1 (MQ=255)
                                                                                                       ccGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAggcg  >  2:3106548/1‑101 (MQ=255)
                                                                                                       ccGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAggcg  <  1:4016853/101‑1 (MQ=255)
                                                                                                        |                                                                                                      
GGGCCGGGCATATCGACACCGATCCGGCGAAAATCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAA.GGAGCCTGTTTTATCTTGCCTCAAA.TGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCG  >  NC_014034/243510‑243709

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 28 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG