Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 124,859 G→C 5.6% A615P (GCT→CCT) ‡ aceE → pyruvate dehydrogenase, decarboxylase component E1, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913124,8590GC 5.6% 418.5 / 12.8 126A615P (GCT→CCT) ‡aceEpyruvate dehydrogenase, decarboxylase component E1, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base G (60/59);  new base C (4/3);  total (64/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.62e-01

Only 100 of 130 total aligned reads displayed.

TTCTGCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTA  >  NC_000913/124710‑125004
                                                                                                                                                     |                                                                                                                                                 
ttCTGCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGggcgg                                                                                                                                                   >  1:559099/1‑150 (MQ=255)
ttCTGCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGggcgg                                                                                                                                                   <  1:195841/150‑1 (MQ=255)
  cTGCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCt                                                                                                                                                 >  2:684293/1‑150 (MQ=255)
    gCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTgg                                                                                                                                               <  1:1483377/150‑1 (MQ=255)
    gCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTgg                                                                                                                                               <  2:823319/150‑1 (MQ=255)
     cAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGTCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGc                                                                                                                                              <  1:156819/150‑1 (MQ=255)
     cAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGc                                                                                                                                              <  2:308451/150‑1 (MQ=255)
      aGGACGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCCACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCGTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCg                                                                                                                                             <  2:948869/150‑1 (MQ=255)
      aGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCg                                                                                                                                             <  2:689635/150‑1 (MQ=255)
       ggAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGa                                                                                                                                            >  1:1352462/1‑150 (MQ=255)
        gAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGAc                                                                                                                                           <  2:905747/150‑1 (MQ=255)
         aaGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGAcc                                                                                                                                          <  2:1449238/150‑1 (MQ=255)
          aGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGAccc                                                                                                                                         >  2:702524/1‑149 (MQ=255)
          aGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCa                                                                                                                                         >  2:1140056/1‑150 (MQ=255)
           gggATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATGTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGAGCTGGGCGGCTGGCGACCAg                                                                                                                                        >  2:1091181/1‑150 (MQ=255)
           gggATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAg                                                                                                                                        <  2:515604/150‑1 (MQ=255)
              aTCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGCTCCCGTTCTACCTCTATTACTCGATGATCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCaa                                                                                                                                     >  1:94035/1‑150 (MQ=255)
               tCAACGAGCTGGGCGCGGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGTCGCCAGCCGCCCAGCaca                                                                                                                                    <  2:556718/150‑3 (MQ=255)
                         gggCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGGGGCtt                                                                                                                          >  1:924817/1‑150 (MQ=255)
                           gcgcAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTcc                                                                                                                        >  2:1401113/1‑150 (MQ=255)
                             gcAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTg                                                                                                                      >  2:61944/1‑150 (MQ=255)
                              cAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGa                                                                                                                     <  2:19362/150‑1 (MQ=255)
                                ggTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATc                                                                                                                   >  2:208616/1‑150 (MQ=255)
                                ggTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATc                                                                                                                   >  1:1159568/1‑150 (MQ=255)
                                   tgttCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATcggc                                                                                                                <  2:878422/150‑1 (MQ=255)
                                            ctggcAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTcc                                                                                                       <  2:823890/150‑1 (MQ=255)
                                            ctggcAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTcc                                                                                                       <  1:684293/150‑1 (MQ=255)
                                             tggcAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCg                                                                                                      <  1:619201/150‑1 (MQ=255)
                                             gggcAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTGCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCg                                                                                                      <  1:18152/149‑1 (MQ=255)
                                              ggcAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCgg                                                                                                     >  1:1627004/1‑150 (MQ=255)
                                                cAGCGGCGACCTCTTACAGCGCCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTc                                                                                                   <  1:1641959/150‑1 (MQ=255)
                                                cAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTc                                                                                                   >  2:311681/1‑150 (MQ=255)
                                                cAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGAGCCCGTTCTACATCGATGACTCGATGTTCGGCTTCCAGCGTATTGGCGCTCTGTGCTGGGCTGCTGGCGACCAGCAAGCGCGTGGCTTCCTGCTCGGCGGTACTTCCGGTc                                                                                                   >  1:1196695/1‑150 (MQ=255)
                                                  gcggcgACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGt                                                                                                 <  1:1754614/150‑1 (MQ=255)
                                                   cggcgACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTa                                                                                                <  1:109709/150‑1 (MQ=255)
                                                   cggcgACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGGGGCTTCCTGATCGGCGGTACTTCCGGTCGTa                                                                                                <  1:992231/150‑1 (MQ=255)
                                                    ggcgACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTac                                                                                               <  1:1674837/150‑1 (MQ=255)
                                                        aCCTCTTACAGCACCAACAATCTGCCGATTATCCCGTTCTACATCTATTACTCGATGTTCGGCTTACAGCGTATTGGCGATATGTGCTGGGCGGCTGGCGTACAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTaccacc                                                                                           <  2:1496692/150‑1 (MQ=255)
                                                           tctTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGCACTTCCGGTCGTACCACCCTg                                                                                        <  1:1015275/150‑1 (MQ=255)
                                                            ctTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGa                                                                                       <  1:1201451/150‑1 (MQ=255)
                                                              tACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAAc                                                                                     <  1:1558837/150‑1 (MQ=255)
                                                              tACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAAc                                                                                     >  2:352295/1‑150 (MQ=255)
                                                               aCAGCACCAACAATCTGCCGATGATCCCGTTCTACATCGATTTCTCGATGTTCGGGTTCCAGCGGATAGGCGATCTGTGCTGTGCGGCTGGCGAGCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCTGTCGTACCACCCTGACCg                                                                                    >  2:92765/1‑150 (MQ=255)
                                                                cAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACgg                                                                                   >  1:644772/1‑150 (MQ=255)
                                                                  gCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCg                                                                                 >  1:1212773/1‑150 (MQ=255)
                                                                  gCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCg                                                                                 >  2:35993/1‑150 (MQ=255)
                                                                  gCACCAACAATCTGCCGATGATCCCGTGCTACATCTATTACTCGATGTTCGGCTGCCCGCGTATTGGCGATCTGTGCGGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCg                                                                                 <  1:1140913/150‑1 (MQ=255)
                                                                   cACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGa                                                                                <  1:526976/150‑1 (MQ=255)
                                                                   cACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCGGGCGACCAGCAAGCGCGGGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGa                                                                                <  1:1148200/150‑1 (MQ=255)
                                                                     ccAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAg                                                                              >  1:1781825/1‑150 (MQ=255)
                                                                     ccAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAg                                                                              >  2:1490528/1‑150 (MQ=255)
                                                                     ccAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAg                                                                              >  1:1009630/1‑150 (MQ=255)
                                                                      caacaaTCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAgg                                                                             <  1:1464236/150‑1 (MQ=255)
                                                                      caacaaTCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAgg                                                                             >  2:899041/1‑150 (MQ=255)
                                                                        acaaTCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTc                                                                           <  2:1162773/150‑1 (MQ=255)
                                                                           aTCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTgc                                                                        <  1:1720080/150‑1 (MQ=255)
                                                                                ccGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCAc                                                                   >  1:4966/1‑150 (MQ=255)
                                                                                 cGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGTCGACCAGCAAGCGCGAGGCTTACTGATCGGCGTTACTTCCGCTCGTACCACCCTGATCGTAGAAGGTCTGCAGCACg                                                                  >  1:1105843/1‑150 (MQ=255)
                                                                                    tgatCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAg                                                               >  2:76692/1‑150 (MQ=255)
                                                                                        cccGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATgg                                                           <  2:210766/150‑1 (MQ=255)
                                                                                         ccGTTCTACATCTATTACTCGATGTTCGGCTTCCCGCGTAGTGGCTATCTGTGCTGGGCGGCGGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGt                                                          <  2:101784/150‑1 (MQ=255)
                                                                                         ccGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTAGTACCACCCTGAACGGAGAAGGTCTGCAGCACGAAGATGGt                                                          >  2:1673690/1‑150 (MQ=255)
                                                                                          cGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTc                                                         <  2:1011567/150‑1 (MQ=255)
                                                                                             tCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTcaca                                                      >  1:678075/1‑150 (MQ=255)
                                                                                               tACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCTGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGc                                                    >  1:980792/1‑150 (MQ=255)
                                                                                               tACATCTATTACTCGATGTTCGGCATCCTGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGc                                                    <  1:47097/150‑1 (MQ=255)
                                                                                                aCATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCAGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGcc                                                   >  1:411776/1‑150 (MQ=255)
                                                                                                   tCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCcaca                                                >  1:1371735/1‑150 (MQ=255)
                                                                                                      aTTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTc                                             >  1:966036/1‑150 (MQ=255)
                                                                                                      aTTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCCGCGGTACTTCCGGGCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTc                                             <  1:1237584/150‑1 (MQ=255)
                                                                                                       ttACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCa                                            <  1:836508/150‑1 (MQ=255)
                                                                                                        tACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAg                                           <  1:892351/150‑1 (MQ=255)
                                                                                                        tACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAg                                           <  2:1352462/150‑1 (MQ=255)
                                                                                                        tACTCGATGTTCGGCTTCCAGCGTATTGGCGAGCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAg                                           <  2:1295529/150‑1 (MQ=255)
                                                                                                         aCTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGt                                          >  2:1005121/1‑150 (MQ=255)
                                                                                                         aCTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGt                                          <  1:1221953/150‑1 (MQ=255)
                                                                                                         aCTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGt                                          >  2:1309467/1‑150 (MQ=255)
                                                                                                             gATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCt                                      <  1:962743/150‑1 (MQ=255)
                                                                                                                    ggCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGCCTATc                               >  1:118696/1‑150 (MQ=255)
                                                                                                                    ggCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATc                               >  1:255970/1‑150 (MQ=255)
                                                                                                                    ggCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATc                               >  1:553245/1‑150 (MQ=255)
                                                                                                                        tCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTCAAGCACGAAGATGGACACAGCCACATTCAGTCGCTGACTATCCCGa                           >  1:1121791/1‑150 (MQ=255)
                                                                                                                         ccAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACCGCCACATTCAGTCGCTGACTATCCCGaa                          >  1:373701/1‑150 (MQ=255)
                                                                                                                          catcGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAAc                         <  2:1682406/147‑1 (MQ=255)
                                                                                                                          cAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAAc                         <  1:352295/150‑1 (MQ=255)
                                                                                                                                aTTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTAtc                   >  1:48523/1‑150 (MQ=255)
                                                                                                                                      gATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTAc             >  1:1351124/1‑150 (MQ=255)
                                                                                                                                       aTCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACg            <  1:1083459/150‑1 (MQ=255)
                                                                                                                                         cTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGACCTGTATCTCTTACGAc          >  1:1067538/1‑150 (MQ=255)
                                                                                                                                          cttgCTGGTCGCCAGCCGCCCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCGCCGCCACCTTCAGTCGCTGACTCTCCCGAACTGTATCTCTTACGAcc         >  1:109954/3‑150 (MQ=255)
                                                                                                                                          cttgCTGGTCGCCAGCCGCCCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGAcc         >  1:729886/3‑150 (MQ=255)
                                                                                                                                           gtgCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGAccc        >  2:818868/1‑150 (MQ=255)
                                                                                                                                           gtgCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGAccc        >  1:1649547/1‑150 (MQ=255)
                                                                                                                                           gtgCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGAccc        >  1:1657131/1‑150 (MQ=255)
                                                                                                                                              cTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgc     <  2:1709905/150‑1 (MQ=255)
                                                                                                                                              cTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgc     >  1:810686/1‑150 (MQ=255)
                                                                                                                                                gggCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgctt   >  1:855402/1‑150 (MQ=255)
                                                                                                                                                gggCGGCTGGCGACCAGCAAGCGCGGGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgctt   <  2:924817/150‑1 (MQ=255)
                                                                                                                                                 ggctgcTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgctta  <  1:642038/150‑1 (MQ=255)
                                                                                                                                                 ggcggcTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgctta  <  1:76692/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
TTCTGCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTA  >  NC_000913/124710‑125004

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: