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 F1 I1 R1 154 26.6 2069536 90.2% 1866721 108.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 3,515,314 Δ1 bp intergenic (+253/‑570) tRNA‑Pro (CGG) → / → pilB tRNA‑Pro/Type IV pilus assembly protein PilB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604793,515,3140G.100.0% 115.8 / NA 26intergenic (+253/‑570)tRNA‑Pro (CGG)/pilBtRNA‑Pro/Type IV pilus assembly protein PilB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (13/13);  total (13/13)

GCCGCGGAAAACAGCTGGCTACGACCCACTACTGCATGGGGCGCGAGCTTGCCGTCGCGCACCAGCCGCCGGAAGGTGGCAATTGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCGTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCTTCATCAGGTCT  >  AM260479/3515186‑3515440
                                                                                                                                |                                                                                                                               
gCCGCGGAAAACAGCTGGCTACGACCCACTACTGCATGGGGCGCGAGCTTGCCGTCGCGCACCAGCCGCCGGAAGGTGGCAATTGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCTTCcgc                                                                                                                           <  1:926443/134‑1 (MQ=255)
    cgGAAAACAGCTGGCTACGACCCACTACTGCATGGGGCGCGAGCTTGCCGTCGCGCACCAGCCGCCGGAAGGTGGCAATTGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCTTCCGCGAgg                                                                                                                       >  2:748043/1‑134 (MQ=255)
          aCAGCTGGCTACGACCCACTACTGCATGGGGCGCGAGCTTGCCGTCGCGCACCAGCCGCCGGAAGGTGGCAATTGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCTTCCGCGAGGTGGCCa                                                                                                                 <  1:873126/134‑1 (MQ=255)
                 gCTACGACCCACTACTGCATGGGGCGCGAGCTTGCCGTCGCGCACCAGCCGCCGGAAGGTGGCAATTGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCTTCCGCGAGGTGGCCAAAGACt                                                                                                           <  1:201957/133‑1 (MQ=255)
                  cTACGACCCACTACTGCATGGGGCGCGAGCTTGCCGTCGCGCACCAGCCGCCGGAAGGTGGCAATTGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCTTCCGCGAGGTGGCCAAAGACTTc                                                                                                         >  2:975139/1‑134 (MQ=255)
                                      gggCGCGAGCTTGCCGTCGCGCACCAGCCGCCGGAAGGTGGCAATTGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGt                                                                                     >  2:654824/1‑134 (MQ=255)
                                           cgAGCTTGCCGTCGCGCACCAGCCGCCGGAAGGTGGCAATTGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCTTCcgc                                                                                                                           <  1:652800/91‑1 (MQ=255)
                                           cgAGCTTGCCGTCGCGCACCAGCCGCCGGAAGGTGGCAATTGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCTTCcgc                                                                                                                           >  2:652800/1‑91 (MQ=255)
                                           cgAGCTTGCCGTCGCGCACCAGCCGCCGGAAGGTGGCAATTGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTc                                                                                          >  1:270256/1‑124 (MQ=255)
                                           cgAGCTTGCCGTCGCGCACCAGCCGCCGGAAGGTGGCAATTGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTc                                                                                          <  2:270256/124‑1 (MQ=255)
                                             aGCTTGCCGTCGCGCACCAGCCGCCGGAAGGTGGCAATTGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGa                                                                               >  2:180635/1‑133 (MQ=255)
                                                ttGCCGTCGCGCACCAGCCGCCGGAAGGTGGCAATTGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCTTCcg                                                                                                                            >  1:873400/1‑85 (MQ=255)
                                                ttGCCGTCGCGCACCAGCCGCCGGAAGGTGGCAATTGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCTTCcg                                                                                                                            <  2:873400/85‑1 (MQ=255)
                                                  gCCGTCGCGCACCAGCCGCCGGAAGGTGGCAATTGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGc                                                                           <  1:975139/132‑1 (MQ=255)
                                                                    ccgGAAGGTGGCAATTGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTg                                                        >  2:628315/1‑133 (MQ=255)
                                                                          ggTGGCAATTGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCTTCCGCGAGGTGGCCaa                                                                                                                >  1:916160/1‑71 (MQ=255)
                                                                          ggTGGCAATTGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCTTCCGCGAGGTGGCCaa                                                                                                                <  2:916160/71‑1 (MQ=255)
                                                                          ggTGGCAATTGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACt                                                 <  2:115685/134‑1 (MQ=255)
                                                                          ggTGGCAATTGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACt                                                 <  2:438950/134‑1 (MQ=255)
                                                                           gTGGCAATTGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACtt                                                <  1:654824/134‑1 (MQ=255)
                                                                                   tGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCt                                        >  1:785699/1‑134 (MQ=255)
                                                                                       aCTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTg                                    >  2:324330/1‑134 (MQ=255)
                                                                                                         gccTCCGCGGCGGTGAAATCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGAc                  <  1:324330/134‑1 (MQ=255)
                                                                                                                    ggTGAAATCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTtcat         >  2:857222/1‑132 (MQ=255)
                                                                                                                        aaaTCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGt    <  1:748043/133‑1 (MQ=255)
                                                                                                                         aaTCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGtct  >  1:621606/1‑134 (MQ=255)
                                                                                                                                |                                                                                                                               
