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 F1 I0 R1 2892 75.3 5698161 66.1% 3766484 59.5

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,482 T→G 23.7% G372G (GGA→GGC)  arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,4820TG23.7% 189.4 / 91.0 135G372G (GGA→GGC) arsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base T (49/54);  new base G (21/11);  total (70/65)
Fisher's exact test for biased strand distribution p-value = 1.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.40e-01

Only 100 of 136 total aligned reads displayed.

ATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCAT  >  minE/2595414‑2595552
                                                                    |                                                                      
atatTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTcc                                                                      >  1:3141656/1‑71 (MQ=255)
 tatTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTcc                                                                      >  1:3336725/1‑70 (MQ=255)
 tatTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTcc                                                                      >  1:4974575/1‑70 (MQ=255)
   ttcttcTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCaa                                                                   <  1:5141489/71‑1 (MQ=255)
    tcttctGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTccc                                                                     >  1:1699749/1‑68 (MQ=255)
     cttctGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCaaa                                                                  <  1:4285372/70‑1 (MQ=255)
       tctGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCg                                                               >  1:855357/1‑71 (MQ=255)
          gCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAg                                                            >  1:4940678/1‑71 (MQ=255)
          gCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAg                                                            >  1:2687148/1‑71 (MQ=255)
           cGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGc                                                              <  1:2964023/68‑1 (MQ=255)
           cGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGc                                                              <  1:2793191/68‑1 (MQ=255)
            gAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATc                                                                <  1:3178037/65‑1 (MQ=255)
            gAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGc                                                           >  1:3121949/1‑70 (MQ=255)
             aaaGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCa                                                         <  1:1807966/71‑1 (MQ=255)
             aaaGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCa                                                         <  1:1218542/71‑1 (MQ=255)
              aaGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCaa                                                        <  1:4062936/71‑1 (MQ=255)
              aaGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCaa                                                        <  1:1582911/71‑1 (MQ=255)
              aaGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCaa                                                        <  1:3251813/71‑1 (MQ=255)
               aGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAAt                                                       >  1:3615097/1‑71 (MQ=255)
                 tACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATca                                                     <  1:2751941/71‑1 (MQ=255)
                 tACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATca                                                     <  1:2745415/71‑1 (MQ=255)
                 tACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATca                                                     >  1:3863892/1‑71 (MQ=255)
                  acatGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATca                                                     >  1:1514912/4‑70 (MQ=255)
                    gTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACAt                                                  <  1:4778729/71‑1 (MQ=255)
                    gTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACAt                                                  <  1:3790476/71‑1 (MQ=255)
                    atGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACAt                                                  <  1:3814260/70‑1 (MQ=255)
                     tGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACAt                                                  >  1:665170/1‑70 (MQ=255)
                        aGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTgg                                               >  1:2543930/1‑70 (MQ=255)
                           cAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATa                                           <  1:1972875/71‑1 (MQ=255)
                              agTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAAc                                        <  1:3884190/71‑1 (MQ=255)
                                tagcgtagcCAGGCTACCAATTGGGGTAATTTTCGGTCCCa                                                                    <  1:698790/41‑1 (MQ=255)
                                tagcgtagcCAGGCTACCAATTGGGGTAATTTTCGGTCCCa                                                                    <  1:892936/41‑1 (MQ=255)
                                     ttgcCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATcac                                                    <  1:5693837/50‑1 (MQ=255)
                                     ttgcCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCAt                                     >  1:308917/3‑67 (MQ=255)
                                     ttgcCAGACTGCCAACCGGGGTGATTTTCGGGCCTAAATCACAGCCAATcac                                                    <  1:5169098/50‑1 (MQ=255)
                                     tagcCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATcac                                                    <  1:3109399/52‑1 (MQ=255)
                                     tagcCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCAt                                     <  1:5261264/67‑1 (MQ=255)
                                      agcCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAACCGCAGCCAATCACATTGGCATAAACCATCGctt                                >  1:1749143/1‑71 (MQ=255)
                                       gcCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATa                                           <  1:4551350/59‑1 (MQ=255)
                                       gcCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTg                                                >  1:830136/1‑54 (MQ=255)
                                       gcCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGcctc                               >  1:4111743/1‑71 (MQ=255)
                                       gcCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGcctc                               >  1:138197/1‑71 (MQ=255)
                                        ccAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAAcc                                       <  1:1278617/62‑1 (MQ=255)
                                         cAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATaaa                                         >  1:725078/1‑59 (MQ=255)
                                          aGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCttt                            <  1:1136169/71‑1 (MQ=255)
                                            gCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATa                                           <  1:5274453/54‑1 (MQ=255)
                                            gCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATa                                           <  1:5621723/54‑1 (MQ=255)
                                            gCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGa                          >  1:3076353/1‑71 (MQ=255)
                                             ctgccAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGcctct                              >  1:445603/4‑66 (MQ=255)
                                             ctgccAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGcctct                              >  1:3569147/4‑66 (MQ=255)
                                             ctgccAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGcctct                              >  1:2284202/4‑66 (MQ=255)
                                             ctgccAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGcctct                              >  1:5456389/4‑66 (MQ=255)
                                              tgccAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATca                                                     <  1:5408309/40‑1 (MQ=16)
                                              tgccAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATca                                                     <  1:4793352/40‑1 (MQ=16)
                                              tgccAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATca                                                     <  1:5516054/40‑1 (MQ=16)
                                                ccAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCATTCACATTGGCATAAACCATCGCTTCTTTGa                          >  1:3344975/1‑67 (MQ=255)
                                                ccAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATc                                                      >  1:1156533/1‑39 (MQ=255)
                                                ccAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATc                                    <  1:4385994/57‑1 (MQ=255)
                                                ccAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAAc                      >  1:3451052/1‑71 (MQ=255)
                                                ccAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGacg                        <  1:530882/69‑3 (MQ=255)
                                                 cAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATcaca                                                   <  1:2562623/41‑1 (MQ=255)
                                                  aaTCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCa                                      >  1:1377362/1‑53 (MQ=255)
                                                  aaTCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGc                    <  1:2905130/71‑1 (MQ=255)
                                                    ttGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAAcc                                       <  1:5342575/50‑1 (MQ=255)
                                                    ttGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGc                    >  1:3748577/1‑69 (MQ=255)
                                                    ttGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGc                    >  1:3177689/1‑69 (MQ=255)
                                                    tcggggTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCa                                      >  1:2381800/3‑51 (MQ=255)
                                                     tGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGc                    >  1:4094631/1‑68 (MQ=255)
                                                      ggggTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAgt                >  1:4442622/1‑70 (MQ=255)
                                                      ggggTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAg                 <  1:3745095/70‑1 (MQ=255)
                                                      ggggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATaaa                                         >  1:1519117/1‑46 (MQ=255)
                                                      ggggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAAc                                        <  1:2011736/47‑1 (MQ=255)
                                                      ggggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAAc                                        <  1:4203771/47‑1 (MQ=255)
                                                      ggggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGcttctt                             >  1:2215367/1‑58 (MQ=255)
                                                      ggggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATaa                       >  1:623043/1‑64 (MQ=255)
                                                      ggggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATaa                       >  1:1869377/1‑64 (MQ=255)
                                                      ggggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATaa                       >  1:2880212/1‑64 (MQ=255)
                                                      ggggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGc                    >  1:4738066/1‑67 (MQ=255)
                                                      ggggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGc                    >  1:4763680/1‑67 (MQ=255)
                                                       gggTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAgtc               >  1:3301054/1‑69 (MQ=255)
                                                       gggTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAgtc               >  1:5121328/1‑69 (MQ=255)
                                                       gggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGa                <  1:720008/70‑1 (MQ=255)
                                                       gggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGa                <  1:2410889/70‑1 (MQ=255)
                                                        ggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGAt                         <  1:1616517/60‑1 (MQ=255)
                                                         gTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCa                                      <  1:5212239/46‑1 (MQ=255)
                                                         gTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTg                           <  1:4444327/57‑1 (MQ=255)
                                                          tAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGcc            <  1:848510/71‑1 (MQ=255)
                                                          tAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGcc            >  1:4639318/1‑71 (MQ=255)
                                                           aaTTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATa                                           <  1:5461221/39‑1 (MQ=255)
                                                           aaTTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCg           >  1:1316500/1‑71 (MQ=255)
                                                            aTTTTCGGTCCCAAATCTCAGCCAATCACATTGGCATAAACAATCTATTCTTTGATAACGCCAGATGcc            >  1:5660729/1‑69 (MQ=255)
                                                            aTTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGc                    <  1:3192559/61‑1 (MQ=255)
                                                             ttttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCttt                            <  1:4251608/52‑1 (MQ=255)
                                                             ttttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGt          <  1:2214288/70‑1 (MQ=255)
                                                             ttttCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGta         >  1:3302211/1‑70 (MQ=255)
                                                              tttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATa                        <  1:1283530/55‑1 (MQ=255)
                                                              tttCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGtac        >  1:4252284/1‑69 (MQ=255)
                                                              tttCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGtac        >  1:947643/1‑69 (MQ=255)
                                                               ttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgct       >  1:1945590/1‑71 (MQ=255)
                                                                   gTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGcc                   >  1:1162773/1‑55 (MQ=255)
                                                                    |                                                                      
ATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCAT  >  minE/2595414‑2595552

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