Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,639,586 +G 2.9% coding (1271/1605 nt) mccB → methylcrotonoyl‑CoA carboxylase subunit beta

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,639,5861.→G 2.9% 7.7 ‑4.5 ‑0.0173coding (1271/1605 nt)mccBmethylcrotonoyl‑CoA carboxylase subunit beta
Reads supporting (aligned to +/- strand):  new base (5/0):  ref base (67/101):  total (72/101)
Fisher's exact test strand distribution p-value = 1.15e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.95e-01

Only 100 of 222 total aligned reads displayed.

CGGCATGGGGCGAAGATGGTGACGGCCGTCGCCACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTCCCGCATCTCGGTGATGGGCGGCGAACAGGCGGCGGGGGTGCTTGCCACGGTCAAG  >  NC_014034/1639486‑1639686
                                                                                                     |                                                                                                      
cGGCATGGGGCGAAGATGGTGACGGCCGTCGCCACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGc                                                                                                         <  1:3677365/101‑1 (MQ=255)
 ggCATGGGGCGAAGATGGTGACGGCCGTCGCCACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGc.g                                                                                                       <  2:1378829/101‑1 (MQ=255)
   cATGGGGCGAAGATGGTGACGGCCGTCGCCACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGc.gg.g                                                                                                    >  2:2519164/1‑101 (MQ=255)
    aTGGGGCGAAGATGGTGACGGCCGTCGCCACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTAAGGCATGGc.gg.gg                                                                                                   >  2:2714867/1‑100 (MQ=255)
     tGGGGCGAAGATGGTGACGGCCGTCGCCACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGc.gg.ggg                                                                                                  >  1:980882/1‑99 (MQ=255)
     tGGGGCGAAGATGGTGACGGCCGTCGCCACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGc.gg.gcc                                                                                                  <  1:2266085/101‑1 (MQ=255)
      ggggCGAAGATGGTGACGGCCGTCGCCACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGc.gg.ggcc                                                                                                 >  1:3377837/1‑98 (MQ=255)
      ggggCGAAGATGGTGACGGCCGTCGCCACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGc.gg.gccg                                                                                                 <  1:4451053/101‑1 (MQ=255)
       gggCGAAGATGGTGACGGCCGTCGCCACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGc.gg.ggcgg                                                                                                >  2:2868554/1‑101 (MQ=255)
        ggCGAAGATGGTGACGGCCGTCGCCACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGc.gg.gccggg                                                                                               <  2:1139844/101‑1 (MQ=255)
         gCGAAGATGGTGACGGCCGTCGCCACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGc.gg.ggagcgc                                                                                              >  2:3961269/1‑95 (MQ=255)
          cGAAGATGGTGACGGCCGTCGCCACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGc.gg.ggcggggc                                                                                             >  1:1752502/1‑99 (MQ=255)
            aaGATGGTGACGGCCGTCGCCACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGCGgg.gccgggcca                                                                                            >  1:3203273/1‑100 (MQ=255)
            aaGATGGTGACGGCCGTCGCCACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GGCGGGCCtc                                                                                           >  2:550488/1‑100 (MQ=255)
             aGATGGTGACGGCCGTCGCCACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTAt                                                                                          <  2:1906073/101‑1 (MQ=255)
               aTGGTGACGGCCGTCGCCACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGCGgg.gccgggcccttt                                                                                         >  1:120820/1‑97 (MQ=255)
               aTGGTGACGGCCGTCGCCACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.gg.ggccggccccttc                                                                                        >  2:4739335/1‑96 (MQ=255)
               aTGGTGACGGCCGTCGCCACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTc                                                                                        <  1:4122368/101‑1 (MQ=255)
               aTGGTGACGGCCGTCGCCACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTc                                                                                        <  1:2522515/101‑1 (MQ=255)
                tGGGGACGGCCGGCGCCACGCCGAACGTGCCGAAGATCCCGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCg                                                                                       <  1:1385815/101‑1 (MQ=255)
                 ggTGACGGCCGTCGCCACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGc                                                                                      <  1:609816/101‑1 (MQ=255)
                  gTGACGGCCGTCGCCACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTAtctgcc                                                                                     <  2:4019498/101‑6 (MQ=255)
                  gTGACGGCCGTCGCCACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGcc                                                                                     <  1:216705/101‑1 (MQ=255)
                     aCGGCTGTCGCCACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCg                                                                                  <  2:1041485/101‑1 (MQ=255)
                     aCGGCCGTCGCCACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGTCCCg                                                                                  <  2:1969872/101‑1 (MQ=255)
                       ggCCGTCGCCACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGt                                                                                <  2:2959038/101‑1 (MQ=255)
                        gCCGTCGCCACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GGCGGCGGGCCAATT.GCCCCgctc                                                                               >  2:2158045/1‑98 (MQ=255)
                          cGTCGCCACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGCGGG.GCGGGGCCCTTTCCCCCCGGTTc                                                                              >  2:3774809/1‑101 (MQ=255)
                          cGTCGCCACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGCGGG.GCCGGGCCCTTTGCCCCCGGtcc                                                                              >  1:3305931/1‑99 (MQ=255)
                           gTCGCCACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GGCGGGGCCCTTTGccccccggccc                                                                            >  1:3253395/1‑94 (MQ=255)
                            tCGCCACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCtg                                                                           <  1:2643165/101‑1 (MQ=255)
                             cGCCACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GGCGGGGCTATTCGCCCCCGGTCCcgt                                                                          >  2:3543145/1‑101 (MQ=255)
                             cGCCACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GGCGGGGCCATTTGCCCCCGGTCCtgt                                                                          >  1:1444465/1‑101 (MQ=255)
                               ccACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GGCGGGCCCATTCTCCCCCCGTCCTGTgg                                                                        >  2:4369387/1‑101 (MQ=255)
                               ccACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTgg                                                                        <  2:1617044/101‑1 (MQ=255)
                                 acgacgAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGAc                                                                      <  2:3359286/101‑1 (MQ=255)
                                   gacgaACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCt                                                                    <  1:4293813/101‑1 (MQ=255)
                                   gacgaACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCt                                                                    <  1:288427/101‑1 (MQ=255)
                                   gacgaACGTGCCGAAGATCACGCTGCTGGGGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGCGGG.GCGGGGCCTTTTCCCCCCGGTCCCTGGGGCcc                                                                     >  2:2246509/1‑101 (MQ=255)
                                    acgaACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.g                                                                  <  1:1877018/101‑1 (MQ=255)
                                      gaACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGGCTATTTGCCCCCGGTCCTGTGGGCCC.GGc                                                                >  2:2496549/1‑101 (MQ=255)
                                        aCGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGccc                                                              <  2:278916/101‑1 (MQ=255)
                                            gCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACt                                                          <  2:980805/101‑1 (MQ=255)
                                            gCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACt                                                          <  1:2756719/101‑1 (MQ=255)
                                             ccGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTc                                                         <  2:3172141/101‑1 (MQ=255)
                                             ccGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTc                                                         <  1:502736/101‑1 (MQ=255)
                                              cGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGGCTATTCTCCCCCGGTCCTGTGGGCCC.GGGCCCACTcc                                                        >  2:200060/1‑101 (MQ=255)
                                              cGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGGCCCTTCGCCCCCGCTCCTGTGGACCTCGGccccctcc                                                         >  2:1734432/1‑96 (MQ=255)
                                              cGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTcc                                                        <  2:901292/101‑1 (MQ=255)
                                              cGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTcc                                                        <  1:242006/101‑1 (MQ=255)
                                               gAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCCGGGCTATTTCCCCCCGGTCCCCTGGGCCC.GGCCCCACTccc                                                       >  1:4170801/1‑101 (MQ=255)
                                                   aTCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTCCCGCAt                                                   <  1:924356/101‑1 (MQ=255)
                                                   aTCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTCCCGCAt                                                   <  2:1011248/101‑1 (MQ=255)
                                                   aTCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTCCCGCAt                                                   <  2:716327/101‑1 (MQ=255)
                                                      aCGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTCCCGCAtctc                                                <  2:1851098/101‑1 (MQ=255)
                                                       cGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTCCCGCATCTCg                                               <  1:2509967/101‑1 (MQ=255)
                                                         ctgctgGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTCCCGCATCTCGGt                                             <  1:1801615/101‑1 (MQ=255)
                                                          tgctgGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTTCCCCCCGGTCCTGTGGGCCC.GGCCCCACTCCCCCCTCTCgggg                                            >  1:1644992/1‑99 (MQ=255)
                                                          tgctgGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTCCCGCATCTCGGTg                                            <  2:3636691/101‑1 (MQ=255)
                                                          tgctgGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTCCCGCATCTCGGTg                                            <  2:4217463/101‑1 (MQ=255)
                                                           gctgGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GGCGGGGCTATTCCCCCCGGGTTCTGTGGACCT.GGGCCCACTCCCCCATCTCggcgg                                           >  2:1302496/1‑98 (MQ=255)
                                                            cgatctggcggcTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTCCCGCATCTCGGTGAt                                          <  1:792487/95‑1 (MQ=255)
                                                             tggtggGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTCCCGCATCTCGGTGATg                                         <  2:3106418/101‑1 (MQ=255)
                                                             tggtggGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCCGGTCCTGTGGGCCT.GGCCCAACTCCCCCCTCTCGGGGATg                                         >  2:3486284/1‑101 (MQ=255)
                                                              ggtggGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GGCGGGGCTATTTGCCCCCGTTCCTGTGGACCC.CGCCCAACTCCCGCCTCTCCGTGATgg                                        >  2:4641455/1‑101 (MQ=255)
                                                              ggtggGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GGCGGGGCTATTCGCCCCCGGTCCTGTGGGCCT.GGCCCCACTCCCCCATCTCGGGGATgg                                        >  1:3996846/1‑101 (MQ=255)
                                                              ggtggGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTCCCGCATCTCGGTGATgg                                        <  2:624384/101‑1 (MQ=255)
                                                                tggGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCCGGGCTATTCGCCCCGGGTCCTTGGGACCC.GGCCCCACTCCCCCACCTCCGGGAGgggg                                      >  2:4521766/1‑100 (MQ=255)
                                                                 gggCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGGCCATTTGCCCCCGCTCCTGTGGGCCC.GGGCCCACTCCCCCATTTTGGTGAAGggcg                                     >  2:1096639/1‑101 (MQ=255)
                                                                 gggCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCCCCCCCGGTCCCGTGGACCT.GGCCCCACTCCCCCCTCCCCGGGATggggg                                     >  1:2370304/1‑99 (MQ=255)
                                                                   gcggcTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTCCCGCATCTCGGTGATGggcggc                                   <  1:302436/101‑1 (MQ=255)
                                                                   gcggaTCGTTCGGGGCGGGCAATTTCGGCATGGC.GG.GCCGGGCCTATTCGCCCCTGTTCCTGTGGACCT.GGCCCAACTCCCGCATCTCGGTGATGggcggc                                   <  2:2354394/101‑1 (MQ=255)
                                                                       cTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTCCCGCATCTCGGTGATGGGCGGCGAAc                               <  1:2821046/101‑1 (MQ=255)
                                                                       cTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTCCCGCATCTCGGTGATGGGCGGCGAAc                               <  2:34688/101‑1 (MQ=255)
                                                                        tCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GGCGGGGCTATTTTCCCCGGGTCCTGTGGGCCT.GGCCCCACTCCCGCACCTCGGGGATGggggggggaaa                              >  2:1360730/1‑95 (MQ=255)
                                                                        tCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTCCCGCATCTCGGTGATGGGCGGCGAACa                              <  2:4279536/101‑1 (MQ=255)
                                                                         cGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCCGGCCTATTTGCCCCCGGTTCTGTGGACCC.GGGCCAACCCCCGCAACTCGGCGATGGGGGGGGGACAg                             >  1:4648182/1‑101 (MQ=255)
                                                                             cGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGGCTATTCGCCCCCGTTCCTGTGGACCT.GGCCCAACTCCCGCCTCTCGGGGATGGGGGGCGGACAggcgg                         >  1:1026705/1‑101 (MQ=255)
                                                                             cGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGGCTATTCGCCCCCGTTCCTGTGGACCT.GGCCCAACTCCCGCATCTCGGTGATGGGCGGCGAACAggcgg                         >  2:3730548/1‑101 (MQ=255)
                                                                             cGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTCCCGCATCTCGGTGATGGGCGGCGAACAggcgg                         <  1:3775044/101‑1 (MQ=255)
                                                                              ggggCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTCCCGCATCTCGGTGATGGGCGGCGAACAggcggc                        <  1:3868563/101‑1 (MQ=255)
                                                                              ggggCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCCGTTCCTGTGGACCT.CGCCCAACTCCCCCCTCTCGGTGATGGGCGGGGAACAagcggg                        >  2:1982827/1‑100 (MQ=255)
                                                                                ggcgggcAATTACGGCATGGC.GG.GCCGGGGCTATTTGCCCCCGGTCCTGTGGGCCT.GGCCCAACTCCCGCATCTCGGTGATGGGCGGGGGACAGGcggggg                      >  1:2475236/1‑101 (MQ=255)
                                                                                  cgggcAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGGCCAACTCCCGCATCTCGGTGATGGGCGGCGAACAGGCGGCgggg                    >  1:407974/1‑101 (MQ=255)
                                                                                    ggcAATTACGGCATGGC.GG.GCCGGGGCTATTCGCCCCCGTTCCTGTGGACCT.GGCCCCACTCCCCCATCTCGGGGATGGGGGGCGAACAGGCGGGgggggg                  >  2:689613/1‑100 (MQ=255)
                                                                                    ggcAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTCCCGCATCTCGGTGATGGGCGGCGAACAGGCGGCGGGGGt                  <  1:3598654/101‑1 (MQ=255)
                                                                                      cAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGAGCT.GGCCCAACTCCCGCATCTCGGTGATGGGCGGCGCACAGGCGGGGGGGGGGc                >  1:740192/1‑101 (MQ=255)
                                                                                      cAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTCCCGCATCTCGGTGATGGGCGGCGAACAGGCGGCGGGGGTGc                <  2:1852826/101‑1 (MQ=255)
                                                                                        aTTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTCCCGCATCTCGGTGATGGGCGGCGAACAGGCGGCGGGGGTGCtt              <  2:3248876/101‑1 (MQ=255)
                                                                                        aTTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCCGGTCCTGTGGACCT.GGCCCAACTCCCGCATCTCGGTGATGGGCGGCGAACAGGGGGCggggggggtt              >  2:425097/1‑96 (MQ=255)
                                                                                          tACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGGTCCTGTGGGCCT.GGCCCAACTCCCGCATCTCGGTGATGGGCGGCGGACAGGGGGCggggggggctgc            >  2:14696/1‑94 (MQ=255)
                                                                                           aCGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTCCCGCATCTCGGTGATGGGCGGCGAACAGGCGGCGGGGGTGCTTGcc           <  1:2040395/101‑1 (MQ=255)
                                                                                              gCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTCCCGCATCTCGGTGATGGGCGGCGAACAGGCGGCGGGGGTGCTTGCCACg        <  1:662504/101‑1 (MQ=255)
                                                                                              gCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTCCCGCATCTCGGTGATGGGCGGCGAACAGGCGGCGGGGGGGGTTGCCACg        >  1:3244817/1‑101 (MQ=255)
                                                                                                aTGGG.GG.GGCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTCCCGCATCTCGGTGATGGGCGGCGAACAGGCGGCGGGGGTGCTTGCCACGGt      <  1:1370221/101‑1 (MQ=255)
                                                                                                aTGGC.GG.GCCGGGGCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTCCCGCATCTCGGTGATGGGGGGCGAACAGGCGGCGGGGGGGGTTGCCACggg      >  2:4314918/1‑100 (MQ=255)
                                                                                                aTGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTCCCGCATCTCGGTGATGGGCGGCGAACAGGCGGCGGGGGTGCTTGCCACGGt      <  1:820682/101‑1 (MQ=255)
                                                                                                 tGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTCCCGCATCTCGGTGATGGGCGGCGAACAGGCGGCGGGGGTGCTTGCCACGGTc     <  1:1587010/101‑1 (MQ=255)
                                                                                                   gC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTCCCGCATCTCGGTGATGGGCGGCGAACAGGCGGCGGGGGTGCTTGCCACGGTCaa   <  1:4520652/101‑1 (MQ=255)
                                                                                                    c.gg.gccgggcCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTCCCGCATCTCGGTGATGGGCGGCGAACAGGCGGCGGGGGGGCTTGCCACGGTCAAg  >  2:276461/1‑101 (MQ=255)
                                                                                                     |                                                                                                      
CGGCATGGGGCGAAGATGGTGACGGCCGTCGCCACGACGAACGTGCCGAAGATCACGCTGCTGGTGGGCGGCTCGTTCGGGGCGGGCAATTACGGCATGGC.GG.GCCGGGCCTATTCGCCCCGGTTCCTGTGGACCT.GGCCCAACTCCCGCATCTCGGTGATGGGCGGCGAACAGGCGGCGGGGGTGCTTGCCACGGTCAAG  >  NC_014034/1639486‑1639686

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 28 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG