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 I205 R1 220 17.5 961422 97.1% 933540 86.5

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,852,078 T→G T467P (ACG→CCG)  yidK ← solute:sodium symporter family transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,852,0780TG100.0% 80.4 / NA 24T467P (ACG→CCG) yidKsolute:sodium symporter family transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (11/13);  total (11/13)

CCGCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGTACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGATTATGCCAATCCCCATCGC  >  NZ_CP009273/3851991‑3852162
                                                                                       |                                                                                    
ccGCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGAc                                                                                    <  1:131752/90‑1 (MQ=255)
ccGCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGAc                                                                                    <  2:404186/90‑1 (MQ=255)
ccGCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGAc                                                                                    <  1:207727/90‑1 (MQ=255)
     aaCGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGcc                                                                               <  1:437834/90‑1 (MQ=255)
               tGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATaga                                                                     >  2:235506/1‑90 (MQ=255)
                 aaGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAgg                                                                   <  1:364135/90‑1 (MQ=255)
                     tGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAAt                                                               >  2:465654/1‑90 (MQ=255)
                      gAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATg                                                              >  1:442620/1‑90 (MQ=255)
                        aCGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGa                                                            >  1:474801/1‑90 (MQ=255)
                           gCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGt                                                         >  1:271763/1‑90 (MQ=255)
                                gcgcgcGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAAc                                                    <  2:371068/90‑1 (MQ=255)
                                gcgcgcGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAAc                                                    <  2:26416/90‑1 (MQ=255)
                                          tGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGAt                                          <  1:368225/90‑1 (MQ=255)
                                          tGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGAt                                          <  1:370961/90‑1 (MQ=255)
                                                    gATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGt                                >  1:423358/1‑90 (MQ=255)
                                                    gATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGt                                >  1:76875/1‑90 (MQ=255)
                                                     aTCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTg                               <  1:166245/90‑1 (MQ=255)
                                                                cACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGAt                    <  1:123900/90‑1 (MQ=255)
                                                                 aCGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGAtt                   <  2:423358/90‑1 (MQ=255)
                                                                       aCGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGc                        <  1:450849/79‑1 (MQ=255)
                                                                       aCGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGc                        >  2:450849/1‑79 (MQ=255)
                                                                            gATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGATTATGCCAATccc        >  1:182619/1‑90 (MQ=255)
                                                                                  gAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGATTATGCCAATCCCCATCGc  >  1:4637/1‑90 (MQ=255)
                                                                                  gAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGATTATGCCAATCCCCATCGc  >  1:357414/1‑90 (MQ=255)
                                                                                       |                                                                                    
CCGCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGTACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGATTATGCCAATCCCCATCGC  >  NZ_CP009273/3851991‑3852162

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ 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
GTTTCATGTCCACCGCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGTACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGATTATGCCAATCCCCATCGC  >  NZ_CP009273/3851979‑3852162
                                                                                                   |                                                                                    
GTTTCATGTCCACCGCAAACGCATCTTTGAAGGTGAATGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGG                                                                                      >  SRR3722087.104401/1‑100 (MQ=60)
            CCGCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCAC                                                                          <  SRR3722087.132985/100‑1 (MQ=60)
            CCGCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCAC                                                                          <  SRR3722087.209540/100‑1 (MQ=60)
                 AACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGA                                                                     <  SRR3722087.442952/100‑1 (MQ=60)
                        CTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATG                                                              >  SRR3722087.447793/1‑100 (MQ=60)
                          TTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGA                                                            >  SRR3722087.480323/1‑100 (MQ=60)
                             AAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGT                                                         >  SRR3722087.274435/1‑100 (MQ=60)
                             AAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGT                                                         <  SRR3722087.368264/100‑1 (MQ=60)
                                                      TGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGT                                <  SRR3722087.372416/100‑1 (MQ=60)
                                                      TGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGT                                >  SRR3722087.77586/1‑100 (MQ=60)
                                                      TGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGT                                >  SRR3722087.428292/1‑100 (MQ=60)
                                                      TGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGT                                <  SRR3722087.375195/100‑1 (MQ=60)
                                                                 ATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGA                     <  SRR3722087.167760/100‑1 (MQ=60)
                                                                        GCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGATTATGCC              <  SRR3722087.456108/100‑1 (MQ=60)
                                                                            CACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGATTATGCCAATC          <  SRR3722087.125071/100‑1 (MQ=60)
                                                                              CGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGATTATGCCAATCCC        >  SRR3722087.184237/1‑100 (MQ=60)
                                                                                    CGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGATTATGCCAATCCCCATCGC  >  SRR3722087.361447/1‑100 (MQ=60)
                                                                                    CGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGATTATGCCAATCCCCATCGC  >  SRR3722087.4681/1‑100 (MQ=60)
                                                                                                   |                                                                                    
GTTTCATGTCCACCGCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGTACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGATTATGCCAATCCCCATCGC  >  NZ_CP009273/3851979‑3852162

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