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
A2 F1 I0 R1 4577 57.4 3436498 77.7% 2670158 63.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,902,939 C→T 41.8% A34A (GCG→GCA)  recB ← exonuclease V (RecBCD complex), beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,902,9390CT41.8% 33.5 / 61.0 86A34A (GCG→GCA) recBexonuclease V (RecBCD complex), beta subunit
Reads supporting (aligned to +/- strand):  ref base C (32/18);  new base T (26/10);  total (58/28)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01

GGCGGGGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTAAG  >  minE/1902884‑1902999
                                                       |                                                            
ggcgggGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAgg                                                 <  1:988879/69‑1 (MQ=255)
ggcgggGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAgg                                                 <  1:78656/69‑1 (MQ=255)
ggcgggGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAgg                                                 <  1:553775/69‑1 (MQ=255)
ggcgggGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAgg                                                 <  1:4962/69‑1 (MQ=255)
ggcgggGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAgg                                                 <  1:1375592/69‑1 (MQ=255)
ggcgggGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAgg                                                 <  1:2791531/69‑1 (MQ=255)
ggcgggGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAgg                                                 <  1:2509135/69‑1 (MQ=255)
                ggAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGcc                                    >  1:424370/1‑66 (MQ=255)
                ggAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGcc                                    >  1:343303/1‑66 (MQ=255)
                ggAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGcc                                    >  1:3356915/1‑66 (MQ=255)
                ggAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGcc                                    >  1:3327195/1‑66 (MQ=255)
                ggAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGcc                                    >  1:3003246/1‑66 (MQ=255)
                ggAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGcc                                    >  1:199557/1‑66 (MQ=255)
                   aCCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga                              >  1:2405765/1‑69 (MQ=255)
                   aCCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga                              >  1:2476041/1‑69 (MQ=255)
                   aCCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga                              >  1:2347071/1‑69 (MQ=255)
                   aCCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga                              >  1:1055290/1‑69 (MQ=255)
                   aCCGCCTAGCCCAAGTAACAGGCGCAAATAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga                              >  1:2511078/1‑69 (MQ=255)
                   aCCGCCTAGCCCAAGTAACAGGCGCAAATAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga                              >  1:1465924/1‑69 (MQ=255)
                   aCCGCCTAGCCCAAGTAACAGGCGCAAATAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga                              >  1:1718677/1‑69 (MQ=255)
                   aCCGCCTAGCCCAAGTAACAGGCGCAAATAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga                              >  1:1783011/1‑69 (MQ=255)
                   aCCGCCTAGCCCAAGTAACAGGCGCAAATAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga                              >  1:311617/1‑69 (MQ=255)
                   aCCGCCGAGCCCAAGTAACAGGCGCAAATAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga                              >  1:436901/1‑69 (MQ=255)
                     cGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCagag                             <  1:2555926/68‑1 (MQ=255)
                              cAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATc                    >  1:2200833/1‑68 (MQ=255)
                              cAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCa                   >  1:2942689/1‑69 (MQ=255)
                              cAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCa                   >  1:901665/1‑69 (MQ=255)
                              cAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCa                   >  1:2756831/1‑69 (MQ=255)
                              cAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCa                   >  1:931329/1‑69 (MQ=255)
                              cAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCa                   >  1:1037088/1‑69 (MQ=255)
                              cAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCa                   >  1:2500947/1‑69 (MQ=255)
                              cAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCa                   >  1:236655/1‑69 (MQ=255)
                              cAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCa                   >  1:164080/1‑69 (MQ=255)
                              cAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCa                   >  1:2157157/1‑69 (MQ=255)
                                aGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTc                        <  1:379912/62‑1 (MQ=255)
                                aGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTc                        <  1:312606/62‑1 (MQ=255)
                                aGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTc                        <  1:2557048/62‑1 (MQ=255)
                                aGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGATTc                        <  1:2649884/62‑1 (MQ=255)
                                  tAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGgcg               <  1:2090219/69‑1 (MQ=255)
                                  tAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGgcg               <  1:1952624/69‑1 (MQ=255)
                                  tAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGgcg               <  1:645660/69‑1 (MQ=255)
                                  tAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGgcg               <  1:2483563/69‑1 (MQ=255)
                                       ggCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCAc          <  1:1527602/69‑1 (MQ=255)
                                       ggCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCAc          <  1:2012254/69‑1 (MQ=255)
                                         cgcAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCAccc        >  1:3378495/1‑69 (MQ=255)
                                         cgcAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCAccc        >  1:3236411/1‑69 (MQ=255)
                                         cgcAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCAccc        >  1:1161508/1‑69 (MQ=255)
                                         cgcAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCAccc        >  1:3017733/1‑69 (MQ=255)
                                         cgcAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCAccc        >  1:1597890/1‑69 (MQ=255)
                                         cgcAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCAccc        >  1:2036836/1‑69 (MQ=255)
                                         cgcAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCAccc        >  1:979653/1‑69 (MQ=255)
                                         cgcAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCCGCAGAGGCTTCAATCAGGCGCTCAccc        >  1:3185470/1‑69 (MQ=255)
                                         cgcAAATAGAGAGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAAt                     >  1:1435927/1‑56 (MQ=255)
                                              aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga                              <  1:1474722/42‑1 (MQ=255)
                                              aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga                              <  1:1289084/42‑1 (MQ=255)
                                              aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga                              <  1:366484/42‑1 (MQ=255)
                                              aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga                              <  1:2228451/42‑1 (MQ=255)
                                              aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga                              <  1:2268229/42‑1 (MQ=255)
                                              aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga                              <  1:1706986/42‑1 (MQ=255)
                                              aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga                              <  1:2866754/42‑1 (MQ=255)
                                              aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga                              <  1:624881/42‑1 (MQ=255)
                                              aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga                              <  1:1812502/42‑1 (MQ=255)
                                              aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa   >  1:906316/1‑69 (MQ=255)
                                              aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa   >  1:192487/1‑69 (MQ=255)
                                              aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa   >  1:166284/1‑69 (MQ=255)
                                              aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa   >  1:1209308/1‑69 (MQ=255)
                                              aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa   >  1:1613812/1‑69 (MQ=255)
                                              aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa   >  1:1229520/1‑69 (MQ=255)
                                              aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa   >  1:627866/1‑69 (MQ=255)
                                              aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa   >  1:1944118/1‑69 (MQ=255)
                                              aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa   >  1:3338341/1‑69 (MQ=255)
                                              aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa   >  1:1396920/1‑69 (MQ=255)
                                              aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa   >  1:3302193/1‑69 (MQ=255)
                                              aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa   >  1:3301179/1‑69 (MQ=255)
                                              aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa   >  1:3117603/1‑69 (MQ=255)
                                              aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa   >  1:3077913/1‑69 (MQ=255)
                                              aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa   >  1:1975495/1‑69 (MQ=255)
                                              aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa   >  1:1035632/1‑69 (MQ=255)
                                              aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa   >  1:2402818/1‑69 (MQ=255)
                                              aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa   >  1:2259970/1‑69 (MQ=255)
                                              aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa   >  1:2172295/1‑69 (MQ=255)
                                              aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGATTCAATCAGGCGCTCACCCTGTaa   >  1:2122219/1‑69 (MQ=255)
                                              aTAGAGCGCTGCAATCGTAAAGGTCTTGCCTGTGCCGGCaga                              <  1:148723/42‑1 (MQ=255)
                                               tAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTAAg  >  1:1209134/1‑69 (MQ=255)
                                               tAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTAAg  >  1:2043541/1‑69 (MQ=255)
                                               tAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTAAg  >  1:222217/1‑69 (MQ=255)
                                                       |                                                            
GGCGGGGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTAAG  >  minE/1902884‑1902999

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑

GATK/CNVnator alignment

N/A