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 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% 126.0 / NA 37E531G (GAA→GGA) ygfKputative selenate reductase subunit YgfK
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/18);  total (19/18)

ATCACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGAAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGATTGGTGCA  >  NZ_CP009273/3010929‑3011092
                                                                                 |                                                                                  
aTCCCCAGCGCAAATACAACTGTCCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGtt                                                                            <  2:436810/90‑1 (MQ=255)
  cACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTgg                                                                          >  2:744809/1‑90 (MQ=255)
  cACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTgg                                                                          >  1:712863/1‑90 (MQ=255)
  cACCAGTGCCAATACAACTGCACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGCAAAAAGTCGCGCTGGGAAAAgggagg                                                                          >  1:301513/1‑86 (MQ=255)
    ccAGTGCCAATACAACTGTCCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGa                                                                        <  2:668554/90‑1 (MQ=255)
    ccAGTGCAAATAAAACTGTCCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGa                                                                        <  2:249549/90‑1 (MQ=255)
     cggTGCCAATACAACTGCACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGAt                                                                       <  1:222388/88‑1 (MQ=255)
     cAGTCCAAATAAAACTGTCCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGAt                                                                       <  2:155121/90‑1 (MQ=255)
     cAGTCCAAAAACACCTGTCCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGAt                                                                       <  2:487503/90‑1 (MQ=255)
       gTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGa                                                                     >  1:415500/1‑90 (MQ=255)
       gTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGa                                                                     >  2:519952/1‑90 (MQ=255)
          caaaTACAACTGTCCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAata                                                                  <  2:727625/88‑1 (MQ=255)
            aaTACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAAtata                                                                >  1:577328/1‑90 (MQ=255)
                caccTGTCCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCa                                                            <  1:528880/87‑1 (MQ=255)
                aaccTGTCCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCa                                                            <  2:311962/87‑1 (MQ=255)
                 aaCTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCaa                                                           >  1:713615/1‑90 (MQ=255)
                  aCTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAAc                                                          >  1:236435/1‑90 (MQ=255)
                     gTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCt                                                       >  1:671441/1‑90 (MQ=255)
                       aCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGTGAAAAGGTTGGGATGAATATAAGCAACGCTgg                                                     >  1:441009/1‑90 (MQ=255)
                          cGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGcac                                                  <  1:748154/90‑1 (MQ=255)
                          cGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGcac                                                  >  1:656629/1‑90 (MQ=255)
                          cGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGcac                                                  >  2:392027/1‑90 (MQ=255)
                           gCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGcaca                                                 <  1:440668/90‑1 (MQ=255)
                             cTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACaaa                                               >  1:48590/1‑90 (MQ=255)
                              tGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAAc                                              <  2:541242/90‑1 (MQ=255)
                                   tACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGcc                                         <  1:644913/90‑1 (MQ=255)
                                         aGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAGCCAGCCGGTTCt                                   >  1:748013/1‑90 (MQ=255)
                                         aGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCt                                   >  1:650512/1‑90 (MQ=255)
                                          gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGGTGAATATAAGCAAGGCTGGCACAAACCAGCCGGTTCTg                                  >  2:154923/1‑90 (MQ=255)
                                          gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg                                  >  1:96024/1‑90 (MQ=255)
                                                    tatCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGAt                                                                       <  2:565834/43‑1 (MQ=255)
                                                    tatCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGAt                                                                       >  1:565834/1‑43 (MQ=255)
                                                       ccGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATcc                     <  2:498285/90‑1 (MQ=255)
                                                                    aaggTCGCGCTGGGAAACGGTTGGGATGAATACAAGCAGCGCTGGCGCAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGAtt        <  1:427547/87‑1 (MQ=255)
                                                                    aaaGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAGCCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGAtt        <  2:748013/90‑1 (MQ=255)
                                                                    aaaGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGAtt        <  1:744809/90‑1 (MQ=255)
                                                                          gcgcTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGATTGGTGCa  >  2:551185/1‑90 (MQ=255)
                                                                                 |                                                                                  
ATCACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGAAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGATTGGTGCA  >  NZ_CP009273/3010929‑3011092

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 27 ≤ 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
TATTTGCGATCACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGAAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGATTGGTGCAGGTC  >  NZ_CP009273/3010921‑3011096
                                                                                         |                                                                                      
TATTTGCGATCACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGG                                                                              >  SRR3722073.727269/1‑100 (MQ=60)
TATTTGCGATCACCAGTGCCAATACAACTGCACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGCAAAAAGTCGCGCTGGGAAAAGGgagg                                                                              >  SRR3722073.307887/1‑96 (MQ=60)
     GCGATCACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGA                                                                         >  SRR3722073.424377/1‑100 (MQ=60)
          CACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATA                                                                    >  SRR3722073.589383/1‑100 (MQ=60)
             gggacaaaaaaaAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGC                                                                 <  SRR3722073.392215/88‑1 (MQ=60)
             CGGTGCCAATACAACTGCACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGC                                                                 <  SRR3722073.226997/100‑1 (MQ=60)
               GTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAA                                                               >  SRR3722073.728035/1‑100 (MQ=60)
                TGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAAC                                                              >  SRR3722073.241380/1‑100 (MQ=60)
                   CAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCT                                                           >  SRR3722073.685176/1‑100 (MQ=60)
                     ATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGTGAAAAGGTTGGGATGAATATAAGCAACGCTGG                                                         >  SRR3722073.450433/1‑100 (MQ=60)
                        CACCTGTCCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCAC                                                      <  SRR3722073.539950/100‑1 (MQ=60)
                        CAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCAC                                                      >  SRR3722073.670097/1‑100 (MQ=60)
                           CTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAA                                                   >  SRR3722073.49733/1‑100 (MQ=60)
                                  CGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCG                                            <  SRR3722073.763322/100‑1 (MQ=60)
                                   GCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGG                                           <  SRR3722073.450086/100‑1 (MQ=60)
                                       GGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAGCCAGCCGGTTCT                                       >  SRR3722073.763179/1‑100 (MQ=60)
                                       GGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCT                                       >  SRR3722073.663884/1‑100 (MQ=60)
                                        GATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTG                                      >  SRR3722073.98138/1‑100 (MQ=60)
                                           TACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTT                                   <  SRR3722073.658202/100‑1 (MQ=60)
                                                  GTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCctgtctcttatacacatctgacgctgccgacgatatg                            >  SRR3722073.577646/1‑63 (MQ=60)
                                                                            AAGGTCGCGCTGGGAAACGGTTGGGATGAATACAAGCAGCGCTGGCGCAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGATTGGTGCAGGTC  <  SRR3722073.436701/100‑1 (MQ=60)
                                                                            AAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGATTGGTGCAGGTC  <  SRR3722073.759896/100‑1 (MQ=60)
                                                                                         |                                                                                      
TATTTGCGATCACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGAAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGATTGGTGCAGGTC  >  NZ_CP009273/3010921‑3011096

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