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 I228 R1 214 26.5 1448192 97.1% 1406194 86.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,854,126 T→C M121V (ATG→GTG)  ydjK ← MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,854,1260TC100.0% 79.6 / NA 25M121V (ATG→GTG) ydjKMFS transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (11/14);  total (11/14)

TGACCACGTTCCACGATTTCTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCATCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAATGGCCTCATAGA  >  NZ_CP009273/1854040‑1854209
                                                                                      |                                                                                   
tGACCACGTTCCACGATTTCTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCcaccac                                                                                  >  1:473205/1‑90 (MQ=255)
tGACCACGTTCCACGATTTCTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCcaccac                                                                                  <  1:672939/90‑1 (MQ=255)
          ccACGATTTCTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACggcag                                                                        >  2:488177/1‑90 (MQ=255)
                               tatTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTc                                                   <  2:528894/90‑1 (MQ=255)
                                  tCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGt                                                <  1:243366/90‑1 (MQ=255)
                                   cGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTg                                               >  2:486425/1‑90 (MQ=255)
                                   cGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTg                                               >  2:330085/1‑90 (MQ=255)
                                   cGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTg                                               >  2:559843/1‑90 (MQ=255)
                                         aaCCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAg                                         >  2:552954/1‑90 (MQ=255)
                                          aCCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGc                                        <  1:486425/90‑1 (MQ=255)
                                                     gCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCa                             >  2:345103/1‑90 (MQ=255)
                                                     gCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCa                             >  1:648209/1‑90 (MQ=255)
                                                       gTAACCAGTAAAGCTCCCCGCCCATCCCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATc                           <  1:716563/90‑1 (MQ=255)
                                                       gTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATc                           <  1:241399/90‑1 (MQ=255)
                                                            cAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGc                      <  2:479586/90‑1 (MQ=255)
                                                              gTAAAGCTCCCAGCCCAACTCCCACCACAAAACggc                                                                          >  2:403334/1‑36 (MQ=39)
                                                              gTAAAGCTCCCAGCCCAACTCCCACCACAAAACggc                                                                          <  1:403334/36‑1 (MQ=39)
                                                                aaaGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAata                  <  1:559433/90‑1 (MQ=255)
                                                                   gCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGa               >  2:473865/1‑90 (MQ=255)
                                                                     tCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAAt             <  1:501674/90‑1 (MQ=255)
                                                                         aGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAATGGcc         <  2:186459/90‑1 (MQ=255)
                                                                         aGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAATGGcc         <  2:287281/90‑1 (MQ=255)
                                                                         aGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAATGGcc         <  1:670714/90‑1 (MQ=255)
                                                                         aGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAATGGcc         <  1:219911/90‑1 (MQ=255)
                                                                              aaCTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAATGGCCTCATa    >  2:702325/1‑90 (MQ=255)
                                                                               aCTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAATGGCCTCATag   >  2:165166/1‑90 (MQ=255)
                                                                                cTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAATGGCCTCATaga  <  1:488177/90‑1 (MQ=255)
                                                                                      |                                                                                   
TGACCACGTTCCACGATTTCTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCATCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAATGGCCTCATAGA  >  NZ_CP009273/1854040‑1854209

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 32 ≤ 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
AAACCCGACTTGACCACGTTCCACGATTTCTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCATCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAATGGCCTCATAGAGGATGAACGC  >  NZ_CP009273/1854030‑1854219
                                                                                                |                                                                                             
AAACCCGACTTGACCACGTTCCACGATTTCTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCAC                                                                                            <  SRR3722112.373656/100‑1 (MQ=60)
AAACCCGACTTGACCACGTTCCACGATTTCTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCAC                                                                                            >  SRR3722112.479963/1‑100 (MQ=60)
          TGACCACGTTCCACGATTTCTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAG                                                                                  <  SRR3722112.683449/100‑1 (MQ=60)
                  cacattgacattcgtcggcagcgtcagatgtgtataagagacagAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAG                                                                          <  SRR3722112.408659/56‑1 (MQ=60)
                                            TCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCAC                                                <  SRR3722112.246251/100‑1 (MQ=60)
                                                    ACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACC                                        <  SRR3722112.493431/100‑1 (MQ=60)
                                                     CCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCA                                       >  SRR3722112.658310/1‑100 (MQ=60)
                                                                 GTAACCAGTAAAGCTCCCCGCCCATCCCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATAC                           <  SRR3722112.727748/100‑1 (MQ=60)
                                                                 GTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATAT                           <  SRR3722112.244265/100‑1 (MQ=60)
                                                                  gatgtgtataagagacagGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATG                          <  SRR3722112.544229/82‑1 (MQ=60)
                                                                          AAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAATGGCCT                  <  SRR3722112.567919/100‑1 (MQ=60)
                                                                               TCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAATGGCCTCATAG             <  SRR3722112.508947/100‑1 (MQ=60)
                                                                                   AGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAATGGCCTCATAGAGGA         <  SRR3722112.222560/100‑1 (MQ=60)
                                                                                   AGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAATGGCCTCATAGAGGA         <  SRR3722112.681187/100‑1 (MQ=60)
                                                                                          CTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAATGGCCTCATAGAGGATGAACGC  <  SRR3722112.495217/100‑1 (MQ=60)
                                                                                          CTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAATGGCCTCATAGAGGATGAACGC  >  SRR3722112.60353/1‑100 (MQ=60)
                                                                                                |                                                                                             
AAACCCGACTTGACCACGTTCCACGATTTCTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCATCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAATGGCCTCATAGAGGATGAACGC  >  NZ_CP009273/1854030‑1854219

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