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 2,795,792 A→T 40.3% A94A (GCA→GCT)  dgkA → diacylglycerol kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,795,7920AT40.3% 65.2 / 116.8 119A94A (GCA→GCT) dgkAdiacylglycerol kinase
Reads supporting (aligned to +/- strand):  ref base A (30/41);  new base T (29/19);  total (59/60)
Fisher's exact test for biased strand distribution p-value = 6.28e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.73e-01

Only 100 of 119 total aligned reads displayed.

TAGCGCCATCGAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTAC  >  minE/2795729‑2795845
                                                               |                                                     
tAGCGCCATCGAAGCTGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCTAAAg                                                   <  1:2957704/68‑1 (MQ=255)
tAGCGCCATCGAAGCTGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCTAAAg                                                   <  1:2495264/68‑1 (MQ=255)
tAGCGCCATCGAAGCTGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCTAAAg                                                   <  1:3069381/68‑1 (MQ=255)
tAGCGCCATCGAAGCTGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCTAAAg                                                   <  1:586126/68‑1 (MQ=255)
tAGCGCCATCGAAGCTGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCTAAAg                                                   <  1:1512248/68‑1 (MQ=255)
tAGCGCCATCGAAGCTGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCTAAAg                                                   <  1:2117609/68‑1 (MQ=255)
tAGCGCCATCGAAGCTGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCTAAAg                                                   <  1:2091501/68‑1 (MQ=255)
tAGCGCCATCGAAGCTGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCTAAAg                                                   <  1:1810603/68‑1 (MQ=255)
tAGCGCCATCGAAGCTGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCTAAAg                                                   <  1:793641/68‑1 (MQ=255)
tAGCGCCATCGAAGCTGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCTAAAg                                                   <  1:1268530/68‑1 (MQ=255)
tAGCGCCATCGAAGCTGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCTAAAg                                                   <  1:1209668/68‑1 (MQ=255)
tAGCGCCATCGAAGCTGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCTAAAg                                                   <  1:2807841/68‑1 (MQ=255)
tAGCGCCATCGAAGCTGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCTAAAg                                                   <  1:3185270/68‑1 (MQ=255)
tAGCGCCATCGAAGCTGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCTAAAg                                                   <  1:2017959/68‑1 (MQ=255)
tAGCGCCATCGAAGCTGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACACGCTAAAg                                                   <  1:548848/68‑1 (MQ=255)
tAGCGCCATCGAAGCTGTGGGTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCTAAAg                                                   <  1:2166783/68‑1 (MQ=255)
tAGCGCCATCGAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGa                                                  <  1:3224604/69‑1 (MQ=255)
 aGCGCCATCGAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGGGCAAAAGat                                                 >  1:1743336/1‑69 (MQ=255)
 aGCGCCATCGAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGat                                                 >  1:642285/1‑69 (MQ=255)
 aGCGCCATCGAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGat                                                 >  1:3260999/1‑69 (MQ=255)
 aGCGCCATCGAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGat                                                 >  1:3010857/1‑69 (MQ=255)
 aGCGCCATCGAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGat                                                 >  1:1328533/1‑69 (MQ=255)
          gAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAg                                                   <  1:1736464/58‑1 (MQ=255)
          gAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAg                                                   <  1:2375253/58‑1 (MQ=255)
          gAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAg                                                   <  1:2325062/58‑1 (MQ=255)
          gAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAg                                                   <  1:1082414/58‑1 (MQ=255)
          gAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATcc                                        >  1:2057328/1‑69 (MQ=255)
          gAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATcc                                        >  1:3133073/1‑69 (MQ=255)
          gAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATcc                                        >  1:3367935/1‑69 (MQ=255)
          gAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATcc                                        >  1:1630753/1‑69 (MQ=255)
          gAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATcc                                        >  1:2569417/1‑69 (MQ=255)
          gAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATcc                                        >  1:927535/1‑69 (MQ=255)
          gAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATcc                                        >  1:1031045/1‑69 (MQ=255)
                  ggTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGt                                >  1:2431518/1‑69 (MQ=255)
                  ggTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGt                                >  1:1387072/1‑69 (MQ=255)
                  ggTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGt                                >  1:2232713/1‑69 (MQ=255)
                   gTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTg                               >  1:414339/1‑69 (MQ=255)
                   gTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTg                               >  1:1911691/1‑69 (MQ=255)
                     tGACCTAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCt                             <  1:2032651/69‑1 (MQ=255)
                     tGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCt                             <  1:2905930/69‑1 (MQ=255)
                     tGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCt                             <  1:2298835/69‑1 (MQ=255)
                     tGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCt                             <  1:692907/69‑1 (MQ=255)
                     tGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCt                             <  1:675164/69‑1 (MQ=255)
                     tGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCt                             <  1:1867340/69‑1 (MQ=255)
                     tGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCt                             <  1:428067/69‑1 (MQ=255)
                     tGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCt                             <  1:3394617/69‑1 (MQ=255)
                     tGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCt                             <  1:2176530/69‑1 (MQ=255)
                         cGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGAtt                         <  1:2596272/69‑1 (MQ=255)
                         cGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGAtt                         <  1:2702086/69‑1 (MQ=255)
                         cGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGAtt                         <  1:2716873/69‑1 (MQ=255)
                          gAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGAtt                         <  1:1669261/68‑1 (MQ=255)
                               ggCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCAtt                   <  1:3255607/69‑1 (MQ=255)
                               ggCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCAtt                   <  1:2667484/69‑1 (MQ=255)
                                  tctGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTAtc                <  1:1586865/69‑1 (MQ=255)
                                  tctGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTAtc                <  1:593750/69‑1 (MQ=255)
                                      aaTACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCAtt                   >  1:912675/1‑62 (MQ=255)
                                      aaTACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCAtt                   >  1:1068364/1‑62 (MQ=255)
                                      aaTACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCAtt                   >  1:2909098/1‑62 (MQ=255)
                                      aaTACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCAtt                   >  1:3045541/1‑62 (MQ=255)
                                      aaTACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCAtt                   >  1:3115381/1‑62 (MQ=255)
                                      aaTACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCAtt                   >  1:2726505/1‑62 (MQ=255)
                                      aaTACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCAtt                   >  1:2113389/1‑62 (MQ=255)
                                            aTGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTg       <  1:2594014/68‑1 (MQ=255)
                                            aTGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTg       <  1:3228104/68‑1 (MQ=255)
                                            aTGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTg       <  1:3232082/68‑1 (MQ=255)
                                            aTGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTg       <  1:2430836/68‑1 (MQ=255)
                                            aTGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTg       <  1:290007/68‑1 (MQ=255)
                                            aTGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTg       <  1:3000537/68‑1 (MQ=255)
                                            aTGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTg       <  1:207964/68‑1 (MQ=255)
                                            aTGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTg       <  1:128099/68‑1 (MQ=255)
                                            aTGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTg       <  1:823060/68‑1 (MQ=255)
                                            aTGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTg       <  1:347769/68‑1 (MQ=255)
                                             tGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTg       <  1:289025/67‑1 (MQ=255)
                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:1502563/1‑69 (MQ=255)
                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:1996311/1‑69 (MQ=255)
                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:1925583/1‑69 (MQ=255)
                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:547318/1‑69 (MQ=255)
                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:1565471/1‑69 (MQ=255)
                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:1521982/1‑69 (MQ=255)
                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:2116988/1‑69 (MQ=255)
                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:1492644/1‑69 (MQ=255)
                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:1367059/1‑69 (MQ=255)
                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:1331839/1‑69 (MQ=255)
                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:783812/1‑69 (MQ=255)
                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:1118890/1‑69 (MQ=255)
                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:3363858/1‑69 (MQ=255)
                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:3259407/1‑69 (MQ=255)
                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:231508/1‑69 (MQ=255)
                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:2359348/1‑69 (MQ=255)
                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:2476793/1‑69 (MQ=255)
                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:2548072/1‑69 (MQ=255)
                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:2582087/1‑69 (MQ=255)
                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:2990289/1‑69 (MQ=255)
                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:2780125/1‑69 (MQ=255)
                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:2904703/1‑69 (MQ=255)
                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTAATTGCCATTATCGTCGCCGTGAtt    >  1:388389/1‑69 (MQ=255)
                                              gAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGAtt    <  1:1155618/69‑1 (MQ=255)
                                               aGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTa   >  1:3146831/1‑69 (MQ=255)
                                               aGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTa   >  1:2480845/1‑69 (MQ=255)
                                                gCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCCGTGCTGATTGCCATTATCGTCGCCGTGATTAc  >  1:2957406/1‑69 (MQ=255)
                                                               |                                                     
TAGCGCCATCGAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTAC  >  minE/2795729‑2795845

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