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 502,414 T→C D42G (GAT→GGT)  ybaK ← Cys‑tRNA(Pro)/Cys‑tRNA(Cys) deacylase YbaK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273502,4140TC96.2% 79.0 / ‑3.3 26D42G (GAT→GGT) ybaKCys‑tRNA(Pro)/Cys‑tRNA(Cys) deacylase YbaK
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (13/12);  total (14/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.47e-01

CACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGATCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTAT  >  NZ_CP009273/502356‑502486
                                                          |                                                                        
cACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTc                                           >  1:165838/1‑90 (MQ=255)
cACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTc                                           >  2:1031670/1‑90 (MQ=255)
           gtTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAAt                                <  2:435925/90‑1 (MQ=255)
            tttCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAAtt                               >  1:117201/1‑90 (MQ=255)
            tttCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAAtt                               >  1:1043555/1‑90 (MQ=255)
            tttCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAAtt                               >  2:1100208/1‑90 (MQ=255)
            tttCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAAtt                               >  2:888216/1‑90 (MQ=255)
                 tatCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGttt                          <  1:930839/90‑1 (MQ=255)
                   tCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGtt                           <  1:460132/87‑1 (MQ=255)
                   tCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGtt                           <  1:373571/87‑1 (MQ=255)
                   tCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGtt                           >  2:373571/1‑87 (MQ=255)
                   tCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGtt                           >  2:460132/1‑87 (MQ=255)
                   tCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCa                        >  1:991755/1‑90 (MQ=255)
                     aCCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGc                      >  2:1081791/1‑90 (MQ=255)
                      ccGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGcc                     <  2:73515/90‑1 (MQ=255)
                        gTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTc                                                                  >  1:823443/1‑43 (MQ=255)
                        gTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTc                                                                  <  2:823443/43‑1 (MQ=255)
                          tCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGAt                 >  1:49152/1‑90 (MQ=255)
                            aCTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCg               <  2:1121354/90‑1 (MQ=255)
                             cTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGt              >  2:93609/1‑90 (MQ=255)
                                   ccaGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGt        <  2:372797/90‑1 (MQ=255)
                                    cagcagCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTa       <  2:387775/90‑1 (MQ=255)
                                       cagcGTTTTGTAGACCTGATCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGt    >  1:363812/1‑90 (MQ=255)
                                         gcGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTAt  <  2:311236/90‑1 (MQ=255)
                                         gcGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTAt  <  2:85536/90‑1 (MQ=255)
                                         gcGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTAt  <  2:953874/90‑1 (MQ=255)
                                                          |                                                                        
CACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGATCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTAT  >  NZ_CP009273/502356‑502486

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ACTTTTTTAAGATCCAGTTGACCTGCGACCGGCGTAACGGCCACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGATCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTATGGATTTGAAACGAAATCTTGTTTTT  >  NZ_CP009273/502315‑502511
                                                                                                   |                                                                                                 
ACTTTTTTAAGATCCAGTTGACCTGCGACCGGCGTAACGGCCACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGAC                                                                                                   >  SRR3722089.967890/1‑100 (MQ=60)
           GTCCAGTTGACCTGCGACCGGCGTAACGGCCACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAA                                                                                        <  SRR3722089.417090/100‑1 (MQ=60)
                               GCGTAACGGCCACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTC                                                                    >  SRR3722089.167915/1‑100 (MQ=60)
                                           CGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATT                                                        >  SRR3722089.1059183/1‑100 (MQ=60)
                                           CGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATT                                                        >  SRR3722089.118625/1‑100 (MQ=60)
                                                  GTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCA                                                 >  SRR3722089.1006469/1‑100 (MQ=60)
                                                       TCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTctgtctcttatacacatctgacgctgccgacgaccat                                            >  SRR3722089.835038/1‑63 (MQ=60)
                                                         ATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGAT                                          <  SRR3722089.378153/100‑1 (MQ=60)
                                                         ATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGAT                                          <  SRR3722089.465628/100‑1 (MQ=60)
                                                         ATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGAT                                          >  SRR3722089.49763/1‑100 (MQ=60)
                                                          TATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATC                                         <  SRR3722089.944438/100‑1 (MQ=60)
                                                                      CTGCCACCAGCAGCGTTTTGTAGACCTGATCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGT                             >  SRR3722089.368294/1‑100 (MQ=60)
                                                                                       TTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTATctgtctcttatacac            >  SRR3722089.650204/1‑85 (MQ=60)
                                                                                       TTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTATGGATTTGAAACGAAA            >  SRR3722089.462442/1‑100 (MQ=60)
                                                                                       TTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTATGGATTTGAAACGAAA            >  SRR3722089.50355/1‑100 (MQ=60)
                                                                                                 GACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCtgtctcttatacacatctgacgctgccgacgaccatctt  >  SRR3722089.388964/1‑61 (MQ=60)
                                                                                                   |                                                                                                 
ACTTTTTTAAGATCCAGTTGACCTGCGACCGGCGTAACGGCCACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGATCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTATGGATTTGAAACGAAATCTTGTTTTT  >  NZ_CP009273/502315‑502511

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