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,579,168 T→C E339G (GAA→GGA)  ggt ← gamma‑glutamyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,579,1680TC100.0% 107.1 / NA 31E339G (GAA→GGA) ggtgamma‑glutamyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/12);  total (19/12)

TGATCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATTCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGCCGTATTTCTTCATAT  >  NZ_CP009273/3579087‑3579252
                                                                                 |                                                                                    
tGATCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCgg                                                                              >  1:269949/1‑90 (MQ=255)
     ggCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCggc                                                                         <  2:194868/90‑1 (MQ=255)
      gCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCggcg                                                                        >  2:423518/1‑90 (MQ=255)
       cAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCggcgt                                                                       >  1:209246/1‑90 (MQ=255)
       cAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCggcgt                                                                       <  2:43769/90‑1 (MQ=255)
       cAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCggcgt                                                                       <  2:225369/90‑1 (MQ=255)
                tGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTAtt                                                              >  2:547979/1‑90 (MQ=255)
                tGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTAtt                                                              >  2:437203/1‑90 (MQ=255)
                tGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTAtt                                                              >  2:307284/1‑90 (MQ=255)
                tGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTAtt                                                              >  2:103092/1‑90 (MQ=255)
                     tAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTcc                                                         >  1:75483/1‑90 (MQ=255)
                           ttATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCt                                                   >  1:121449/1‑90 (MQ=255)
                                ggTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCAt                                              >  2:70630/1‑90 (MQ=255)
                                ggTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCAt                                              >  2:109537/1‑90 (MQ=255)
                                 gTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATg                                             >  1:15149/1‑90 (MQ=255)
                                  tCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCAt                                              <  1:743407/88‑1 (MQ=255)
                                  tCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCAt                                              >  2:743407/1‑88 (MQ=255)
                                    aGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGAtt                                          <  1:70630/90‑1 (MQ=255)
                                     gcgcCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGAttt                                         >  1:44028/1‑90 (MQ=255)
                                           gCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTgcatcg                                   >  2:335456/1‑90 (MQ=255)
                                           gCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTgcatcg                                   >  2:640959/1‑90 (MQ=255)
                                                   aCTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGgcg                           <  1:547979/90‑1 (MQ=255)
                                                   aCTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGgcg                           <  1:386640/90‑1 (MQ=255)
                                                   aCTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGgcg                           >  1:24977/1‑90 (MQ=255)
                                                           aaaaTCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCaaa                   <  1:423518/90‑1 (MQ=255)
                                                           aaaaTCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCaaa                   <  2:44028/90‑1 (MQ=255)
                                                             aaTCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGc                 <  2:67014/90‑1 (MQ=255)
                                                              aTCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGcc                <  2:15149/90‑1 (MQ=255)
                                                              aTCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGcc                <  2:4328/90‑1 (MQ=255)
                                                                       gtcAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGCCGTATttctt       >  1:87640/3‑90 (MQ=255)
                                                                            gATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGCCGTATTTCTTCatat  >  1:276785/1‑90 (MQ=255)
                                                                                 |                                                                                    
TGATCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATTCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGCCGTATTTCTTCATAT  >  NZ_CP009273/3579087‑3579252

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 31 ≤ 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
GCTTTATTGATATCAATTTGATCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATTCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGCCGTATTTCTTCATATCGAAGTTTTCCAG  >  NZ_CP009273/3579069‑3579265
                                                                                                   |                                                                                                 
GCTTTATTGATATCAGTTTGATCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGGACTTTGAAAAAATCCGGGTCGCCAAGATATC                                                                                                   >  SRR3722073.635796/1‑100 (MQ=60)
 CTTTATTGATATCAATTTGATCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCC                                                                                                  <  SRR3722073.392473/100‑1 (MQ=60)
 CTTTATTGATATCAATTTGATCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCC                                                                                                  <  SRR3722073.242756/100‑1 (MQ=60)
 CTTTATTGATATCAACTTGATCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCC                                                                                                  <  SRR3722073.310016/100‑1 (MQ=60)
        GATATCAATTTGATCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGG                                                                                           >  SRR3722073.275658/1‑100 (MQ=60)
               ATTTGATCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGT                                                                                    >  SRR3722073.213532/1‑100 (MQ=60)
                             AGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCC                                                                      >  SRR3722073.77228/1‑100 (MQ=60)
                                   GGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCT                                                                >  SRR3722073.124028/1‑100 (MQ=60)
                                         GGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATG                                                          >  SRR3722073.15535/1‑100 (MQ=60)
                                             TTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTT                                                      >  SRR3722073.45083/1‑100 (MQ=60)
                                                  GGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATC                                                 <  SRR3722073.758471/100‑1 (MQ=60)
                                                      AGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCAT                                             <  SRR3722073.72276/100‑1 (MQ=60)
                                                           CTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCG                                        >  SRR3722073.25608/1‑100 (MQ=60)
                                                                     ACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGC                              <  SRR3722073.394872/100‑1 (MQ=60)
                                                                     ACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGC                              <  SRR3722073.559420/100‑1 (MQ=60)
                                                                             AAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGCCGTATTTC                      <  SRR3722073.432583/100‑1 (MQ=60)
                                                                               AATCCTGGTGGTCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGCCGTATTTCTT                    >  SRR3722073.89605/1‑100 (MQ=60)
                                                                                    GGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGCCGTATTTCTTCATAT               >  SRR3722073.282644/1‑100 (MQ=60)
                                                                                            AAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGCCGTATTTCTTCctgtctcttata       >  SRR3722073.150984/1‑88 (MQ=60)
                                                                                                 ATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGCCGTATTTCTTCATATCGAAGTTTTCCAG  >  SRR3722073.652415/1‑100 (MQ=60)
                                                                                                   |                                                                                                 
GCTTTATTGATATCAATTTGATCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATTCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGCCGTATTTCTTCATATCGAAGTTTTCCAG  >  NZ_CP009273/3579069‑3579265

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