GCCGCGGAAAACAGCTGGCTACGACCCACTACTGCATGGGGCGCGAGCTTGCCGTCGCGCACCAGCCGCCGGAAGGTGGCAATTGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCGTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCTTCATCAGGTCT  >  AM260479/3515186‑3515440

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GGCCGCGGAAAACAGCTGGCTACGACCCACTACTGCATGGGGCGCGAGCTTGCCGTCGCGCACCAGCCGCCGGAAGGTGGCAATTGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCGTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCTTCATCAGGTCTTCTGCGGTCATGAT  >  AM260479/3515185‑3515454
                                                                                                                                 |                                                                                                                                             
gtgggctcggagatgtgtataagagacagCTACTGCATGGGGCGCGAGCTTGCCGTCGCGCACCAGCCGCCGGAAGGTGGCAATTGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCTTCCGCGAGGTGGCCAAAGAC                                                                                                                          <  GWNJ‑0478:712:GW2002102894th:2:2111:11203:5337/121‑1 (MQ=60)
GGCCGCGGAAAACAGCTGGCTACGACCCACTACTGCATGGGGCGCGAGCTTGCCGTCGCGCACCAGCCGCCGGAAGGTGGCAATTGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCTTCCGCGAGGTGGCCAAAGAC                                                                                                                          <  GWNJ‑0478:712:GW2002102894th:2:2211:8030:93857/150‑1 (MQ=60)
          AACAGCTGGCTACGACCCACTACTGCATGGGGCGCGAGCTTGCCGTCGCGCACCAGCCGCCGGAAGGTGGCAATTGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGC                                                                                                                <  GWNJ‑0478:712:GW2002102894th:2:2206:18081:65659/150‑1 (MQ=60)
                 GGCTACGACCCACTACTGCATGGGGCGCGAGCTTGCCGTCGCGCACCAGCCGCCGGAAGGTGGCAATTGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAG                                                                                                          <  GWNJ‑0478:712:GW2002102894th:2:1112:15034:89672/149‑1 (MQ=60)
                             CTACTGCATGGGGCGCGAGCTTGCCGTCGCGCACCAGCCGCCGGAAGGTGGCAATTGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTG                                                                                              >  GWNJ‑0478:712:GW2002102894th:2:1203:10475:6502/1‑149 (MQ=60)
                                  GCATGGGGCGCGAGCTTGCCGTCGCGCACCAGCCGCCGGAAGGTGGCAATTGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCTTCCGCGAGGTGGCCAAAGACTgtctcttatacacatctccgagcccacgagacg                                                                                        >  GWNJ‑0478:712:GW2002102894th:2:2206:13409:67498/1‑117 (MQ=60)
                                                  TGCCGTCGCGCACCAGCCGCCGGAAGGTGGCAATTGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGG                                                                          <  GWNJ‑0478:712:GW2002102894th:2:2213:5294:73824/148‑1 (MQ=60)
                                                            CACCAGCCGCCGGAAGGTGGCAATTGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGctgtctcttatacacatctccgagcccacgagacggagctacatctcgt                                                              >  GWNJ‑0478:712:GW2002102894th:2:2211:16977:34220/1‑101 (MQ=60)
                                                                     CCGGAAGGTGGCAATTGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTC                                                     >  GWNJ‑0478:712:GW2002102894th:2:2201:13429:89440/1‑150 (MQ=60)
                                                                           GGTGGCAATTGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGC                                               <  GWNJ‑0478:712:GW2002102894th:2:2111:19901:16535/150‑1 (MQ=60)
                                                                                                         CGCCTCCGCGGCGGTGAAATCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTC                 <  GWNJ‑0478:712:GW2002102894th:2:1204:1531:74015/150‑1 (MQ=60)
                                                                                                           CCTCCGCGGCGGTGAAATCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTT               >  GWNJ‑0478:712:GW2002102894th:2:2108:19417:10690/1‑150 (MQ=60)
                                                                                                                     GGTGAAATCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCAT     >  GWNJ‑0478:712:GW2002102894th:2:2109:3495:77182/1‑150 (MQ=60)
                                                                                                                        GAAATCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGAT  >  GWNJ‑0478:712:GW2002102894th:2:1102:21040:42971/1‑150 (MQ=60)
                                                                                                                        GAAATCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGA   <  GWNJ‑0478:712:GW2002102894th:2:2114:18840:54498/149‑1 (MQ=60)
                                                                                                                                 |                                                                                                                                             
GGCCGCGGAAAACAGCTGGCTACGACCCACTACTGCATGGGGCGCGAGCTTGCCGTCGCGCACCAGCCGCCGGAAGGTGGCAATTGAAACTTCGAGGTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCGTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCTTCATCAGGTCTTCTGCGGTCATGAT  >  AM260479/3515185‑3515454

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