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% 60.4 144intergenic (+29/‑55)betB/betAbetaine‑aldehyde dehydrogenase/choline dehydrogenase
Reads supporting (aligned to +/- strand):  new base (22/100):  ref base (18/0):  total (42/102)
Fisher's exact test strand distribution p-value = 5.32e-12
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.96e-12

Only 100 of 194 total aligned reads displayed.

GCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGCCGGTCCCCCCCCGG.ACCCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCGGGCGG  >  NC_014034/985590‑985785
                                                                                                      |                                                                                               
gCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCGCCGG.Ac                                                                                                 >  2:2513724/1‑101 (MQ=255)
   gcgATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGGCCCCCCCGG.Tcccc                                                                                              >  1:1927414/1‑101 (MQ=255)
    cgATCGTGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.Accccc                                                                                             <  1:949449/101‑1 (MQ=255)
     gATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCG.Acccccc                                                                                            >  2:4249401/1‑101 (MQ=255)
      aTCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.Accccccg                                                                                           <  2:5518169/101‑2 (MQ=255)
      aTCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.Accccccg                                                                                           <  1:5771196/101‑2 (MQ=255)
       tCGAGCACTATTCGCTGTTGACGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACTCAGGCC.GGTGGGTCCCCCCCCGG.Tccccccgg                                                                                          <  2:4813100/101‑3 (MQ=255)
       tCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.Accccccgg                                                                                          <  2:4765500/101‑3 (MQ=255)
       tCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.Accccccgg                                                                                          <  1:5260741/101‑3 (MQ=255)
       tCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.Accccccgg                                                                                          <  1:75089/101‑3 (MQ=255)
       tCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.Accccccgg                                                                                          <  1:2714339/101‑3 (MQ=255)
        cGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.Accccccggc                                                                                         <  2:1351746/101‑4 (MQ=255)
        cGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGGCCCCCCCCG.Accccccccc                                                                                         >  1:4693251/1‑99 (MQ=255)
        cGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCG.GCCCCCCCgc                                                                                         >  1:3357499/1‑100 (MQ=255)
        cGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCC.GCCCCCCCgc                                                                                         >  1:4429772/1‑100 (MQ=255)
         gAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCCG.GCCCCCCCgcg                                                                                        >  2:3667130/1‑99 (MQ=255)
         gAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCG.ACCCCCCCgcg                                                                                        >  1:1832200/1‑99 (MQ=255)
          aGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGt                                                                                      <  1:1670752/101‑1 (MQ=255)
           gCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCGG.CCCCCCCCgcgga                                                                                      >  2:4250265/1‑97 (MQ=255)
            cACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.Accccccggcggat                                                                                     >  2:5792222/1‑94 (MQ=255)
            cACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGTAt                                                                                    <  1:5108114/101‑1 (MQ=255)
            cACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGGCCCCCCCCC.GCCCCCCCgcgcgt                                                                                     >  1:5316002/1‑96 (MQ=255)
            cACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCGG.Accccccggcggat                                                                                     >  2:3805289/1‑94 (MQ=255)
             aCTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTGCCCCCCCGC.GCCCCCCCgcgcatt                                                                                    >  2:2292097/1‑95 (MQ=255)
              cTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGTAttt                                                                                  <  1:5082362/101‑1 (MQ=255)
              cTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGTAttt                                                                                  <  1:832033/101‑1 (MQ=255)
               tataCGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGGGGCCCCCCCCG.CCCCCCCCCGCgcttt                                                                                   >  1:2118867/1‑97 (MQ=255)
               tataCGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGACCATATTGAAACGAAGGCCGGGGGGGGACCCCCCCGGCACCCCCGCgggtttt                                                                                    >  1:4312271/1‑96 (MQ=255)
                ataCGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCCGGGGGGGGCCCCCCCCCG.Ccccccccccgcttttg                                                                                  >  1:3082185/1‑92 (MQ=255)
                   cGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCG.A.CCCCCCGGCGGATTTGCGcc                                                                             >  2:3235768/1‑101 (MQ=255)
                      aGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAg                                                                          <  2:1603706/101‑1 (MQ=255)
                      aGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGGCCCCCCCCG.GCCCCCCC.GCGGATTTGCGCCAAg                                                                          >  1:3604843/1‑101 (MQ=255)
                      aGGTGAAGGCCGTCGATGTCGGCATGGGGCCGGGCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAg                                                                          <  2:1916966/101‑1 (MQ=255)
                        gTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGaa                                                                        <  1:1507264/101‑1 (MQ=255)
                         tGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.GCCCCCCC.GCGGATTTGCGCCCAGaaa                                                                       >  1:892380/1‑101 (MQ=255)
                         ggAAGGCCGTCTATGTCGGCAGGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGaaa                                                                       <  2:1799087/100‑1 (MQ=255)
                           aaGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCG.GCCCCCCC.GCGGATTTTGGCCCACaaaaa                                                                     >  2:2295853/1‑101 (MQ=255)
                           aaGGCCGTCTATGTCGGCATGGGGCCGGTCAAGAGCCCCTATTGAATCGCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGaaaaa                                                                     <  2:4782019/101‑1 (MQ=37)
                           aaGGCCGTCTATGTCGGCATGGGGACGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGaaaaa                                                                     <  1:3713559/101‑1 (MQ=38)
                            aGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGTGGGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGaaaaaa                                                                    <  2:4026840/101‑1 (MQ=255)
                             ggCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAAc                                                                   <  1:1937014/101‑1 (MQ=255)
                              gCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACa                                                                  <  1:1136163/101‑1 (MQ=255)
                               ccGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAg                                                                 <  1:3235968/101‑1 (MQ=255)
                               ccGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGACCCCCCCGG..GACCCCCGCCGTATTTGAGCCAAGAAAAAACaa                                                                 >  2:1306918/1‑100 (MQ=25)
                                 gTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG..GCCCCCCGGCGGATTTGCGGCAAGAAAAAACAggg                                                               >  2:3453396/1‑101 (MQ=255)
                                 gTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCGG.ACCCCCCC.GCGTATTTGCGCCCAGAAAAAAAAggg                                                               >  2:5181725/1‑101 (MQ=255)
                                  tCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGt                                                              <  2:1742172/101‑1 (MQ=255)
                                     aTGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCa                                                           <  1:2604771/101‑1 (MQ=255)
                                     aTGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCG.GCCCCCCCGCGGCTTTTGGCCCAGAAAAAAAAAGGGttg                                                            >  1:2183292/1‑100 (MQ=255)
                                     aTGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCC.GCCCCCCCGGCGCATTTG.GCCCAGaaaaaaaaagcggca                                                           >  1:1378761/1‑92 (MQ=255)
                                     aTGGCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCa                                                           <  2:2871045/101‑1 (MQ=38)
                                      tGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAg                                                          <  2:660681/101‑1 (MQ=255)
                                       gTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAgg                                                         <  1:1312761/101‑1 (MQ=255)
                                       gTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAgg                                                         <  2:719436/101‑1 (MQ=255)
                                        tCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGc                                                        <  1:555580/101‑1 (MQ=255)
                                           gctgggggCCGGTCGAGAGCCCCTATTGAATCGCAGGCC.GGGGGGTCCCCCCCCGG..TCCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCaga                                                     <  1:938127/97‑1 (MQ=38)
                                           gCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCaga                                                     <  1:704830/101‑1 (MQ=255)
                                            cATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCG.ACCCCCCC.GCGGATTTGCGCCCAGAAAAAACAGGGGTCAGGCagag                                                    >  1:2886947/1‑101 (MQ=255)
                                              tGGGGCCGGTCGAGAGCCGCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAggg                                                  <  2:1927306/101‑1 (MQ=38)
                                              tGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGGCCCCCCCCG.GCCCCCCCGCCCGTTTTGGCCCAAAAAAAAAAAGGG.TCAGGGAGAggg                                                  >  2:1988392/1‑101 (MQ=255)
                                                  gCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCGG.A.CCCCCCGGCGGATTTGCGCCAAGAAAAAACAGGGTTAAGGCAGAGGGTggc                                              >  2:3879514/1‑101 (MQ=255)
                                                    cGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCG.CCCCCCCCGCCGATTTTGC.CCCAGAAAAAAAAGGGGTCAAGCAGAGGGAcgcgc                                            >  2:814598/1‑101 (MQ=255)
                                                     ggTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTcgcgcg                                           <  2:1059480/101‑1 (MQ=255)
                                                     ggTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGCATTTCCGCCACGAAAAACAAGGGTGCGGGCGGGggggggcggg                                           >  1:5589741/1‑94 (MQ=255)
                                                     ggTCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCGG.GCCCCCCC.GCGTATTTTTGGCCAGAAAAAAAAAGGGTCAGGCAGAGGGGggcgcg                                           >  1:1615422/1‑101 (MQ=255)
                                                       tCGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCgat                                         <  1:3488737/101‑1 (MQ=255)
                                                       gcGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGGGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCgat                                         <  1:3667358/100‑1 (MQ=38)
                                                       gcGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCgat                                         <  1:2769326/100‑1 (MQ=39)
                                                        cGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCgatg                                        <  2:2118750/101‑1 (MQ=255)
                                                        cGAGAGCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCgatg                                        <  2:4692970/101‑1 (MQ=255)
                                                          agagCCCCTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCgatgga                                      <  2:3166436/101‑1 (MQ=255)
                                                            agCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCGG.ACCCCCCC.GCCAATTTGCGCCAAAAAAAAAAAGGGCACAGGCAGGGGGGCGCGCgttgggag                                    >  1:3722339/1‑98 (MQ=255)
                                                            agCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCG.A.CCCCCCGGCGGATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGGCGCGCgatggatg                                    >  2:1518612/1‑101 (MQ=255)
                                                             gCCCCTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCCG.ACCCCCCC.GCGGATTTGCGCCCAGAAAAAAAAGGGGTCAGGCAGAGGGGCGCGGgatgggtgg                                   >  2:3038040/1‑100 (MQ=255)
                                                               cccTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCgg                                 <  2:5723171/101‑1 (MQ=255)
                                                                ccTATTGAAACGCAGGCC.GGGGGGGCCCCCCCCGG.ACCCCCCC.GCGGCCTTTCGCCCAGAAAAAACAGGGGTCAGGCAGAGGGGCGCGCGCTGGATGCggg                                >  1:2570237/1‑100 (MQ=255)
                                                                 cTATTGAAACGCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGAt                               <  1:4442636/101‑1 (MQ=255)
                                                                   aTTGAACCGC.GGCC.CGGGGGGCCCCCCCCGG..GCCCCCCGGCGTTATTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTAt                            <  2:5589449/101‑1 (MQ=255)
                                                                       aaaCGCAGGCC.GGGGGGGCCCCCCCCGG.ACCCCCCC.GCGGATTTTCGCCAAGAAAAAAAAGGGGTCAGGCAGAGGGGCGGGCGATGGGTGGGGATTAtgtg                         >  2:1249947/1‑101 (MQ=255)
                                                                         aCGCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGAt                       <  1:1734386/101‑1 (MQ=255)
                                                                          cGCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATc                      <  2:3604615/101‑1 (MQ=255)
                                                                           gCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCg                     <  2:4429500/101‑1 (MQ=255)
                                                                           gCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCg                     <  1:5464107/101‑1 (MQ=255)
                                                                           gCAGGCC.GGGGGGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCg                     <  1:4250525/101‑1 (MQ=255)
                                                                              ggcc.ggggGGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTg                  >  1:2458608/1‑101 (MQ=255)
                                                                              ggcc.ggggGGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTg                  <  2:4561970/101‑1 (MQ=255)
                                                                               gcc.gggggGTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTgg                 <  2:4056662/101‑1 (MQ=255)
                                                                                 gccggggtgccccccccCG.GGCCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTgg                 <  1:39366/92‑1 (MQ=38)
                                                                                 ccggggggTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTggg                <  1:4376902/95‑1 (MQ=255)
                                                                                  cggggggTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTgggg               <  2:3716922/96‑1 (MQ=255)
                                                                                  cggggggTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTgggg               <  1:2295979/96‑1 (MQ=255)
                                                                                   ggggggTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCGGAGGGCCGCGCGATGGATGCGGATTATGTGATCGTTGGGGc              >  1:3837291/5‑101 (MQ=38)
                                                                                        gTCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAAAAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGGAGGcc         >  2:3327072/1‑100 (MQ=39)
                                                                                         tCCCCCCCCGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTc        <  1:3307221/101‑1 (MQ=255)
                                                                                         tCCCCCCCCAG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTc        <  1:5366042/101‑1 (MQ=255)
                                                                                         gccccccccGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGCTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTc        <  2:2953406/100‑1 (MQ=255)
                                                                                             cccccGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCgggg    >  1:951266/1‑100 (MQ=255)
                                                                                             cccccGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCGggc    <  2:4337812/101‑1 (MQ=255)
                                                                                              ccccGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCGggcg   >  2:1428884/1‑101 (MQ=255)
                                                                                               cccGG.A.CCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCGggcgg  >  2:1084681/1‑101 (MQ=255)
                                                                                                      |                                                                                               
GCCGCGATCGAGCACTATACGCAGGTGAAGGCCGTCTATGTCGGCATGGGGCCGGTCGAGAGCCCCTATTGAAACGCAGGCC.GGCCGGTCCCCCCCCGG.ACCCCCCCGGCGTATTTGCGCCAAGAAAAAACAGGGTTCAGGCAGAGGGTCGCGCGATGGATGCGGATTATGTGATCGTTGGGGCAGGCTCGGGCGG  >  NC_014034/985590‑985785

Base quality scores:  ATCG  < 3 ≤  ATCG  < 12 ≤  ATCG  < 23 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG