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 I227 R1 223 21.3 1167702 97.0% 1132670 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% 66.7 / NA 21M121V (ATG→GTG) ydjKMFS transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (10/11);  total (10/11)

TCCACGATTTCTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCATCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAATGGCC  >  NZ_CP009273/1854049‑1854202
                                                                             |                                                                            
tCCACGATTTCTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGAAAAGCTCCCAGCCCAACTCCCACCACAAAACggca                                                                  <  2:76531/90‑1 (MQ=255)
 ccACGATTTCTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACggcag                                                                 <  1:186655/90‑1 (MQ=255)
       tttCTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATa                                                           >  2:350650/1‑90 (MQ=255)
         tCTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAg                                                         >  2:361000/1‑90 (MQ=255)
          cTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAgg                                                        <  1:433589/90‑1 (MQ=255)
            aCCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGaa                                                      <  2:203869/90‑1 (MQ=255)
              cGGGTATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAAt                                                    <  2:233517/90‑1 (MQ=255)
                    tatatTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCatat                                              >  1:252768/1‑90 (MQ=255)
                          cGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTg                                        >  2:361611/1‑90 (MQ=255)
                          cGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTg                                        >  2:415521/1‑90 (MQ=255)
                                     gCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCg                             >  1:436905/1‑90 (MQ=255)
                                     gCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCg                             >  1:579289/1‑90 (MQ=255)
                                         aCAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACaa                         <  1:348817/90‑1 (MQ=255)
                                              gTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATc                    <  1:350650/90‑1 (MQ=255)
                                               tAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCg                   <  2:367651/90‑1 (MQ=255)
                                                  ccAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAgg                <  1:24886/90‑1 (MQ=255)
                                                    aGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCa              >  2:356477/1‑90 (MQ=255)
                                                       aaaGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAata           <  1:5828/90‑1 (MQ=255)
                                                               cAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAATGGc   >  2:280669/1‑90 (MQ=255)
                                                               cAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAATGGc   >  2:516199/1‑90 (MQ=255)
                                                                aGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAATGGcc  <  1:356477/90‑1 (MQ=255)
                                                                             |                                                                            
TCCACGATTTCTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCATCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAATGGCC  >  NZ_CP009273/1854049‑1854202

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ 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
GAAACCCGACTTGACCACGTTCCACGATTTCTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCATCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAATGGCCTCATAGAGGAT  >  NZ_CP009273/1854029‑1854213
                                                                                                 |                                                                                       
GAAACCCGACTTGACCACGTTCCACGATTTCTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCA                                                                                       <  SRR3722111.254751/100‑1 (MQ=60)
                     gagatctacacattgacattcgtcggcagcgtcagatgtgtataagagacagGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGA                                                                  <  SRR3722111.398796/48‑1 (MQ=60)
                     CCACGATTTCTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGA                                                                  <  SRR3722111.188984/100‑1 (MQ=60)
                              accgagatctacacattgacattcgtcggcagcgtcagatgtgtataagagacagGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATAT                                                         <  SRR3722111.518696/45‑1 (MQ=60)
                              CTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATAT                                                         >  SRR3722111.255742/1‑100 (MQ=60)
                              CTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATAT                                                         <  SRR3722111.440066/100‑1 (MQ=60)
                                               GGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCG                                        >  SRR3722111.443448/1‑100 (MQ=60)
                                               GGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCG                                        >  SRR3722111.588375/1‑100 (MQ=60)
                                                             ACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGC                          <  SRR3722111.353535/100‑1 (MQ=60)
                                                                  GTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATAT                     <  SRR3722111.355412/100‑1 (MQ=60)
                                                                      CCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAAT                 <  SRR3722111.25201/100‑1 (MQ=60)
                                                                           AAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAATGGCCT            <  SRR3722111.5897/100‑1 (MQ=60)
                                                                                    AGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAATGGCCTCATAGAGGA   <  SRR3722111.361379/100‑1 (MQ=60)
                                                                                     GCCCAACTCCCACCACAAAACGTCTGGCAATTAGTAAATCCATATTCGGTGAGAAACCACCGACAACCATCGGGGCAATATGAATGGCCTgataacggat  >  SRR3722111.400081/1‑90 (MQ=60)
                                                                                                 |                                                                                       
GAAACCCGACTTGACCACGTTCCACGATTTCTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCATCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAATGGCCTCATAGAGGAT  >  NZ_CP009273/1854029‑1854213

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