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 I199 R1 220 57.5 3163680 97.0% 3068769 86.1

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,7310AG97.6% 124.7 / ‑3.6 41D40G (GAT→GGT) entFenterobactin non‑ribosomal peptide synthetase EntF
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (25/15);  total (25/16)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

CGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGATTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGATAACGGCGAA  >  NZ_CP009273/609648‑609816
                                                                                   |                                                                                     
cGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGcc                                                                                 >  1:518998/1‑90 (MQ=255)
cGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGcc                                                                                 >  2:114261/1‑90 (MQ=255)
cGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGcc                                                                                 <  2:1465265/90‑1 (MQ=255)
  gCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCAt                                                                               >  1:1218746/1‑90 (MQ=255)
  gCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCAt                                                                               >  1:13290/1‑90 (MQ=255)
       tGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTg                                                                          >  1:1560397/1‑90 (MQ=255)
       tGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTg                                                                          >  1:475721/1‑90 (MQ=255)
       tGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTg                                                                          >  2:1442568/1‑90 (MQ=255)
             gCAGAAAAACTGTCAGAATTACCCACCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCcgc                                                                    >  1:306652/1‑90 (MQ=255)
               aGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCcgcgc                                                                  <  2:331031/90‑1 (MQ=255)
               aGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCcgcgc                                                                  <  2:363599/90‑1 (MQ=255)
                  aaaaCTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCggt                                                               >  2:374451/1‑90 (MQ=255)
                  aaaaCTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCggt                                                               >  2:1028580/1‑90 (MQ=255)
                    aaCTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCggtgg                                                             <  1:667987/90‑1 (MQ=255)
                                 aCCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAgc                                                <  2:778363/90‑1 (MQ=255)
                                   ccTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAgcgc                                              <  1:255151/90‑1 (MQ=255)
                                     tCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCaa                                            <  2:627638/90‑1 (MQ=255)
                                      ccGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAg                                           >  2:932371/1‑90 (MQ=255)
                                      ccGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAg                                           >  1:631660/1‑90 (MQ=255)
                                       cGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGc                                          <  2:1245099/90‑1 (MQ=255)
                                            ggAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATa                                     >  2:1458539/1‑90 (MQ=255)
                                             gAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGATTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATAc                                    <  1:772927/90‑1 (MQ=255)
                                                cGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCt                                 >  1:895882/1‑90 (MQ=255)
                                                cGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCt                                 >  2:1317289/1‑90 (MQ=255)
                                                 gTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGAc                                                    <  1:1522679/70‑1 (MQ=255)
                                                 gTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGAc                                                    >  2:1522679/1‑70 (MQ=255)
                                                 gTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTg                                >  2:313493/1‑90 (MQ=255)
                                                  tGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGc                               >  1:1285479/1‑90 (MQ=255)
                                                     cgcATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGa                            <  1:1442568/90‑1 (MQ=255)
                                                       cATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATg                          <  2:1573299/90‑1 (MQ=255)
                                                         ttACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCg                        >  2:1260299/1‑90 (MQ=255)
                                                              ttGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTa                   >  1:1270263/1‑90 (MQ=255)
                                                               tGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTAc                  <  2:217716/90‑1 (MQ=255)
                                                               tGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCCCAAGCAGATACGCTGCGGATGCGTTTTAc                  <  1:684262/90‑1 (MQ=255)
                                                                gAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACg                 >  2:1548912/1‑90 (MQ=255)
                                                                 aGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACgg                <  1:1564745/90‑1 (MQ=255)
                                                                 aGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACgg                <  2:1357517/90‑1 (MQ=255)
                                                                     aaCCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGa            >  2:1290049/1‑90 (MQ=255)
                                                                     aaCCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGa            >  1:1410874/1‑90 (MQ=255)
                                                                     aaCCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGa            >  1:478343/1‑90 (MQ=255)
                                                                              ggttggttCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGATAACGGCGa   <  1:1557634/90‑1 (MQ=255)
                                                                              ggttggttCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGATAACGGCGa   <  2:1293391/90‑1 (MQ=255)
                                                                              ggttggttCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGATAACGGCGa   <  1:236016/90‑1 (MQ=255)
                                                                              ggttggttCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGATAACGGCGa   <  1:1290049/90‑1 (MQ=255)
                                                                              ggttggttCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGATAACGGCGa   <  1:67021/90‑1 (MQ=255)
                                                                               gTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGATAACGGCGaa  >  2:1330265/1‑90 (MQ=255)
                                                                                   |                                                                                     
CGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGATTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGATAACGGCGAA  >  NZ_CP009273/609648‑609816

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ 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
GGTCGCCGCACAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGATTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGATAACGGCGAAGTCTGGCAG  >  NZ_CP009273/609633‑609825
                                                                                                  |                                                                                              
GGTCGCCGCACAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGT                                                                                               <  SRR3722078.1202372/100‑1 (MQ=60)
     CCGCACAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCC                                                                                          >  SRR3722078.525141/1‑100 (MQ=60)
       GCACAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCAT                                                                                        >  SRR3722078.1237917/1‑100 (MQ=60)
       GCACAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCAT                                                                                        >  SRR3722078.13451/1‑100 (MQ=60)
            GCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTG                                                                                   >  SRR3722078.481365/1‑100 (MQ=60)
            GCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTG                                                                                   >  SRR3722078.1585791/1‑100 (MQ=60)
                  CATCTGGATGGCAGAAAAACTGTCAGAATTACCCACCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGC                                                                             >  SRR3722078.310396/1‑100 (MQ=60)
                                   AACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACT                                                            <  SRR3722078.675626/100‑1 (MQ=60)
                                           GAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAG                                                    >  SRR3722078.638927/1‑100 (MQ=60)
                                            taagagacagCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGC                                                   <  SRR3722078.1547500/90‑1 (MQ=60)
                                                  CCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATAC                                             <  SRR3722078.258275/100‑1 (MQ=60)
                                                     CCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCT                                          >  SRR3722078.907916/1‑100 (MQ=60)
                                                       GCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGC                                        >  SRR3722078.1305934/1‑100 (MQ=60)
                                                            GAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGATTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATG                                   <  SRR3722078.782909/100‑1 (MQ=60)
                                                                   GCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTA                            >  SRR3722078.1290415/1‑100 (MQ=60)
                                                                    CGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTAC                           <  SRR3722078.1466061/100‑1 (MQ=60)
                                                                          ACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGA                     >  SRR3722078.1433765/1‑100 (MQ=60)
                                                                          ACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGA                     >  SRR3722078.484017/1‑100 (MQ=60)
                                                                              TGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCCCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGATAAC                 <  SRR3722078.692268/100‑1 (MQ=60)
                                                                                AGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGATAACGG               <  SRR3722078.1590214/100‑1 (MQ=60)
                                                                                             GGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGATAACGGCGAAGTCTGGCAG  <  SRR3722078.238890/100‑1 (MQ=60)
                                                                                             GGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGATAACGGCGAAGTCTGGCAG  <  SRR3722078.1582987/100‑1 (MQ=60)
                                                                                             GGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGATAACGGCGAAGTCTGGCAG  <  SRR3722078.409427/100‑1 (MQ=60)
                                                                                             GGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGATAACGGCGAAGTCTGGCAG  <  SRR3722078.1310584/100‑1 (MQ=60)
                                                                                             GGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGATAACGGCGAAGTCTGGCAG  <  SRR3722078.617394/100‑1 (MQ=60)
                                                                                             GGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGATAACGGCGAAGTCTGGCAG  <  SRR3722078.67803/100‑1 (MQ=60)
                                                                                                  |                                                                                              
GGTCGCCGCACAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGATTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGATAACGGCGAAGTCTGGCAG  >  NZ_CP009273/609633‑609825

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