Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,607,916 Δ1 bp 36.1% coding (17/1320 nt) ECB_01510 ← putative tail component of prophage

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,607,9160T→.36.1% 117.5 ‑0.2 0.083coding (17/1320 nt)ECB_01510putative tail component of prophage
Reads supporting (aligned to +/- strand):  new base (12/18):  ref base (24/29):  total (36/47)
Fisher's exact test strand distribution p-value = 8.18e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GCAAGATGCGGCGTGGCGGTTATCTGAAAATTGATCTTACCCAGCAATAGTGGACACGCGGCT  >  REL606/1607886‑1607948
                              |                                
gCAAGATGCGGCGTGGCGGTTATCTGAAAATTGAt                              <  1:18543752/35‑1 (MQ=40)
 cAAGATGCGGCGTGGCGGTTATCTGAAAATTGATc                             <  1:9624194/35‑1 (MQ=255)
 cAAGATGCGGCGTGGCGGTTATCTGAAAATGATCt                            <  1:14916390/35‑1 (MQ=1)
 cAAGATGCGGCGTGGCGGTTATCTGAAAATGATCt                            <  1:754947/35‑1 (MQ=1)
   aGATGCGGCGTGGCGGTTATCTGAAAATTGATCtt                           <  1:7995069/35‑1 (MQ=255)
   aGATGCGGCGTGGCGGTTATCTGAAAATGATCTTa                          <  1:8408131/35‑1 (MQ=11)
    gATGCGGCGTGGCGGTTATCTGAAAATTGATCTTa                          <  1:2768753/35‑1 (MQ=255)
    gATGCGGCGTGGCGGTTATCTGAAAATGATCTTAc                         <  1:86773/35‑1 (MQ=14)
     aTGCGGCGTGGCGGTTATCTGAAAATTGATCTTAc                         <  1:8360444/35‑1 (MQ=255)
     aTGCGGCGTGGCGGTTATCTGAAAATGATCTTAcc                        <  1:2896476/35‑1 (MQ=14)
      tGCGGCGTGGCGGTTATCTGAAAATTGATCTTAcc                        <  1:18587753/35‑1 (MQ=255)
      tGCGGCGTGGCGGTTATCTGAAAATGATCTTAccc                       <  1:12037829/35‑1 (MQ=17)
        cggcgTGGCGGTTATCTGAAAATTGATCTTACCCa                      <  1:6467694/35‑1 (MQ=255)
        cggcgTGGCGGTTATCTGAAAATTGATCTTACCCa                      <  1:2141460/35‑1 (MQ=255)
        cggcgTGGCGGTTATCTGAAAATGATCTTACCCan                     <  1:19116429/35‑2 (MQ=17)
         ggcgTGGCGGTTATCTGAAAATTGATCTTACCCAg                     <  1:714062/35‑1 (MQ=255)
         ggcgTGGCGGTTATCTGAAAATTGATCTTACCCAg                     <  1:1317060/35‑1 (MQ=255)
         ggcgTGGCGGTTATCTGAAAATTGATCTTACCCAg                     <  1:12679187/35‑1 (MQ=255)
         ggcgTGGCGGTTATCTGAAAATGATCTTACACAGc                    <  1:5642447/35‑1 (MQ=17)
          gcgTGGCGGTTATCTGAAAATGATCTTACCCAGCa                   <  1:9666645/35‑1 (MQ=21)
           cgTGGCGGTTATCTGAAAATTGATCTTACCCAGCa                   <  1:3003507/35‑1 (MQ=255)
            gTGGCGGTTATCTGAAAATTGATCTTACCCAGCaa                  <  1:19548761/35‑1 (MQ=255)
            gTGGCGGTTATCTGAAAATGATCTTACCCAGCAAt                 <  1:17770994/35‑1 (MQ=25)
             tGGCGGTTATCTGAAAATGATCTTACCCAGCAATa                <  1:20254620/35‑1 (MQ=25)
              ggCGGTTATCTGAAAATTGATCTTACCCAGCAATa                <  1:11042134/35‑1 (MQ=255)
               gCGGTTATCTGAAAATTGATCTTACCCAGCAATAg               <  1:9667480/35‑1 (MQ=255)
                cGGTTATCTGAAAATTGATCTTACCCAGCAATAGt              <  1:17161084/35‑1 (MQ=255)
                cGGTTATCTGAAAATTGATCTTACCCAGCAATAGt              >  1:14916442/1‑35 (MQ=255)
                cGGTTATCTGAAAATGATCTTACCCAGCAATAgag             >  1:5821899/1‑33 (MQ=18)
                cGGTTATCTGAAAATGATCTTACCCAGCAATAGTg             <  1:16924156/35‑1 (MQ=18)
                cGGTTATCTGAAAATGATCTTACCCAGCAATAGTg             >  1:7848567/1‑35 (MQ=18)
                cGGTTATCTGAAAATGATCTTACCCAGCAATAGTg             >  1:1416904/1‑35 (MQ=18)
                cGGTTATCTGAAAATGATCTTACCCAGCAATAGTg             <  1:1055423/35‑1 (MQ=18)
                cGGTTATCTGAAAATGATCTTACCCAGCAATAGTg             <  1:15675205/35‑1 (MQ=18)
                cGGTTATCTCAAAATTGATCTTACCCAGCAATAGt              >  1:7002073/1‑35 (MQ=255)
                 ggTTATCTGAAACTTGATCTTACCCAGCAATAGTg             >  1:7035198/1‑35 (MQ=255)
                 ggTTATCTGAAAATTGATCTTACCCAGCAATAGTg             <  1:16699404/35‑1 (MQ=255)
                 ggTTATCTGAAAATTGATCTTACCCAGCAATAGTg             >  1:18313293/1‑35 (MQ=255)
                 ggTTATCTGAAAATTGATCTTACCCAGCAATAGTg             <  1:15824653/35‑1 (MQ=255)
                 ggTTATCTGAAAATTGATCTTACCCAGCAATAGTg             <  1:12597631/35‑1 (MQ=255)
                 ggTTATCTGAAAATGATCTTACCCAGCAATAGTgg            >  1:17095253/1‑35 (MQ=18)
                 ggTTATCTGAAAATGATCTTACCCAGCAATAGTgg            >  1:14116230/1‑35 (MQ=18)
                 ggTTATCTGAAAATGATCTTACCCAGCAATAGTgg            <  1:13235403/35‑1 (MQ=18)
                  gTTATCTGAAAATTGATCTTACCCAGCAATAGTgg            >  1:10002320/1‑35 (MQ=255)
                  gTTATCTGAAAATTGATCTTACCCAGCAATAGTgg            <  1:10293477/35‑1 (MQ=255)
                  gTTATCTGAAAATTGATCTTACCCAGCAATAGTgg            >  1:12110953/1‑35 (MQ=255)
                   ttATCTGAAAATTGATCTTACCCAGCAATAGTGGa           >  1:412842/1‑35 (MQ=255)
                   ttATCTGAAAATGATCTTACCCAGCAATAGTGGac          >  1:1563574/1‑35 (MQ=17)
                   ttATCTGAAAATGATCTTACCCAGCAATAGTGGac          <  1:13984342/35‑1 (MQ=17)
                   ttATCTGAAAATGACCTTACCCAGCAATAGTGGac          >  1:10055934/1‑35 (MQ=17)
                    tATCTGAAAATTGATCTTACCCAGCAATAGTGGac          >  1:6074937/1‑35 (MQ=255)
                    tATCTGAAAATGATCTTACCCAGCAATAGTGGaca         >  1:16349039/1‑35 (MQ=14)
                     aTCTGAAAATTGATCTTACCCAGCAATAGTGGaca         >  1:9057435/1‑35 (MQ=255)
                     aTCTGAAAATTGATCTTACCCAGCAATAGTGGaca         >  1:9075744/1‑35 (MQ=255)
                     aTCTGAAAATGATCTTACCCAGCAATAGTGGacac        >  1:10802390/1‑35 (MQ=14)
                     aTCTGAAAATGATCTTACCCAGCAATAGTGGacac        <  1:13692983/35‑1 (MQ=14)
                     aTCTGAAAATGATCTTACCCAGCAATAGTGGacac        <  1:18397026/35‑1 (MQ=14)
                     aTCTGAAAATGATCTTACCCAGCAATAGCGGacac        >  1:268202/1‑35 (MQ=14)
                      tCTGAAAATTGATCTTACCCAGCAATAGTGGacac        <  1:11688175/35‑1 (MQ=255)
                      tCTGAAAATTGATCTTACCCAGCAATAGTGGacac        >  1:2765446/1‑35 (MQ=255)
                      tCTGAAAATTGATCTTACCCAGCAATAGTGGacac        >  1:3948482/1‑35 (MQ=255)
                      tCTGAAAATGATCTTACCCAGCAATAGTGGACAcg       >  1:9066961/1‑35 (MQ=11)
                       cTGAAAATTGATCTTACCCATCAATAGTGGACAcg       >  1:11018612/1‑35 (MQ=255)
                       cTGAAAATTGATCTTACCCAGCAATAGTGGACAcg       >  1:11498553/1‑35 (MQ=255)
                       cTGAAAATTGATCTTACCCAGCAATAGTGGACAcg       >  1:8248664/1‑35 (MQ=255)
                       cTGAAAATTGATCTTACCCAGCAATAGTGGACAcg       >  1:8969902/1‑35 (MQ=255)
                       cTGAAAATGATCTTACCCAGCAATAGTGGACAcgc      >  1:919346/1‑35 (MQ=11)
                        tGAAAATTGATCTTACCCAGCAATAGTGGACAcgc      <  1:19131340/35‑1 (MQ=40)
                        tGAAAATTGATCTTACCCAGCAATAGTGGACAcgc      <  1:5220493/35‑1 (MQ=40)
                         gAAAATTGATCTTACCCAGCAATAGTGGACAcgcg     <  1:1842676/35‑1 (MQ=38)
                         gAAAATTGATCTTACCCAGCAATAGTGGACAcgcg     >  1:16296802/1‑35 (MQ=38)
                         gAAAATTGATCTTACCCAGCAATAGTGGACAcgcg     >  1:14766655/1‑35 (MQ=38)
                         gAAAATTGATCTTACCCAGCAATAGTGGACAcgcg     >  1:12641877/1‑35 (MQ=38)
                         gAAAATTATCTTACCCAGCAATAGTGGACACGCgg    >  1:14274543/1‑35 (MQ=11)
                          aaaaTTGATCTTACCCAGCAATAGTGGACACGCgg    >  1:9981010/1‑35 (MQ=38)
                          aaaaTTGATCTTACCCAGCAATAGTGGACACGCgg    <  1:13102508/35‑1 (MQ=38)
                           aaaTTGATCTTACCCAGCAATAGTGGACACGCGGc   <  1:16706066/35‑1 (MQ=35)
                           aaaTTGATCTTACCCAGCAATAGTGGACACGCGGc   <  1:6358580/35‑1 (MQ=35)
                            aaTTGATCTTACCCAGCAATAGTGGACACGCGGCt  >  1:12352185/1‑35 (MQ=32)
                            aaTTGATCTTACCCAGCAATAGTGGACACGCGGCt  >  1:1304264/1‑35 (MQ=32)
                            aaTTGATCTTACCCAGCAATAGTGGACACGCGGCt  <  1:2655292/35‑1 (MQ=32)
                            aaTTGATCTTACCCAGCAATAGTGGACACGCGGCt  >  1:14665205/1‑35 (MQ=32)
                            aaTTGATCTTACACAGCAATAGTGGACACGCGGCt  <  1:11380142/35‑1 (MQ=12)
                              |                                
GCAAGATGCGGCGTGGCGGTTATCTGAAAATTGATCTTACCCAGCAATAGTGGACACGCGGCT  >  REL606/1607886‑1607948

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ 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.