Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,091,351 C→T G168D (GGC→GAC)  rhaD ← rhamnulose‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,091,3510CT100.0% 1326.5 / NA 368G168D (GGC→GAC) rhaDrhamnulose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (157/211);  total (157/211)

Only 100 of 386 total aligned reads displayed.

CGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGCCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATC  >  NC_000913/4091262‑4091440
                                                                                         |                                                                                         
 ggTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTc                                                                                       >  1:903730/1‑93 (MQ=255)
      ccTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCcaca                                                                                  <  1:2730515/93‑1 (MQ=255)
       cTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAg                                                                                 >  2:5515920/1‑93 (MQ=255)
       cTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAg                                                                                 >  2:2035378/1‑93 (MQ=255)
        tGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGt                                                                                >  1:5254756/1‑93 (MQ=255)
        tGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGt                                                                                >  2:1978722/1‑93 (MQ=255)
         ggCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGtt                                                                               <  2:7825008/93‑1 (MQ=255)
         ggCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGtt                                                                               >  1:7839127/1‑93 (MQ=255)
         ggCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGtt                                                                               >  2:1752407/1‑93 (MQ=255)
         ggCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGtt                                                                               >  1:2981192/1‑93 (MQ=255)
         ggCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGtt                                                                               <  2:3396118/93‑1 (MQ=255)
          gCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCATAGTTg                                                                              >  2:7746428/1‑93 (MQ=255)
          gCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCATAGTTg                                                                              >  2:4271683/1‑93 (MQ=255)
           ccGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTgg                                                                             <  2:3081951/93‑1 (MQ=255)
            cGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGc                                                                            <  1:6256373/93‑1 (MQ=255)
            cGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGc                                                                            >  1:6600868/1‑93 (MQ=255)
             gATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCATAGTTGGCg                                                                           >  1:6947875/1‑93 (MQ=255)
             gATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCg                                                                           >  1:4898675/1‑93 (MQ=255)
              aTTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGa                                                                          <  2:6444710/93‑1 (MQ=255)
              aTTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGa                                                                          <  1:7666966/93‑1 (MQ=255)
              aTTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGa                                                                          >  1:3857422/1‑93 (MQ=255)
              aTTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGa                                                                          <  2:3106419/93‑1 (MQ=255)
              aTTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGa                                                                          <  2:3200573/93‑1 (MQ=255)
               tttCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAg                                                                         >  2:3271682/1‑93 (MQ=255)
                ttCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAgg                                                                        >  2:1482566/1‑92 (MQ=255)
                 tCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTg                                                                       <  1:5319357/93‑1 (MQ=255)
                  cGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGa                                                                      <  1:3995188/93‑1 (MQ=255)
                    tCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAg                                                                    >  2:378219/1‑93 (MQ=255)
                    tCCGTGCCGGGCAACATACACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGCTGGCGAGTGAAg                                                                    >  2:7551583/1‑93 (MQ=255)
                    gccGTGCCGGGCACCATCCACGGCAAAATGCCAGCGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAg                                                                    <  1:4287401/92‑1 (MQ=255)
                     ccGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGa                                                                   <  1:968397/93‑1 (MQ=255)
                      cGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGAc                                                                  <  2:7311092/93‑1 (MQ=255)
                      cGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGAc                                                                  <  2:5614206/93‑1 (MQ=255)
                       gTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGAcc                                                                 >  2:1361601/1‑93 (MQ=255)
                           cGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCgg                                                             <  2:683834/93‑1 (MQ=255)
                           cGGCCACCATCCACGGCAAAATGCCACCGCCATCCGGGAATACCACCAGACACGCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCgg                                                             <  1:1569237/93‑1 (MQ=255)
                              tcACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACTCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGgtgt                                                          <  1:3442046/92‑1 (MQ=255)
                              gCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGgtgt                                                          >  1:1222139/1‑93 (MQ=255)
                              gCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGgtgt                                                          <  2:1069612/93‑1 (MQ=255)
                              gCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGgtgt                                                          <  1:7414399/93‑1 (MQ=255)
                               cACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCATAGTTGGCGAGTGAAGACCGCGGTGTc                                                         >  1:5422127/1‑93 (MQ=255)
                                aCCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCg                                                        >  2:5303739/1‑93 (MQ=255)
                                 ccATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGt                                                       <  2:2714307/93‑1 (MQ=255)
                                 ccATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGt                                                       >  2:5098388/1‑93 (MQ=255)
                                  cATCCGCGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGtt                                                      <  1:2864583/93‑1 (MQ=255)
                                  cATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGtt                                                      <  2:7392028/93‑1 (MQ=255)
                                  cATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGtt                                                      <  2:706430/93‑1 (MQ=255)
                                  cATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGtt                                                      <  2:1667192/93‑1 (MQ=255)
                                   aTCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGttt                                                     <  2:6231261/93‑1 (MQ=255)
                                    tCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGtttt                                                    >  2:2449692/1‑93 (MQ=255)
                                       aCGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCaa                                                 <  2:4580413/93‑1 (MQ=255)
                                       aCGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCaa                                                 <  1:471562/93‑1 (MQ=255)
                                         ggCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGt                                               >  2:3520044/1‑93 (MQ=255)
                                          gCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTa                                              <  2:2985216/93‑1 (MQ=255)
                                          gCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTa                                              <  2:3710809/93‑1 (MQ=255)
                                           cAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTAc                                             <  1:6307240/93‑1 (MQ=255)
                                            aaaaTGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACa                                            >  1:220078/1‑93 (MQ=255)
                                            aaaaTGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACa                                            >  1:4137303/1‑93 (MQ=255)
                                             aaaTGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACAt                                           >  2:1927163/1‑93 (MQ=255)
                                              aaTGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATa                                          >  2:3996287/1‑93 (MQ=255)
                                               aTGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAg                                         <  2:1436491/93‑1 (MQ=255)
                                                 gCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGt                                       <  2:6722803/93‑1 (MQ=255)
                                                   cAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGa                                     <  1:3266634/93‑1 (MQ=255)
                                                   cAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGa                                     <  1:7494881/93‑1 (MQ=255)
                                                    aaCGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAg                                    >  1:5595816/1‑93 (MQ=255)
                                                    aaCGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAg                                    <  1:2121803/93‑1 (MQ=255)
                                                    aaCGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAg                                    <  1:2424714/93‑1 (MQ=255)
                                                    aaCGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAg                                    <  2:3870358/93‑1 (MQ=255)
                                                      cGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTAGTTTTCAAGTACATAGGTGAggg                                  <  1:3687845/93‑1 (MQ=255)
                                                         cATCCGGGAATACAACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGa                               <  2:762171/93‑1 (MQ=255)
                                                          aTCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGAt                              >  1:6465120/1‑93 (MQ=255)
                                                          aTCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGAt                              <  1:5570807/93‑1 (MQ=255)
                                                          aTCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGAt                              >  1:23613/1‑93 (MQ=255)
                                                           accGGGAATACCCCCAGACCCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATc                             <  1:6741344/92‑1 (MQ=255)
                                                            ccGGGAATCCCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCa                            <  2:2243561/93‑1 (MQ=255)
                                                            ccGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCa                            <  1:7232986/93‑1 (MQ=255)
                                                            ccGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCa                            <  1:7114808/93‑1 (MQ=255)
                                                            ccGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCa                            <  1:5683131/93‑1 (MQ=255)
                                                              gggAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAgg                          <  1:6057216/93‑1 (MQ=255)
                                                              gggAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAgg                          >  1:6663508/1‑93 (MQ=255)
                                                               ggAATACCCCCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGt                         <  2:5762866/93‑1 (MQ=255)
                                                                gAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGtt                        <  2:2664960/93‑1 (MQ=255)
                                                                 aagaccaccAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtg                       <  1:1856717/90‑1 (MQ=255)
                                                                  aTACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtgg                      <  2:4431931/93‑1 (MQ=255)
                                                                   tcccaccaGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtggt                     <  1:5044944/91‑1 (MQ=255)
                                                                   tACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtggt                     >  2:602661/1‑93 (MQ=255)
                                                                     ccaccaGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtggtgg                   <  1:1418457/93‑1 (MQ=255)
                                                                        ccaGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcgt                <  1:1000682/93‑1 (MQ=255)
                                                                         caGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCGGGTTGgtggcgtg               <  2:2520453/93‑1 (MQ=255)
                                                                            acacTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGa                                     <  2:4818917/68‑1 (MQ=255)
                                                                            acacTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGa                                     >  1:4818918/1‑68 (MQ=255)
                                                                            acacTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCa            <  1:7570016/93‑1 (MQ=255)
                                                                              acTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGt          <  2:130058/93‑1 (MQ=255)
                                                                               cTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTg         >  1:5781769/1‑93 (MQ=255)
                                                                                tCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGc        >  2:4498134/1‑93 (MQ=255)
                                                                                 cGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCa       >  1:5854806/1‑93 (MQ=255)
                                                                                   gTGCTGTCTTCCCATAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATg     >  1:5993732/1‑93 (MQ=255)
                                                                                    tgctgtCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGa    <  1:2548877/93‑1 (MQ=255)
                                                                                     gctgtCTTCCCATAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGAt   <  1:5183713/93‑1 (MQ=255)
                                                                                      ctgtCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATc  >  2:5870139/1‑93 (MQ=255)
                                                                                         |                                                                                         
CGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGCCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATC  >  NC_000913/4091262‑4091440

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