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
A6 F21 I0 R1 2169 63.0 3952800 76.2% 3012033 62.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 609,783 G→C 19.2% I157M (ATC→ATG)  pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE609,7830GC19.2% 214.2 / 44.0 99I157M (ATC→ATG) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (47/33);  new base C (14/5);  total (61/38)
Fisher's exact test for biased strand distribution p-value = 2.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.19e-01

Only 100 of 103 total aligned reads displayed.

CCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGA  >  minE/609715‑609853
                                                                    |                                                                      
ccGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGAt                                                                      <  1:1036060/71‑1 (MQ=255)
 cGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATg                                                                     >  1:1591625/1‑71 (MQ=255)
  gatgatACGGCCACGGCCATATGGATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGt                                                                    >  1:3618208/1‑71 (MQ=255)
  gatgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGt                                                                    >  1:2490979/1‑71 (MQ=255)
  gataatGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCAtat                                                                    <  1:3795941/71‑3 (MQ=255)
   atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTc                                                                   >  1:8738/1‑71 (MQ=255)
   atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAAGAtag                                                                    >  1:1777735/1‑68 (MQ=255)
      atACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCgg                                                                <  1:2110042/71‑1 (MQ=255)
      atACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCg                                                                 <  1:264934/70‑1 (MQ=255)
       tACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCg                                                                 >  1:403082/1‑69 (MQ=255)
       tACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCg                                                                 >  1:1819363/1‑69 (MQ=255)
          ggccacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAg                                                              <  1:3154099/69‑1 (MQ=255)
            ccacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAgtgt                                                           >  1:2180042/1‑70 (MQ=255)
            ccacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAgt                                                             >  1:2812195/1‑68 (MQ=255)
             cacggccaTATGCATCTGTCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaa                                                         <  1:5535/71‑1 (MQ=255)
             cacggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaa                                                         <  1:3006292/71‑1 (MQ=255)
              acggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaaa                                                        >  1:1205796/1‑71 (MQ=255)
               cggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa                                                          >  1:2006983/1‑68 (MQ=255)
               cggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa                                                          <  1:1759126/68‑1 (MQ=255)
               cggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaaa                                                        >  1:2393813/1‑70 (MQ=255)
               cggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaaa                                                        >  1:2389073/1‑70 (MQ=255)
               cggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaaa                                                        >  1:2683014/1‑70 (MQ=255)
               cggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaaa                                                        >  1:3015164/1‑70 (MQ=255)
                tgccaGATGCATATGGCAGACAGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACg                                                      >  1:2894211/2‑71 (MQ=255)
                ggccaTATGCATCTGGCAGACAGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACg                                                      <  1:1154638/71‑1 (MQ=255)
                      atGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGa                                                               >  1:2681569/1‑56 (MQ=255)
                      atGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGa                                                               >  1:3667396/1‑56 (MQ=255)
                       tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                                                        >  1:89121/1‑62 (MQ=255)
                       tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                               <  1:1557305/71‑1 (MQ=255)
                       tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                               <  1:3863764/71‑1 (MQ=255)
                       tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                               <  1:3713503/71‑1 (MQ=255)
                        gCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGa                                                               >  1:2082483/1‑54 (MQ=255)
                        ccATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCggt                                                               >  1:1971649/2‑53 (MQ=16)
                         cATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACg                                             <  1:3294142/71‑1 (MQ=255)
                         cATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaa                                                         >  1:2722572/1‑59 (MQ=255)
                         cATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaa                                                         >  1:3389602/1‑59 (MQ=255)
                          aTCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATcaa                                                  <  1:2573065/65‑3 (MQ=255)
                          aTCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGc                                            >  1:934271/1‑71 (MQ=255)
                          aTCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAg                                                              >  1:84688/1‑53 (MQ=255)
                          aTCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGa                                                  >  1:3472631/1‑65 (MQ=255)
                          aTCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGa                                                  >  1:622401/1‑65 (MQ=255)
                            cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                                                        >  1:430826/1‑57 (MQ=255)
                            cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                                                        >  1:75521/1‑57 (MQ=255)
                            cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                          >  1:3202735/1‑71 (MQ=255)
                            cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                          >  1:683712/1‑71 (MQ=255)
                            cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                          >  1:3124160/1‑71 (MQ=255)
                            cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                          >  1:1200076/1‑71 (MQ=255)
                              ggCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGa                                                  <  1:3232434/61‑1 (MQ=255)
                               gCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGa                                                  <  1:2694772/60‑1 (MQ=255)
                                 aGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa                                                          <  1:1952353/50‑1 (MQ=255)
                                  aaCCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGtta                                    >  1:3413941/2‑70 (MQ=255)
                                   aCCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                                 <  1:540759/57‑1 (MQ=255)
                                   aCCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgt                                   <  1:794064/71‑1 (MQ=255)
                                   aCCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgt                                   >  1:1962760/1‑71 (MQ=255)
                                    ccGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAAc                                                       >  1:3332771/1‑50 (MQ=255)
                                    ccGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgtg                                  <  1:1080095/71‑1 (MQ=255)
                                     cGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGa                                 <  1:472274/71‑1 (MQ=255)
                                      ggTCAGAACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCACTGGTtatggg                                >  1:2884392/1‑69 (MQ=255)
                                       gTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGa                                 >  1:2894459/1‑69 (MQ=255)
                                        tCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGGCGAACACGCCCTGGTTGTGAGtt                              <  1:1754112/71‑1 (MQ=255)
                                         cAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGttt                             >  1:10591/1‑71 (MQ=255)
                                          aGAACTCCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGtttt                            <  1:2585645/71‑1 (MQ=255)
                                          aGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTTTAAACGTCGAACACGCCCTGGGTGTGAGtttt                            <  1:3339867/71‑1 (MQ=255)
                                          aGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCg                                                   <  1:3287562/48‑1 (MQ=255)
                                          aGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGttt                             <  1:858471/70‑1 (MQ=255)
                                          aGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGttt                             <  1:186373/70‑1 (MQ=255)
                                          aGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGttt                             <  1:2490538/70‑1 (MQ=255)
                                            aaCACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTa                           <  1:816445/70‑1 (MQ=255)
                                             acgccAGATTTACGGCATCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                         >  1:1815777/4‑71 (MQ=255)
                                              cacCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACgg                        >  1:1129451/1‑71 (MQ=255)
                                              cacCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACgg                        >  1:437384/1‑71 (MQ=255)
                                               gccAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGttt                             <  1:1631026/64‑1 (MQ=255)
                                               gccAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                         >  1:2061190/2‑69 (MQ=255)
                                                 cAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTAt                     >  1:1531336/1‑71 (MQ=255)
                                                  aGATTTACGGCAGCGCAGCATATCCGGTGAGTAAAc                                                       >  1:3579797/1‑36 (MQ=12)
                                                  aGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGtt                                     >  1:1452631/1‑54 (MQ=255)
                                                   gATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                         >  1:1326924/1‑65 (MQ=255)
                                                   gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACg                                             >  1:1028313/1‑45 (MQ=255)
                                                   gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGt                       >  1:1099139/1‑67 (MQ=255)
                                                   gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTAtt                    >  1:593872/1‑70 (MQ=255)
                                                   gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTc                   <  1:1668607/71‑1 (MQ=255)
                                                      ttACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTg                                        >  1:2056657/1‑47 (MQ=255)
                                                      ttACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGt                >  1:3031529/1‑71 (MQ=255)
                                                      ttACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGt                >  1:1708068/1‑71 (MQ=255)
                                                       tACGGCAACGCAGGATGTCCGGAGTGTAAACGTCg                                                   <  1:627472/35‑1 (MQ=255)
                                                       tACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGt                <  1:2080340/70‑1 (MQ=255)
                                                        aCGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTg               <  1:2314118/70‑1 (MQ=255)
                                                        aCGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTg               <  1:3194530/70‑1 (MQ=255)
                                                         cGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGa              <  1:1607697/70‑1 (MQ=255)
                                                             aaCGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtt         >  1:1050678/1‑71 (MQ=255)
                                                              aCGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                          >  1:1502410/1‑37 (MQ=255)
                                                              aCGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                          >  1:649819/1‑37 (MQ=255)
                                                              aCGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGa           >  1:3117557/1‑68 (MQ=255)
                                                                gCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAg                                <  1:82200/45‑1 (MQ=255)
                                                                gCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAttttt      <  1:1930890/71‑1 (MQ=255)
                                                                 cAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtttttt     >  1:1831887/1‑71 (MQ=255)
                                                                  aGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAttttttt    >  1:1797344/1‑71 (MQ=255)
                                                                   ggATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTg   <  1:1372368/71‑1 (MQ=255)
                                                                    gATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGa  <  1:2937464/71‑1 (MQ=255)
                                                                    gATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGa  <  1:2074140/71‑1 (MQ=255)
                                                                    |                                                                      
CCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGA  >  minE/609715‑609853

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