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
A2 F5 I181 R1 159 32.7 1279082 96.9% 1239430 121.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 923,620 A→G V344V (GTT→GTC cydC ← cysteine/glutathione ABC transporter ATP‑binding protein/permease CydC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273923,6200AG100.0% 54.2 / NA 16V344V (GTT→GTCcydCcysteine/glutathione ABC transporter ATP‑binding protein/permease CydC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (11/5);  total (11/5)

ATGCGCGGGTCAGCTGTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGC  >  NZ_CP009273/923484‑923754
                                                                                                                                        |                                                                                                                                      
aTGCGCGGGTCAGCTGTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGtgaactggac                                                                                                                                      <  2:639416/139‑1 (MQ=255)
 tGCGCGGGTCAGCTGTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGGACa                                                                                                                                     >  2:432761/1‑139 (MQ=255)
        gTCAGCTGTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTACGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGGACATCCCGt                                                                                                                               >  2:320698/1‑138 (MQ=255)
                ttgtAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGGACATCCCGTAACGTCAGc                                                                                                                      >  2:59367/1‑139 (MQ=255)
                          gtTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCTGGATAAGTGAACTGGACATCCCGTAACGTCAGCGAAACGCGAt                                                                                                            >  2:150433/1‑139 (MQ=255)
                          gtTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGGACATCCCGTAACGTCAGCGAAACGCGAt                                                                                                            >  1:388401/1‑139 (MQ=255)
                                   ccGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGGACATCCCGTAACGTCAGCGAAACGCGATCGGCAacac                                                                                                   >  2:584659/1‑139 (MQ=255)
                                          cGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGGACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTg                                                                                            <  1:186417/139‑1 (MQ=255)
                                                 cccGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGGACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTAt                                                                                       >  1:309301/1‑137 (MQ=255)
                                                 cccGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGGACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTAt                                                                                       >  1:22267/1‑137 (MQ=255)
                                                                         ggCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGGACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGa                                                             <  1:584659/139‑1 (MQ=255)
                                                                                     aaGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCTGGATAAGTGAACTGGACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGt                                                   <  1:150433/137‑1 (MQ=255)
                                                                                     aaGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGGACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGt                                                   >  2:210303/1‑137 (MQ=255)
                                                                                       gagaAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGGACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCa                                                 >  1:421655/1‑137 (MQ=255)
                                                                                            aTCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGGACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGAtac                                          >  2:632976/1‑139 (MQ=255)
                                                                                                          cTGTTGAGATTGCTCCGGAGAAGTGAACTGGACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCa                            <  1:210303/139‑1 (MQ=255)
                                                                                                                         cGGATAAGTGAACTGGACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGa             <  2:421655/139‑1 (MQ=255)
                                                                                                                                    actggacATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGc  <  1:432761/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                      
ATGCGCGGGTCAGCTGTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGC  >  NZ_CP009273/923484‑923754

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTG  >  NZ_CP009273/923499‑923765
                                                                                                                         |                                                                                                                                                 
GTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGGACATCCCGTAACGTCAGCGAAACGCGATC                                                                                                                      >  SRR3722135.393826/1‑151 (MQ=60)
                       CATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGGACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATC                                                                                                 >  SRR3722135.22559/1‑149 (MQ=60)
                       CATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGGACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATC                                                                                                 >  SRR3722135.313729/1‑149 (MQ=60)
                          CCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGGACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAA                                                                                            <  SRR3722135.189037/151‑1 (MQ=60)
                                                         CGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGGACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTC                                                             <  SRR3722135.592335/151‑1 (MQ=60)
                                                             GTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGGACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAG                                                           >  SRR3722135.427438/1‑149 (MQ=60)
                                                                     GAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCTGGATAAGTGAACTGGACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGT                                                   <  SRR3722135.152520/149‑1 (MQ=60)
                                                                                          CCTGTTGAGATTGCTCCGGAGAAGTGAACTGGACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCC                            <  SRR3722135.213283/151‑1 (MQ=60)
                                                                                                                    AACTGGACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTG  <  SRR3722135.438666/151‑1 (MQ=60)
                                                                                                                         |                                                                                                                                                 
GTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTG  >  NZ_CP009273/923499‑923765

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