Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 985,671 +G 3.1% intergenic (+10/‑74) betB → / → betA betaine‑aldehyde dehydrogenase/choline dehydrogenase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034985,6711.→G 3.1% 2.7 ‑3.3 ‑0.096intergenic (+10/‑74)betB/betAbetaine‑aldehyde dehydrogenase/choline dehydrogenase
Reads supporting (aligned to +/- strand):  new base (3/0):  ref base (30/63):  total (33/63)
Fisher's exact test strand distribution p-value = 3.82e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.89e-01

Only 100 of 122 total aligned reads displayed.

GGGGCGGGAAAATTCCAAGGCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGCCGGTCCCCCCCCGGACCCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTG  >  NC_014034/985571‑985769
                                                                                                     |                                                                                                   
ggggCGGGAAAATTCCAAGGCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAggcc                                                                                                      <  1:4614607/101‑1 (MQ=255)
ggggCGGGAAAATTCCAAGGCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAggcc                                                                                                      <  1:2373274/101‑1 (MQ=255)
 gggCGGGAAAATTCCAAGGCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAggcc.g                                                                                                    <  1:3590109/101‑1 (MQ=255)
 gggCGGGAAAATTCCAAGGCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAggcc.g                                                                                                    <  2:1014717/101‑1 (MQ=255)
   gCGGGAAAATTCCAAGGCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAggcc.ggg                                                                                                  <  1:4254760/101‑2 (MQ=255)
   gCGGGAAAATTCCAAGGCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAggcc.ggg                                                                                                  <  2:2121248/101‑2 (MQ=255)
    cGGGAAAATTCCAAGGCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAggcc.gggg                                                                                                 >  1:2915059/1‑99 (MQ=255)
    cGGGAAAATTCCAAGGCCGCGATCGAGCACTATACGCAGGAGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAggcc.gggg                                                                                                 <  1:2232985/101‑3 (MQ=255)
     gggAAAATTCCAAGGCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAggcc.ggggg                                                                                                <  2:902478/101‑4 (MQ=255)
     gggAAAATTCCAAGGCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAggcc.ggggg                                                                                                >  1:1757201/1‑98 (MQ=255)
      ggAAAATTCCAAGGCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAggcc.gggggg                                                                                               <  1:1242927/101‑5 (MQ=255)
       gAAAATTCCAAGGCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAggcc.ggggggt                                                                                              <  1:3589751/101‑6 (MQ=255)
       gAAAATTCCAAGGCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAggcc.ggggggt                                                                                              <  2:18761/101‑6 (MQ=255)
       gAAAATTCCAAGGCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAggcc.ggggggt                                                                                              >  1:2552630/1‑96 (MQ=255)
       gAAAATTCCAAGGCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAggcc.ggggggg                                                                                              >  1:1987105/1‑96 (MQ=255)
       gAAAATTCCAAGGCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCCATTGAACCGCAggcc.ggggggg                                                                                              >  2:3743997/1‑96 (MQ=255)
        aaaaTTCCAAGGCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTc                                                                                             <  2:3629277/101‑1 (MQ=255)
         aaaTTCCAAGGCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAggcccgggggggg                                                                                             >  1:2372188/1‑92 (MQ=255)
          aaTTCCAAGGCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAggcc.gggggggccc                                                                                           >  1:3358791/1‑93 (MQ=255)
           aTTCCAAGGCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTcccc                                                                                          <  2:4331674/101‑1 (MQ=255)
           aTTCCAAGGCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTcccc                                                                                          <  2:4039394/101‑1 (MQ=255)
             tCCAAGGCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCG.GGGGGGTgccccc                                                                                        >  2:1527292/1‑101 (MQ=255)
              ccAAGGCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTccccccc                                                                                       <  1:2699838/101‑1 (MQ=255)
              ccAAGGCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTccccccc                                                                                       <  1:4417154/101‑1 (MQ=255)
              ccAAGGCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGcgccccc                                                                                       >  1:3449989/1‑101 (MQ=255)
               cAAGGCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTcccccccc                                                                                      <  1:2031830/101‑1 (MQ=255)
                aaGGCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTccccccccc                                                                                     >  2:3028821/1‑100 (MQ=255)
                aaGGCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCg                                                                                     <  1:2143906/101‑1 (MQ=255)
                 aGGCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCgg                                                                                    <  1:987306/101‑1 (MQ=255)
                 aGGCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCgg                                                                                    >  2:283308/1‑101 (MQ=255)
                  ggCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGa                                                                                   <  1:460286/101‑1 (MQ=255)
                   gCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGAc                                                                                  <  2:2886657/101‑1 (MQ=255)
                        gATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGAcccccc                                                                             <  2:1619538/101‑1 (MQ=255)
                        gATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGAcccccc                                                                             <  1:4381320/101‑1 (MQ=255)
                         aTCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGAccccccg                                                                            <  2:3580353/101‑2 (MQ=255)
                         aTCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCGGCccccccg                                                                            >  2:511316/1‑100 (MQ=255)
                          tCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGAccccccgg                                                                           <  2:813315/101‑3 (MQ=255)
                           caaGCACTATACGCAGGTGAAGGCCGCCTACGTCGGCATGGGCCCGGTCGAGTGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.CCCCCCCGGc                                                                         <  2:288696/99‑1 (MQ=255)
                           cGAGCAGTATACGCAGGTGAAGGCCGTCTAGGCCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGAccccccggc                                                                          <  2:2582465/101‑4 (MQ=255)
                           cGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGAccccccggc                                                                          <  2:3922479/101‑4 (MQ=255)
                           cGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCGAccccccggc                                                                          >  1:669769/1‑98 (MQ=255)
                            gAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGAccccccggcg                                                                         <  2:2326559/101‑5 (MQ=255)
                             aGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCGG.GCCCCCCGGCGt                                                                       >  1:4157615/1‑101 (MQ=255)
                             aGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCCGGCCCCCCgggg                                                                        >  1:2303203/1‑99 (MQ=255)
                                 cTATACGCAGGTGAAGGCCGTGTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCGGCCCCCCCGCCgactt                                                                    >  1:3433068/1‑97 (MQ=255)
                                 cTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTAttt                                                                   <  2:2372026/101‑1 (MQ=255)
                                    taCGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTgcg                                                                <  1:1224806/101‑1 (MQ=255)
                                      cGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGcc                                                              <  1:651505/101‑1 (MQ=255)
                                       gCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCa                                                             <  2:3449735/101‑1 (MQ=255)
                                        cAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCGCCCCCCGGA.CCCCCCGGCGTATTTGCGCCaa                                                            >  2:603143/1‑101 (MQ=38)
                                         aGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAg                                                           <  1:1624914/101‑1 (MQ=255)
                                          ggTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGa                                                          <  1:1258204/101‑1 (MQ=255)
                                          ggTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGa                                                          <  2:734370/101‑1 (MQ=255)
                                            tGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGaaa                                                        <  2:2626279/101‑1 (MQ=255)
                                              aaGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGaaaaa                                                      <  2:1520331/101‑1 (MQ=255)
                                               aGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGaaaaaa                                                     <  2:1165628/101‑1 (MQ=255)
                                                ggCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.CCCCCCCGGGGTATTTGCGCCCCGAAAAAAc                                                    >  2:1342317/1‑101 (MQ=255)
                                                    gTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.Aggg                                               >  1:600479/1‑101 (MQ=255)
                                                    gTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCCGACCCCCCC.GCGGATTTGCGCCCAGAAAAAAC.Aggg                                               >  1:2194666/1‑101 (MQ=255)
                                                      cTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGtt                                             <  2:3079937/101‑1 (MQ=255)
                                                       tATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTc                                            <  1:2249841/101‑1 (MQ=255)
                                                       tATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTc                                            <  1:238193/101‑1 (MQ=255)
                                                        aTGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCa                                           <  2:3103621/101‑1 (MQ=255)
                                                         tGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAg                                          <  2:4060997/101‑1 (MQ=255)
                                                         tGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAg                                          <  1:3715996/101‑1 (MQ=255)
                                                         tGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGGCCCCCCCCG.CCCCCCCCCCGGATTTTCCCCCAGAAAAAAA.AAGGGTCAg                                          >  1:407133/1‑101 (MQ=255)
                                                         tGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCGGCCCCCCC.GCGGATTTTTGCCCAAAAAAAACAAGGGTTcg                                           >  2:3724915/1‑100 (MQ=255)
                                                          gTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGGGCCCCCCCCCCCCCCCCCCCCGGGTTTTGCCCCAAAAAAAAAA.Aggggcgg                                           >  2:920036/1‑97 (MQ=255)
                                                          gTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCGGACCCCCCC.GCGGATTTTCGCCAAGAAAAAAA.Agggggccgg                                         >  1:3349066/1‑95 (MQ=255)
                                                            cGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCa                                       <  1:1988949/101‑1 (MQ=255)
                                                             ggCATTGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCag                                      <  2:1366010/101‑1 (MQ=37)
                                                                aTGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAgg                                   <  2:1757082/101‑1 (MQ=255)
                                                                aTGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCGGCCCCCCC.GCGGATTTGGGCCCAAAAAAAAA.AGGGGTCAGGCGGAgg                                   >  1:3664188/1‑101 (MQ=255)
                                                                 tGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGGGTATTTGGGCCAAGAAAAAAA.AGGGGTCCGGCAGAggg                                  >  2:2810594/1‑101 (MQ=25)
                                                                  ggggCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGGGCCCCCCCCGGCCCCCCCCGG.G.ATTTTCGCCCAGAAAAAAA.AAGGGTCAGGGCGAgggg                                 >  1:1096846/1‑100 (MQ=255)
                                                                    ggCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTcg                               <  1:1342412/101‑1 (MQ=255)
                                                                     gCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCCAGAAAAAAC.AGGGTTCAGGCAGAGGGGcgc                              >  1:2783558/1‑101 (MQ=37)
                                                                     gCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCCGA.CCCCCCGGCGTATTTGCGCCCAGAAAAAAC.AGGGTTCAGGCAGAGGGGcgc                              >  2:3038723/1‑101 (MQ=37)
                                                                      ccGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTcgcg                             <  1:4051929/101‑1 (MQ=255)
                                                                       cGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTcgcgc                            <  2:2476182/101‑1 (MQ=255)
                                                                        ggTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTcgcgcg                           <  2:3358542/101‑1 (MQ=255)
                                                                         gTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGGcggggga                          >  2:597675/1‑96 (MQ=38)
                                                                          tCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCgat                         <  2:211928/101‑1 (MQ=255)
                                                                           cGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCgatg                        <  1:1090858/101‑1 (MQ=255)
                                                                            gagaGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCgatgg                       <  2:4495059/101‑1 (MQ=255)
                                                                            gagaGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCgatgg                       >  2:2952846/1‑101 (MQ=39)
                                                                             agagCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCgatggg                      >  2:1544431/1‑100 (MQ=37)
                                                                              gagCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCCGA.CCCCCCGGCGAATTTTGGCCAAGAAAAAAC.AGGGGTTAGGCCGGGGGGcgggcggggggt                     >  2:4562779/1‑95 (MQ=255)
                                                                               agCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCGGACCCCCCGGCCGT.TTTGCCCCCAAAAAAAAC.AGGGTTCAAGGACAGGGTCCCGGgagggatg                    >  1:2088681/1‑101 (MQ=255)
                                                                                gCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGc                   <  1:4562506/101‑1 (MQ=255)
                                                                                    cTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGAt               <  1:283340/101‑1 (MQ=255)
                                                                                     tATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGAtt              <  1:1577123/101‑1 (MQ=255)
                                                                                         gAAACGCGGGCC.GG.GGGGCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTAtgtg         <  2:3427923/101‑1 (MQ=37)
                                                                                           aaCGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGa        <  1:76440/101‑1 (MQ=255)
                                                                                            aCGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGAt       <  1:2736031/101‑1 (MQ=255)
                                                                                             cGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATc      <  1:1605852/101‑1 (MQ=255)
                                                                                                aGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGtt   <  2:3713897/101‑1 (MQ=255)
                                                                                                aGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGtt   <  2:3669124/101‑1 (MQ=255)
                                                                                                aGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGGCGCGCGATGGATGCGGATTATGGGAGCGtt   >  1:626465/1‑101 (MQ=25)
                                                                                                 ggcc.ggggGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTg  >  1:3836414/1‑101 (MQ=255)
                                                                                                     |                                                                                                   
GGGGCGGGAAAATTCCAAGGCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGCCGGTCCCCCCCCGGACCCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTG  >  NC_014034/985571‑985769

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