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 I200 R1 216 30.8 1683622 97.2% 1636480 87.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.97.3% 144.3 / ‑1.4 37coding (1253/2187 nt)espL4type III secretion system effector EspL4
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (9/27);  total (10/27)
Fisher's exact test for biased strand distribution p-value = 2.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-01

TACAACATTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTAT  >  NZ_CP009273/4211078‑4211248
                                                                                     |                                                                                     
tACAACATTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGAttttt                                                                                   >  1:407556/1‑90 (MQ=255)
    acaTTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGAtttttttcg                                                                               <  1:691398/90‑3 (MQ=255)
          ttATCCTTACGTGATATCAGTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTa                                                                        <  2:384414/90‑1 (MQ=255)
           tATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTaa                                                                       <  1:817498/90‑1 (MQ=255)
           tATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTaa                                                                       <  2:501645/90‑1 (MQ=255)
            aTCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAAt                                                                      <  1:244524/90‑1 (MQ=255)
             tCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAGTATGCATGAGATTTTTTTCGAGAGTAATc                                                                     <  1:689151/90‑1 (MQ=255)
             tCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATc                                                                     <  2:166667/90‑1 (MQ=255)
             tCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATc                                                                     <  1:98688/90‑1 (MQ=255)
             tCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATc                                                                     <  1:525344/90‑1 (MQ=255)
               cTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCt                                                                   >  2:536703/1‑90 (MQ=255)
                ttACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCtt                                                                  >  2:144390/1‑90 (MQ=255)
                 tACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTc                                                                 >  2:77431/1‑90 (MQ=255)
                   cGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCat                                                               <  1:374833/90‑1 (MQ=255)
                          tCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGt                                                        >  1:19446/1‑90 (MQ=255)
                            gCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGttt                                                      >  2:160087/1‑90 (MQ=255)
                                aaaaTAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGa                                                  <  1:429347/90‑1 (MQ=255)
                                   ataatCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACaa                                               <  1:144390/90‑1 (MQ=255)
                                    taatCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAAt                                              <  2:636561/90‑1 (MQ=255)
                                     aatCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATg                                             <  2:258412/90‑1 (MQ=255)
                                     aatCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATg                                             <  1:77431/90‑1 (MQ=255)
                                     aatCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATg                                             <  2:657845/90‑1 (MQ=255)
                                      atCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATg                                             >  1:173395/1‑90 (MQ=255)
                                       tCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGaa                                           >  1:612582/1‑90 (MQ=255)
                                           gAAAAACCATTTTTATCTTTTGCCTCCAGATTATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGa                                       >  1:29315/1‑90 (MQ=255)
                                              aaaCCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGaaaa                                    <  1:370442/90‑1 (MQ=255)
                                               aaCCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAAt                                   <  1:554682/90‑1 (MQ=255)
                                               aaCCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAAt                                   <  2:258149/90‑1 (MQ=255)
                                               aaCCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAAt                                   <  2:398846/90‑1 (MQ=255)
                                                            ttttGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGata                      <  1:132868/90‑1 (MQ=255)
                                                            ttttGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGata                      <  2:603218/90‑1 (MQ=255)
                                                                   tCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCCCAACGATATCCGCAt               <  2:831690/90‑1 (MQ=255)
                                                                        aaTATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCa          <  2:405635/90‑1 (MQ=255)
                                                                         atatGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCAt         <  2:647999/90‑1 (MQ=255)
                                                                            tGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATa                                  >  2:425368/1‑62 (MQ=255)
                                                                            tGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATa                                  <  1:425368/62‑1 (MQ=255)
                                                                                tGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTtat  <  2:642891/90‑1 (MQ=255)
                                                                                tGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTtat  >  2:472940/1‑90 (MQ=255)
                                                                                tGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTtat  <  1:71604/90‑1 (MQ=255)
                                                                                     |                                                                                     
TACAACATTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTAT  >  NZ_CP009273/4211078‑4211248

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ACATTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATAGCCAAATATAAACC  >  NZ_CP009273/4211082‑4211263
                                                                                        |                                                                                             
ACATTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCC                                                                                   <  SRR3722079.698669/100‑1 (MQ=60)
       TATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATAT                                                                            <  SRR3722079.826257/100‑1 (MQ=60)
        ATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATC                                                                           <  SRR3722079.246583/100‑1 (MQ=60)
         TCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAGTATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCC                                                                          <  SRR3722079.696398/100‑1 (MQ=60)
         TCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCC                                                                          <  SRR3722079.99468/100‑1 (MQ=60)
         TCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCC                                                                          <  SRR3722079.530217/100‑1 (MQ=60)
            TTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGT                                                                       >  SRR3722079.19582/1‑100 (MQ=60)
               CGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTC                                                                    <  SRR3722079.377868/100‑1 (MQ=60)
                        GCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATG                                                            >  SRR3722079.174903/1‑100 (MQ=60)
                         CTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAA                                                          >  SRR3722079.618778/1‑100 (MQ=60)
                            AAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTG                                                       <  SRR3722079.433063/100‑1 (MQ=60)
                             AAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGATTATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGA                                                      >  SRR3722079.29531/1‑100 (MQ=60)
                               ATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAA                                                    <  SRR3722079.145630/100‑1 (MQ=60)
                                 AATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAAT                                                  <  SRR3722079.78064/100‑1 (MQ=60)
                                          AAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAAC                                         <  SRR3722079.373413/100‑1 (MQ=60)
                                           AACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACG                                        <  SRR3722079.560022/100‑1 (MQ=60)
                                            gatgtgtataagagacagCTCCAGAATATGTATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGA                                       <  SRR3722079.429040/82‑1 (MQ=60)
                                                        TTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTAC                           <  SRR3722079.133976/100‑1 (MQ=60)
                                                                   GAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATA                >  SRR3722079.61128/1‑100 (MQ=60)
                                                                   GAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATA                >  SRR3722079.336460/1‑100 (MQ=60)
                                                                   GAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATA                >  SRR3722079.770533/1‑100 (MQ=60)
                                                                            TGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATAGCCAAATAT       <  SRR3722079.72202/100‑1 (MQ=60)
                                                                                  TTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATAGCCAAATATAAACC  <  SRR3722079.114361/100‑1 (MQ=60)
                                                                                  TTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATAGCCAAATATAAACC  <  SRR3722079.277434/100‑1 (MQ=60)
                                                                                        |                                                                                             
ACATTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATAGCCAAATATAAACC  >  NZ_CP009273/4211082‑4211263

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