Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,339,313 A→C 100% Y51S (TAT→TCT)  fis → DNA‑binding protein Fis

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,339,3130A→C100.0% 573.0 160Y51S (TAT→TCT) fisDNA‑binding protein Fis
Reads supporting (aligned to +/- strand):  new base (80/80):  ref base (0/0):  total (80/80)

Only 100 of 165 total aligned reads displayed.

TCAACTGAATGGTCAGGATGTGAATGACCTCTATGAGCTGGTACTGGCTGAAGTAGAACAGCCC  >  REL606/3339281‑3339344
                                |                               
tCAACTGAATGGTCAGGATGTGAATGACCTCTCTg                               >  1:2052917/1‑35 (MQ=255)
tCAACTGAATGGTCAGGATGTGAATGACCTCTCTg                               >  1:15971494/1‑35 (MQ=255)
tCAACTGAATGGTCAGGATGTGAATGACCTCTCTg                               <  1:329357/35‑1 (MQ=255)
 cAACTGAATGGTCAGGATGTGAATGACCTCTCTGa                              <  1:16468759/35‑1 (MQ=255)
  aaCTGAATGGTCAGGATGTGAATGACCTCTCTGAg                             >  1:12953258/1‑35 (MQ=255)
  aaCTGAATGGTCAGGATGTGAATGACCTCTCTGAg                             <  1:18888820/35‑1 (MQ=255)
   aCTGAATGGTCAGGATGTGAATGACCTCTCTGAGc                            >  1:10457483/1‑35 (MQ=255)
   aCTGAATGGTCAGGATGTGAATGACCTCTCTGAGc                            >  1:20357937/1‑35 (MQ=255)
   aCTGAATGGTCAGGATGTGAATGACCTCTCTGAGc                            <  1:16981629/35‑1 (MQ=255)
   aCTGAATGGTCAGGATGTGAATGACCTCTCTGAGc                            >  1:11178817/1‑35 (MQ=255)
    cTGAATGGTCAGGATGTGAATGACCTCTCTGAGCt                           <  1:13405553/35‑1 (MQ=255)
    cTGAATGGTCAGGATGTGAATGACCTCTCTGAGCt                           >  1:1704448/1‑35 (MQ=255)
    cTGAATGGTCAGGATGTGAATGACCTCTCTGAGCt                           >  1:7935813/1‑35 (MQ=255)
    cTGAATGGTCAGGATGTGAATGACCTCTCTGAGCt                           <  1:15045201/35‑1 (MQ=255)
    cTGAATGGTCAGGATGTGAATGACCTCTCTGAGCt                           <  1:9661798/35‑1 (MQ=255)
     tGAATGGTCAGGATGTGAATGACCTCTCTGAGCTg                          >  1:18182716/1‑35 (MQ=255)
     tGAATGGTCAGGATGTGAATGACCTCTCTGAGCTg                          <  1:19193144/35‑1 (MQ=255)
     tGAATGGTCAGGATGTGAATGACCTCTCTGAGCTg                          >  1:3031606/1‑35 (MQ=255)
      gAATGGTCAGGATGTGAATGACCTCTCTGAGCTgg                         >  1:15802772/1‑35 (MQ=255)
      gAATGGTCAGGATGTGAATGACCTCTCTGAGCTgg                         <  1:2535272/35‑1 (MQ=255)
      gAATGGTCAGGATGTGAATGACCTCTCTGAGCTgg                         >  1:3232214/1‑35 (MQ=255)
      gAATGGTCAGGATGTGAATGACCTCTCTGAGCTgg                         >  1:6917628/1‑35 (MQ=255)
      gAATGGTCAGGATGTGAATGACCTCTCTGAGCTgg                         >  1:848824/1‑35 (MQ=255)
       aaTTGTCAGGATGTGAATGACCTCTCTGAGCTGGt                        >  1:4079872/1‑35 (MQ=37)
       aaTGGTCAGGATGTGAATGACCTCTCTGAGCTGGt                        <  1:20047129/35‑1 (MQ=255)
       aaTGGTCAGGATGTGAATGACCTCTCTGAGCTGGt                        <  1:3705589/35‑1 (MQ=255)
       aaTGGTCAGGATGTGAATGACCTCTCTGAGCTGGt                        >  1:10865486/1‑35 (MQ=255)
       aaTGGTCAGGATGCGAATGACCTCTCTGAGCTGGt                        >  1:4010365/1‑35 (MQ=37)
        acggTCAGGATGTGAATGACCTCTCTGAGCTGGTa                       <  1:16272547/33‑1 (MQ=38)
        aTGGTCAGGATGTGAATGACCTCTCTGAGCTGGTn                       >  1:2945317/1‑35 (MQ=255)
        aTGGTCAGGATGTGAATGACCTCTCTGAGCTGGTa                       >  1:7085116/1‑35 (MQ=255)
        aTGGTCAGGATGTGAATGACCTCTCTGAGCTGGTa                       >  1:9860811/1‑35 (MQ=255)
        aTGGTCAGGATGTGAATGACCTCTCTGAGCTGGTa                       <  1:10147506/35‑1 (MQ=255)
        aTGGTCAGGATGTGAATGACCTCTCTGAGCTGGTa                       <  1:12099864/35‑1 (MQ=255)
        aTGGTCAGGATGTGAATGACCTATCTGAGCTGGTa                       <  1:15279358/35‑1 (MQ=37)
         tGGTCAGGATGTGAATGACCTCTCTGAGCTGGTAc                      >  1:15666167/1‑35 (MQ=255)
         tGGTCAGGATGTGAATGACCTCTCTGAGCTGGTAc                      <  1:8651212/35‑1 (MQ=255)
         tGGTCAGGATGTGAATGACCTCTCTGAGCTGGTAc                      >  1:1452670/1‑35 (MQ=255)
         tGGTCAGGATGTGAATGACCTCTCTGAGCTGGTAc                      <  1:11059/35‑1 (MQ=255)
          tgTCAGGGTGTGAATGACCTCTCTGAGCTGGTACt                     <  1:18549284/34‑1 (MQ=25)
          ggTCAGGATGTGAATGACCTCTCTGAGCTGGTACt                     <  1:13056898/35‑1 (MQ=255)
          ggTCAGGATGTGAATGACCTCTCTGAGCTGGTACt                     >  1:13842071/1‑35 (MQ=255)
          ggTCAGGATGTGAATGACCTCTCTGAGCTGGTACt                     <  1:4212227/35‑1 (MQ=255)
          ggTCAGGATGTGAATGACCTCTCTGAGCTGGTACt                     <  1:5990483/35‑1 (MQ=255)
          ggTCAGGATGTGAATGACCTCTCTGAGCTGGTACt                     >  1:10492702/1‑35 (MQ=255)
           gTCAGGATGTGAATGACCTCTCTGAGCTGGTACTg                    <  1:7061397/35‑1 (MQ=255)
           gTCAGGATGTGAATGACCTCTCTGAGCTGGTACTg                    <  1:1779990/35‑1 (MQ=255)
           gTCAGGATGTGAATGACCTCTCTGAGCTGGTACTg                    <  1:13954059/35‑1 (MQ=255)
           gTCAGGATGTGAATGACCTCTCTGAGCTGGTACTg                    >  1:14881451/1‑35 (MQ=255)
            tCAGGATGTGAATGACCTCTCTGAGCTGGTACTgg                   >  1:8168420/1‑35 (MQ=255)
            tCAGGATGTGAATGACCTCTCTGAGCTGGTACTgg                   <  1:12686414/35‑1 (MQ=255)
             cAGGATGTGAATGACCTCTCTGAGCTGGTACTGGc                  <  1:7737873/35‑1 (MQ=255)
             cAGGATGTGAATGACCTCTCTGAGCTGGTACTGGc                  >  1:18723667/1‑35 (MQ=255)
             cAGGATGTGAATGACCTCTCTGAGCTGGTACTGGc                  >  1:13735060/1‑35 (MQ=255)
             cAGGATGTGAATGACCTCTCTGAGCTGGTACTGGc                  >  1:15801953/1‑35 (MQ=255)
              aGGATGTGAATGACCTCTCTGAGCTGGTACTGGCt                 <  1:10661647/35‑1 (MQ=255)
                gATGTGAATGACCTCTCTGAGCTGGTACTGGCTGa               >  1:14369791/1‑35 (MQ=255)
                 aTGTGAATGACCTCTCTGAGCTGGTACTGGCTgga              >  1:11324912/1‑33 (MQ=37)
                 aTGTGAATGACCTCTCTGAGCTGGTACTGGCTgga              >  1:1538750/1‑33 (MQ=37)
                 aTGTGAATGACCTCTCTGAGCTGGTACTGGCTGaa              <  1:16278474/35‑1 (MQ=255)
                 aTGTGAATGACCTCTCTGAGCTGGTACTGGCTGaa              <  1:6692415/35‑1 (MQ=255)
                 aTGTGAATGACCTCTCTGAGCTGGTACTGGCTGaa              >  1:14127658/1‑35 (MQ=255)
                 aTGTGAATGACCTCTCTGAGCTGGTACTGGCTGaa              >  1:4979723/1‑35 (MQ=255)
                  tgtgAATGACCTCTCTGAGCTGGTACTGGCTGGAg             >  1:3049442/1‑35 (MQ=37)
                  tgtgAATGACCTCTCTGAGCTGGTACTGGCTGGAg             >  1:14780381/1‑35 (MQ=37)
                  tgtgAATGACCTCTCTGAGCTGGTACTGGCTGGAg             >  1:12867857/1‑35 (MQ=37)
                  tgtgAATGACCTCTCTGAGCTGGTACTGGCTGGAg             >  1:17762611/1‑35 (MQ=37)
                  tgtgAATGACCTCTCTGAGCTGGTACTGGCTGAAg             <  1:3078521/35‑1 (MQ=255)
                  tgtgAATGACCTCTCTGAGCTGGTACTGGCTGAAg             <  1:1937196/35‑1 (MQ=255)
                  tgtgAATGACCTCTCTGAGCTGGTACTGGCTGAAg             <  1:7632882/35‑1 (MQ=255)
                  tgtgAATGACCTCTCTGAGCTGGTACTGGCTGAAg             <  1:4869934/35‑1 (MQ=255)
                   gtgAATGACCTCTCTGAGCTGGTACTGGCTGAAGt            <  1:8789066/35‑1 (MQ=255)
                   gtgAATGACCTCTCTGAGCTGGTACTGGCTGAAGt            <  1:16383026/35‑1 (MQ=255)
                    tgAATGACCTCTCTGAGCTGGTACTGGCTGGAgga           >  1:298497/1‑33 (MQ=25)
                    tgAATGACCTCTCTGAGCTGGTACTGGCTGAAgga           >  1:5403858/1‑33 (MQ=37)
                    tgAATGACCTCTCTGAGCTGGTACTGGCTGAAGTa           <  1:11595934/35‑1 (MQ=255)
                    tgAATGACCTCTCTGAGCTGGTACTGGCTGAAGTa           <  1:5277704/35‑1 (MQ=255)
                     gAATGACCTCTCTGAGCTGGTACTGGCTGAAGTAg          <  1:9507296/35‑1 (MQ=255)
                      aaTGACCTCTCTGAGCTGGTACTGGCTGGAGTAGa         >  1:1798277/1‑35 (MQ=37)
                      aaTGACCTCTCTGAGCTGGTACTGGCTGGAGAAGa         >  1:7569711/1‑35 (MQ=21)
                      aaTGACCTCTCTGAGCTGGTACTGGCTGAAGTAGa         >  1:3819194/1‑35 (MQ=255)
                      aaTGACCTCTCTGAGCTGGTACTGGCTGAAGTAGa         <  1:10916961/35‑1 (MQ=255)
                      aaTGACCTCTCTGAGCTGGTACTGGCTGAAGCAGa         >  1:18697358/1‑35 (MQ=37)
                       aTGACCTCTCTGAGCTGGTACTGGCTGGAGGAGaa        >  1:19401988/1‑35 (MQ=25)
                       aTGACCTCTCTGAGCTGGTACTGGCTGGAGCAGaa        >  1:14000015/1‑35 (MQ=25)
                       aTGACCTCTCTGAGCTGGTACTGGCTGAAGTAGaa        >  1:10158490/1‑35 (MQ=255)
                       aTGACCTCTCTGAGCTGGTACTGGCTGAAGTAGaa        <  1:14149762/35‑1 (MQ=255)
                        tGACCTCTCTGAGCTGGTACTGGCTGAAGTAGAAc       <  1:12005721/35‑1 (MQ=255)
                        tGACCTCTCTGAGCTGGTACTGGCTGAAGTAGAAc       <  1:8754362/35‑1 (MQ=255)
                        tGACCTCTCTGAGCTGGTACTGGCTGAAGTAGAAc       <  1:9265077/35‑1 (MQ=255)
                         gACCTCTCTGAGCTGGTACTGGCTGAAGTAGAACa      <  1:3593862/35‑1 (MQ=255)
                         gACCTCTCTGAGCTGGTACTGGCTGAAGGAGAACa      >  1:12733038/1‑35 (MQ=37)
                          aCCTCTCTGAGCTGGTACTGGCTGAAGTAGAACAg     <  1:11844454/35‑1 (MQ=255)
                          aCCTCTCTGAGCTGGTACTGGCTGAAGTAGAACAg     <  1:3418199/35‑1 (MQ=255)
                          aCCTCTCTGAGCTGGTACTGGCTGAAGTAGAACAg     <  1:8986957/35‑1 (MQ=255)
                           ccTCTCTGAGCTGGTACTGGCTGGAGTAGAAAAGc    >  1:9005651/1‑35 (MQ=25)
                           ccTCTCTGAGCTGGTACTGGCTGAAGTAGAACAGc    <  1:8095345/35‑1 (MQ=255)
                            ctctCTGAGCTGGTACTGGCTGGAGTAGAACAGcc   >  1:4784331/1‑35 (MQ=37)
                             tctCTGAGCTGGTACTGGCTGAAGTAGAACAGccc  <  1:16593384/35‑1 (MQ=255)
                             tctCTGAGCTGGTACTGGCTGAAGTAGAACAGccc  <  1:2460984/35‑1 (MQ=255)
                                |                               
TCAACTGAATGGTCAGGATGTGAATGACCTCTATGAGCTGGTACTGGCTGAAGTAGAACAGCCC  >  REL606/3339281‑3339344

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.