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 465,464 T→C S380P (TCG→CCG)  mdlA → SmdA family multidrug ABC transporter permease/ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273465,4640TC100.0% 86.2 / NA 26S380P (TCG→CCG) mdlASmdA family multidrug ABC transporter permease/ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (10/16);  total (10/16)

GGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACT  >  NZ_CP009273/465382‑465540
                                                                                  |                                                                            
ggAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTCCCCTGTTGCCGCTctg                                                                       <  2:148869/90‑3 (MQ=255)
  aaaaCGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTc                                                                     <  2:516610/90‑1 (MQ=255)
                             gTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAggggg                                          >  1:751615/1‑90 (MQ=255)
                                 gATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAGGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGatat                                      <  1:142315/90‑1 (MQ=255)
                                 gATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGatat                                      >  2:708736/1‑90 (MQ=255)
                                 gATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGatat                                      <  2:159523/90‑1 (MQ=255)
                                       gggTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCtt                                >  1:281720/1‑90 (MQ=255)
                                       gggTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCtt                                >  2:170158/1‑90 (MQ=255)
                                        ggTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCTTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCttt                               <  2:721719/90‑1 (MQ=255)
                                         gTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTc                              <  1:337103/90‑1 (MQ=255)
                                         gTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTc                              <  1:128960/90‑1 (MQ=255)
                                         gTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTc                              <  1:280725/90‑1 (MQ=255)
                                         gTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTc                              <  2:346668/90‑1 (MQ=255)
                                           aTCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCAt                            >  1:647903/1‑90 (MQ=255)
                                           aTCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCAt                            <  2:712817/90‑1 (MQ=255)
                                           aTCTGCGCGCGGAAAGGATACCGCCAAAGAACCCTGTTGCCGAACATTCCGCGCCATTTCGACGTCAGCGCGGGGGATATCCGCTTTCAt                            >  1:653459/1‑90 (MQ=255)
                                                 gggCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATAtt                      >  1:470966/1‑90 (MQ=255)
                                                 gggCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATAtt                      <  2:537798/90‑1 (MQ=255)
                                                      gACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTc                                                    >  2:94326/1‑55 (MQ=255)
                                                      gACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTc                                                    <  1:94326/55‑1 (MQ=255)
                                                          ggTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACg             <  2:593345/90‑1 (MQ=255)
                                                            ttCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGaa           >  1:715943/1‑90 (MQ=255)
                                                            ttCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGaa           <  1:67128/90‑1 (MQ=255)
                                                            ttCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGaa           <  1:662035/90‑1 (MQ=255)
                                                                  cAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACa     <  2:67355/90‑1 (MQ=255)
                                                                     aaGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACt  >  1:70271/1‑90 (MQ=255)
                                                                                  |                                                                            
GGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACT  >  NZ_CP009273/465382‑465540

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 29 ≤ 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
CCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTC  >  NZ_CP009273/465374‑465541
                                                                                          |                                                                             
CCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTC                                                                      <  SRR3722073.329008/100‑1 (MQ=60)
                           agcgtcagatgtgtataagagacagTCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGG                                           <  SRR3722073.96414/75‑1 (MQ=60)
                           CTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGG                                           >  SRR3722073.766834/1‑100 (MQ=60)
                                     GTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTT                                 >  SRR3722073.287674/1‑100 (MQ=60)
                                         GATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAGGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCAT                             <  SRR3722073.145267/100‑1 (MQ=60)
                                         GATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCAT                             >  SRR3722073.661242/1‑100 (MQ=60)
                                               GGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATT                       >  SRR3722073.480991/1‑100 (MQ=60)
                                                 GTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCC                     <  SRR3722073.131661/100‑1 (MQ=60)
                                                 GTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCC                     <  SRR3722073.286661/100‑1 (MQ=60)
                                                 GTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCC                     <  SRR3722073.344258/100‑1 (MQ=60)
                                                          GGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAA            >  SRR3722073.730413/1‑100 (MQ=60)
                                                                   GTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACT   >  SRR3722073.71908/1‑100 (MQ=60)
                                                                    TTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTC  <  SRR3722073.675605/100‑1 (MQ=60)
                                                                    TTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTC  <  SRR3722073.68700/100‑1 (MQ=60)
                                                                                          |                                                                             
CCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTC  >  NZ_CP009273/465374‑465541

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