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 I205 R1 220 17.5 961422 97.1% 933540 86.5

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 886,189 A→G G7G (GGT→GGC grxA ← glutaredoxin 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273886,1890AG100.0% 83.7 / NA 25G7G (GGT→GGCgrxAglutaredoxin 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/7);  total (18/7)

TACTGATACTGAAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGT  >  NZ_CP009273/886106‑886274
                                                                                   |                                                                                     
tactgatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCaaaa                                                                                 >  2:451018/1‑90 (MQ=255)
tactgatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCaaaa                                                                                 >  2:191497/1‑90 (MQ=255)
tactgatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCa                                                                                    >  1:387325/1‑87 (MQ=255)
tactgatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCa                                                                                    <  2:387325/87‑1 (MQ=255)
  ctgatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCa                                                                                    <  1:217133/85‑1 (MQ=255)
  ctgatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCa                                                                                    >  2:217133/1‑85 (MQ=255)
       actgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACgg                                                                          >  1:75416/1‑90 (MQ=255)
                  tcGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCattatt                                                               >  2:361723/1‑90 (MQ=255)
                                              gATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAAt                                   >  2:74796/1‑90 (MQ=255)
                                              gATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAAt                                   >  1:34531/1‑90 (MQ=255)
                                                  ttttGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATc                               <  1:412011/90‑1 (MQ=255)
                                                               caGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGcc                                         <  1:203855/67‑1 (MQ=255)
                                                               caGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGcc                                         >  2:203855/1‑67 (MQ=255)
                                                                      ggCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTg           <  1:74796/90‑1 (MQ=255)
                                                                       gCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGa          >  1:231344/1‑90 (MQ=255)
                                                                         aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGaaa        >  2:38430/1‑90 (MQ=255)
                                                                         aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGaaa        >  2:442041/1‑90 (MQ=255)
                                                                         aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGaaa        >  2:93513/1‑90 (MQ=255)
                                                                         aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGaaa        >  2:112068/1‑90 (MQ=255)
                                                                         aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGaaa        >  1:397067/1‑90 (MQ=255)
                                                                         aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGaaa        >  1:294337/1‑90 (MQ=255)
                                                                           cccGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTAtt                <  1:84192/80‑1 (MQ=255)
                                                                           cccGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTAtt                >  2:84192/1‑80 (MQ=255)
                                                                           cccGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGt      <  1:442041/90‑1 (MQ=255)
                                                                               aaCGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGt  >  1:106040/1‑90 (MQ=255)
                                                                                   |                                                                                     
TACTGATACTGAAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGT  >  NZ_CP009273/886106‑886274

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AATATCTACATACTGATACTGAAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACC  >  NZ_CP009273/886096‑886280
                                                                                             |                                                                                           
AATATCTACATACTGATACTGAAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAA                                                                                       >  SRR3722087.391797/1‑100 (MQ=60)
       ACATACTGATACTGAAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGG                                                                                <  SRR3722087.219048/100‑1 (MQ=60)
       ACATACTGATACTGAAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGG                                                                                >  SRR3722087.76106/1‑100 (MQ=60)
                                              TTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAAT                                         >  SRR3722087.34807/1‑100 (MQ=60)
                                                  gtataagagacagTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATC                                     <  SRR3722087.205640/87‑1 (MQ=60)
                                                            TTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAAT                           <  SRR3722087.416808/100‑1 (MQ=60)
                                                                       CACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGA                >  SRR3722087.233467/1‑100 (MQ=60)
                                                                         CAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAA              >  SRR3722087.297339/1‑100 (MQ=60)
                                                                         CAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAA              >  SRR3722087.401655/1‑100 (MQ=60)
                                                                           GTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGT            <  SRR3722087.84981/100‑1 (MQ=60)
                                                                               GGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGT        >  SRR3722087.107058/1‑100 (MQ=60)
                                                                                GGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTT       <  SRR3722087.75481/100‑1 (MQ=60)
                                                                                     CCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACC  <  SRR3722087.447207/100‑1 (MQ=60)
                                                                                             |                                                                                           
AATATCTACATACTGATACTGAAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACC  >  NZ_CP009273/886096‑886280

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