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,015,838 G→C 24.8% P281P (CCG→CCC)  metC → cystathionine beta‑lyase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,015,8380GC24.8% 200.8 / 67.6 125P281P (CCG→CCC) metCcystathionine beta‑lyase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (14/80);  new base C (19/12);  total (33/92)
Fisher's exact test for biased strand distribution p-value = 1.57e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 125 total aligned reads displayed.

AGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAA  >  minE/2015772‑2015897
                                                                  |                                                           
aGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCa                                                           <  1:2229207/69‑1 (MQ=255)
aGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCa                                                           >  1:1704759/1‑69 (MQ=255)
         tttGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGcgcg                                                  <  1:60189/69‑1 (MQ=255)
              gTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGtt                                              <  1:500730/68‑1 (MQ=255)
              gTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGtt                                              <  1:443798/68‑1 (MQ=255)
              gTCAACATCATGAAAGCAGTATGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGtt                                              <  1:3399390/68‑1 (MQ=255)
                cAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAac                                           <  1:1086066/69‑1 (MQ=255)
                cAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAac                                           <  1:2915920/69‑1 (MQ=255)
                cAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAac                                           <  1:1208876/69‑1 (MQ=255)
                cAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAac                                           <  1:3042189/69‑1 (MQ=255)
                cAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAac                                           <  1:3059581/69‑1 (MQ=255)
                cAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAac                                           <  1:2509861/69‑1 (MQ=255)
                cAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAac                                           <  1:136949/69‑1 (MQ=255)
                cAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAac                                           <  1:1630664/69‑1 (MQ=255)
                cAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAac                                           <  1:2361095/69‑1 (MQ=255)
                cAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCCAGTTAac                                           <  1:1038301/69‑1 (MQ=255)
                   catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAaccac                                        <  1:1519952/69‑1 (MQ=255)
                   catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAaccac                                        <  1:708550/69‑1 (MQ=255)
                   catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAaccac                                        <  1:2424879/69‑1 (MQ=255)
                   catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAaccac                                        <  1:3399263/69‑1 (MQ=255)
                   catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAaccac                                        <  1:3321296/69‑1 (MQ=255)
                   catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAaccac                                        <  1:1236079/69‑1 (MQ=255)
                   catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAaccac                                        <  1:990212/69‑1 (MQ=255)
                   catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAacca                                         <  1:3429090/68‑1 (MQ=255)
                   catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAacca                                         <  1:359847/68‑1 (MQ=255)
                   catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAacca                                         <  1:2551652/68‑1 (MQ=255)
                   catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAacca                                         <  1:2176321/68‑1 (MQ=255)
                   catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAacca                                         <  1:2923755/68‑1 (MQ=255)
                   catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAacca                                         <  1:488700/68‑1 (MQ=255)
                       atGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTg                                    >  1:408373/1‑69 (MQ=255)
                         gAAAGCAGTCTGAAAGTGGCTGAATGGCTTGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc                                   >  1:1294104/1‑68 (MQ=255)
                         gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGc                                   <  1:2971152/68‑1 (MQ=255)
                         gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGc                                   <  1:628064/68‑1 (MQ=255)
                         gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGc                                   <  1:2121903/68‑1 (MQ=255)
                         gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGc                                   <  1:1635444/68‑1 (MQ=255)
                         gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGc                                   <  1:3183957/68‑1 (MQ=255)
                         gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGc                                   <  1:1414412/68‑1 (MQ=255)
                         gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGc                                   <  1:2479816/68‑1 (MQ=255)
                         gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGtt                                              >  1:2487544/1‑57 (MQ=255)
                         gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTg                                    >  1:2303020/1‑67 (MQ=255)
                         gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc                                   >  1:1094467/1‑68 (MQ=255)
                         gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc                                   >  1:283834/1‑68 (MQ=255)
                         gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc                                   >  1:3011618/1‑68 (MQ=255)
                         gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc                                   >  1:2083284/1‑68 (MQ=255)
                         gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc                                   >  1:1305925/1‑68 (MQ=255)
                         gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc                                   >  1:1340347/1‑68 (MQ=255)
                         gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc                                   >  1:134188/1‑68 (MQ=255)
                         gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc                                   >  1:2484602/1‑68 (MQ=255)
                         gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc                                   >  1:758152/1‑68 (MQ=255)
                         gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc                                   >  1:2315036/1‑68 (MQ=255)
                         gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc                                   >  1:1655359/1‑68 (MQ=255)
                         gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc                                   >  1:1692589/1‑68 (MQ=255)
                         gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc                                   >  1:2203668/1‑68 (MQ=255)
                                  cTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAacca                                         <  1:1854868/53‑1 (MQ=255)
                                   tGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGcgcg                                                  >  1:836475/1‑43 (MQ=255)
                                       aGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTaa                    <  1:1491898/69‑1 (MQ=255)
                                       aGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTaa                    <  1:3318333/69‑1 (MQ=255)
                                       aGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTaa                    <  1:2706342/69‑1 (MQ=255)
                                       aGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTaa                    <  1:546135/69‑1 (MQ=255)
                                        gTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTaa                    <  1:1060870/68‑1 (MQ=255)
                                        gTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTaa                    <  1:2437436/68‑1 (MQ=255)
                                          ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt                                     <  1:2931938/49‑1 (MQ=255)
                                          ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt                                     <  1:1162778/49‑1 (MQ=255)
                                          ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt                                     <  1:982790/49‑1 (MQ=255)
                                          ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt                                     <  1:745524/49‑1 (MQ=255)
                                          ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt                                     <  1:737135/49‑1 (MQ=255)
                                          ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt                                     <  1:637286/49‑1 (MQ=255)
                                          ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt                                     <  1:2266201/49‑1 (MQ=255)
                                          ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt                                     <  1:322791/49‑1 (MQ=255)
                                          ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt                                     <  1:3157357/49‑1 (MQ=255)
                                          ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt                                     <  1:2928959/49‑1 (MQ=255)
                                          ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt                                     <  1:3069843/49‑1 (MQ=255)
                                          ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt                                     <  1:2970741/49‑1 (MQ=255)
                                          ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt                                     <  1:2911939/49‑1 (MQ=255)
                                          ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt                                     <  1:2751695/49‑1 (MQ=255)
                                          ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt                                     <  1:273948/49‑1 (MQ=255)
                                                   gctggcAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAAtt        >  1:3300975/1‑69 (MQ=255)
                                                   gctggcAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAAtt        >  1:3131526/1‑69 (MQ=255)
                                                   gctggcAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGCCACGaa          <  1:446231/67‑1 (MQ=255)
                                                   gctggcAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGCCACGaa          <  1:746635/67‑1 (MQ=255)
                                                   gctggcAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGCCACGaa          <  1:1587551/67‑1 (MQ=255)
                                                   gctggcAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGCCACGaa          <  1:561828/67‑1 (MQ=255)
                                                   gctggcAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGCCACGaa          <  1:2038803/67‑1 (MQ=255)
                                                   gctggcAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGCCACGaa          <  1:2216875/67‑1 (MQ=255)
                                                   gctggcAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGCCACGaa          <  1:230316/67‑1 (MQ=255)
                                                   gctggcAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGCCACGaa          <  1:2641345/67‑1 (MQ=255)
                                                   gctggcAGAACATCCCCAAGTAGCGCGAGTTAAACACCCTGCTCTGCCTGGCAGTAAAGGCCACGaa          <  1:2519435/67‑1 (MQ=255)
                                                   gctggcAGAACATCCCCAAGTAGCGCGACTTAACCACCCTGCTCTGCCTGGCAGTAAAGGCCACGaa          <  1:2296054/67‑1 (MQ=255)
                                                         aGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGaa  <  1:1247180/69‑1 (MQ=255)
                                                         aGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGaa  <  1:1889297/69‑1 (MQ=255)
                                                         aGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGaa  <  1:963995/69‑1 (MQ=255)
                                                         aGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGaa  <  1:1867486/69‑1 (MQ=255)
                                                         aGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGaa  <  1:1754557/69‑1 (MQ=255)
                                                         aGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGaa  <  1:3040452/69‑1 (MQ=255)
                                                         aGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGaa  <  1:157391/69‑1 (MQ=255)
                                                         aGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGa   <  1:2009794/68‑1 (MQ=255)
                                                            aCATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGaa  >  1:444014/1‑66 (MQ=255)
                                                            aCATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGaa  >  1:2231118/1‑66 (MQ=255)
                                                            aCATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGaa  >  1:2502209/1‑66 (MQ=255)
                                                            aCATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGaa  >  1:2524745/1‑66 (MQ=255)
                                                                  |                                                           
AGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAA  >  minE/2015772‑2015897

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

GATK/CNVnator alignment

N/A