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% 39.7 162intergenic (+29/‑55)betB/betAbetaine‑aldehyde dehydrogenase/choline dehydrogenase
Reads supporting (aligned to +/- strand):  new base (36/110):  ref base (14/0):  total (50/112)
Fisher's exact test strand distribution p-value = 2.04e-08
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.26e-10

Only 100 of 197 total aligned reads displayed.

GCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGCC.GGTCCCCCCCCGGACCCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCGGGCGG  >  NC_014034/985593‑985785
                                                                                                   |                                                                                               
  gATCGAGCCCTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGAcccccc                                                                                            <  1:787082/101‑1 (MQ=255)
  gATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGAcccccc                                                                                            <  2:3385823/101‑1 (MQ=255)
   aTCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCCGCccccccc                                                                                           >  2:3158725/1‑101 (MQ=255)
   aTCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGGCCCCCCCCGGGccccccg                                                                                           >  2:3127732/1‑100 (MQ=255)
   aTCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGGCCCCCCCCCGCgcccccc                                                                                           >  2:1336781/1‑101 (MQ=255)
    tCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGGCCCCCCCCCGCCCCCCCCg                                                                                          >  1:2381828/1‑101 (MQ=255)
       aGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGG.GGGCCCCCCCCGGCccccccccgg                                                                                        >  1:4440886/1‑98 (MQ=255)
       aGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGt                                                                                      <  1:569008/101‑1 (MQ=255)
        gCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTa                                                                                     <  1:1839797/101‑1 (MQ=255)
          aCTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTAtt                                                                                   <  1:3333585/101‑1 (MQ=255)
          aCTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTAta                                                                                   <  1:1183924/101‑2 (MQ=255)
           cTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTAttt                                                                                  <  1:1303879/101‑1 (MQ=255)
           cTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTAttt                                                                                  <  1:884021/101‑1 (MQ=255)
           cTATACGCAGGTGAAGGCCATCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTAttt                                                                                  <  1:2844416/101‑1 (MQ=255)
             gtgcGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGGGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTgc                                                                                <  2:18950/98‑1 (MQ=255)
             ataCGCAGGTGAAGGCCGTCTATGTGGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTgc                                                                                <  1:927550/101‑1 (MQ=255)
             ataCGCAGGTGAAGGCCGTCTATGTCGTCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGGCCCCCCCCCGCCCCCCCGGGCGATTTTg                                                                                 >  1:3108011/1‑101 (MQ=255)
             ataCGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTgc                                                                                <  1:123158/101‑1 (MQ=255)
             ataCGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGGCCCCCCCCCGA.CCCCCCGGCGTAttttg                                                                                >  2:2736830/1‑99 (MQ=255)
                 gCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCa                                                                            <  2:4834885/101‑1 (MQ=255)
                   aGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAg                                                                          <  1:3543675/101‑1 (MQ=255)
                   aGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAg                                                                          <  1:2551694/101‑1 (MQ=255)
                   aGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAg                                                                          <  2:110038/101‑1 (MQ=255)
                   aGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAg                                                                          <  2:1162603/101‑1 (MQ=255)
                     gTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGaa                                                                        <  1:3133819/101‑1 (MQ=255)
                     gTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGaa                                                                        <  2:1532447/101‑1 (MQ=255)
                      tGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGaaa                                                                       <  1:4163443/101‑1 (MQ=255)
                       gAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGaaaa                                                                      <  2:4099211/101‑1 (MQ=255)
                        aaGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGaaaaa                                                                     <  1:4959766/101‑1 (MQ=255)
                          ggCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAc                                                                   <  2:3231764/101‑1 (MQ=255)
                          ggCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAc                                                                   <  2:1450978/101‑1 (MQ=255)
                            ccGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAg                                                                 <  1:4050763/101‑1 (MQ=255)
                            ccGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGGCCCCCCCCGGA.CCCCCCGGCGGATTTGCGCCAAGAAAAAACAg                                                                 >  2:535820/1‑101 (MQ=38)
                              gTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAggg                                                               <  1:1844705/101‑1 (MQ=255)
                               tCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGGCCCCCCCCGGACCCCCCC.GCGGATTTTCCCCAAGAAAAAAAAgggg                                                              >  1:678733/1‑100 (MQ=255)
                                cTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGtt                                                             <  2:1021677/101‑1 (MQ=255)
                                cTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGtt                                                             <  2:1508307/101‑1 (MQ=255)
                                cTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGtt                                                             <  2:159731/101‑1 (MQ=255)
                                 tATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGG.GGGCCCCCCCCCCCCCCCCCCCCCGGTTTTGCCCCAAAAAAAAAAAggga                                                              >  1:1087090/1‑100 (MQ=255)
                                 tATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTc                                                            <  2:4760612/101‑1 (MQ=255)
                                 tATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTc                                                            <  1:4179484/101‑1 (MQ=255)
                                 tATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGGCCCCCCCCGGACCCCCCGGGCG.CTTTGGGGCAAGAAAAAACAGGGTTc                                                            >  1:1482065/1‑101 (MQ=255)
                                 tATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGGCCCCCCCCCGACCCCCCC.GCGGGTTTTCGCCCAGAAAAAACAGGGGTc                                                            >  2:3087374/1‑101 (MQ=255)
                                     tCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACTCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGc                                                        <  1:2204518/101‑1 (MQ=38)
                                       ggCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCCGGACCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCa                                                       <  1:4439608/101‑1 (MQ=38)
                                       ggCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGGCCCCCCCCCGGCCCCCCC.GCGGATTTGCGCCCAGAAAAAAAAAGGGTCAGGCag                                                      >  1:376019/1‑101 (MQ=255)
                                          aTGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAgg                                                   <  1:2251752/101‑1 (MQ=255)
                                              ggCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTcg                                               <  2:3981265/101‑1 (MQ=255)
                                                ccGGTCGAGAGCCCTTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTcgcg                                             <  2:4520477/101‑1 (MQ=37)
                                                ccGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTTCGCCAAGAAAAAACAGGGGTTAGGCCGAGGGGcgcg                                             >  2:2155164/1‑101 (MQ=25)
                                                  ggTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTcgcgcg                                           <  1:716825/101‑1 (MQ=255)
                                                  ggTCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTcgcgcg                                           <  2:678644/101‑1 (MQ=255)
                                                    tCGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCgat                                         <  2:4002441/101‑1 (MQ=255)
                                                     cGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCgatg                                        <  2:2053540/101‑1 (MQ=255)
                                                     cGAGAGCCCCTATTGAAACGCAGGCC.GGGG.GGGCCCCCCCCCGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAAAAGGGGTCAGGCAGAGGGGCGCGCgatg                                        >  2:257431/1‑101 (MQ=25)
                                                      gagaTCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCgatgg                                       <  2:1190955/101‑1 (MQ=37)
                                                      gagaGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCgatgg                                       <  2:4180797/101‑1 (MQ=255)
                                                      gagaGCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCgatgg                                       <  1:4877184/101‑1 (MQ=255)
                                                      gagaGCCCCTATTGAAACGCAGGCC.GGGG.GGGGCCCCCCCCCCGCCCCCCCGCGGGTTTGGGCCAAAAAAAAAAAGGGGTAAGGCGGCgggggcgggggtg                                        >  2:1805769/1‑91 (MQ=255)
                                                       agagCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCgatgga                                      <  1:3944900/101‑1 (MQ=255)
                                                       agagCCCCTATTGAAACGAAGGCC.GGGG.GGCCCCCCCCCGGACCCCCCC.GCGTAATTTGGCCAAGAAAAAACAAGGGTCAGGCAATGGGAcccgcgcgagg                                      >  2:825882/1‑96 (MQ=255)
                                                        gagCCCCTATTGAAACGCAGGCC.GGGG.GGGTGCCCCCCCGCCCCCCCCGCGGTATTT.TGCCCAGAAAAAAAAAGGTTCAGGCCGAGGGGcgcgcggggggt                                     >  1:459006/1‑95 (MQ=255)
                                                        gagCCCCTATTGAAACGCAGGCC.GGGG.GGGCCCCCCCCCGACCCCCCC.GCGTATTTGCGCCCAGAAAAAACAGGGGTCCGGCAGAGGGGCGCGGgatgggt                                     >  2:2983284/1‑99 (MQ=255)
                                                         agCCCCTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCgatggatg                                    <  1:1279622/101‑1 (MQ=255)
                                                         agCCCCTATTGAAACGCAGGCC.GGGG.GGGCCCCCCCCGGGCCCCCCCGG.GGATCTGCGCCCAGAAAAAACAAGGTTTAGGCAGAGGGTCGCGCgatggatg                                    >  2:2859977/1‑101 (MQ=255)
                                                           ccccTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCg                                  <  2:3499434/101‑1 (MQ=255)
                                                             ccTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGa                                <  1:922301/101‑1 (MQ=255)
                                                              cTATTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGACCCCCC..GCGGATTTTCGCCAAAAAAAAAAAAGGGTCAGGCAGAGGGTCGGGCGGGGGAGGCGGGtt                              >  2:1690565/1‑101 (MQ=255)
                                                                aTTGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTa                             <  1:2167912/101‑1 (MQ=255)
                                                                 ttGAAACGCAGGCC.GGGGGGGGCCCCCCCCGCGCCCCCCCGGGGTTTTTG.GCCCAAAAAAAAAAGGGGTCGGGGGAAGGGGGGGGCGGTGGAGGCGGGTTa                             >  2:3529968/1‑101 (MQ=255)
                                                                 ttGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTAt                            <  1:4695798/101‑1 (MQ=255)
                                                                 ttGAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTAt                            <  1:4406685/101‑1 (MQ=255)
                                                                 ttGAAACGCAGGCC.GGGG.GGGCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTAt                            <  1:1052255/101‑1 (MQ=38)
                                                                 ttGAAACGCAGGCC.GGGG.GGGCCCCCCCCCGGCCCCCCC.GCGGATTTTCGGCAAAAAAAAACAGGGGTCAGGCAGAGGGGCGCGCGAGGGGTGCGGGTTAt                            >  1:3036160/1‑101 (MQ=255)
                                                                   gAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTAtgt                          <  2:3052747/101‑1 (MQ=255)
                                                                   gAAACGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGGCGCGCGATGGATGCGGATTAtgt                          >  2:2922174/1‑101 (MQ=37)
                                                                    aaaCGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTAtgtg                         <  2:2542883/101‑1 (MQ=255)
                                                                    aaaCGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTAtgtg                         <  1:3868471/101‑1 (MQ=255)
                                                                     aaCGCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGa                        <  1:1695913/101‑1 (MQ=255)
                                                                      aCGCAGGCT.GGGG.GCTCCCCCCCCGGA.CCCCCCTGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGAt                       <  1:4358992/101‑1 (MQ=38)
                                                                       ctctggcc.ggtg.GGACCCCCCCCGG.TCCCCCCGGCTTATGTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGTATTATGTGATc                      <  2:126297/97‑1 (MQ=255)
                                                                        gCAGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCg                     <  2:1888814/101‑1 (MQ=255)
                                                                          aGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGtt                   >  2:446731/1‑101 (MQ=255)
                                                                          aGGCC.GGGG.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGtt                   <  2:3758769/101‑1 (MQ=255)
                                                                           ggcc.gggg.GGTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTg                  <  2:1048385/101‑1 (MQ=255)
                                                                                tccgatctCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTgggg               <  2:2632723/94‑1 (MQ=38)
                                                                                 ggggggTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGc              >  1:4748658/5‑101 (MQ=255)
                                                                                 ggggggTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGc              >  2:740767/5‑101 (MQ=255)
                                                                                 ggggggTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGGCGCGCGATGGATGCGGATTATGTGATCGTTGGGGc              >  1:196922/5‑101 (MQ=38)
                                                                                  gggggTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCa             >  1:825397/4‑101 (MQ=255)
                                                                                  gggggTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCa             >  1:515782/4‑101 (MQ=255)
                                                                                  gggggTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCa             >  2:2419153/4‑101 (MQ=255)
                                                                                   ggggcccccccccGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAg            <  1:1905759/96‑1 (MQ=255)
                                                                                    gggTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAgg           >  2:3946293/2‑101 (MQ=255)
                                                                                    gggTCCCCCCCCGG.GCCCCCCGGCGTATTTTCGCCAAGAAAAAAAAGGGGTCAGGCAGAGGGTTGCGCGAAGGAAGCGGATTAAGTGATCGTTggggatgg           >  2:4667025/2‑97 (MQ=25)
                                                                                     ggTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGc          >  1:4978931/1‑101 (MQ=255)
                                                                                     ggTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGc          >  2:221639/1‑101 (MQ=255)
                                                                                     ggTCCCCCCCCGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGGCGCGCGAGGGATGCGGATTATGTGATCGTTGGGGGAGGc          >  2:2763754/1‑101 (MQ=38)
                                                                                        ccccccccGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTggggagatcgga       >  1:2632908/1‑93 (MQ=38)
                                                                                        ccccccccGGA.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGATCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCg       >  1:2074261/1‑101 (MQ=255)
                                                                                                   |                                                                                               
GCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGCC.GGTCCCCCCCCGGACCCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCGGGCGG  >  NC_014034/985593‑985785

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 23 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG