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
A16 F60 I1 R1 14 59.5 2012980 98.7% 1986811 139.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,517,267 A→G W151R (TGG→CGG)  fecR ← regulator for fec operon, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,517,2670AG100.0% 99.6 / NA 28W151R (TGG→CGG) fecRregulator for fec operon, periplasmic
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (14/14);  total (14/14)

CTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCG  >  NC_000913/4517136‑4517396
                                                                                                                                   |                                                                                                                                 
cTGTGTGAAATTATCGTGCTGGTGGTCGGCAAATTCTTTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCGGAGCCGGa                                                                                                                           <  1:467406/140‑1 (MQ=255)
cTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCGGAGCCGGa                                                                                                                           <  2:197167/140‑1 (MQ=255)
cTGCGTGAAATTATCCAGCTGGCGGACGGTAAATTTTGTACCTAAAGCTGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCGGAGCCGGa                                                                                                                           <  2:465480/140‑1 (MQ=255)
    gTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGCTAGTAATGGCGCTTGCACCGTACCGGAGCCGGACGGt                                                                                                                       >  1:469281/1‑140 (MQ=255)
    gTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCATGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCGGAGCCGGACGGt                                                                                                                       >  1:736850/1‑140 (MQ=255)
                                  tCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCGGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATccgc                                                                                         <  1:641469/140‑1 (MQ=255)
                                       cccTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCGGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCgcgc                                                                                    <  2:591609/140‑1 (MQ=255)
                                       cccTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCGGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCgcgc                                                                                    <  2:681963/140‑1 (MQ=255)
                                       cccTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCGGAGCCGGACGGTGCGCTGATGAGCTTCAAAAGGCACATCCGCCgcgc                                                                                    >  1:304368/1‑140 (MQ=255)
                                                 gTGAGCTGCCCCTGACCAGTCATGACACGGCGGAGACGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCGGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTAtt                                                                          <  2:469281/140‑1 (MQ=255)
                                                    aGCTGGCCCTGACGGGGCAGGAAACGGAAGCGGGGTTGCAGGGCATCTGTCGCGGTAGTAATGGCGATTTCACCGTACCGGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAg                                                                       <  2:458993/140‑1 (MQ=255)
                                                       tGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCGGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGt                                                                    <  1:789321/140‑1 (MQ=255)
                                                              gACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCGGAGCCGGACGGTGCGCTGATGGGCATCAAATCGCTCCTACGCCGCGCTTTGGCCATTCAGCGTGAGCAgg                                                             >  2:454006/1‑140 (MQ=255)
                                                              gACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCGGAGCCGGACGGTGCGCTGATGCGCATCAAAACACACATCCGCCGGGCTTTGGGTATTCAGCGGGAGCAgg                                                             >  1:321340/1‑140 (MQ=255)
                                                                cGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCGGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGa                                                           >  1:623827/1‑140 (MQ=255)
                                                                 gggTCAGGACACGGAAGGGGCGTTGCAGGGCTTCTTTCGCGGTAGTAATGGCGATTTCACCGTACCGGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAg                                                          <  2:736850/140‑1 (MQ=255)
                                                                 gggTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCGGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAg                                                          <  2:41573/140‑1 (MQ=255)
                                                                   gTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCGGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGcc                                                        >  2:182945/1‑140 (MQ=255)
                                                                       ggACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCGGAGCCGGACGGCGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCt                                                    >  2:509747/1‑140 (MQ=255)
                                                                             ggAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCGGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGt                                              >  1:702879/1‑140 (MQ=255)
                                                                                 ggggCGTTGCAGGGCATCTTTCACGTTAGTAATGGCGATTTCACCGTACCGGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCt                                          <  1:454006/140‑1 (MQ=255)
                                                                                  gggCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCGGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTg                                         >  1:846434/1‑140 (MQ=255)
                                                                                    gCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCGGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGAc                                       >  2:906644/1‑140 (MQ=255)
                                                                                             ggCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCGGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGTGCCTTCGTCCAGTTGCTGACGGCTGACtg                              >  1:338066/1‑140 (MQ=255)
                                                                                                            agtaATGGCGATTTCACCGTACCGGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACt                               <  2:120359/124‑1 (MQ=255)
                                                                                                            agtaATGGCGATTTCACCGTACCGGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACt                               >  1:120359/1‑124 (MQ=255)
                                                                                                                       tttCACCGTACCGGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGcc    >  1:988124/1‑140 (MQ=255)
                                                                                                                       tttCACCGTACCGGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGcc    >  1:353083/1‑140 (MQ=255)
                                                                                                                         tCACCGTACCGGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCg  <  2:623827/140‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                 
CTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCG  >  NC_000913/4517136‑4517396

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

GATK/CNVnator alignment

N/A