Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 826,038 Δ1 bp 4.0% coding (410/672 nt) thiN → thiamine pyrophosphokinase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034826,0380C→. 4.0% 14.2 ‑0.3 ‑0.5174coding (410/672 nt)thiNthiamine pyrophosphokinase
Reads supporting (aligned to +/- strand):  new base (4/3):  ref base (82/85):  total (86/88)
Fisher's exact test strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.21e-01

Only 100 of 218 total aligned reads displayed.

GGGCTGCGGCTTGATCACACGCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGC.TTCCGGGCGCTCGACCGGGCTGCGCTGGCCGATC  >  NC_014034/825938‑826138
                                                                                                     |                                                                                                     
 ggCTGCGGCTTGATCACACGCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.cc                                                                                                      <  2:1929704/101‑1 (MQ=255)
  gCTGCGGCTTGATCACACGCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.ccc                                                                                                     <  1:3714895/101‑1 (MQ=255)
    tGCGGCTTGATCACACGCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.ccccc                                                                                                   <  1:1088405/101‑1 (MQ=255)
    tGCGGCTTGATCACACGCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.ccccc                                                                                                   <  1:870451/101‑1 (MQ=255)
    tGCGGCTTGATCACACGCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.ccccc                                                                                                   <  2:355713/101‑1 (MQ=255)
     gCGGCTTGATCACACGCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.cccccc                                                                                                  >  1:3904200/1‑101 (MQ=255)
        gCTTGATCACCCGCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.ccccccccc                                                                                               <  1:534834/101‑1 (MQ=255)
        gCTTGATCACACGCTGGCGGCGCTGACCGAACTGACCTCGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.ccccccccc                                                                                               <  1:3305800/101‑1 (MQ=255)
        gCTTGATCACACGCGGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.ccccccccc                                                                                               <  2:1746825/101‑1 (MQ=255)
         cTTGATCACACGCTGGCGGCGCTGACCGAACTGACCGCGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCg                                                                                              <  2:569572/101‑1 (MQ=255)
         cTTGATCACACGCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCg                                                                                              <  1:2464428/101‑1 (MQ=255)
          ttGATCACACGCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCgg                                                                                             >  2:2436222/1‑101 (MQ=255)
          ttGATCACACGCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCgg                                                                                             >  1:4371510/1‑101 (MQ=255)
          ttGATCACACGCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCgg                                                                                             <  1:2590522/101‑1 (MQ=255)
           tGATCACACGCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGc                                                                                            >  1:3972919/1‑101 (MQ=255)
            gATCACACGCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGcg                                                                                           >  1:1236452/1‑100 (MQ=255)
            gATCACACGCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCt                                                                                           <  1:2786497/101‑1 (MQ=255)
             aTCACACGCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTg                                                                                          >  1:2737042/1‑101 (MQ=255)
               cacccacTGGCGGCGCTGACCGAACTGGCCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGgc                                                                                        <  2:3893730/95‑1 (MQ=255)
                 cacGCTGGCGGCGCTGACCGAACTGACCACTCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGgcgc                                                                                      >  1:1663704/1‑101 (MQ=255)
                  ccGCTGGCGGCGCTGACCGAACTGGCCACGCGGCCCGGGCAGAGCGTGATCCTGACCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCt                                                                                     <  2:3956622/100‑1 (MQ=255)
                    gCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTgg                                                                                   >  2:2979283/1‑101 (MQ=255)
                    gCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTgg                                                                                   >  1:2390470/1‑101 (MQ=255)
                     cTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGa                                                                                  <  1:2893281/101‑1 (MQ=255)
                       ggcggcGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCCGGATc                                                                                >  2:2618870/1‑101 (MQ=255)
                        gcggcgCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCt                                                                               <  1:3111288/101‑1 (MQ=255)
                         cggcgCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTg                                                                              >  2:4065904/1‑101 (MQ=255)
                         cggcgCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTg                                                                              <  2:628095/101‑1 (MQ=255)
                         cggcgCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTg                                                                              >  1:1834428/1‑101 (MQ=255)
                           gcgcTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGcc                                                                            >  2:1341274/1‑101 (MQ=255)
                            cgcTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCCCGCTGGCGCTGGATCTGcc                                                                            >  1:3827842/1‑101 (MQ=255)
                            cgcTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGccc                                                                           <  2:3528575/101‑1 (MQ=255)
                             gcTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCt                                                                          <  1:3490860/101‑1 (MQ=255)
                              cTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCtt                                                                         <  2:2293668/101‑1 (MQ=255)
                              cTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCtt                                                                         <  2:4730051/101‑1 (MQ=255)
                                gACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATTTTTCGCG.CCCCCCCCCCGCTGGCGCTGGATCCGCCCTTgg                                                                       >  1:4687396/1‑101 (MQ=255)
                                 aCCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTggg                                                                      <  1:3357778/101‑1 (MQ=255)
                                 aCCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTggg                                                                      <  1:176652/101‑1 (MQ=255)
                                  ccGAACTGGCCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGc                                                                     <  1:1414248/101‑1 (MQ=255)
                                  ccGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGc                                                                     <  2:147352/101‑1 (MQ=255)
                                   cGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCa                                                                    <  2:1747874/101‑1 (MQ=255)
                                    gAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCAc                                                                   >  1:1304615/1‑101 (MQ=255)
                                    gAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCAc                                                                   <  2:1504158/101‑1 (MQ=255)
                                    gAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTGGGGCAc                                                                   <  2:191056/101‑1 (MQ=255)
                                        tGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCgg                                                               <  2:3712880/101‑1 (MQ=255)
                                         gACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGc                                                              >  1:1528325/1‑101 (MQ=255)
                                          cccACCCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCt                                                             <  1:4514508/100‑1 (MQ=255)
                                          aCCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCt                                                             <  2:4687140/101‑1 (MQ=255)
                                             aCGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTc                                                          >  2:3314314/1‑101 (MQ=255)
                                              cgcgGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCg                                                         <  2:1645926/101‑1 (MQ=255)
                                                cgGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCt                                                       <  2:4371263/101‑1 (MQ=255)
                                                 ggCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTg                                                      >  2:4296403/1‑101 (MQ=255)
                                                 ggCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTg                                                      >  1:2261639/1‑101 (MQ=255)
                                                    ccGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCAt                                                   <  2:3257758/101‑1 (MQ=255)
                                                    ccGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG..CCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATc                                                  >  2:2745108/1‑101 (MQ=255)
                                                     cGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCAGGATCTGGATCGGCACTTGGGCACGCGGCGGTCGCtcaatc                                                  >  2:1140736/1‑96 (MQ=255)
                                                      gggCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATcc                                                 <  2:3856061/101‑1 (MQ=255)
                                                      gggCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATcc                                                 <  1:319993/101‑1 (MQ=255)
                                                       ggCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCg                                                <  2:1435735/101‑1 (MQ=255)
                                                       ggCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCCGCGGGCGCTGGATCTGCCCTTGGGCACGAGGCGGTCGCGGCGTCCg                                                >  2:4266374/1‑101 (MQ=255)
                                                        gCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATcatt                                               <  2:669290/101‑4 (MQ=255)
                                                         cAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGAt                                              <  2:268784/101‑1 (MQ=255)
                                                           gATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATgg                                            <  1:252223/101‑1 (MQ=255)
                                                           gATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCACTGGATCTGCCCTTGGGCACGCGGCTGGCGCTGCATCCGATgg                                            >  2:3322380/1‑101 (MQ=255)
                                                              cGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGAc                                         >  1:2447384/1‑101 (MQ=255)
                                                              cGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCCGGTGGGCCTGGATGTGCCCCTGGGCACGGGGCCGTCCCTGGATCCCATGGGAc                                         >  2:4635752/1‑101 (MQ=255)
                                                              cGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG..CCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGAcc                                        >  2:1709823/1‑101 (MQ=255)
                                                               gTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGAcc                                        <  2:454922/101‑1 (MQ=255)
                                                                  aTCCTGATCGCCGAGGACGAGGTGATCTTTCGCG..CCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGC.t                                   >  2:4153852/1‑101 (MQ=255)
                                                                   tCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGC.t                                   <  2:1240639/101‑1 (MQ=255)
                                                                           gCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGC.TTCCcgggg                           >  2:2467669/1‑99 (MQ=255)
                                                                             cGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTTGCGCTGGATCTGTCCTTCGGCACGCCGCTGTCGCTGGGTCCGATGGGACCGGC.GTCCGGGGGCt                         >  2:2941045/1‑101 (MQ=255)
                                                                              gAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTTCCGGgggcc                         >  2:3733311/1‑100 (MQ=255)
                                                                              gAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGC.TTCCGGGCGCTc                        >  1:4654360/1‑101 (MQ=255)
                                                                               gggACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGC.TTCCGGGCGCTCg                       <  2:1580102/100‑1 (MQ=255)
                                                                               aGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGC.TTTCGGGGGCCCg                       >  2:1188666/1‑101 (MQ=255)
                                                                                ggACGAGGTGATCTTTCGCG..CCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGC.TTCCGGGCGCTCGAc                     <  1:4154098/101‑1 (MQ=255)
                                                                                 gACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGC.TTCCGGGCGCTCGAc                     <  1:710684/101‑1 (MQ=255)
                                                                                 gACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGC.TTCCGGGCGCCCGAc                     >  2:201214/1‑101 (MQ=255)
                                                                                  aCGAGGTGATCTTTCGCG.CCCCCCCCCGGCCGGCGCCGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCCTCCGGTGGGACCGGG.TTTCCGgcgccccccc                    >  1:3343483/1‑95 (MQ=255)
                                                                                    gAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGC.TTTCTGGGGCTTGAcccg                  >  1:4214546/1‑99 (MQ=255)
                                                                                     aGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGC.TTCCGGGCGCTCGACCggg                 <  1:2978362/101‑1 (MQ=255)
                                                                                         gATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGC.TTCCGGGGGCTCGCCCGGGCCgc             >  1:559831/1‑101 (MQ=255)
                                                                                         gATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGC.TTCCGGGCGCTCGACCGGGCTgc             <  1:1191860/101‑1 (MQ=255)
                                                                                            cTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGC.TTCCGGGCGCTCGACCGGGCTGCGCt          <  1:4356100/101‑1 (MQ=255)
                                                                                            cTTTCGCG..CCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGC.TTCCGGGCGCTCGACCGGGCTgcgccg         >  1:3982392/1‑99 (MQ=255)
                                                                                             tttCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGC.TTCTGGGCGCTCGACCGGGCTGCGCTg         >  1:4465888/1‑101 (MQ=255)
                                                                                             tttCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGC.TTCCGGGGGCTCGACCGGGCTGCCCTg         >  2:1842420/1‑101 (MQ=255)
                                                                                             tttCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGC.TTCCGGGCGCTCGACCGGGCTGCGCTg         <  2:3772949/101‑1 (MQ=255)
                                                                                             tttCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGC.TTCCGGGCGCTCCCCCGGGCTgcgccg         >  2:3208457/1‑99 (MQ=255)
                                                                                             tttCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGC.TTCCCGGCCCTCGGCCGGGCTggggtg         >  2:2184232/1‑96 (MQ=255)
                                                                                               tCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGC.TTCCGGGCGCTCGACCGGGCTGCGCTGGc       <  1:4253697/101‑1 (MQ=255)
                                                                                                cgcg.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGC.TTCCGGGCGCTCGACCGGGCTGCGCTGGcc      <  1:189522/101‑1 (MQ=255)
                                                                                                cgcg.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGC.TTCCGGGCGCTCGACCGGGCGGCGCTGGcc      >  1:1415294/1‑101 (MQ=255)
                                                                                                cgcg.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGC.TTCCGGGCGCTCGACCCGGCTGCGCTGGcc      >  2:2552807/1‑101 (MQ=255)
                                                                                                cgcg.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGC.TTCCCGGGGCTCGACCGGGCTGCGCTGGcc      >  2:3153171/1‑101 (MQ=255)
                                                                                                   g.cCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGC.TTCCGGGCGCTCGCCCGGGGTGCGCGGGCCGa  >  2:433011/1‑100 (MQ=255)
                                                                                                   g.cCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGC.TTCCGGGCGCTCGACCGGGCTGCGCTGGCCGAt   <  2:2380135/101‑1 (MQ=255)
                                                                                                    ccccccccccGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGC.TTCCGGGCGCTCGACCGGGCTGCGCTGGCCGAt   <  2:3972677/100‑1 (MQ=255)
                                                                                                     cccccccccGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGC.TTCCGGGCGCTCGACCGGGCTGCGCTGGCCGATc  <  2:1377447/101‑1 (MQ=255)
                                                                                                     |                                                                                                     
GGGCTGCGGCTTGATCACACGCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGC.TTCCGGGCGCTCGACCGGGCTGCGCTGGCCGATC  >  NC_014034/825938‑826138

Base quality scores:  ATCG  < 3 ≤  ATCG  < 19 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG