Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,578,097 A→G T148A (ACC→GCC)  tktB → transketolase 2 isozyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,578,0970AG99.7% 1033.6 / NA 298T148A (ACC→GCC) tktBtransketolase 2 isozyme
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (143/154);  total (143/155)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 311 total aligned reads displayed.

GGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGCTG  >  NC_000913/2578007‑2578186
                                                                                          |                                                                                         
ggTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGcc                                                                                         <  2:3526225/93‑1 (MQ=255)
  tttGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTa                                                                                       >  2:2313693/1‑93 (MQ=255)
    tGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTAtg                                                                                     <  2:3118514/93‑1 (MQ=255)
     ggCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTAtgt                                                                                    <  1:4473809/93‑1 (MQ=255)
       cGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATgtgt                                                                                  <  2:6321176/93‑1 (MQ=255)
        gAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGtt                                                                                 <  2:2990283/93‑1 (MQ=255)
        gAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGtt                                                                                 <  1:522728/93‑1 (MQ=255)
           cGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTAt                                                                              <  1:1949572/93‑1 (MQ=255)
           cGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTAt                                                                              >  2:581295/1‑93 (MQ=255)
            gCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATg                                                                             >  1:1562576/1‑93 (MQ=255)
              cGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATggg                                                                           <  1:7928052/93‑1 (MQ=255)
              cGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATggg                                                                           <  2:3751425/93‑1 (MQ=255)
               gTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGc                                                                          <  2:2157287/93‑1 (MQ=255)
                 cGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCCGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCTCCTATGCGTTTATGGGCGa                                                                        >  1:4199215/1‑93 (MQ=255)
                  gggCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGAc                                                                       >  1:7506319/1‑93 (MQ=255)
                  gggCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGAc                                                                       <  2:2698643/93‑1 (MQ=255)
                  gggCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGAc                                                                       <  2:7849522/93‑1 (MQ=255)
                   ggCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACg                                                                      >  2:4598532/1‑93 (MQ=255)
                    gCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACgg                                                                     >  2:1686928/1‑93 (MQ=255)
                      tGGCGATAGCAGAGCGTACACTGGCTGCGCAGTTTAACCAGCCAGACAATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCTACGGCt                                                                   >  2:6737489/1‑93 (MQ=255)
                      tGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCt                                                                   <  1:4805382/93‑1 (MQ=255)
                      tGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCt                                                                   >  1:8009581/1‑93 (MQ=255)
                      tGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCt                                                                   >  1:1195248/1‑93 (MQ=255)
                         cGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCCATCACTTCGCCTATGTGTTTATGGGCGACGGCTGcc                                                                <  1:6286324/93‑1 (MQ=255)
                          gATAGCAGGGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCt                                                               <  2:7610787/93‑1 (MQ=255)
                           aTAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTg                                                              >  1:7114893/1‑93 (MQ=255)
                            tAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGa                                                             >  2:4836483/1‑93 (MQ=255)
                            tAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGa                                                             >  2:5175746/1‑93 (MQ=255)
                            tAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGa                                                             <  2:1235751/93‑1 (MQ=255)
                             aGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGAt                                                            >  1:5655690/1‑93 (MQ=255)
                              gCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATg                                                           >  1:330614/1‑93 (MQ=255)
                              gCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGGTTATGGGCGACGGCTGCCTGATg                                                           >  2:4736375/1‑93 (MQ=255)
                                agagCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGa                                                         <  1:8649603/93‑1 (MQ=255)
                                 gagCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGaa                                                        <  1:3936015/93‑1 (MQ=255)
                                   gCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAgg                                                      <  1:5456412/93‑1 (MQ=255)
                                   gCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAgg                                                      >  1:6920681/1‑93 (MQ=255)
                                     gTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTa                                                    <  2:8254534/93‑1 (MQ=255)
                                     gTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTa                                                    >  1:8499047/1‑93 (MQ=255)
                                      tACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTAt                                                   >  2:5686289/1‑93 (MQ=255)
                                      tACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTAt                                                   <  1:2599864/93‑1 (MQ=255)
                                      tACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTAt                                                   >  1:597497/1‑93 (MQ=255)
                                      tACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTAt                                                   >  2:6159472/1‑93 (MQ=255)
                                        cacTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTAttt                                                 >  1:5740080/1‑93 (MQ=255)
                                        cacTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTAttt                                                 >  1:4579853/1‑93 (MQ=255)
                                                gcgcAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGaa                                         <  2:8773848/93‑1 (MQ=255)
                                                gcgcAGTTTAACCAGACAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGaa                                         <  1:7813061/93‑1 (MQ=255)
                                                 cgcAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAg                                        <  1:5169117/93‑1 (MQ=255)
                                                 cgcAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAg                                        >  2:859814/1‑93 (MQ=255)
                                                  gcAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGt                                       >  1:6435448/1‑93 (MQ=255)
                                                   cAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTc                                      <  1:5041376/93‑1 (MQ=255)
                                                   cAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTc                                      <  2:4156609/93‑1 (MQ=255)
                                                    aGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCt                                     <  1:4940115/93‑1 (MQ=255)
                                                    aGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCt                                     >  1:1007143/1‑93 (MQ=255)
                                                     gTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTg                                    >  2:2012867/1‑93 (MQ=255)
                                                      tttAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGt                                   <  2:1017539/93‑1 (MQ=255)
                                                       ttAACCAGCCAGACGATGAGATCGTCGATCACTTCGCCTATGTCTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGtt                                  <  2:4263331/93‑1 (MQ=255)
                                                       ttAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGtt                                  <  2:8202517/93‑1 (MQ=255)
                                                        tAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGGAGGTATTTCCCACGAAGTCTGTTc                                 <  1:5851209/93‑1 (MQ=255)
                                                         aaCCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCg                                >  1:8072079/1‑93 (MQ=255)
                                                          aCCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGc                               <  1:6346356/93‑1 (MQ=255)
                                                           ccagccagACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCt                              <  1:3078655/93‑1 (MQ=255)
                                                            cagccagACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTg                             <  2:3270184/93‑1 (MQ=255)
                                                             agccagACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTgg                            <  1:2292145/93‑1 (MQ=255)
                                                              gccagACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTggc                           <  2:1461659/93‑1 (MQ=255)
                                                              gccagACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTggc                           >  2:2265704/1‑93 (MQ=255)
                                                               ccagACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTggca                          >  2:4337548/1‑93 (MQ=255)
                                                               ccagACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTggca                          >  1:7027225/1‑93 (MQ=255)
                                                                 cgACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTggcagg                        <  1:6673467/92‑1 (MQ=255)
                                                                 agACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTggcagg                        <  2:7343321/93‑1 (MQ=255)
                                                                  gACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTggtaggc                       >  2:6600216/1‑93 (MQ=255)
                                                                  gACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTggcaggc                       <  2:345349/93‑1 (MQ=255)
                                                                  gACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTggcaggc                       <  2:1095663/93‑1 (MQ=255)
                                                                  gACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTggcaggc                       >  2:5569722/1‑93 (MQ=255)
                                                                  gACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTggcaggc                       >  2:5287217/1‑93 (MQ=255)
                                                                   aCCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTggcaggca                      <  1:3981392/93‑1 (MQ=255)
                                                                   aCCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTggcaggca                      >  2:1299748/1‑93 (MQ=255)
                                                                       tGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCt                  >  1:1461195/1‑93 (MQ=255)
                                                                       tGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCt                  <  2:587639/93‑1 (MQ=255)
                                                                       tGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCt                  >  1:3529883/1‑93 (MQ=255)
                                                                        gagaTCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTg                 >  1:4309564/1‑93 (MQ=255)
                                                                         agaTCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTgg                <  2:4510295/93‑1 (MQ=255)
                                                                          gaTCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTggg               <  1:2977335/93‑1 (MQ=255)
                                                                           aTCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGa              >  1:7625475/1‑93 (MQ=255)
                                                                           aTCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGATATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGa              <  2:5410528/93‑1 (MQ=255)
                                                                              gtcgATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTg           <  1:4108893/93‑1 (MQ=255)
                                                                              gtcgATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTg           >  1:2892534/1‑93 (MQ=255)
                                                                              gtcgATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTg           >  1:6815936/1‑93 (MQ=255)
                                                                               tcgATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTgg          <  1:3811242/93‑1 (MQ=255)
                                                                                 gATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGc        >  1:1793821/1‑93 (MQ=255)
                                                                                  aTCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCa       <  1:942963/93‑1 (MQ=255)
                                                                                  aTCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCa       >  2:1919076/1‑93 (MQ=255)
                                                                                  aTCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCa       >  2:2380015/1‑93 (MQ=255)
                                                                                  aTCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCa       >  1:7866262/1‑93 (MQ=255)
                                                                                    cACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAg     >  2:4085701/1‑93 (MQ=255)
                                                                                     aCTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGc    >  1:7952572/1‑93 (MQ=255)
                                                                                      cTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGCt   <  2:3577446/93‑1 (MQ=255)
                                                                                       ttCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGCTg  <  2:499396/93‑1 (MQ=255)
                                                                                       ttCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGCTg  >  1:1696786/1‑93 (MQ=255)
                                                                                       ttCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGCTg  >  1:6639445/1‑93 (MQ=255)
                                                                                       ttCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATCTCCCACGAAGCCTGTTCGCTGGCAGACACGCTGGGACTGGGCAAGCTg  >  1:2486728/1‑93 (MQ=37)
                                                                                          |                                                                                         
GGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGCTG  >  NC_000913/2578007‑2578186

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: