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
A1 F3 I96 R1 172 34.7 2005794 96.2% 1929573 85.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,306,692 A→T I115N (ATC→AAC)  phnM ← alpha‑D‑ribose 1‑methylphosphonate 5‑triphosphate diphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,306,6920AT100.0% 97.1 / NA 27I115N (ATC→AAC) phnMalpha‑D‑ribose 1‑methylphosphonate 5‑triphosphate diphosphatase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (13/14);  total (13/14)

GTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGCT  >  NZ_CP009273/4306605‑4306777
                                                                                       |                                                                                     
gtggtgTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGttc                                                                                     >  2:850380/1‑90 (MQ=255)
   gtgtGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGttcatc                                                                                  <  1:567969/90‑1 (MQ=255)
        aTGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTtctc                                                                             >  2:671553/1‑90 (MQ=255)
          gCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCAtct                                                                                 >  2:927810/1‑84 (MQ=255)
          gCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCAtct                                                                                 <  1:927810/84‑1 (MQ=255)
           cGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAg                                                                          <  1:850380/90‑1 (MQ=255)
                      agcgcagATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGcc                                                               <  2:93442/90‑1 (MQ=255)
                                 aGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGccgccg                                                    >  2:691518/1‑90 (MQ=255)
                                 aGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGccgccg                                                    >  2:729854/1‑90 (MQ=255)
                                     ggTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTcgc                                                <  2:341964/90‑1 (MQ=255)
                                      gTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTcgcg                                               >  2:480798/1‑90 (MQ=255)
                                      gTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTcgcg                                               >  1:444210/1‑90 (MQ=255)
                                      gTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTcgcg                                               >  1:697375/1‑90 (MQ=255)
                                      gTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTcgcg                                               >  2:872183/1‑90 (MQ=255)
                                      gTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTcgcg                                               >  1:862692/1‑90 (MQ=255)
                                        gCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCa                                             <  1:800006/90‑1 (MQ=255)
                                         cTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCAc                                            <  1:325015/90‑1 (MQ=255)
                                           cGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGt                                          >  1:198005/1‑90 (MQ=255)
                                                cgGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCa                                     <  1:405935/90‑1 (MQ=255)
                                                cgGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCa                                     <  2:290494/90‑1 (MQ=255)
                                                cgGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCa                                     >  2:982394/1‑90 (MQ=255)
                                                  gTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAAt                                   >  1:297751/1‑90 (MQ=255)
                                                        gccgcGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCAc                             <  1:858371/90‑1 (MQ=255)
                                                          cgcgTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACgg                           <  1:691518/90‑1 (MQ=255)
                                                            cgTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCa                         <  1:575792/90‑1 (MQ=255)
                                                               ttCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATcc                      <  2:239013/90‑1 (MQ=255)
                                                                     gTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACg                <  1:729854/90‑1 (MQ=255)
                                                                             gATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATcc                      >  2:321156/1‑76 (MQ=255)
                                                                             gATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATcc                      <  1:321156/76‑1 (MQ=255)
                                                                                   gttgttCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGCt  >  1:894489/1‑90 (MQ=255)
                                                                                       |                                                                                     
GTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGCT  >  NZ_CP009273/4306605‑4306777

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGCTCGCCACCATC  >  NZ_CP009273/4306595‑4306787
                                                                                                 |                                                                                               
caccacgctgccgctgtttgaaaaactggtGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTC                                                                                               <  SRR3722236.303248/70‑1 (MQ=60)
     GCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTT                                                                                          <  SRR3722236.120746/100‑1 (MQ=60)
      tcgtcggcagcgtcagatgtgtataagagaCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTC                                                                                         <  SRR3722236.139174/70‑1 (MQ=60)
             GTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGAT                                                                                  <  SRR3722236.576711/100‑1 (MQ=60)
              TGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATT                                                                                 <  SRR3722236.944193/100‑1 (MQ=60)
                     CGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGC                                                                          <  SRR3722236.865257/100‑1 (MQ=60)
                                      GATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCG                                                         >  SRR3722236.877822/1‑100 (MQ=60)
                                      GATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCG                                                         >  SRR3722236.450578/1‑100 (MQ=60)
                                      GATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCG                                                         >  SRR3722236.709132/1‑100 (MQ=60)
                                           AGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGT                                                    >  SRR3722236.200953/1‑100 (MQ=60)
                                                  GCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAAT                                             >  SRR3722236.302059/1‑100 (MQ=60)
                                                  GCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAAT                                             <  SRR3722236.813894/100‑1 (MQ=60)
                                                   CTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATT                                            <  SRR3722236.329658/100‑1 (MQ=60)
                                                          CGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGG                                     <  SRR3722236.411615/100‑1 (MQ=60)
                                                                  GCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGT                             <  SRR3722236.873415/100‑1 (MQ=60)
                                                                    CGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTAC                           <  SRR3722236.703158/100‑1 (MQ=60)
                                                                      CGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGG                         <  SRR3722236.584665/100‑1 (MQ=60)
                                                                         acagGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGG                      <  SRR3722236.325750/96‑1 (MQ=60)
                                                                               GTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGC                <  SRR3722236.742258/100‑1 (MQ=60)
                                                                                   CTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGCT            >  SRR3722236.910240/1‑100 (MQ=60)
                                                                                           GCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACctgtctcttatacacatctgacg    >  SRR3722236.817640/1‑77 (MQ=60)
                                                                                             GTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGCTCGCCACCATC  >  SRR3722236.769216/1‑100 (MQ=60)
                                                                                                 |                                                                                               
CAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGCTCGCCACCATC  >  NZ_CP009273/4306595‑4306787

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: