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% 48.2 119intergenic (+29/‑55)betB/betAbetaine‑aldehyde dehydrogenase/choline dehydrogenase
Reads supporting (aligned to +/- strand):  new base (29/72):  ref base (16/0):  total (45/74)
Fisher's exact test strand distribution p-value = 3.21e-08
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.43e-11

Only 100 of 136 total aligned reads displayed.

CGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGCCGGTCCCCCCCCGGACCCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCGGGCGGC  >  NC_014034/985594‑985786
                                                                                                 |                                                                                                 
cgATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCGGAccccc                                                                                               >  2:1958063/1‑101 (MQ=255)
  aTCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGAccccccg                                                                                             <  1:1368220/101‑2 (MQ=255)
  aTCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGAccccccg                                                                                             <  1:2562304/101‑2 (MQ=255)
  aTCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCGGccccccc                                                                                             >  1:384138/1‑101 (MQ=255)
    cGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGGGCCCCTATTGAAACGCAGGCC.GGGGGGGCACCCCCCGGACCCCCCCg.g                                                                                          >  2:2842636/1‑101 (MQ=255)
    cGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCGCCCCCCGGACCCCCCCg.g                                                                                          >  2:2983096/1‑101 (MQ=255)
      aGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGt                                                                                       <  2:1288110/101‑1 (MQ=255)
       gCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTa                                                                                      <  1:3385592/101‑1 (MQ=255)
        cACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTAt                                                                                     <  1:4257107/101‑1 (MQ=255)
        cACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTAt                                                                                     <  1:1058821/101‑1 (MQ=255)
        cACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCGGCGCCCCCgcggct                                                                                       >  2:3313670/1‑96 (MQ=255)
        cACTAAACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTAt                                                                                     <  1:916778/101‑1 (MQ=255)
          cTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCGGCCCCCCCgcgggttt                                                                                     >  1:3578754/1‑94 (MQ=255)
           tataCGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGGGCCCCCCCCCCGCCCCCCCgcggtttt                                                                                     >  2:1610297/1‑94 (MQ=255)
            ataCGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCGGCCCCCCCg.gggcttttg                                                                                  >  1:2902401/1‑93 (MQ=255)
             taCGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTgcg                                                                                <  2:3693707/101‑1 (MQ=255)
             taCGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGT.CCCCCCCGGA.CCCCCCG.GCGTATTTgcgc                                                                               >  2:4027046/1‑101 (MQ=255)
                gCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCa                                                                             <  2:2039701/101‑1 (MQ=255)
                  aGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAg                                                                           <  1:3343275/101‑1 (MQ=255)
                  aGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGGCCCCCCCGG.CCCCCCCG.TCGGCTTTGGGCCAAg                                                                           >  2:560924/1‑101 (MQ=255)
                    gTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGaa                                                                         <  1:2372869/101‑1 (MQ=255)
                    gTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGaa                                                                         <  2:2790149/101‑1 (MQ=255)
                     tGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGaaa                                                                        <  1:3716347/101‑1 (MQ=255)
                     tGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCGA.CCCCCCG.GCGTATTTGCGCCCAGaaa                                                                        >  2:2094068/1‑101 (MQ=255)
                        aGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGaaaaaa                                                                     <  2:2813518/101‑1 (MQ=255)
                         ggCGGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAAc                                                                    <  2:1004248/101‑1 (MQ=38)
                          gCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACa                                                                   <  1:1958192/101‑1 (MQ=255)
                            cGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAgg                                                                 <  1:3004092/101‑1 (MQ=255)
                             gTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAggg                                                                <  1:3540504/101‑1 (MQ=255)
                              tCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCG.CCCCCCCC.GCGGATTTGCGCCCAGaaaaaaaaaggg                                                               >  1:3995430/1‑95 (MQ=255)
                               cTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAggggt                                                              >  2:2741883/1‑99 (MQ=37)
                               cTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGtt                                                              <  1:3313902/101‑1 (MQ=255)
                                tATGTCGGCGTGGGGCCGGTCGAGAGCCCCTCTTGTAACGCGGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTc                                                             <  1:1720676/101‑1 (MQ=25)
                                tATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTc                                                             <  2:403690/101‑1 (MQ=255)
                                tATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTc                                                             <  2:34441/101‑1 (MQ=255)
                                 aTGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCa                                                            <  1:2983317/101‑1 (MQ=255)
                                 aTGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCa                                                            >  2:2126863/1‑101 (MQ=255)
                                  tGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAg                                                           <  1:3344106/101‑1 (MQ=255)
                                   gTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAgg                                                          <  1:697844/101‑1 (MQ=255)
                                    tCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGc                                                         <  1:1727079/101‑1 (MQ=255)
                                     cGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCa                                                        <  1:3481901/101‑1 (MQ=255)
                                     cGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCa                                                        <  1:159440/101‑1 (MQ=255)
                                     cGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCa                                                        <  1:2184280/101‑1 (MQ=255)
                                      ggCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCag                                                       <  1:469571/101‑1 (MQ=255)
                                      ggCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCACCGACACCCCC..GCGGGTTTGCGCACCCAAAAAAAAAGGATCTGGCag                                                       >  2:1778932/1‑101 (MQ=255)
                                      ggCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCGGA.CCCCCCG.GCGGCTTTGCGCCCAGAAAAAACAGGGGTCAGGCag                                                       >  1:2728838/1‑101 (MQ=255)
                                       gCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCaga                                                      <  2:2127411/101‑1 (MQ=255)
                                        cATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCagag                                                     <  2:1094771/101‑1 (MQ=255)
                                        cATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCagag                                                     <  2:41595/101‑1 (MQ=255)
                                        cATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCGGCCCCCCC..GCGGATTTTCGCCCAGAAAAAAAAAGGGTCAGGCagag                                                     >  2:3289667/1‑101 (MQ=255)
                                         aTGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAgg                                                    <  1:2890261/101‑1 (MQ=255)
                                         aTGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCGGA.CCCCCCG.GCGGATTTGCGCCCAGAAAAAACAGGGGTCAGGCAGAgg                                                    >  2:698046/1‑101 (MQ=255)
                                           ggggCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCGACCCCCCCGCGCGT.TTTG.GCCCACAAAAAAAAAGGGTGAGGCCGAgggg                                                  >  2:3226572/1‑100 (MQ=255)
                                           ggggCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCGACCCCCCC..GCGTATTTGCGCCCAGAAAAAACAGGGGTCAGGCAGCgggg                                                  >  2:428562/1‑100 (MQ=255)
                                           ggggCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGAGCCCCCCCCCGCCCCCCCG.GCGGCTTATCGGCAGGAAAAAAAAGGGCTCCggaaagggg                                                   >  1:3212370/1‑94 (MQ=255)
                                            gggCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTc                                                 <  1:986243/101‑1 (MQ=255)
                                              gCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTcgc                                               <  1:2327359/101‑1 (MQ=255)
                                                cGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGGTCAGGCAGAGGGGcgcgc                                             >  2:4193734/1‑101 (MQ=37)
                                                 ggTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGACTCCCCC..GCGCATTTTCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTcgggcc                                            >  2:4399941/1‑100 (MQ=255)
                                                 ggTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTcgcgcg                                            <  1:3169540/101‑1 (MQ=255)
                                                 ggTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCGGCCCCCCC..GCGGATTTGCGCCAAGAAAAAAAAAGGGTCAGGCAGAGGGGcgcgcg                                            >  1:2621818/1‑101 (MQ=255)
                                                  gTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGCTCCCCCCCGG.CCCCCCCC.CCTTGATTGAGCCAAGAAAAAAAAGGTGTCTCGTAGGGGGTggcgcgg                                           >  1:1411019/1‑100 (MQ=255)
                                                   tCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCgat                                          <  2:952938/101‑1 (MQ=255)
                                                    cGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCgatg                                         <  2:384099/101‑1 (MQ=255)
                                                    cGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCgatg                                         <  1:1072619/101‑1 (MQ=255)
                                                    cGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCgatg                                         <  1:374078/101‑1 (MQ=255)
                                                    cGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGCCTTTCCGCCAAGAAAAAACAGGGATCAGGCAGCGGGGCGCGCgaag                                         >  2:2800890/1‑99 (MQ=255)
                                                     gagaGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCgatgg                                        <  2:3804686/101‑1 (MQ=255)
                                                      agagCCCCTATTGAAACGCAGGCC.GGGGGGT.CCCCCCCGG.CCCCCCCG.GCGGATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGGCGCGCgatggag                                      >  1:3603049/1‑100 (MQ=255)
                                                        gatctcTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCgatggatg                                     <  1:1926289/96‑1 (MQ=37)
                                                        agCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGGTCAGGCAGAGGGGCGCGCgatgggtg                                     >  1:575129/1‑101 (MQ=37)
                                                         gCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGc                                    <  2:2212920/101‑1 (MQ=255)
                                                         gCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGc                                    <  1:3050180/101‑1 (MQ=255)
                                                          ccccTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGGCGCGCGATGGATGCg                                   >  1:2950780/1‑101 (MQ=37)
                                                           cccTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCgg                                  <  1:1657175/101‑1 (MQ=255)
                                                           cccTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCCGACCCCCCC..GCGGATTTGCGCCAAGAAAAAACAGGGGTCAGGGAGAGGGGCGCGCGATGGGTGCgg                                  >  2:2378971/1‑101 (MQ=255)
                                                           cccTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCGGCCCCCCC..GCGGATTTTCGCCCAGAAAAAAAAGGGGTCAGGCAGAGGGGCGGGGGATGGGTGCgg                                  >  1:1166729/1‑101 (MQ=255)
                                                            ccTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGa                                 <  1:2094205/101‑1 (MQ=255)
                                                            ccTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGa                                 <  1:2742073/101‑1 (MQ=255)
                                                              tATTGAAACGCAGGCC.GGGGTGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGAtt                               <  2:3175800/101‑1 (MQ=38)
                                                               aTTGAAACGCAGGCC.GGGGGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTa                              <  1:428602/101‑1 (MQ=255)
                                                                          ggcc.ggggGGTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGGGATCGTTg                   >  2:1355478/1‑101 (MQ=38)
                                                                             ccggggggTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTggg                 <  2:1149829/95‑1 (MQ=255)
                                                                              cggggggTCCCCCCCCCGGACCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTggg                 <  2:2148668/96‑1 (MQ=38)
                                                                               ggggggTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGGGGATTATGTGATCGGTGGGGc               >  2:4068593/5‑101 (MQ=38)
                                                                               ggggggTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGc               <  2:2455459/97‑1 (MQ=255)
                                                                               ggggggTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGCTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGc               >  1:1415661/5‑101 (MQ=38)
                                                                                gggggTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCa              >  2:3869841/4‑101 (MQ=255)
                                                                                 ggggTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAg             <  1:759026/99‑1 (MQ=255)
                                                                                 ggggTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAg             >  2:1969990/3‑101 (MQ=255)
                                                                                  gggTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAgg            <  2:1360409/100‑1 (MQ=255)
                                                                                   ggTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGc           >  2:1486475/1‑101 (MQ=255)
                                                                                   gccccccccccGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGc           <  1:2949088/98‑1 (MQ=255)
                                                                                    gTCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCt          <  1:3184162/101‑1 (MQ=255)
                                                                                     tCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTc         >  1:2927396/1‑101 (MQ=255)
                                                                                     tCCCCCCCCGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTc         <  1:1610403/101‑1 (MQ=255)
                                                                                      ccccccccGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCg        >  1:1706364/1‑101 (MQ=255)
                                                                                       cccccccGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCgg       <  1:2127002/101‑1 (MQ=255)
                                                                                           cccGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCGggcgg   >  2:4356639/1‑101 (MQ=255)
                                                                                           cccGGA.CCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCGggcgg   >  1:762297/1‑101 (MQ=255)
                                                                                                 |                                                                                                 
CGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGCCGGTCCCCCCCCGGACCCCCCCG.GCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCGGGCGGC  >  NC_014034/985594‑985786

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