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,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% 81.1 / NA 24T2A (ACC→GCC) hemAglutamyl‑tRNA reductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/6);  total (18/6)

GGGCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCA  >  NZ_CP009273/1259086‑1259258
                                                                                       |                                                                                     
gggCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGcc                                                                                     >  2:547723/1‑90 (MQ=255)
gggCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGcc                                                                                     >  2:111670/1‑90 (MQ=255)
                 gATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTAt                                                                    >  1:694640/1‑90 (MQ=255)
                 gATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTAt                                                                    >  1:110251/1‑90 (MQ=255)
                   tCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCa                                                                  <  1:18019/90‑1 (MQ=255)
                   tCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCa                                                                  >  1:324696/1‑90 (MQ=255)
                   tCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCa                                                                  >  2:445294/1‑90 (MQ=255)
                   tCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCa                                                                  >  1:752844/1‑90 (MQ=255)
                                   tgatgCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCAcc                                                  <  1:547723/90‑1 (MQ=255)
                                     atgCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTg                                                >  2:135407/1‑90 (MQ=255)
                                     atgCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTg                                                >  1:97938/1‑90 (MQ=255)
                                       gCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTa                                              >  1:190552/1‑90 (MQ=255)
                                        cAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCCCCTGTAt                                             >  1:743140/1‑90 (MQ=255)
                                         aaGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATc                                            <  2:190552/90‑1 (MQ=255)
                                          aGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCg                                           >  1:23358/1‑90 (MQ=255)
                                                        aTCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACgtgt                             >  2:13777/1‑90 (MQ=255)
                                                             cGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCAGAAAACGGCACCTGTATCGCTGCGAGAACCTGTATCGt                        >  1:666005/1‑90 (MQ=255)
                                                              gTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGtt                       >  2:523202/1‑90 (MQ=255)
                                                                   tattatTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGc                  <  1:676608/90‑1 (MQ=255)
                                                                            ccGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGc         >  1:583931/1‑90 (MQ=255)
                                                                            ccGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGc         >  1:511265/1‑90 (MQ=255)
                                                                              gCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCTCTGCGGGAACAGGTATCGTTTTCGCCGGATAAGCTc       >  1:335416/1‑90 (MQ=255)
                                                                                   cATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCa  <  2:639784/90‑1 (MQ=255)
                                                                                   cATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCa  <  2:752844/90‑1 (MQ=255)
                                                                                       |                                                                                     
GGGCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCA  >  NZ_CP009273/1259086‑1259258

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 30 ≤ 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
CTCCAAAAGGGGGCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTC  >  NZ_CP009273/1259076‑1259253
                                                                                                 |                                                                                
tctacacactgcatatcgtcggcagcgtcagatgtgtataagagacagTGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCC                                                                                <  SRR3722073.433968/52‑1 (MQ=60)
 CCCAAACGGGGGCGCTCTCTTTTATTGATATTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCC                                                                               <  SRR3722073.582900/100‑1 (MQ=60)
          cggcgaccaccgagatctacacactgcatatcgtcggcagcgtcagatgtgtataagagacaGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCAC                                                                      <  SRR3722073.347493/38‑1 (MQ=60)
                 CTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTAT                                                               >  SRR3722073.112642/1‑100 (MQ=60)
                 CTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTAT                                                               >  SRR3722073.708760/1‑100 (MQ=60)
                   CTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCA                                                             >  SRR3722073.768080/1‑100 (MQ=60)
                   CTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCA                                                             >  SRR3722073.331593/1‑100 (MQ=60)
                             TCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAAC                                                   <  SRR3722073.18474/100‑1 (MQ=60)
                                     ATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTG                                           >  SRR3722073.100085/1‑100 (MQ=60)
                                       CCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTA                                         >  SRR3722073.194328/1‑100 (MQ=60)
                                        CTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCCCCTGTAT                                        >  SRR3722073.758199/1‑100 (MQ=60)
                                          GTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCG                                      >  SRR3722073.23952/1‑100 (MQ=60)
                                             TGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTG                                   <  SRR3722073.559161/100‑1 (MQ=60)
                                                             ACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCAGAAAACGGCACCTGTATCGCTGCGAGAACCTGTATCGT                   >  SRR3722073.679647/1‑100 (MQ=60)
                                                                            GTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGC    >  SRR3722073.522040/1‑100 (MQ=60)
                                                                            GTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGC    >  SRR3722073.596119/1‑100 (MQ=60)
                                                                             TATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCT   <  SRR3722073.690443/100‑1 (MQ=60)
                                                                              ATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCTCTGCGGGAACAGGTATCGTTTTCGCCGGATAAGCTC  >  SRR3722073.342533/1‑100 (MQ=60)
                                                                                                 |                                                                                
CTCCAAAAGGGGGCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTC  >  NZ_CP009273/1259076‑1259253

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