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% 17.7 104intergenic (+29/‑55)betB/betAbetaine‑aldehyde dehydrogenase/choline dehydrogenase
Reads supporting (aligned to +/- strand):  new base (24/71):  ref base (6/1):  total (31/73)
Fisher's exact test strand distribution p-value = 2.43e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.89e-06

Only 100 of 144 total aligned reads displayed.

CCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGCCGGTCCCCCCCCGGACCCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCGGGCGG  >  NC_014034/985591‑985785
                                                                                                    |                                                                                                
ccGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGAcc                                                                                                 <  2:1714632/101‑1 (MQ=255)
ccGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGAcc                                                                                                 <  2:3759324/101‑1 (MQ=255)
  gcgATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGAcccc                                                                                               <  1:3180661/101‑1 (MQ=255)
  gcgATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCGGAcccc                                                                                               >  1:1065622/1‑101 (MQ=255)
   cgATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGTGTCCCCCCCCGGAccccc                                                                                              >  1:150483/1‑101 (MQ=255)
   cgATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGAccccc                                                                                              <  2:4084470/101‑1 (MQ=255)
    gATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGAcccccc                                                                                             <  2:2464632/101‑1 (MQ=255)
    gATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGAcccccc                                                                                             <  1:2769424/101‑1 (MQ=255)
    gATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCGGcccccc                                                                                             >  2:2408924/1‑101 (MQ=255)
     aTCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGAccccccg                                                                                            <  2:1602422/101‑2 (MQ=255)
     aTCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGAccccccg                                                                                            <  1:286531/101‑2 (MQ=255)
     aTCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGAccccccg                                                                                            <  1:2733170/101‑2 (MQ=255)
     aTCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCGCccccccc                                                                                            >  1:2563804/1‑101 (MQ=255)
      tCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGAccccccgg                                                                                           <  2:3728294/101‑3 (MQ=255)
        gAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGAccccccggcg                                                                                         <  1:3046448/101‑5 (MQ=255)
         aGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGt                                                                                      <  2:2559127/101‑1 (MQ=255)
          gCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTa                                                                                     <  1:1142011/101‑1 (MQ=255)
            aCTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTAtt                                                                                   <  2:873763/101‑1 (MQ=255)
            aCTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCCGccccccccccggtt                                                                                     >  1:3930302/1‑94 (MQ=255)
             cTATACGCAGGTGAAGGCCGTCTATGTCGGCATTGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTAttt                                                                                  <  1:2265546/101‑1 (MQ=255)
              tataCGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGTGTCCCCCCCCGGA.CCCCCCG.GCGTATTTg                                                                                 <  2:150459/101‑1 (MQ=255)
              tataCGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGGCCCCCCCCGCCCCCCCCG.GCggtttt                                                                                  >  1:2461147/1‑96 (MQ=255)
              tataCGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCGGA.CCCCCCG.CCGTATTTg                                                                                 >  1:398600/1‑101 (MQ=255)
                taCGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTgcg                                                                               <  1:3361640/101‑1 (MQ=255)
                 aCGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGACCCCCCC..GCGTATTTgcgc                                                                              >  2:591491/1‑101 (MQ=255)
                 aCGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTgcgc                                                                              <  2:3595537/101‑1 (MQ=255)
                   gCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGACCCCCCGG.GGGGATTTggcca                                                                             >  2:1848604/1‑97 (MQ=255)
                   gCAGGGGAAGCCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCa                                                                            <  2:1881932/101‑1 (MQ=255)
                    cAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCaa                                                                           <  1:2274779/101‑1 (MQ=255)
                     tggTGAAGGCTGTCTAAGTCGGCATGGGGCCGGTCGAAAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAg                                                                          <  1:2441080/100‑1 (MQ=255)
                     aGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAg                                                                          <  2:1999101/101‑1 (MQ=255)
                     aGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAg                                                                          <  2:2131632/101‑1 (MQ=255)
                     aGGTGAAGGCCGTCTATCTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAg                                                                          <  2:3127340/101‑1 (MQ=255)
                         gAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGaaaa                                                                      <  2:737603/101‑1 (MQ=255)
                          aaGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGaaaaa                                                                     <  1:2564140/101‑1 (MQ=255)
                           aGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCGA.CCCCCCG.GCGGATTTGCGCCAAGaaaaaa                                                                    >  1:2055596/1‑101 (MQ=255)
                            ggCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAAc                                                                   <  2:3433675/101‑1 (MQ=255)
                             gCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACa                                                                  <  2:646831/101‑1 (MQ=255)
                             gCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACa                                                                  <  2:2563613/101‑1 (MQ=255)
                               cGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAAAAgg                                                                >  2:983925/1‑101 (MQ=37)
                                  cTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGtt                                                             <  2:697815/101‑1 (MQ=255)
                                   tATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGGGGCCCCCCCCCGCCCCCCCC.GCGTTTTTG.GCCAAAAAAAAAAAAGGta                                                             >  1:3849/1‑100 (MQ=255)
                                   tATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCGACCCCCCCG.CCG.ACTTGCGCCCAGAAAAAACAAGGGTc                                                            >  2:1145413/1‑101 (MQ=255)
                                    aTGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCa                                                           <  2:25810/101‑1 (MQ=255)
                                    aTGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCa                                                           <  1:2616132/101‑1 (MQ=255)
                                     tGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAg                                                          <  1:3078362/101‑1 (MQ=255)
                                      gTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCaag                                                         >  2:4191584/1‑99 (MQ=39)
                                      gTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAgg                                                         <  2:395026/101‑1 (MQ=255)
                                       tCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGc                                                        <  1:786654/101‑1 (MQ=255)
                                       tCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGc                                                        <  1:3596746/101‑1 (MQ=255)
                                        cGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCa                                                       <  2:2256516/101‑1 (MQ=255)
                                        cGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCa                                                       <  2:2167026/101‑1 (MQ=255)
                                         ggCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCGAGCCCCCC..GCGGTGTTGTGCCGAGAAAAAACAAGGGTCAAGCag                                                      >  2:3826844/1‑101 (MQ=255)
                                          gCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCaga                                                     <  1:2889875/101‑1 (MQ=255)
                                            aTGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCGA.CCCCCCG.GCGGATTTGCGCCCAGAAAAAACAGGGTTCAGGCAGAgg                                                   >  1:657398/1‑101 (MQ=255)
                                               gggCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAACACAGGGTTCAGGCAGAGGGTc                                                >  2:1209114/1‑101 (MQ=38)
                                                  ccGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGCGCCCCCCG.CCGT.TTTGCCCCCAGACAAAACAGGGGTGAGACAGAggggccgg                                             >  2:4223402/1‑96 (MQ=255)
                                                   cGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCGGACCCCCCC..GCGGATTTCGGCCAAAAAAAAAAAGGGTGCAGGGAGAGGGGggcgc                                            >  1:3495027/1‑101 (MQ=255)
                                                    ggTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTGCCCCCCCGGA.CCCCCCG.GCGGCTTTGCGCCAACAAAAACAAAGGTAAAGGCAGGAGGGcgcgcg                                           >  2:3842886/1‑101 (MQ=255)
                                                     gTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCga                                          <  1:2189271/101‑1 (MQ=255)
                                                     gTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCCGGCCCCCCCGTGCAT.CTTGCCCCACGAAAAAAACGGGTT.AGGCAGAggggggggacg                                          >  2:4184007/1‑93 (MQ=255)
                                                       cGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCgatg                                        <  2:3741614/101‑1 (MQ=255)
                                                         agagCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGGTCAGGCAGAGGGGCGCGCgatgga                                      >  1:2751257/1‑101 (MQ=37)
                                                         agagCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCG.GCCCCCCG.GCGGGTTTGCGCCAAGAAAACAGAGGGTGCAGGGAGAGGGGCGCGCgaaggg                                      >  2:3971968/1‑100 (MQ=255)
                                                          gagCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCgatggat                                     <  1:591555/101‑1 (MQ=255)
                                                          gagCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCgatggat                                     <  1:4151861/101‑1 (MQ=255)
                                                             ccccTATTGAAGCGCAGGCC.GGGGGGGCTCCCCCCGG.CCCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCAATGGATGCg                                  <  2:3296833/101‑1 (MQ=37)
                                                             ccccTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCg                                  <  2:2055458/101‑1 (MQ=255)
                                                             ccccTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCg                                  <  1:266617/101‑1 (MQ=255)
                                                              cccTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCgg                                 <  1:1527722/101‑1 (MQ=255)
                                                              cccTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.GCCCCCCG.GCGTATTTTGGGCAAGAAAAAACAAGGCTCAGGGCGAGGGGCGCGGGACGGATGCgg                                 >  1:3347371/1‑101 (MQ=255)
                                                                cTATTGAAACGCAGGCC.GGGGGGTCCCCCCNCGG.CCCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGNGCGATGGATGCGGAt                               <  1:744529/101‑1 (MQ=38)
                                                                  aTTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTa                             <  1:4221067/101‑1 (MQ=255)
                                                                    tGAAACGCAGGCC.GGGGGGGCCCCCCCCCG.GCCCCCCG.GCGGATTTTCGCCAAGAAAAAAAAGGGGTCAGGCCGGGGGGCGGGCGGTGGGGGCGGGTTAtg                           >  1:1579393/1‑101 (MQ=255)
                                                                     gAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGGCGCGCGATGGATGCGGATTAtgt                          >  2:2622346/1‑101 (MQ=37)
                                                                      aaaCGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTAtgtg                         <  2:1489757/101‑1 (MQ=255)
                                                                       aaCGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGa                        <  2:657321/101‑1 (MQ=255)
                                                                        agccAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGAt                       <  1:2103846/98‑1 (MQ=38)
                                                                         cGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATc                      <  2:963200/101‑1 (MQ=255)
                                                                         cGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATc                      <  1:2365132/101‑1 (MQ=255)
                                                                          gCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCg                     <  1:3141772/101‑1 (MQ=255)
                                                                          gCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCg                     <  1:1209208/101‑1 (MQ=255)
                                                                           cAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTAATTTCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGCCGCGCGATGGATGCGGATTATGTGATCGt                    >  2:450944/1‑101 (MQ=38)
                                                                             ggcc.ggggGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTg                  <  2:2268261/101‑1 (MQ=255)
                                                                                 cgggggggccccccccGGA.CCCCCCG.GCGGATTTGCGCCAAGAAAAAACAGGGGTCAGGCAGAGGGGCGCGCGATGGGTGCGGGTTATGTGATCGTTgggg               >  2:539145/9‑101 (MQ=255)
                                                                                 cggggggTCCCCCCCCCGA.CCCCCCG.GCGTATTTGCGCCCAGAAAAAACAAGGGTCAGGCCGAGGGGCGCGCGGTGGATGCGGATTATGGTTTCGTTtggg               >  2:661276/6‑101 (MQ=16)
                                                                                  ggggggTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGCTTATGGGATCGTTGGGGc              >  2:3437574/5‑101 (MQ=38)
                                                                                  ggggggTCCCCCCCCGGA.CCCCCCG.GCGCATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTTGGGc              >  1:1978034/5‑101 (MQ=38)
                                                                                    ggggTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAg            <  1:984016/99‑1 (MQ=255)
                                                                                     gggTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAgg           <  1:3469143/100‑1 (MQ=255)
                                                                                     gggTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAgg           <  1:904878/100‑1 (MQ=255)
                                                                                      ggTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGc          <  1:852369/101‑1 (MQ=255)
                                                                                       gTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCt         <  1:3077192/101‑1 (MQ=255)
                                                                                        tCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTc        <  2:889700/101‑1 (MQ=255)
                                                                                          cccccccGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTccg      >  2:3365429/1‑99 (MQ=255)
                                                                                           ccccccGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCggg     <  1:1576112/101‑1 (MQ=255)
                                                                                            cccccGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCGggc    <  2:2460963/101‑1 (MQ=255)
                                                                                             ccccGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCGggcg   >  2:1513137/1‑101 (MQ=255)
                                                                                              cccGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCGggcgg  >  2:1031966/1‑101 (MQ=255)
                                                                                              cccGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCGggcgg  >  2:3554767/1‑101 (MQ=255)
                                                                                                    |                                                                                                
CCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGCCGGTCCCCCCCCGGACCCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCGGGCGG  >  NC_014034/985591‑985785

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG