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 501,960 G→C 40.0% L139L (CTG→CTC)  modA → molybdate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE501,9600GC40.0% 31.5 / 57.3 60L139L (CTG→CTC) modAmolybdate transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (18/18);  new base C (14/10);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

GACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAAAGAAGCACTGCAAA  >  minE/501916‑502018
                                            |                                                          
gACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCg                                                  <  1:2178921/55‑1 (MQ=255)
gACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCg                                                  <  1:1151964/55‑1 (MQ=255)
gACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCg                                                  <  1:1653938/55‑1 (MQ=255)
gACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGtt                                    >  1:1220561/1‑69 (MQ=255)
gACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGtt                                    >  1:1912252/1‑69 (MQ=255)
gACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGtt                                    <  1:2088448/69‑1 (MQ=255)
gACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGtt                                    <  1:1075403/69‑1 (MQ=255)
gACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGtt                                    >  1:2619027/1‑69 (MQ=255)
gACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGtt                                    >  1:2970461/1‑69 (MQ=255)
gACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGtt                                    <  1:3113121/69‑1 (MQ=255)
gACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGtt                                    <  1:121057/69‑1 (MQ=255)
gACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGtt                                    >  1:1193379/1‑69 (MQ=255)
gACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGtt                                    >  1:32705/1‑69 (MQ=255)
gACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGtt                                    <  1:3420608/69‑1 (MQ=255)
gACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGtt                                    <  1:115067/69‑1 (MQ=255)
gACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGtt                                    <  1:1079247/69‑1 (MQ=255)
              cTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGAGCCGGAACATGtt                                    <  1:2994035/55‑1 (MQ=255)
                      cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGa                                           >  1:3425117/3‑40 (MQ=37)
                      cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                >  1:3119218/3‑67 (MQ=255)
                      cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                >  1:215007/3‑67 (MQ=255)
                      cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                >  1:2758511/3‑67 (MQ=255)
                      cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                >  1:2715456/3‑67 (MQ=255)
                      cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                >  1:445187/3‑67 (MQ=255)
                      cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                >  1:485498/3‑67 (MQ=255)
                      cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                >  1:1681129/3‑67 (MQ=255)
                      cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                >  1:1762549/3‑67 (MQ=255)
                      cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                >  1:1730154/3‑67 (MQ=255)
                      cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                >  1:1721745/3‑67 (MQ=255)
                      cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                >  1:1109510/3‑67 (MQ=255)
                      cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                >  1:1186386/3‑67 (MQ=255)
                      cgctgctgAAAGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                >  1:1395921/3‑67 (MQ=255)
                      cACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCaaa                >  1:1101447/1‑67 (MQ=255)
                      cACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCaaa                >  1:2756262/1‑67 (MQ=255)
                      cACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCaaa                >  1:1595590/1‑67 (MQ=255)
                      cACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCaaa                >  1:3196363/1‑67 (MQ=255)
                      cACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCaaa                >  1:1699066/1‑67 (MQ=255)
                      cACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCaaa                >  1:604149/1‑67 (MQ=255)
                      cACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCaaa                >  1:761543/1‑67 (MQ=255)
                      cACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCaaa                >  1:1778036/1‑67 (MQ=255)
                       aCTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAaaga             >  1:1885744/1‑69 (MQ=255)
                       aCTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAaaga             >  1:150431/1‑69 (MQ=255)
                       aCTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAaaga             >  1:1622763/1‑69 (MQ=255)
                       aCTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAaaga             >  1:2077611/1‑69 (MQ=255)
                         tgctgAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAaagaag           <  1:281908/69‑1 (MQ=255)
                         tgctgAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAaagaag           <  1:2418655/69‑1 (MQ=255)
                         tgctgAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAaagaac           <  1:379301/69‑2 (MQ=255)
                         tgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTTCTGGCATTTATGCGAaagaag           <  1:2386204/69‑1 (MQ=255)
                         tgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCGAaagaag           <  1:2568018/69‑1 (MQ=255)
                         tgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCGAaagaag           <  1:2532615/69‑1 (MQ=255)
                         tgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCGAaagaag           <  1:2451247/69‑1 (MQ=255)
                         tgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCGAaagaag           <  1:1232763/69‑1 (MQ=255)
                         tgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCGAaagaag           <  1:1205917/69‑1 (MQ=255)
                         tgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCGAaagaag           <  1:1976817/69‑1 (MQ=255)
                         tgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCGAaagaag           <  1:2347321/69‑1 (MQ=255)
                         tgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCGAaagaag           <  1:1105023/69‑1 (MQ=255)
                         tgcagAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCGAaagaag           <  1:261820/69‑1 (MQ=255)
                                  gctgTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAAAGAAGCACTGCaaa  <  1:3384187/66‑1 (MQ=255)
                                  gCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAAAGAAGCACTGCaaa  <  1:2247919/69‑1 (MQ=255)
                                  gCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAAAGAAGCACTGCaaa  <  1:1167438/69‑1 (MQ=255)
                                  gCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAAAGAAGCACTGCaaa  <  1:915316/69‑1 (MQ=255)
                                            |                                                          
GACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAAAGAAGCACTGCAAA  >  minE/501916‑502018

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

GATK/CNVnator alignment

N/A