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,8420CT99.8% 1530.3 / ‑5.8 414intergenic (+48/‑78)rbsB/rbsKD‑ribose periplasmic binding protein/ribokinase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (191/222);  total (192/222)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

Only 100 of 420 total aligned reads displayed.

GGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCA  >  NC_000913/3934752‑3934931
                                                                                          |                                                                                         
ggCTAAGTATCCGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAg                                                                                         >  2:5202377/1‑93 (MQ=255)
ggCTAAGTATCCGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAg                                                                                         >  2:7347361/1‑93 (MQ=255)
 gCTAAGTATCCGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGa                                                                                        <  1:1975534/93‑1 (MQ=255)
  cTAAGTATCCGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGAt                                                                                       <  2:2920306/93‑1 (MQ=255)
   tAAGTATCCGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATg                                                                                      <  2:1858443/93‑1 (MQ=255)
   tAAGTATCCGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATg                                                                                      >  1:1260007/1‑93 (MQ=255)
    aaGTATCCGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGa                                                                                     >  1:7078756/1‑93 (MQ=255)
    aaGTATCCGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGa                                                                                     <  2:5262527/93‑1 (MQ=255)
     aGTATCCGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATgga                                                                                    <  2:4375759/93‑3 (MQ=255)
     aGTATCCGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGaa                                                                                    >  2:5488042/1‑93 (MQ=255)
      gTATCCGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAAc                                                                                   <  1:6930351/93‑1 (MQ=255)
      gTATCCGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAAc                                                                                   >  2:7181320/1‑93 (MQ=255)
      gTATCCGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAAc                                                                                   <  2:2324550/93‑1 (MQ=255)
      gTATCCGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAAc                                                                                   <  2:6863719/93‑1 (MQ=255)
       tATCCGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACg                                                                                  >  2:6552399/1‑93 (MQ=255)
        aTCCGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACgt                                                                                 <  2:150172/93‑1 (MQ=255)
         tCCGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACgtg                                                                                <  2:3129255/93‑1 (MQ=255)
          ccGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACgtgt                                                                               <  2:1608301/93‑1 (MQ=255)
           cGGTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTa                                                                              >  1:1763849/1‑93 (MQ=255)
             gTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATCAATACGTCATCGATAGATGAACGTGTAat                                                                            >  2:179362/1‑93 (MQ=255)
             gTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAat                                                                            <  1:1005116/93‑1 (MQ=255)
             gTTGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAat                                                                            >  1:7202181/1‑93 (MQ=255)
              ttGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAata                                                                           <  2:1127476/93‑1 (MQ=255)
              ttGATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAata                                                                           <  2:4842264/93‑1 (MQ=255)
                gATCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAAtata                                                                         <  2:6176864/93‑1 (MQ=255)
                 aTCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATaa                                                                        <  2:7533581/93‑1 (MQ=255)
                 aTCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATaa                                                                        >  1:5703129/1‑93 (MQ=255)
                  tCTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATaaa                                                                       >  1:6369560/1‑93 (MQ=255)
                   cTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAg                                                                      >  2:2443724/1‑93 (MQ=255)
                   cTGGAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAg                                                                      <  1:86208/93‑1 (MQ=255)
                    tggaaCTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGa                                                                     >  1:4185355/4‑93 (MQ=255)
                    tggaaCTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGa                                                                     >  2:5037199/4‑93 (MQ=255)
                    tggaaCTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGa                                                                     <  2:2241552/90‑1 (MQ=255)
                    tggaaCTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGa                                                                     <  2:4334465/90‑1 (MQ=255)
                      gatcTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGaaa                                                                   <  1:5930351/90‑1 (MQ=255)
                      gaaCTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGaaa                                                                   >  1:7765254/2‑93 (MQ=255)
                      gaaCTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGaaa                                                                   <  1:2926984/92‑1 (MQ=255)
                       aaCTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGaaaa                                                                  <  1:4253649/93‑1 (MQ=255)
                       aaCTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGaaaa                                                                  <  2:4791180/93‑1 (MQ=255)
                       aaCTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGaaaa                                                                  <  1:2777134/93‑1 (MQ=255)
                       aaCTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGaaaa                                                                  <  1:2371763/93‑1 (MQ=255)
                       aaCTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGaaaa                                                                  <  1:2304915/93‑1 (MQ=255)
                          tGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCa                                                               <  2:422440/93‑1 (MQ=255)
                          tGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCa                                                               <  1:7247707/93‑1 (MQ=255)
                          tGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCa                                                               >  1:6750612/1‑93 (MQ=255)
                          tGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCa                                                               >  1:6137746/1‑93 (MQ=255)
                          tGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCa                                                               >  1:5338787/1‑93 (MQ=255)
                           ggTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAg                                                              <  2:7698355/93‑1 (MQ=255)
                             ttgttAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAggg                                                            <  1:7651464/93‑1 (MQ=255)
                             ttgttAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAggg                                                            >  2:3555467/1‑93 (MQ=255)
                             ttgttAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAggg                                                            <  1:1069802/93‑1 (MQ=255)
                             ttgttAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAACATAAAGAAAAGCAggg                                                            <  1:3402677/93‑1 (MQ=255)
                             ttgttAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAggg                                                            >  1:1404096/1‑93 (MQ=255)
                               gttAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCa                                                          >  1:4745301/1‑93 (MQ=255)
                                ttAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCAc                                                         <  2:7858097/93‑1 (MQ=255)
                                 tAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCAcg                                                        <  2:5580059/93‑1 (MQ=255)
                                 tAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCAcg                                                        <  2:3806557/93‑1 (MQ=255)
                                  aaGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCAcgc                                                       <  1:3777110/93‑1 (MQ=255)
                                   aGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCAcgcg                                                      >  1:4994272/1‑93 (MQ=255)
                                   aGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCAcgcc                                                      >  1:3915055/1‑92 (MQ=255)
                                    gCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACgcgc                                                     >  1:3600011/1‑93 (MQ=255)
                                    gCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACgcgc                                                     >  2:5831026/1‑93 (MQ=255)
                                         agtTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCAccc                                                <  2:5731518/93‑1 (MQ=255)
                                          gtTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCt                                               >  1:620539/1‑93 (MQ=255)
                                           ttttAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTa                                              <  1:5048751/93‑1 (MQ=255)
                                            tttAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTaa                                             >  2:1008232/1‑93 (MQ=255)
                                            tttAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCCaa                                             >  2:3264382/1‑93 (MQ=255)
                                                 tCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACgg                                        <  1:4477758/93‑1 (MQ=255)
                                                  cAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGt                                       <  2:3957810/93‑1 (MQ=255)
                                                    ggTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTgg                                     <  2:858557/93‑1 (MQ=255)
                                                     gTTGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGACCGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGgc                                    >  1:5414870/1‑93 (MQ=255)
                                                      ttGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTTACACGGTGgcg                                   >  2:4153303/1‑93 (MQ=255)
                                                      ttGTATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTCACACGGGGgcc                                   >  2:2052104/1‑92 (MQ=255)
                                                         tATGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGGGGCGCAt                                >  1:7283808/1‑93 (MQ=255)
                                                          aTGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCAtt                               <  2:6421511/93‑1 (MQ=255)
                                                           tGACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCAttt                              <  2:749041/93‑1 (MQ=255)
                                                            gACCTGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCAtttt                             >  1:6555580/1‑93 (MQ=255)
                                                               cttaTGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATg                          <  2:6795679/90‑1 (MQ=255)
                                                                tGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGGCACCCTAACACGGGGGGGGCTTTTTTgg                         >  1:1956532/1‑93 (MQ=255)
                                                                tGATGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATgg                         >  2:5926544/1‑93 (MQ=255)
                                                                  aTGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGAc                       <  1:3569543/93‑1 (MQ=255)
                                                                   tGGTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACa                      <  1:3227258/93‑1 (MQ=255)
                                                                     gTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATc                    <  2:6349495/93‑1 (MQ=255)
                                                                     gTGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATc                    <  1:211743/93‑1 (MQ=255)
                                                                      tGACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATcc                   <  1:2510293/93‑1 (MQ=255)
                                                                       gACATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATccc                  <  2:1346307/93‑1 (MQ=255)
                                                                         cATAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCCTTTTTTTGACATCCCGa                >  2:6217308/1‑93 (MQ=255)
                                                                          aTAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGaa               >  2:2146172/1‑93 (MQ=255)
                                                                           tAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAat              >  1:7435393/1‑93 (MQ=255)
                                                                           tAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAat              <  1:733959/93‑1 (MQ=255)
                                                                           tAAATACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAat              >  1:3512245/1‑93 (MQ=255)
                                                                            aaaTACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAata             <  1:5507335/93‑1 (MQ=255)
                                                                            aaaTACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAata             <  2:2667763/93‑1 (MQ=255)
                                                                             aaTACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAatat            <  2:2525824/93‑1 (MQ=255)
                                                                               tACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGc          >  1:4055856/1‑93 (MQ=255)
                                                                               tACGTCATCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGc          >  2:4015071/1‑93 (MQ=255)
                                                                                     aTCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGGTCCTCCCTAACACGGTGGCTCATTTTATGGACATCCCGAATATGCAAAACg    <  2:676243/93‑1 (MQ=255)
                                                                                      tCGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTCACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACgc   >  2:3459722/1‑93 (MQ=255)
                                                                                       cGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACgca  <  1:6039550/93‑1 (MQ=255)
                                                                                       cGATAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACgca  <  2:5156701/93‑1 (MQ=255)
                                                                                          |                                                                                         
GGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCA  >  NC_000913/3934752‑3934931

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