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 F2 I175 R1 221 18.5 1047116 95.1% 995807 84.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,616,816 T→C V363A (GTC→GCC)  radA → DNA repair protein RadA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,616,8160TC100.0% 105.5 / NA 30V363A (GTC→GCC) radADNA repair protein RadA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/15);  total (15/15)

AAACCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGC  >  NZ_CP009273/4616733‑4616902
                                                                                   |                                                                                      
aaaCCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGCcggcgg                                                                                  >  2:344459/1‑90 (MQ=255)
aaaCCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGCcggcgg                                                                                  >  1:462507/1‑90 (MQ=255)
    cGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGCCGGCGGCGTg                                                                              <  2:400314/90‑1 (MQ=255)
         ggCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGCCGGCGGCGTGAAGGt                                                                         <  1:240839/90‑1 (MQ=255)
         ggCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGCCGGCGGCGTGAAGGt                                                                         <  2:275032/90‑1 (MQ=255)
                   ctggctgTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGCCGGCGGCGTGAAGGTAACCGAAAcc                                                               >  2:70417/1‑90 (MQ=255)
                              gCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGCCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACtt                                                    <  1:344459/90‑1 (MQ=255)
                               cACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGCCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTa                                                   <  2:9541/90‑1 (MQ=255)
                               cACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGCCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTa                                                   <  1:363176/90‑1 (MQ=255)
                                   gTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGCcggc                                                                                    <  2:468857/53‑1 (MQ=255)
                                   gTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGCcggc                                                                                    >  1:468857/1‑53 (MQ=255)
                                               tGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGCCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGa                                   <  1:99111/90‑1 (MQ=255)
                                                 cAAATGGCCGATCAGGATGTGTTTGTGAACGTGGCCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATg                                 <  2:237221/90‑1 (MQ=255)
                                                     tGGCCGATCAGGATGTGTTTGTGAACGTGGCCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGttt                             >  1:200301/1‑90 (MQ=255)
                                                     tGGCCGATCAGGATGTGTTTGTGAACGTGGCCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGttt                             >  1:424410/1‑90 (MQ=255)
                                                      ggCCGATCAGGATGTGTTTGTGAACGTGGCCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTc                            <  2:424410/90‑1 (MQ=255)
                                                      ggCCGATCAGGATGTGTTTGTGAACGTGGCCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTc                            >  1:114290/1‑90 (MQ=255)
                                                        ccGATCAGGATGTGTTTGTGAACGTGGCCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCa                          >  2:453723/1‑90 (MQ=255)
                                                          gATCAGGATGTGTTTGTGAACGTGGCCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGc                        >  1:52573/1‑90 (MQ=255)
                                                                gATGTGTTTGTGAACGTGGCCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGtt                                              >  1:369949/1‑62 (MQ=255)
                                                                gATGTGTTTGTGAACGTGGCCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGtt                                              <  2:369949/62‑1 (MQ=255)
                                                                gATGTGTTTGTGAACGTGGCCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTgcgc                  <  2:200301/90‑1 (MQ=255)
                                                                  tgtgTTTGTGAACGTGGCCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGa                >  2:433204/1‑90 (MQ=255)
                                                                    tgtTTGTGAACGTGGCCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACa              >  1:318062/1‑90 (MQ=255)
                                                                     gtTTGTGAACGTGGCCGGCGGCGTGAAGGTAACCg                                                                    <  2:274172/35‑1 (MQ=38)
                                                                     gtTTGTGAACGTGGCCGGCGGCGTGAAGGTAACCg                                                                    >  1:274172/1‑35 (MQ=38)
                                                                     gtTTGTGAACGTGGCCGGCGGCGTGAAGGTAACCGAAAcc                                                               <  2:319431/40‑1 (MQ=39)
                                                                     gtTTGTGAACGTGGCCGGCGGCGTGAAGGTAACCGAAAcc                                                               >  1:319431/1‑40 (MQ=39)
                                                                             aCGTGGCCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGc     >  2:183917/1‑90 (MQ=255)
                                                                                tGGCCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGc  <  2:446521/90‑1 (MQ=255)
                                                                                   |                                                                                      
AAACCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGC  >  NZ_CP009273/4616733‑4616902

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ 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
GGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAG  >  NZ_CP009273/4616723‑4616891
                                                                                             |                                                                           
GGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGCCGGCGG                                                                       <  SRR3722048.119571/100‑1 (MQ=60)
GGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGCCGGCGG                                                                       >  SRR3722048.470130/1‑100 (MQ=60)
                   GGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGCCGGCGGCGTGAAGGTAACCGAAACC                                                    <  SRR3722048.244246/100‑1 (MQ=60)
                                   GTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGCCGGCGGCGTGAAGGctgtctcttatacacatctccgagccc                                    >  SRR3722048.476618/1‑73 (MQ=60)
                                        GCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGCCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTG                               <  SRR3722048.349844/100‑1 (MQ=60)
                                         CACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGCCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGC                              <  SRR3722048.368987/100‑1 (MQ=60)
                                                     GGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGCCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTT                  >  SRR3722048.203105/1‑100 (MQ=60)
                                                     GGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGCCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTT                  >  SRR3722048.431333/1‑100 (MQ=60)
                                                      GTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGCCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTC                 >  SRR3722048.115958/1‑100 (MQ=60)
                                                         TGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGCCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAG              <  SRR3722048.100555/100‑1 (MQ=60)
                                                          GCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGCCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGC             >  SRR3722048.53322/1‑100 (MQ=60)
                                                                GGCCGATCAGGATGTGTTTGTGAACGTGGCCGGCGCCGTGAAGGTAACCGAAACCAGTGCCGACTTAGGGTTACTGCTGGCGctgtctcttatacaaatc       >  SRR3722048.375907/1‑82 (MQ=60)
                                                                    GATCAGGATGTGTTTGTGAACGTGGCCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACA   >  SRR3722048.322897/1‑100 (MQ=60)
                                                                     ATCAGGATGTGTTTGTGAACGTGGCCGGCGGCGTGAAGGTAACCGAAACCAGTGCCtgtctcttatacacatctccgagcccacgagacctctctacatc  >  SRR3722048.278218/1‑56 (MQ=60)
                                                                     ATCAGGATGTGTTTGTGAACGTGGCCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTctgtctcttatacacatctccgagcccacgagacctctct  >  SRR3722048.324295/1‑60 (MQ=60)
                                                                                             |                                                                           
GGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAG  >  NZ_CP009273/4616723‑4616891

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