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 I207 R1 218 45.0 2482936 97.1% 2410930 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,011,010 A→G E531G (GAA→GGA)  ygfK → putative selenate reductase subunit YgfK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,011,0100AG100.0% 162.5 / NA 47E531G (GAA→GGA) ygfKputative selenate reductase subunit YgfK
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/22);  total (25/22)

TTTGCGATCACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGAAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGATT  >  NZ_CP009273/3010923‑3011086
                                                                                        |                                                                            
tttGCGATACCAGCGCCAAACACAACTGTCCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGgaa                                                                            <  2:85369/90‑1 (MQ=255)
       tccccAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTg                                                                     <  1:606606/87‑1 (MQ=255)
          ccAGTCCAAAAAAAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGa                                                                  <  1:634543/90‑1 (MQ=255)
           cAGTGCAAAAAAAACTGTCCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGAt                                                                 <  2:711002/90‑1 (MQ=255)
           cAGTCCAAATACAACTGTCCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGAt                                                                 <  2:576379/90‑1 (MQ=255)
                ccAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAata                                                            <  2:120854/90‑1 (MQ=255)
                        aaCTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGgaaaa                                                                          >  2:868962/1‑69 (MQ=255)
                        aaCTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGgaaaa                                                                          <  1:868962/69‑1 (MQ=255)
                        aaCTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCaa                                                     >  1:951696/1‑90 (MQ=255)
                        aaCTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCaa                                                     >  2:910438/1‑90 (MQ=255)
                        aaCTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGCTGGGAGGAATATAAGCaa                                                     >  1:449727/1‑90 (MQ=255)
                         aCTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAAc                                                    >  1:345494/1‑90 (MQ=255)
                                 cGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGcac                                            >  2:825057/1‑90 (MQ=255)
                                 cGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGcac                                            >  1:802468/1‑90 (MQ=255)
                                 cGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGcac                                            >  1:26686/1‑90 (MQ=255)
                                      ggATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAAcc                                       >  2:81711/1‑90 (MQ=255)
                                          tACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGcc                                   <  2:961863/90‑1 (MQ=255)
                                           aCGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCg                                  <  2:133299/90‑1 (MQ=255)
                                                 gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGgaaaa                                                                          >  1:60962/1‑44 (MQ=255)
                                                 gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGgaaaa                                                                          <  2:60962/44‑1 (MQ=255)
                                                 gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGAt                                                                 >  2:222055/1‑53 (MQ=255)
                                                 gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGAt                                                                 <  1:222055/53‑1 (MQ=255)
                                                 gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAata                                                            >  2:1198558/1‑58 (MQ=255)
                                                 gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAata                                                            <  1:1198558/58‑1 (MQ=255)
                                                 gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAAc                                        <  1:563982/78‑1 (MQ=255)
                                                 gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAAc                                        >  2:563982/1‑78 (MQ=255)
                                                 gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTc                              >  1:949364/1‑88 (MQ=255)
                                                 gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTc                              <  2:949364/88‑1 (MQ=255)
                                                 gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCt                             >  2:346476/1‑89 (MQ=255)
                                                 gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCt                             <  1:346476/89‑1 (MQ=255)
                                                 gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg                            >  1:944698/1‑90 (MQ=255)
                                                 gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg                            >  1:827221/1‑90 (MQ=255)
                                                 gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg                            >  2:272623/1‑90 (MQ=255)
                                                 gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg                            >  1:382802/1‑90 (MQ=255)
                                                   gcgcTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGt                          >  2:1109813/1‑90 (MQ=255)
                                                     gcTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTc                        >  1:748930/1‑90 (MQ=255)
                                                     gcTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTc                        <  2:183174/90‑1 (MQ=255)
                                                      cTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCa                       <  2:1151697/90‑1 (MQ=255)
                                                       tGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCAc                      <  1:66777/90‑1 (MQ=255)
                                                        gAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTg                                                                     <  2:305477/42‑1 (MQ=255)
                                                        gAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTg                                                                     >  1:305477/1‑42 (MQ=255)
                                                         aaTATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGc                    <  2:116657/90‑1 (MQ=255)
                                                         aaTATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGc                    <  2:345494/90‑1 (MQ=255)
                                                           tatCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCa                  >  1:1143943/1‑90 (MQ=255)
                                                             tCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATc                >  1:702163/1‑90 (MQ=255)
                                                             tCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATc                >  1:608037/1‑90 (MQ=255)
                                                                 cgACCTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGt            <  2:586948/90‑1 (MQ=255)
                                                                  gAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGtt           >  2:435582/1‑90 (MQ=255)
                                                                  gAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGtt           >  1:291655/1‑90 (MQ=255)
                                                                   aaCTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTg          <  2:827221/90‑1 (MQ=255)
                                                                      tGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCg       <  2:438414/90‑1 (MQ=255)
                                                                           aaaGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGAtt  <  2:988467/90‑1 (MQ=255)
                                                                                        |                                                                            
TTTGCGATCACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGAAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGATT  >  NZ_CP009273/3010923‑3011086

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ 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
GTCATATTTGCGATCACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGAAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTT  >  NZ_CP009273/3010917‑3011077
                                                                                             |                                                                   
GTCATTTTTGCGATCACCAGTCCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGG                                                               <  SRR3722089.945004/100‑1 (MQ=60)
 TCATATTTGCGACCCCCAGTTCAAAAACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGT                                                              <  SRR3722089.107774/100‑1 (MQ=60)
 TCATATTGGCGACCCCCAGTGCAAAAACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGT                                                              <  SRR3722089.280213/100‑1 (MQ=60)
        aaaagagacggGTGCCAATACAACTGTCCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATG                                                       <  SRR3722089.881425/89‑1 (MQ=60)
           gatcctacgttcaaatccaaccgtcCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAAT                                                    <  SRR3722089.12222/75‑1 (MQ=60)
             TCCCCAGTCCCAAGACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATAT                                                  <  SRR3722089.614133/100‑1 (MQ=60)
                ccagtccaaaaaaAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAG                                               <  SRR3722089.642625/87‑1 (MQ=60)
                 gcagcgtcagatgtgtataagagacaGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGC                                              <  SRR3722089.224873/74‑1 (MQ=60)
                   GTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAA                                            >  SRR3722089.965674/1‑100 (MQ=60)
                   GTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGCTGGGAGGAATATAAGCAA                                            >  SRR3722089.455121/1‑100 (MQ=60)
                    TGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAAC                                           >  SRR3722089.349792/1‑100 (MQ=60)
                      gtcagatgtgtataagagacaGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAAAGC                                         <  SRR3722089.1216782/79‑1 (MQ=60)
                            CAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCAC                                   >  SRR3722089.813603/1‑100 (MQ=60)
                            CAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCAC                                   >  SRR3722089.27009/1‑100 (MQ=60)
                                          AGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTC                     <  SRR3722089.570594/100‑1 (MQ=60)
                                            GATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGctgtctcttatacacatctgacgctgccgacgacca                   >  SRR3722089.61708/1‑64 (MQ=60)
                                            GATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTG                   >  SRR3722089.387467/1‑100 (MQ=60)
                                            GATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTG                   >  SRR3722089.963315/1‑100 (MQ=60)
                                            GATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTG                   >  SRR3722089.958548/1‑100 (MQ=60)
                                            GATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTG                   >  SRR3722089.838889/1‑100 (MQ=60)
                                                ACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTC               >  SRR3722089.758832/1‑100 (MQ=60)
                                                   ACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATctgtctcttatacacatctgacgctgccgacgaccatc            >  SRR3722089.309357/1‑62 (MQ=60)
                                                     AGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCC          <  SRR3722089.350783/100‑1 (MQ=60)
                                                      GTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCA         >  SRR3722089.1161353/1‑100 (MQ=60)
                                                        GCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATC       >  SRR3722089.711315/1‑100 (MQ=60)
                                                        GCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATC       >  SRR3722089.615592/1‑100 (MQ=60)
                                                            TGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGT   <  SRR3722089.67586/100‑1 (MQ=60)
                                                             GAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTT  >  SRR3722089.295385/1‑100 (MQ=60)
                                                                                             |                                                                   
GTCATATTTGCGATCACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGAAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTT  >  NZ_CP009273/3010917‑3011077

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