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 I195 R1 222 27.1 1551108 96.0% 1489063 86.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,903,061 A→T intergenic (‑38/‑121) mgrB ← / → yobH PhoP/PhoQ regulator MgrB/YobH family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,903,0610AT100.0% 107.8 / NA 28intergenic (‑38/‑121)mgrB/yobHPhoP/PhoQ regulator MgrB/YobH family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (15/13);  total (15/13)

GCAAGCCAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTAGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGC  >  NZ_CP009273/1902976‑1903142
                                                                                     |                                                                                 
gcaagcCAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCAt                                                                               <  1:135567/90‑1 (MQ=255)
   agcCAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCg                                                                            <  1:43156/90‑1 (MQ=255)
   agcCAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCg                                                                            <  1:437573/90‑1 (MQ=255)
     ccAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCgtt                                                                          <  2:520302/90‑1 (MQ=255)
     ccAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCgtt                                                                          <  1:274522/90‑1 (MQ=255)
             acGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAAc                                                                  <  2:126717/90‑1 (MQ=255)
             acGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAAc                                                                  <  1:585781/90‑1 (MQ=255)
              cGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACa                                                                 <  1:241130/90‑1 (MQ=255)
               gACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAg                                                                <  2:319052/90‑1 (MQ=255)
                  aaCCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAgcgc                                                             <  2:731785/90‑1 (MQ=255)
                   aCCAGAACGGCCCCTCGAAACTTTTTCCCTCCAATCTCCGTTTCACCTCCCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGcc                                                            <  1:15683/90‑1 (MQ=255)
                    ccAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCt                                                           <  2:500544/90‑1 (MQ=255)
                            aCCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGt                                                   >  1:425386/1‑90 (MQ=255)
                            aCCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGt                                                   >  1:183279/1‑90 (MQ=255)
                                                         cGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCg                      >  1:493902/1‑90 (MQ=255)
                                                         cGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCg                      >  1:501655/1‑90 (MQ=255)
                                                         cGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCg                      >  2:49491/1‑90 (MQ=255)
                                                             tCACCTACCGTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCgcttg                  >  1:137676/1‑90 (MQ=255)
                                                                  tacctTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCa             >  2:439716/1‑90 (MQ=255)
                                                                  tacctTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCa             >  2:99274/1‑90 (MQ=255)
                                                                       tATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCa        >  2:228112/1‑90 (MQ=255)
                                                                       tATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCa        >  2:631927/1‑90 (MQ=255)
                                                                        aTGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCat       >  2:721392/1‑90 (MQ=255)
                                                                             aTTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCat       >  1:63167/1‑85 (MQ=255)
                                                                             aTTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCat       <  2:63167/85‑1 (MQ=255)
                                                                             aTTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGc  >  2:452179/1‑90 (MQ=255)
                                                                             aTTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGc  >  2:306246/1‑90 (MQ=255)
                                                                             aTTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGc  >  1:518288/1‑90 (MQ=255)
                                                                                     |                                                                                 
GCAAGCCAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTAGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGC  >  NZ_CP009273/1902976‑1903142

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ 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
CAGCAAGCAAGCCAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTAGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGC  >  NZ_CP009273/1902970‑1903142
                                                                                           |                                                                                 
CAGCAAGCAAGCCAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGT                                                                           <  SRR3722073.281742/100‑1 (MQ=60)
      GCAAGCCAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTA                                                                     <  SRR3722073.138389/100‑1 (MQ=60)
         AGCCAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAAC                                                                  <  SRR3722073.44190/100‑1 (MQ=60)
         AGCCAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAAC                                                                  <  SRR3722073.446933/100‑1 (MQ=60)
           CCAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAG                                                                <  SRR3722073.280326/100‑1 (MQ=60)
               cgagatctacacgaccgccatcgtcggcagcgtcagatgtgtataagagacagCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCC                                                            <  SRR3722073.782890/47‑1 (MQ=60)
                   ACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAAC                                                        <  SRR3722073.598018/100‑1 (MQ=60)
                    CGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACT                                                       <  SRR3722073.246189/100‑1 (MQ=60)
                        AACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGT                                                   >  SRR3722073.186845/1‑100 (MQ=60)
                        AACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGT                                                   >  SRR3722073.434486/1‑100 (MQ=60)
                         ACCAGAACGGCCCCTCGAAACTTTTTCCCTCCAATCTCCGTTTCACCTCCCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTC                                                  <  SRR3722073.16082/100‑1 (MQ=60)
                                                     CTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCG                      >  SRR3722073.504361/1‑100 (MQ=60)
                                                     CTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCG                      >  SRR3722073.512262/1‑100 (MQ=60)
                                                         AATCTCCGTTTCACCTACCGTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTG                  >  SRR3722073.140539/1‑100 (MQ=60)
                                                                         ACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGC  >  SRR3722073.529192/1‑100 (MQ=60)
                                                                         ACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGC  >  SRR3722073.64649/1‑100 (MQ=60)
                                                                                           |                                                                                 
CAGCAAGCAAGCCAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTAGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGC  >  NZ_CP009273/1902970‑1903142

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