Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F2 I195 R1 222 27.1 1551108 96.0% 1489063 86.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,211,170 (T)8→7 coding (1246/2187 nt) espL4 ← type III secretion system effector EspL4

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,211,1630T.100.0% 125.2 / NA 30coding (1253/2187 nt)espL4type III secretion system effector EspL4
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (14/16);  total (14/16)

TACAACATTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTAT  >  NZ_CP009273/4211078‑4211248
                                                                                     |                                                                                     
tACAACATTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGAttttt                                                                                   >  2:234195/1‑90 (MQ=255)
tACAACATTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGAttttt                                                                                   >  1:694596/1‑90 (MQ=255)
tACAACATTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGAttttt                                                                                   >  1:239293/1‑90 (MQ=255)
       tttttATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCgaga                                                                           <  1:121943/90‑1 (MQ=255)
           tATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTaa                                                                       <  1:15087/90‑1 (MQ=255)
            aTCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAAt                                                                      <  1:347971/90‑1 (MQ=255)
             tCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATc                                                                     <  1:402698/90‑1 (MQ=255)
                 tACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTc                                                                 >  2:681989/1‑90 (MQ=255)
                 tACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTc                                                                 >  2:26703/1‑90 (MQ=255)
                 tACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTc                                                                 >  1:95193/1‑90 (MQ=255)
                   cGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCat                                                               <  2:611069/90‑1 (MQ=255)
                            gCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGttt                                                      >  2:519720/1‑90 (MQ=255)
                            gCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGttt                                                      >  2:306559/1‑90 (MQ=255)
                              cagaaaTAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCt                                                    >  1:650399/4‑90 (MQ=255)
                                 aaaTAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGaa                                                 >  1:20/1‑89 (MQ=255)
                                     aatCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATg                                             <  2:699899/90‑1 (MQ=255)
                                         cAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAAtt                                         <  1:250211/90‑1 (MQ=255)
                                         cAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAAtt                                         <  1:320731/90‑1 (MQ=255)
                                              aaaCCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGaaaa                                    <  1:238662/90‑1 (MQ=255)
                                                aCCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATa                                  <  2:348890/90‑1 (MQ=255)
                                                       ttATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCaa                           <  2:522249/90‑1 (MQ=255)
                                                            ttttGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGata                      <  2:126856/90‑1 (MQ=255)
                                                            ttttGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGata                      <  2:287813/90‑1 (MQ=255)
                                                            ttttGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGata                      >  2:570519/1‑90 (MQ=255)
                                                             tttGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGatat                     >  1:85243/1‑90 (MQ=255)
                                                             tttGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGatat                     >  2:308082/1‑90 (MQ=255)
                                                              ttGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATc                    >  1:282560/1‑90 (MQ=255)
                                                               tGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATcc                   <  1:681989/90‑1 (MQ=255)
                                                                          tatGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATa        <  2:63638/90‑1 (MQ=255)
                                                                           atGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATaa       >  1:432368/1‑90 (MQ=255)
                                                                              cATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAAttt    >  2:499104/1‑90 (MQ=255)
                                                                                tGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTtat  <  2:92960/90‑1 (MQ=255)
                                                                                tGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTtat  >  1:99859/1‑90 (MQ=255)
                                                                                     |                                                                                     
TACAACATTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTAT  >  NZ_CP009273/4211078‑4211248

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATAGCCAAATATAAACC  >  NZ_CP009273/4211085‑4211263
                                                                                     |                                                                                             
TTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTC                                                                                <  SRR3722073.124532/100‑1 (MQ=60)
TTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTC                                                                                >  SRR3722073.97297/1‑100 (MQ=60)
    TATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATAT                                                                            <  SRR3722073.15472/100‑1 (MQ=60)
     ATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATC                                                                           <  SRR3722073.355341/100‑1 (MQ=60)
      TCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCC                                                                          <  SRR3722073.411285/100‑1 (MQ=60)
             gtataagagacagAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCT                                                                   >  SRR3722073.663770/14‑100 (MQ=60)
                ATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGAA                                                                >  SRR3722073.22/1‑100 (MQ=60)
                                  CAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTC                                              <  SRR3722073.255466/100‑1 (MQ=60)
                                  CAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTC                                              <  SRR3722073.327531/100‑1 (MQ=60)
                                       AAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAAC                                         <  SRR3722073.243662/100‑1 (MQ=60)
                                            ATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATAT                                    >  SRR3722073.87165/1‑100 (MQ=60)
                                             TTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATC                                   >  SRR3722073.288541/1‑100 (MQ=60)
                                                        TGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCAT                        <  SRR3722073.695911/100‑1 (MQ=60)
                                                          CCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAA                      >  SRR3722073.441618/1‑100 (MQ=60)
                                                               AGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTAT                 >  SRR3722073.102033/1‑100 (MQ=60)
                                                                  ATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATAGC              >  SRR3722073.656678/1‑100 (MQ=60)
                                                                       CATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATAGCCAAAT         >  SRR3722073.153499/1‑100 (MQ=60)
                                                                              ATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAGAATAGTCTCAACGATATCCGCATTACCATAATTTATAGCCAAATATAAAC   >  SRR3722073.493565/1‑100 (MQ=60)
                                                                              ATTTTTTTCGAGAGAAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATAGCCAAATATAAAC   >  SRR3722073.174205/1‑100 (MQ=60)
                                                                               TTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATAGCCAAATATAAACC  <  SRR3722073.345981/100‑1 (MQ=60)
                                                                               TTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATAGCCAAATATAAACC  <  SRR3722073.513729/100‑1 (MQ=60)
                                                                                     |                                                                                             
TTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATAGCCAAATATAAACC  >  NZ_CP009273/4211085‑4211263

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