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 F3 I122 R1 193 20.9 1228102 94.7% 1163012 84.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,650,010 (T)6→5 coding (91/333 nt) hdeA ← acid‑activated periplasmic chaperone HdeA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,650,0050T.100.0% 136.6 / NA 32coding (96/333 nt)hdeAacid‑activated periplasmic chaperone HdeA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (14/18);  total (14/18)

GCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGACCACCAAGAATAACGCCTAA  >  NZ_CP009273/3649923‑3650088
                                                                                  |                                                                                   
gcgcTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGtttttgtt                                                                              >  1:121129/1‑87 (MQ=255)
 cgcTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGtttttgtta                                                                             <  2:179329/90‑5 (MQ=255)
     tCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTGTTATCAg                                                                        >  1:270369/1‑90 (MQ=255)
     tCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTGTTATCAg                                                                        <  2:509930/90‑1 (MQ=255)
       aGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTGTTATCAgct                                                                      <  1:422332/90‑1 (MQ=255)
                  cTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTGTTATCAGCTGCTTTTTgcgc                                                           >  1:210563/1‑90 (MQ=255)
                   tGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCa                                                          <  2:210563/90‑1 (MQ=255)
                      aGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATcc                                                       <  2:162969/90‑1 (MQ=255)
                            cTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCa                                                 >  2:254908/1‑90 (MQ=255)
                             tGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCAt                                                <  2:547155/90‑1 (MQ=255)
                             tGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCAt                                                <  2:287170/90‑1 (MQ=255)
                             tGGAAGGACTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCAt                                                <  1:24357/90‑1 (MQ=255)
                                  ggATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTc                                           <  2:533521/90‑1 (MQ=255)
                                  ggATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTc                                           <  2:270369/90‑1 (MQ=255)
                                    aTTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTcac                                         <  2:85550/90‑1 (MQ=255)
                                             caGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCag                                >  1:338905/1‑90 (MQ=255)
                                                  aGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCAt                                                <  2:557956/69‑1 (MQ=255)
                                                  aGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCAt                                                >  1:557956/1‑69 (MQ=255)
                                                  aGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAgaagaa                           >  2:117601/1‑90 (MQ=255)
                                                        tCTTCACAGGTCCAGGAGTTGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGac                     <  2:282203/90‑1 (MQ=255)
                                                             acaGGTCCAGGAGTTGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGAccacca                <  1:89275/90‑1 (MQ=255)
                                                             acaGGTCCAGGAGTTGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGAccacca                <  2:192802/90‑1 (MQ=255)
                                                                      ggAGTTGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGACCACCAAGAATAACg       <  2:338905/90‑1 (MQ=255)
                                                                        aGTTGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGACCACCAAGAATAACGcc     >  1:491911/1‑90 (MQ=255)
                                                                          ttGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTca                                          <  2:267354/51‑1 (MQ=255)
                                                                          ttGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTca                                          >  1:267354/1‑51 (MQ=255)
                                                                           tGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTgg                                   >  1:76032/1‑57 (MQ=255)
                                                                           tGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTgg                                   <  2:76032/57‑1 (MQ=255)
                                                                           tGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCaga                               <  2:389168/61‑1 (MQ=255)
                                                                           tGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCaga                               >  1:389168/1‑61 (MQ=255)
                                                                           tGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGACCACCAAGAATAACGCCTaa  >  2:161655/1‑90 (MQ=255)
                                                                           tGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGACCACCAAGAATAACGCCTaa  >  2:299869/1‑90 (MQ=255)
                                                                           tGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGACCACCAAGAATAACGCCTaa  >  1:84212/1‑90 (MQ=255)
                                                                           tGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGACCACCAAGAATAACGCCTaa  >  1:17592/1‑90 (MQ=255)
                                                                                  |                                                                                   
GCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGACCACCAAGAATAACGCCTAA  >  NZ_CP009273/3649923‑3650088

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ 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
TTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGACCACCAAGAATAACGCCTAATACTTTTTTCATCGTAATAT  >  NZ_CP009273/3649913‑3650108
                                                                                                 |                                                                                                  
TTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTGTT                                                                                                 >  SRR3721989.123341/1‑100 (MQ=60)
     GTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTGTTATCAG                                                                                            >  SRR3721989.275187/1‑100 (MQ=60)
                 ccaccgagatctacactccttacgtcgtcggcagcgtcagatgtgtataagagacagGGTCCAGGAGTTGACCGGTTTTTGTTATCAGCTGCTTTTTGCG                                                                                <  SRR3721989.581373/43‑1 (MQ=60)
                 AGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTGTTATCAGCTGCTTTTTGCG                                                                                <  SRR3721989.431053/100‑1 (MQ=60)
                  GCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTGTTATCAGCTGCTTTTTGCGC                                                                               >  SRR3721989.214381/1‑100 (MQ=60)
                                       TGGAAGGACTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAAC                                                          <  SRR3721989.24776/100‑1 (MQ=60)
                                             GATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAG                                                    >  SRR3721989.345386/1‑100 (MQ=60)
                                                  GTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACctgtctcttat                                               >  SRR3721989.570059/1‑89 (MQ=60)
                                                                       ACAGGTCCAGGAGTTGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGACCACCAAGAATAACGC                          <  SRR3721989.90889/100‑1 (MQ=60)
                                                                        CAGGTCCAGGAGTTGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGACCACCAAGAATAACGCC                         >  SRR3721989.502316/1‑100 (MQ=60)
                                                                          GGTCCAGGAGTTGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGctgtctcttatacacatctgacgctgccg                       >  SRR3721989.272099/1‑71 (MQ=60)
                                                                           GTCCAGGAGTTGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATctgtctcttatacacatctgacgctgccgacgacgtaaggagtgta                      >  SRR3721989.158078/1‑54 (MQ=60)
                                                                           GTCCAGGAGTTGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCctgtctcttatacacatctgacg                      >  SRR3721989.77387/1‑77 (MQ=60)
                                                                           GTCCAGGAGTTGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGACCtgtctcttatacacatct                      >  SRR3721989.397019/1‑82 (MQ=60)
                                                                           GTCCAGGAGTTGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGACCACCAAGAATAACGCCTAA                      >  SRR3721989.17900/1‑100 (MQ=60)
                                                                           GTCCAGGAGTTGACCGGTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGACCACCAAGAATAACGCCTAA                      >  SRR3721989.85733/1‑100 (MQ=60)
                                                                                                TTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGACCACCAAGAATAACGCCTAATACTTTTTTCATCGTAATAT  <  SRR3721989.164592/100‑1 (MQ=60)
                                                                                                 |                                                                                                  
TTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGACCACCAAGAATAACGCCTAATACTTTTTTCATCGTAATAT  >  NZ_CP009273/3649913‑3650108

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