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

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034826,0380C→. 3.4% 11.0 ‑1.6 ‑0.2176coding (410/672 nt)thiNthiamine pyrophosphokinase
Reads supporting (aligned to +/- strand):  new base (1/5):  ref base (87/83):  total (88/88)
Fisher's exact test strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.97e-01

Only 100 of 214 total aligned reads displayed.

GGCTGCGGCTTGATCACACGCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTCCGGGCGCT.CGACCGGGCTGCGCTGGCCGATC  >  NC_014034/825939‑826138
                                                                                                   |                                                                                                     
  cTGCGGCTTGATCACACGCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGGCGAGGTGATCTTTCGCGcccc                                                                                                    >  2:4408191/1‑101 (MQ=255)
   tGCGGCTTGATCACACGCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGGGATCTTTCGCGccccc                                                                                                   >  2:1362313/1‑101 (MQ=255)
    gCGGCTTGATCACACGCTGGCGGCGCTGACCGAACTGACCACCCGGCCCGGGCAGATCGGTATCCTGATCGCCCAGGACGAGGTGACCTGTCGCGcccccc                                                                                                  >  2:529745/1‑101 (MQ=255)
       gCTTGATCACACGCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGccccccccc                                                                                               <  1:4198435/101‑1 (MQ=255)
       gCTTGATCACACGCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGccccccccc                                                                                               >  2:3289015/1‑101 (MQ=255)
        cTTGATCACACGCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATTGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCg                                                                                              >  1:1318798/1‑101 (MQ=255)
        cTTGATCACACGCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCg                                                                                              <  1:2700038/101‑1 (MQ=255)
          tGATCACACGCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGc                                                                                            >  1:184794/1‑101 (MQ=255)
           gATCACACGCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGcg                                                                                           >  2:1480705/1‑100 (MQ=255)
             tCACACGCTGGCGGCGCTGACCGAGCTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTgg                                                                                         <  2:3707494/101‑1 (MQ=255)
               acacGCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGgcg                                                                                       <  2:4264804/101‑1 (MQ=255)
               acacGCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGgcg                                                                                       <  1:4486443/101‑1 (MQ=255)
                cacGCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGgcgc                                                                                      >  1:86459/1‑101 (MQ=255)
                  cGCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGgcgccg                                                                                    >  2:3108450/1‑99 (MQ=255)
                   gCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTgg                                                                                   <  1:4945531/101‑1 (MQ=255)
                   gCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTgg                                                                                   >  1:305990/1‑101 (MQ=255)
                    cTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGa                                                                                  <  2:3652147/101‑1 (MQ=255)
                      ggcggcGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATc                                                                                <  2:2171468/101‑1 (MQ=255)
                      ggcggcGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATc                                                                                <  1:783940/101‑1 (MQ=255)
                      ggcggcGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCCGCTGGCGCTGGGTc                                                                                >  1:549593/1‑101 (MQ=255)
                        cggcgCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCGGCTGGCGCTGGATCTGc                                                                             <  1:3832880/101‑1 (MQ=255)
                         ggcgCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGc                                                                             <  2:475047/101‑1 (MQ=255)
                          gcgcTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGcc                                                                            <  2:1089418/101‑1 (MQ=255)
                           cgcTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGccc                                                                           <  2:171117/101‑1 (MQ=255)
                            gcTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCt                                                                          >  2:2494888/1‑101 (MQ=255)
                             cTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCtt                                                                         >  2:238819/1‑101 (MQ=255)
                             cTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCtt                                                                         <  2:4336066/101‑1 (MQ=255)
                              tGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTg                                                                        <  2:120996/101‑1 (MQ=255)
                               gACCGACCTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTgg                                                                       <  1:2982902/101‑1 (MQ=255)
                               gACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTgg                                                                       <  2:3514923/101‑1 (MQ=255)
                               gACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTgg                                                                       <  1:4594506/101‑1 (MQ=255)
                               gACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTgg                                                                       >  2:89424/1‑101 (MQ=255)
                               gACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTgg                                                                       <  1:1409973/101‑1 (MQ=255)
                                aCCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTggg                                                                      <  1:3014581/101‑1 (MQ=255)
                                 ccGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGc                                                                     <  2:870474/101‑1 (MQ=255)
                                 ccGAACTGACCACGCGGCCCGGGCAGATCGGGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGc                                                                     <  1:3770880/101‑1 (MQ=255)
                                  cGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCa                                                                    <  2:3168863/101‑1 (MQ=255)
                                   gAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCAc                                                                   >  2:1189538/1‑101 (MQ=255)
                                      cTGCCCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCAcgcg                                                                <  1:1365067/101‑1 (MQ=255)
                                      cTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGAGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGACCTGCCCATGGGCAcgcg                                                                >  2:2105822/1‑101 (MQ=255)
                                        gACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGGGCTGGATCTGCCCTTGGGCACGGGGc                                                              >  2:501456/1‑101 (MQ=255)
                                        gACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGc                                                              >  1:3635071/1‑101 (MQ=255)
                                         aCCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGGCGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCt                                                             <  2:3998484/101‑1 (MQ=255)
                                         aCCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCt                                                             <  1:1715113/101‑1 (MQ=255)
                                          ccACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTg                                                            <  2:2336173/101‑1 (MQ=255)
                                            aCGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTc                                                          >  1:3428812/1‑101 (MQ=255)
                                                   ccGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCAt                                                   >  1:3338598/1‑101 (MQ=255)
                                                   ccGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATc                                                  <  2:3798621/101‑1 (MQ=255)
                                                    cGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGGGCTGGGTCTGCCCTTGGGCACGCGGCGGTCGCTGCGTc                                                  >  1:3906945/1‑101 (MQ=255)
                                                    cGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATc                                                  <  1:3846053/101‑1 (MQ=255)
                                                     gggCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATcc                                                 <  1:760101/101‑1 (MQ=255)
                                                     gggCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATcc                                                 >  1:1556739/1‑101 (MQ=255)
                                                     gggCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATcc                                                 >  1:1518140/1‑101 (MQ=255)
                                                      ggCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCg                                                >  1:2342799/1‑101 (MQ=255)
                                                        gaGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGAt                                              <  2:4937655/100‑1 (MQ=255)
                                                        cAGATCGTGATCCTGATCGCGGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGAt                                              <  1:3533123/101‑1 (MQ=255)
                                                        cAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATcgagt                                              >  2:3417389/1‑97 (MQ=255)
                                                          gATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATgg                                            <  2:3744954/101‑1 (MQ=255)
                                                           aTCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATggg                                           <  1:1784313/101‑1 (MQ=255)
                                                             cGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGAc                                         >  1:2814405/1‑101 (MQ=255)
                                                              gTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGAcc                                        >  2:4508084/1‑101 (MQ=255)
                                                               tGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCg                                       <  1:2055692/101‑1 (MQ=255)
                                                                gATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGGACGCGGCCGTCCCCGCATCCGTTGGGACCgg                                      >  2:3901749/1‑101 (MQ=255)
                                                                  tCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCt                                    >  1:1487303/1‑101 (MQ=255)
                                                                   ccTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCtt                                   <  1:2610433/101‑1 (MQ=255)
                                                                    cTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTc                                  <  2:303665/101‑1 (MQ=255)
                                                                    cTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTc                                  <  1:3226648/101‑1 (MQ=255)
                                                                    cTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCATTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTc                                  >  2:4237985/1‑101 (MQ=255)
                                                                    cTGATCGCCGAGGACGAGGTGATCTTTCGCG.CCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTcc                                 <  2:2556912/101‑1 (MQ=255)
                                                                       aTCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTcccg                               >  2:2323658/1‑99 (MQ=255)
                                                                       aTCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTCCgg                               <  1:1888731/101‑1 (MQ=255)
                                                                       aTCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTCCgg                               <  2:268821/101‑1 (MQ=255)
                                                                        tCGGCGAGGACGGGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTCCggg                              <  2:2457676/101‑1 (MQ=255)
                                                                        tCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTCCggg                              <  1:2262012/101‑1 (MQ=255)
                                                                            cGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTCCggggggc                          >  2:4842903/1‑97 (MQ=255)
                                                                             gAGGACGAGGTGATCATTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTCCGGGCGCT.c                        >  1:4258911/1‑101 (MQ=255)
                                                                              aGGACGAGGTGATCTTTCGCGCCCGCCCCCGGCTGGCGCTGGATCGGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTCCGGGCGCT.Cg                       <  2:565548/101‑1 (MQ=255)
                                                                              aGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTTCGGGCGCC.Cg                       >  1:4631335/1‑101 (MQ=255)
                                                                              aGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTCCGGGCGCT.Cg                       <  1:668393/101‑1 (MQ=255)
                                                                              aGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTCCGGGCGCT.Cg                       <  1:3289277/101‑1 (MQ=255)
                                                                               ggACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTTCGGGGGCC.Cgc                      >  1:1469230/1‑100 (MQ=255)
                                                                               ggACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTTCCGGCGCC.Cgg                      >  2:74797/1‑100 (MQ=255)
                                                                                 aagAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTCCGGGCGCT.CGAcc                    <  1:1269384/99‑1 (MQ=255)
                                                                                 aCGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTCCGGGCGCT.CGAcc                    <  2:1335040/101‑1 (MQ=255)
                                                                                 aCGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTCCCGGCGCT.ccccc                    >  1:3189391/1‑97 (MQ=255)
                                                                                  cGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTCCGGGCGCT.CGACCg                   <  1:3990194/101‑1 (MQ=255)
                                                                                  cGAGGTGATCTTTCGCG.CCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTCCGGGCGCT.CGGCCgg                  >  1:4530069/1‑101 (MQ=255)
                                                                                    gggTGATCTTTCGCG.CCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTCCGGGCGCT.CGACCGGGc                <  1:4408510/100‑1 (MQ=255)
                                                                                     ggTGATCTTTCGCGCCCCCCCCCGGCTTGCGCTGGAGCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTCCGGGCGCT.AGACCGGGc                <  2:143335/101‑1 (MQ=255)
                                                                                     ggTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTCCTGGCGCT.CGACCGGGc                >  2:3300776/1‑101 (MQ=255)
                                                                                      gTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTCCCGGGGCC.CGCCCcgggg               >  1:3366074/1‑99 (MQ=255)
                                                                                       tGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTCCGGGGGAC.CGACCCGGGTg              >  2:1986338/1‑101 (MQ=255)
                                                                                        gATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTCCGGGCGCT.CGACCGGGCTgc             <  1:3947034/101‑1 (MQ=255)
                                                                                           cTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTCCGGGCGCT.CGACCGGGCTGCGCt          >  1:2462427/1‑101 (MQ=255)
                                                                                             ttCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTCCCGGGGCG.CCGCCGGGccccgccgg        >  2:3264762/1‑93 (MQ=255)
                                                                                               cgcgCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTCCGGGCGCTCCGACC.GGCTGCGCGGGcc      >  2:2876391/1‑101 (MQ=255)
                                                                                               cgcgCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTCCCGGCGCT.CGACCCGGCCGCGCTGGcc      >  1:4762572/1‑101 (MQ=255)
                                                                                                  gcCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTCCGGGCGCT.CGACCGGGCTGCGCTGGCCgct   >  1:980345/1‑99 (MQ=255)
                                                                                                  gcCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTCCGGGCGCT.CGACCGGGCTGCGCTGGCCGac   >  2:2822295/1‑100 (MQ=255)
                                                                                                  gcCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTCCGGGCGCT.CGACCGGGCTGCGCTGGCCGAt   >  2:358186/1‑101 (MQ=255)
                                                                                                   |                                                                                                     
GGCTGCGGCTTGATCACACGCTGGCGGCGCTGACCGAACTGACCACGCGGCCCGGGCAGATCGTGATCCTGATCGCCGAGGACGAGGTGATCTTTCGCGCCCCCCCCCGGCTGGCGCTGGATCTGCCCTTGGGCACGCGGCTGTCGCTGCATCCGATGGGACCGGCTTCCGGGCGCT.CGACCGGGCTGCGCTGGCCGATC  >  NC_014034/825939‑826138

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 24 ≤  ATCG  < 32 ≤  ATCG  < 40 ≤  ATCG