Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,077,905 A→G S584P (TCG→CCG)  tktA ← transketolase 1 isozyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,077,9050AG100.0% 1143.8 / NA 333S584P (TCG→CCG) tktAtransketolase 1 isozyme
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (161/172);  total (161/172)

Only 100 of 339 total aligned reads displayed.

CCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGACATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGA  >  NC_000913/3077815‑3077994
                                                                                          |                                                                                         
ccGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCa                                                                                         <  2:3313548/93‑1 (MQ=255)
 cGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCAt                                                                                        <  1:4704312/93‑1 (MQ=255)
 cGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCAt                                                                                        <  2:5818076/93‑1 (MQ=255)
  gCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATg                                                                                       <  1:3830560/93‑1 (MQ=255)
  gCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATg                                                                                       <  1:6438346/93‑1 (MQ=255)
  gCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATg                                                                                       >  2:155296/1‑93 (MQ=255)
  gCTTCTACAGCAACGCGTCCAGTAACCGCTTTCGGCAGTACAGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATg                                                                                       <  2:5293972/93‑1 (MQ=255)
   cTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATgg                                                                                      >  2:2474958/1‑93 (MQ=255)
   cTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATgg                                                                                      >  2:1185070/1‑93 (MQ=255)
    ttCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGa                                                                                     <  2:4919252/93‑1 (MQ=255)
     tCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGac                                                                                    >  1:3723759/1‑93 (MQ=255)
     tCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGac                                                                                    >  1:7429736/1‑93 (MQ=255)
         cAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGAcacca                                                                                <  1:2304394/93‑1 (MQ=255)
          aGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGAcaccac                                                                               <  1:1753887/93‑1 (MQ=255)
              aCGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCAcgcgc                                                                           >  1:7461697/1‑93 (MQ=255)
                gcgTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACgcgcgc                                                                         >  2:6109863/1‑93 (MQ=255)
                 cgTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCt                                                                        <  1:6067852/93‑1 (MQ=255)
                   tGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCttt                                                                      >  1:6478688/1‑93 (MQ=255)
                    gCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTc                                                                     <  2:882537/93‑1 (MQ=255)
                    gCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTc                                                                     >  1:1932735/1‑93 (MQ=255)
                    gCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTc                                                                     >  2:6009092/1‑93 (MQ=255)
                     cAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCa                                                                    <  1:3659328/93‑1 (MQ=255)
                     cAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCa                                                                    <  2:830691/93‑1 (MQ=255)
                      aGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCAc                                                                   >  1:423049/1‑93 (MQ=255)
                      aGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCAc                                                                   <  2:3803492/93‑1 (MQ=255)
                           ccGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCtt                                                              <  1:2476182/93‑1 (MQ=255)
                             gCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCg                                                            <  1:4097968/93‑1 (MQ=255)
                             gCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCg                                                            >  2:7428004/1‑93 (MQ=255)
                             gCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCg                                                            <  2:2734563/93‑1 (MQ=255)
                               tttCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGc                                                          >  1:4454402/1‑93 (MQ=255)
                               tttCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGc                                                          <  2:4726579/93‑1 (MQ=255)
                                ttCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGca                                                         >  1:7320954/1‑93 (MQ=255)
                                ttCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGca                                                         <  2:1850061/93‑1 (MQ=255)
                                 tCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGcag                                                        <  1:4930360/93‑1 (MQ=255)
                                   ggCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGcagtc                                                      <  2:3407594/93‑1 (MQ=255)
                                     cAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGcagtcag                                                    <  2:2873648/93‑1 (MQ=255)
                                     cAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGcagtcag                                                    <  2:3683152/93‑1 (MQ=255)
                                      aGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGcagtcagt                                                   <  2:6614299/93‑1 (MQ=255)
                                       gTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGtt                                                  <  2:1283714/93‑1 (MQ=255)
                                       gTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGtt                                                  >  1:138478/1‑93 (MQ=255)
                                        tACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGttt                                                 <  2:6429403/93‑1 (MQ=255)
                                        tACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGttt                                                 >  2:2545009/1‑93 (MQ=255)
                                        tACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGttt                                                 >  1:1016077/1‑93 (MQ=255)
                                          cGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGttttt                                               <  2:3225675/93‑1 (MQ=255)
                                           ggATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCAGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTc                                              <  2:4638702/93‑1 (MQ=255)
                                            gATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCg                                             >  2:4006543/1‑93 (MQ=255)
                                              ttCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTa                                           <  2:439012/93‑1 (MQ=255)
                                              ttCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTa                                           <  1:3253354/93‑1 (MQ=255)
                                              ttCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTa                                           <  2:4272751/93‑1 (MQ=255)
                                              ttCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTa                                           <  2:4779807/93‑1 (MQ=255)
                                              ttCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTa                                           >  2:5044958/1‑93 (MQ=255)
                                               tCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAg                                          <  1:1665751/93‑1 (MQ=255)
                                               tCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAg                                          <  1:1500019/93‑1 (MQ=255)
                                                 aCGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGgc                                        >  2:2923968/1‑93 (MQ=255)
                                                 aCGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGgc                                        <  1:5036216/93‑1 (MQ=255)
                                                   ggTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGgcag                                      >  2:7034770/1‑93 (MQ=255)
                                                      aaGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCaa                                   <  1:998110/93‑1 (MQ=255)
                                                       agcagcATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAAc                                  >  1:5243862/1‑93 (MQ=255)
                                                        gcagcaTCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACa                                 <  2:1103643/93‑1 (MQ=255)
                                                         cagcaTCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAg                                <  1:112989/93‑1 (MQ=255)
                                                         cagcaTCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAg                                <  1:5415009/93‑1 (MQ=255)
                                                          agcaTCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGc                               >  2:2415252/1‑93 (MQ=255)
                                                          agcaTCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGc                               <  1:5986721/93‑1 (MQ=255)
                                                          agcaTCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGc                               >  2:6783981/1‑93 (MQ=255)
                                                           gcaTCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGcc                              <  1:2055823/93‑1 (MQ=255)
                                                            caTCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCa                             <  2:2093099/93‑1 (MQ=255)
                                                            caTCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCa                             <  1:2785128/93‑1 (MQ=255)
                                                             aTCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAg                            >  2:189554/1‑93 (MQ=255)
                                                             aTCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAg                            <  2:430392/93‑1 (MQ=255)
                                                             aTCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAg                            <  1:5740113/93‑1 (MQ=255)
                                                              tCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGt                           >  1:5820043/1‑93 (MQ=255)
                                                              tCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGt                           >  2:1991301/1‑93 (MQ=255)
                                                              tCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGCAGGCAGCAACAGCCAGt                           >  2:1478653/1‑93 (MQ=255)
                                                               ccTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGtt                          <  1:2982260/93‑1 (MQ=255)
                                                                cTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTc                         <  2:6337815/93‑1 (MQ=255)
                                                                  gCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCaa                       >  2:1820573/1‑93 (MQ=255)
                                                                    ttGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACt                     >  1:7520736/1‑93 (MQ=255)
                                                                      gTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTc                   <  2:4880686/93‑1 (MQ=255)
                                                                       tCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCt                  <  2:3562293/93‑1 (MQ=255)
                                                                         aaaTGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGa                >  2:7202571/1‑93 (MQ=255)
                                                                          aaTGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCGCGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGaa               >  1:6223756/1‑93 (MQ=255)
                                                                          aaTGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGaa               >  1:7032609/1‑93 (MQ=255)
                                                                           aTGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAAc              <  2:1507867/93‑1 (MQ=255)
                                                                           aTGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAAc              >  2:6690644/1‑93 (MQ=255)
                                                                           aTGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAAc              <  2:7253940/93‑1 (MQ=255)
                                                                            tGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAAcc             <  1:6908931/93‑1 (MQ=255)
                                                                              cGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCgg           <  1:4044302/93‑1 (MQ=255)
                                                                              cGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCgg           <  1:417515/93‑1 (MQ=255)
                                                                               gTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGt          >  1:7122229/1‑93 (MQ=255)
                                                                                tCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTa         <  1:2220957/93‑1 (MQ=255)
                                                                                tCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTa         <  1:1947938/93‑1 (MQ=255)
                                                                                tCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTa         <  1:2499322/93‑1 (MQ=255)
                                                                                tCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTa         <  1:2875580/93‑1 (MQ=255)
                                                                                tCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTa         >  1:6636567/1‑93 (MQ=255)
                                                                                    tAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGa     <  2:5458426/93‑1 (MQ=255)
                                                                                    tAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGa     >  1:3531140/1‑93 (MQ=255)
                                                                                    tAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGa     <  2:6885623/93‑1 (MQ=255)
                                                                                    tAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGa     >  2:3113790/1‑93 (MQ=255)
                                                                                      gacggcATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATg   >  2:3376034/1‑93 (MQ=255)
                                                                                      gacggcATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATg   >  1:6751939/1‑93 (MQ=255)
                                                                                          |                                                                                         
CCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGACATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGA  >  NC_000913/3077815‑3077994

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