Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 985,690 Δ1 bp 100% intergenic (+29/‑55) betB → / → betA betaine‑aldehyde dehydrogenase/choline dehydrogenase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034985,6900C→.100.0% 21.4 96intergenic (+29/‑55)betB/betAbetaine‑aldehyde dehydrogenase/choline dehydrogenase
Reads supporting (aligned to +/- strand):  new base (33/54):  ref base (6/2):  total (40/56)
Fisher's exact test strand distribution p-value = 6.06e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.61e-06

Only 100 of 115 total aligned reads displayed.

GCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGCCGGTCCCCCCCCGGACCCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCGGGCGGC  >  NC_014034/985590‑985786
                                                                                                     |                                                                                                 
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)
        cGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGAccccccggc                                                                                           <  2:3922479/101‑4 (MQ=255)
        cGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCGAccccccggc                                                                                           >  1:669769/1‑98 (MQ=255)
         gAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGAccccccggcg                                                                                          <  2:2326559/101‑5 (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)
                      aGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAg                                                                            <  1:575436/101‑1 (MQ=255)
                      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)
                        gTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGaa                                                                          <  1:2303373/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)
                                  tCTATGTCGGGATGGGGCCGGTCGAGAGCCCCTATTTAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGt                                                               <  2:3626174/101‑1 (MQ=38)
                                   cTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGtt                                                              <  2:3079937/101‑1 (MQ=255)
                                    tATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTc                                                             <  1:238193/101‑1 (MQ=255)
                                    tATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTc                                                             <  1:2249841/101‑1 (MQ=255)
                                     aTGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCa                                                            <  2:2172046/101‑1 (MQ=255)
                                     aTGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCa                                                            <  2:3103621/101‑1 (MQ=255)
                                     aTGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCa                                                            <  2:6538/101‑1 (MQ=255)
                                      tGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAg                                                           <  2:4060997/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.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAgg                                                          <  2:166714/101‑1 (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)
                                                   ccGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGGGTTTGGGCCCAGAAAAAAA.AGGGGTCAGGCAGGGGGAcgcg                                              >  2:2899379/1‑101 (MQ=255)
                                                    cGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTcgcgc                                             <  2:2476182/101‑1 (MQ=255)
                                                     ggTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTcgcgcg                                            <  2:3358542/101‑1 (MQ=255)
                                                     ggTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCGA.CCCCCCGGCGTATTTGCGCCCAGAAAAAAC.AGGGGTAAGGCAGAGGGGcgcgcg                                            >  1:1291336/1‑101 (MQ=25)
                                                      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.GGGGGGGCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGGTCAGGCAGAGGGGCGCGCgatgg                                        >  1:3468247/1‑101 (MQ=37)
                                                          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)
                                                                ccTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGa                                 <  2:4451833/101‑1 (MQ=255)
                                                                  tATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGAtt                               <  1:1577123/101‑1 (MQ=255)
                                                                    ttGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTAt                             <  1:3208303/101‑1 (MQ=255)
                                                                     tGAAACGCAGGCC.GGGGGGGCCCCCCCCGGA.CCCCCCGGCGGATTTGCGCCAAGAAAAAAC.AGGGGTCAGGCAGAGGGGCGCGCGATGGGTGGGGGTTAtg                            >  1:725949/1‑101 (MQ=255)
                                                                      gAAACGCGGGCC.GG.GGGGCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTAtgtg                          <  2:3427923/101‑1 (MQ=37)
                                                                         aCGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGAt                        <  1:2736031/101‑1 (MQ=255)
                                                                          cGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATc                       <  1:1605852/101‑1 (MQ=255)
                                                                           ccAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCg                      <  1:259311/100‑1 (MQ=37)
                                                                            cAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGt                     <  1:4019853/101‑1 (MQ=255)
                                                                             aGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGtt                    <  2:3669124/101‑1 (MQ=255)
                                                                             aGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGtt                    <  2:3713897/101‑1 (MQ=255)
                                                                             aGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGGCGCGCGATGGATGCGGATTATGGGAGCGtt                    >  1:626465/1‑101 (MQ=25)
                                                                              ggcc.ggggGGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTg                   >  1:3836414/1‑101 (MQ=255)
                                                                                  cggggggcccccccccGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTgggg                <  1:4170206/93‑1 (MQ=38)
                                                                                   ggggggTCCCCCCCCGGA.CCCCCCGGGGTATTTGCGCCAAGAAAAAAC.AGGGGTCAGGCAGAGGGGCGCGCGGTGGGTGCGGGTTATGTGATCGTTGGGGc               >  2:1024646/5‑101 (MQ=25)
                                                                                   ggggggTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGc               >  1:4597083/5‑101 (MQ=255)
                                                                                   ggggggTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGc               >  1:2677547/5‑101 (MQ=255)
                                                                                   ggggggTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGGTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTggggg               >  1:4194866/5‑100 (MQ=38)
                                                                                    gggggTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCa              >  1:635126/4‑101 (MQ=255)
                                                                                    gggggTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCa              <  2:2803443/98‑1 (MQ=255)
                                                                                       ggTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGCCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGc           >  2:4564191/1‑101 (MQ=255)
                                                                                        gTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCGGAGGGCCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAgggt          >  1:3874417/1‑99 (MQ=39)
                                                                                        gTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCt          <  2:1234904/101‑1 (MQ=255)
                                                                                         tCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCCATCAAAAAC.AGGGTTCCAGCAGCGGGTCGCGGGATGGAAGCGGGTTGTGTGATCGTTGGGGCAGGCTc         >  2:526287/1‑101 (MQ=255)
                                                                                         tCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTc         <  1:920116/101‑1 (MQ=255)
                                                                                          ccccccccGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGGTCg        >  2:1030757/1‑101 (MQ=255)
                                                                                          ccccccccGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCg        <  1:3372485/101‑1 (MQ=255)
                                                                                          ccccccccGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCg        >  1:4620204/1‑101 (MQ=255)
                                                                                          ccccccccCGGCCCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTc         <  1:2911155/101‑1 (MQ=255)
                                                                                           cccccccGGA.CCCCCCGGCGTATTTTCGCCCAGAAAAAAA.AGGGGTCAGGCAGAGGGGCGCGCGATGGATGCGGGTTATGTGATCGTTGGGGGAGGCTCgg       >  2:710700/1‑101 (MQ=255)
                                                                                           cccccccGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTccg       >  2:2668161/1‑99 (MQ=255)
                                                                                           cccccccGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCgg       >  2:3220647/1‑101 (MQ=255)
                                                                                             cccccGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCgggg     >  1:4024007/1‑100 (MQ=255)
                                                                                             cccccGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCgggg     >  1:1204418/1‑100 (MQ=255)
                                                                                             cccccGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCGggc     <  1:3744272/101‑1 (MQ=255)
                                                                                             cccccGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCGggc     <  2:1607710/101‑1 (MQ=255)
                                                                                               cccGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCGggcgg   >  2:1276587/1‑101 (MQ=255)
                                                                                               cccGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCATTGGGGCAGGCTCGggcgg   >  2:716365/1‑101 (MQ=255)
                                                                                                ccGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCAggcggc  >  2:1946411/1‑101 (MQ=255)
                                                                                                     |                                                                                                 
GCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGCCGGTCCCCCCCCGGACCCCCCCGGCGTATTTGCGCCAAGAAAAAAC.AGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCGGGCGGC  >  NC_014034/985590‑985786

Base quality scores:  ATCG  < 3 ≤  ATCG  < 12 ≤  ATCG  < 23 ≤  ATCG  < 32 ≤  ATCG  < 40 ≤  ATCG