Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,934,842 C→T intergenic (+48/‑78) rbsB → / → rbsK D‑ribose periplasmic binding protein/ribokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,934,8420CT100.0% 832.1 / NA 227intergenic (+48/‑78)rbsB/rbsKD‑ribose periplasmic binding protein/ribokinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (107/120);  total (107/120)

Only 100 of 231 total aligned reads displayed.

GGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCA  >  NC_000913/3934752‑3934931
                                                                                          |                                                                                         
ggCTAAGTATCCGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAg                                                                                         <  2:1489105/93‑1 (MQ=255)
 gCTAAGTATCCGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGa                                                                                        <  2:84759/93‑1 (MQ=255)
 gCTAAGTATCCGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGa                                                                                        <  2:82494/93‑1 (MQ=255)
 gCTAAGTATCCGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGa                                                                                        <  2:3216495/93‑1 (MQ=255)
  cTAAGTATCCGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGAt                                                                                       <  1:1485182/93‑1 (MQ=255)
  cTAAGTATCCGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGAt                                                                                       <  2:1964178/93‑1 (MQ=255)
   tAAGTATCCGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATg                                                                                      <  2:4303568/93‑1 (MQ=255)
   tAAGCATCCGGATGATCTGGAACTGGTTGTTAAGAAATAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATg                                                                                      >  2:1853530/1‑93 (MQ=255)
     aGTATCCGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGaa                                                                                    >  2:3395663/1‑93 (MQ=255)
     aGTATCCGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGaa                                                                                    <  2:444184/93‑1 (MQ=255)
     aGTATCCGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGaa                                                                                    >  1:2132172/1‑93 (MQ=255)
     aGTATCCGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGaa                                                                                    <  2:729499/93‑1 (MQ=255)
      gTATCCGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAAc                                                                                   <  2:1220495/93‑1 (MQ=255)
      gTATCCGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAAc                                                                                   <  2:3943541/93‑1 (MQ=255)
       tATCCGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACg                                                                                  <  2:466784/93‑1 (MQ=255)
        aTCCGGTTGATCTGGGACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACgt                                                                                 >  1:3658881/1‑93 (MQ=255)
        aTCCGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACgt                                                                                 <  1:3021716/93‑1 (MQ=255)
         tCCGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACgtg                                                                                <  2:2737953/93‑1 (MQ=255)
          ccGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACgtgt                                                                               <  2:2775773/93‑1 (MQ=255)
          ccGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACgtgt                                                                               <  1:2727385/93‑1 (MQ=255)
           cGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTa                                                                              >  1:3930114/1‑93 (MQ=255)
            ggTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTaa                                                                             >  2:476698/1‑93 (MQ=255)
             gTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAat                                                                            >  1:4597957/1‑93 (MQ=255)
             gTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAat                                                                            >  2:4203824/1‑93 (MQ=255)
              ttGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAata                                                                           <  2:2395171/93‑1 (MQ=255)
              ttGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAata                                                                           <  1:1208860/93‑1 (MQ=255)
                gATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAAtata                                                                         <  2:3609622/93‑1 (MQ=255)
                 aTCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATaa                                                                        >  2:540434/1‑93 (MQ=255)
                 aTCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATaa                                                                        >  2:3022227/1‑93 (MQ=255)
                  tCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATaaa                                                                       >  1:1329757/1‑93 (MQ=255)
                      gagcTGGTTGTTAAGCAGTAGTTTTAATCAGTTCGTATGCCCTGAGGATGACATAAATACGTGTTCGATAGCTGAACGTGTAATATAAAGaaa                                                                   <  1:3429347/90‑1 (MQ=255)
                       aaCTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGaaaa                                                                  <  2:3169327/93‑1 (MQ=255)
                       aaCTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGaaaa                                                                  <  1:267191/93‑1 (MQ=255)
                       aaCTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCAACGATAGATGAACGTGTAATATAAAGaaaa                                                                  >  1:1788205/1‑93 (MQ=255)
                        aCTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAg                                                                 <  1:3247556/93‑1 (MQ=255)
                         cTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGc                                                                >  1:1310295/1‑93 (MQ=255)
                         cTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGc                                                                <  1:495991/93‑1 (MQ=255)
                          tGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCa                                                               >  1:1897578/1‑93 (MQ=255)
                           ggTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAg                                                              <  2:3896849/93‑1 (MQ=255)
                           ggTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAg                                                              <  2:3785169/93‑1 (MQ=255)
                           ggTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAg                                                              >  2:3202425/1‑93 (MQ=255)
                           ggTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGATGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAg                                                              <  2:2784635/93‑1 (MQ=255)
                            gttgttAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAgg                                                             >  1:987562/1‑93 (MQ=255)
                            gttgttAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAgg                                                             >  2:632832/1‑93 (MQ=255)
                            gttgttAAGCAGTAGTTTTAATCAGGTTATATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAgg                                                             <  1:3835970/93‑1 (MQ=255)
                             ttgttAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAggg                                                            >  2:3463031/1‑93 (MQ=255)
                                ttAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCAc                                                         <  2:1930204/93‑1 (MQ=255)
                                  aaGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCAcgc                                                       >  1:2070968/1‑93 (MQ=255)
                                  aaGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCAcgc                                                       <  2:1834343/93‑1 (MQ=255)
                                      agtagtTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCa                                                   >  1:4657096/1‑93 (MQ=255)
                                      agtagtTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCa                                                   <  1:1883183/93‑1 (MQ=255)
                                       gtagtTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATTAACGTGTGATATAAAGAAAAGAAGGGGAAGCGacac                                                  >  2:1997626/1‑93 (MQ=255)
                                        tagtTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCcccc                                                 >  1:4099604/1‑93 (MQ=255)
                                            tttAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTaa                                             <  1:1130099/93‑1 (MQ=255)
                                            tttAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTaa                                             <  1:649034/93‑1 (MQ=255)
                                             ttAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAac                                            <  1:2500768/93‑1 (MQ=255)
                                              tAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAaca                                           <  2:4104321/93‑1 (MQ=255)
                                              tAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAaca                                           <  2:2882457/93‑1 (MQ=255)
                                                 tCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACgg                                        <  1:3185724/93‑1 (MQ=255)
                                                 tCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACgg                                        <  2:1120353/93‑1 (MQ=255)
                                                  cAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGt                                       <  1:1189337/93‑1 (MQ=255)
                                                     gTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGgc                                    >  1:591159/1‑93 (MQ=255)
                                                        gTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCa                                 >  1:725187/1‑93 (MQ=255)
                                                         tATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCAt                                >  1:3124931/1‑93 (MQ=255)
                                                          aTGACCTGATGTTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCAtt                               <  2:2585386/93‑1 (MQ=255)
                                                          aTGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCAtt                               <  2:1309379/93‑1 (MQ=255)
                                                            gACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCAtttt                             <  2:425506/93‑1 (MQ=255)
                                                             aCCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCAttttt                            >  2:3614658/1‑92 (MQ=255)
                                                             aCCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTa                            <  1:2454475/93‑1 (MQ=255)
                                                             aCCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTa                            >  1:4541288/1‑93 (MQ=255)
                                                             aCCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTa                            <  2:3056237/93‑1 (MQ=255)
                                                             aCCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTa                            >  1:1563862/1‑93 (MQ=255)
                                                             aCCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTa                            <  2:2059279/93‑1 (MQ=255)
                                                              ccTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGGGGCGCATTTTAt                           >  1:929292/1‑93 (MQ=255)
                                                               cTGATGGTGAAATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAATGCGGGGCGCGCGCCCCCCTAACACGGTGGCGCATTTTATg                          <  2:156761/93‑1 (MQ=255)
                                                                 gATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGa                        <  2:2171950/93‑1 (MQ=255)
                                                                 gATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGa                        >  2:1604457/1‑93 (MQ=255)
                                                                   tGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACa                      <  1:1049353/93‑1 (MQ=255)
                                                                    ggTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACAt                     >  2:4631891/1‑93 (MQ=255)
                                                                    ggTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACAt                     <  2:2709913/93‑1 (MQ=255)
                                                                      tGACATAAATACGTCATCGATAGATGAAAGTGGAATATAAAGAAAAGCAGGGCACGGGCCAACCTAACACGGTGGCGCATTTTATGGACATcc                   >  2:414295/1‑93 (MQ=255)
                                                                       gACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATccc                  <  1:2348646/93‑1 (MQ=255)
                                                                       gACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATccc                  <  1:3282697/93‑1 (MQ=255)
                                                                        aCATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCg                 <  2:4527914/93‑1 (MQ=255)
                                                                         cATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGGCACCCTCACACGGTGGGGCATTTTTTGGTCATcccca                >  2:981584/1‑91 (MQ=255)
                                                                         cATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCCAAAACCGTGGCGCATTTTTTGGACATcccca                >  2:3431740/1‑91 (MQ=255)
                                                                         cATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCCACCCCCTAACACCGTTGGCCCTTTTATTGGCATcccta                >  2:1032289/1‑91 (MQ=255)
                                                                          aTAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGACACCCTAACACGGTGGCGCATTTTATGGACATCCCGaa               >  1:3800118/1‑93 (MQ=255)
                                                                           tAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAat              >  1:3492636/1‑93 (MQ=255)
                                                                           tAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAat              <  2:1115433/93‑1 (MQ=255)
                                                                           tAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGAGCATTTTATGGACATCCCGAat              >  2:2244873/1‑93 (MQ=255)
                                                                           tAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGCGCACGCGCCTCCCTAACACGGTGTCGCATTTTATGGACATCCCGAat              <  1:3588169/93‑1 (MQ=255)
                                                                                 cGTCATCGATAGATGAACGTGTAATATAGAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCaa        <  2:2505438/93‑1 (MQ=255)
                                                                                  gTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCaaa       <  2:3663932/93‑1 (MQ=255)
                                                                                   tCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCaaaa      <  1:3391337/93‑1 (MQ=255)
                                                                                   tCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCaaaa      <  2:3241215/93‑1 (MQ=255)
                                                                                   tCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCaaaa      <  1:197217/93‑1 (MQ=255)
                                                                                   tCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCaaaa      <  1:1871874/93‑1 (MQ=255)
                                                                                    cATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGTGCCACCCTAACACGGGGGGGCATTTTATGGACATCCCGAACATGCaaagc     >  1:1273592/1‑91 (MQ=255)
                                                                                       cGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACgca  <  1:2188132/93‑1 (MQ=255)
                                                                                          |                                                                                         
GGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCA  >  NC_000913/3934752‑3934931

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