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 I194 R1 197 32.2 1845014 95.2% 1756453 85.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,259,173 A→G T2A (ACC→GCC)  hemA → glutamyl‑tRNA reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,259,1730AG100.0% 78.7 / NA 23T2A (ACC→GCC) hemAglutamyl‑tRNA reductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/8);  total (15/8)

GGGCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCA  >  NZ_CP009273/1259086‑1259258
                                                                                       |                                                                                     
gggCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGcc                                                                                     >  2:740981/1‑90 (MQ=255)
gggCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGcc                                                                                     >  1:878030/1‑90 (MQ=255)
           tttATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACt                                                                          <  1:643170/90‑1 (MQ=255)
                 gATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTAt                                                                    >  2:3038/1‑90 (MQ=255)
                                  atgatgCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCAc                                                   >  2:579705/1‑90 (MQ=255)
                                     atgCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTg                                                >  2:249043/1‑90 (MQ=255)
                                        cAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTAt                                             >  1:869306/1‑90 (MQ=255)
                                           gCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGc                                          <  1:813536/90‑1 (MQ=255)
                                              gACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTc                                                                         <  1:665837/56‑1 (MQ=255)
                                              gACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTc                                                                         >  2:665837/1‑56 (MQ=255)
                                              gACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGc                                       >  2:857913/1‑90 (MQ=255)
                                                              gTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGtt                       >  1:711107/1‑90 (MQ=255)
                                                                   tattatTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGc                  <  1:348200/90‑1 (MQ=255)
                                                                      tattTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACgg                                                      >  1:304478/1‑51 (MQ=255)
                                                                      tattTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACgg                                                      <  2:304478/51‑1 (MQ=255)
                                                                      tattTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCgaga                                   >  2:908256/1‑70 (MQ=255)
                                                                      tattTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCgaga                                   <  1:908256/70‑1 (MQ=255)
                                                                       attTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGa              >  2:280829/1‑90 (MQ=255)
                                                                        tttCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGAt             <  1:74484/90‑1 (MQ=255)
                                                                           cccGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAg          >  1:691812/1‑90 (MQ=255)
                                                                            ccGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGc         >  2:471907/1‑90 (MQ=255)
                                                                            ccGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGc         >  1:52541/1‑90 (MQ=255)
                                                                                   cATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCa  <  1:92728/90‑1 (MQ=255)
                                                                                       |                                                                                     
GGGCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCA  >  NZ_CP009273/1259086‑1259258

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
CTCCAAAAGGGGGCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCAGGCGCTTGAC  >  NZ_CP009273/1259076‑1259268
                                                                                                 |                                                                                               
CTCCAAAAGGGGGCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCC                                                                                               >  SRR3722072.893152/1‑100 (MQ=60)
                     TTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAAC                                                                          <  SRR3722072.654635/100‑1 (MQ=60)
                      ctcggagatgtgtataagagacagGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACC                                                                         <  SRR3722072.677743/76‑1 (MQ=60)
                                        CTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTAT                                                       >  SRR3722072.884280/1‑100 (MQ=60)
                                                     GCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACG                                          <  SRR3722072.827707/100‑1 (MQ=60)
                                                            taagagacagACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCG                                   <  SRR3722072.923885/90‑1 (MQ=60)
                                                              CCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTT                                 >  SRR3722072.723781/1‑100 (MQ=60)
                                                                      ACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCctgtctcttatacacatctccgagcccac                         >  SRR3722072.309754/1‑71 (MQ=60)
                                                                           GGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAaaacgatacacgttctcgcagcgatac                    <  SRR3722072.846267/100‑28 (MQ=60)
                                                                           GGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAG                    >  SRR3722072.704188/1‑100 (MQ=60)
                                                                            GTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGC                   >  SRR3722072.53494/1‑100 (MQ=60)
                                                                             TATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCT                  <  SRR3722072.354337/100‑1 (MQ=60)
                                                                                  TTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATC             <  SRR3722072.75874/100‑1 (MQ=60)
                                                                                             CATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCAGGCGCTTGAC  <  SRR3722072.94437/100‑1 (MQ=60)
                                                                                                 |                                                                                               
CTCCAAAAGGGGGCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCAGGCGCTTGAC  >  NZ_CP009273/1259076‑1259268

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ 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: