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 609,731 A→G D40G (GAT→GGT)  entF → enterobactin non‑ribosomal peptide synthetase EntF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273609,7310AG100.0% 63.9 / NA 19D40G (GAT→GGT) entFenterobactin non‑ribosomal peptide synthetase EntF
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (9/10);  total (9/10)

AGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGATTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGATAACGGCGA  >  NZ_CP009273/609644‑609815
                                                                                       |                                                                                    
aGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAggttggtt                                                                                    <  1:64124/90‑1 (MQ=255)
  cccGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCg                                                                                  <  2:548891/90‑1 (MQ=255)
    cGGCATCTGGATGGCAGAAAAACTGTCAGAATTGCCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGcc                                                                                <  1:61755/90‑1 (MQ=255)
    cGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGcc                                                                                <  1:156920/90‑1 (MQ=255)
                ggCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCcg                                                                    >  2:544230/1‑90 (MQ=255)
                        aaCTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCggtgg                                                            <  1:542642/90‑1 (MQ=255)
                                      cccTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGcc                                                                                >  2:399152/1‑56 (MQ=255)
                                      cccTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGcc                                                                                <  1:399152/56‑1 (MQ=255)
                                          ccGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAg                                          >  1:648265/1‑90 (MQ=255)
                                          ccGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAg                                          >  2:123600/1‑90 (MQ=255)
                                           cGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGc                                         >  1:582535/1‑90 (MQ=255)
                                                 gAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATAc                                   <  2:189548/90‑1 (MQ=255)
                                                   gCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGc                                 <  1:48178/90‑1 (MQ=255)
                                                      tGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGc                              >  2:679604/1‑90 (MQ=255)
                                                           cATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATg                         <  1:343798/90‑1 (MQ=255)
                                                                  ttGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTa                  >  2:252109/1‑90 (MQ=255)
                                                                   tGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTAc                 <  1:765580/90‑1 (MQ=255)
                                                                     aGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACgg               <  2:449310/90‑1 (MQ=255)
                                                                     aGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACgg               >  1:209305/1‑90 (MQ=255)
                                                                       ttAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGaa             >  2:312074/1‑90 (MQ=255)
                                                                                  ggttggttCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGATAACGGCGa  <  1:139301/90‑1 (MQ=255)
                                                                                  ggttggttCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGATAACGGCGa  <  1:677976/90‑1 (MQ=255)
                                                                                  ggttggttCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGATAACGGCGa  >  1:202076/1‑90 (MQ=255)
                                                                                       |                                                                                    
AGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGATTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGATAACGGCGA  >  NZ_CP009273/609644‑609815

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 32 ≤ 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
TCGCCGCACAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGATTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGATAACGGCGAAGTCTGGCAG  >  NZ_CP009273/609635‑609825
                                                                                                |                                                                                              
TCGCCGCACAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTC                                                                                             <  SRR3722073.772786/100‑1 (MQ=60)
         AGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACT                                                                                    <  SRR3722073.65629/100‑1 (MQ=60)
             gcgtcagatgtgtataagagacagGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCC                                                                                <  SRR3722073.407648/76‑1 (MQ=60)
             CGGCATCTGGATGGCAGAAAAACTGTCAGAATTGCCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCC                                                                                <  SRR3722073.63204/100‑1 (MQ=60)
             CGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCC                                                                                <  SRR3722073.160097/100‑1 (MQ=60)
                                 AACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACT                                                            <  SRR3722073.553981/100‑1 (MQ=60)
                                         GAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAG                                                    >  SRR3722073.661609/1‑100 (MQ=60)
                                          AATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGC                                                   >  SRR3722073.594692/1‑100 (MQ=60)
                                                            GCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCG                                 <  SRR3722073.49311/100‑1 (MQ=60)
                                                                    CATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGG                         >  SRR3722073.213593/1‑100 (MQ=60)
                                                                    CATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGG                         <  SRR3722073.351081/100‑1 (MQ=60)
                                                                            TGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGATAAC                 <  SRR3722073.781035/100‑1 (MQ=60)
                                                                                 TAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGATAACGGCGA            >  SRR3722073.206171/1‑100 (MQ=60)
                                                                                           GGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGATAACGGCGAAGTCTGGCAG  <  SRR3722073.142197/100‑1 (MQ=60)
                                                                                           GGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGATAACGGCGAAGTCTGGCAG  <  SRR3722073.691831/100‑1 (MQ=60)
                                                                                                |                                                                                              
TCGCCGCACAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGATTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGATAACGGCGAAGTCTGGCAG  >  NZ_CP009273/609635‑609825

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