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 I210 R1 226 19.6 1130256 95.4% 1078264 84.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,753,169 T→C S276P (TCA→CCA)  yfjI → YfjI family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,753,1690TC100.0% 87.1 / NA 27S276P (TCA→CCA) yfjIYfjI family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/14);  total (13/14)

ATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATGAGAGCATCATTAAAATTAATGAAAA  >  NZ_CP009273/2753091‑2753245
                                                                              |                                                                            
attggattTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTCCAAATGAACAt                                                                   >  1:474191/1‑90 (MQ=255)
      ttttttGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTCCAAATGAACATTTGCCg                                                             >  1:27225/1‑90 (MQ=255)
                  tGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTCCAAATGAACATTTGCCGGTATTTCACCaa                                                 <  1:21450/90‑1 (MQ=255)
                  tGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTCCAAATGAACATTTGCCGGTATTTCACCaa                                                 <  2:466944/90‑1 (MQ=255)
                  tGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTCCAAATGAACATTTGCCGGTATTTCACCaa                                                 <  2:530711/90‑1 (MQ=255)
                  tGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTCCAAATGAACATTTGCCGGTATTTCACCaa                                                 <  1:305711/90‑1 (MQ=255)
                           tgCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTCCAAATGAACATTTGCCGGTAtt                                                        <  2:21804/74‑1 (MQ=255)
                           tgCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTCCAAATGAACATTTGCCGGTAtt                                                        >  1:21804/1‑74 (MQ=255)
                           tgCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTCCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATg                                        >  2:106415/1‑90 (MQ=255)
                              cAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTCCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGaa                                     <  2:505645/90‑1 (MQ=255)
                                gCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTCCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAAt                                   <  2:279068/90‑1 (MQ=255)
                                    gCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTCCAAATGAACATTTGCCGGTATTTCACCaa                                                 <  1:235396/72‑1 (MQ=255)
                                    gCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTCCAAATGAACATTTGCCGGTATTTCACCaa                                                 >  2:235396/1‑72 (MQ=255)
                                         aaCACAAGGTAACAGAAAAATTTCCAACCCAATTTTTCCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATga                          >  2:73175/1‑90 (MQ=255)
                                             caAGGTAACAGAAAAATTTCCAACCCAATTTTTCCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATGAGAGc                      <  1:127332/90‑1 (MQ=255)
                                             caAGGTAACAGAAAAATTTCCAACCCAATTTTTCCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATGAGAGc                      <  2:94073/90‑1 (MQ=255)
                                               aGGTAACAGAAAAATTTCCAACCCAATTTTTCCAAATGa                                                                       >  1:240502/1‑39 (MQ=39)
                                               aGGTAACAGAAAAATTTCCAACCCAATTTTTCCAAATGa                                                                       <  2:240502/39‑1 (MQ=39)
                                               aGGTAACAGAAAAATTTCCAACCCAATTTTTCCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATGAGAGcat                    >  2:457105/1‑90 (MQ=255)
                                               aGGTAACAGAAAAATTTCCAACCCAATTTTTCCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATGAGAGcat                    >  2:142997/1‑90 (MQ=255)
                                               aGGTAACAGAAAAATTTCCAACCCAATTTTTCCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATGAGAGcat                    >  2:41978/1‑90 (MQ=255)
                                               aGGTAACAGAAAAATTTCCAACCCAATTTTTCCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATGAGAGcat                    <  1:371930/90‑1 (MQ=255)
                                               aGGTAACAGAAAAATTTCCAACCCAATTTTTCCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATGAGAGcat                    >  1:132858/1‑90 (MQ=255)
                                                    aCAGAAAAATTTCCAACCCAATTTTTCCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATGAGAGCATCATTa               >  1:28518/1‑90 (MQ=255)
                                                                 cAACCCAATTTTTCCAAATGAACATTTGCCGGTATTTCACCaa                                                 >  2:215387/1‑43 (MQ=255)
                                                                 cAACCCAATTTTTCCAAATGAACATTTGCCGGTATTTCACCaa                                                 <  1:215387/43‑1 (MQ=255)
                                                                 cAACCCAATTTTTCCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATGAGAGCATCATTAAAATTAATGaaaa  <  2:27225/90‑1 (MQ=255)
                                                                              |                                                                            
ATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATGAGAGCATCATTAAAATTAATGAAAA  >  NZ_CP009273/2753091‑2753245

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATGAGAGCATCATTAAAATT  >  NZ_CP009273/2753081‑2753237
                                                                                        |                                                                    
GGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTCCAAATGAACAT                                                           >  SRR3722092.483597/1‑100 (MQ=60)
      GGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTCCAAATGAACATTTGCCG                                                     >  SRR3722092.27697/1‑100 (MQ=60)
                           GTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTCCAAATGAACATTTGCCGGTATTTCACCAACGTCTgtct                                >  SRR3722092.22167/1‑96 (MQ=60)
                            atggtcgtcggcagcgtcagatgtgtataagagacagGAAAAATTTCCAACCCAATTTTTCCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGG                               <  SRR3722092.218932/63‑1 (MQ=60)
                            agagacagGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTCCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGG                               <  SRR3722092.239239/92‑1 (MQ=60)
                            TGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTCCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGG                               <  SRR3722092.311091/100‑1 (MQ=60)
                            TGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTCCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGC                               <  SRR3722092.21811/100‑1 (MQ=60)
                                               CTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTCCAAATGAACATTTGCCGctgtctcttatacacatctgacgctgccgacgaccatctta            >  SRR3722092.244408/1‑59 (MQ=60)
                                               CTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTCCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATGAGAGCAT            >  SRR3722092.135104/1‑100 (MQ=60)
                                                    ACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTCCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATGAGAGCATCATTA       >  SRR3722092.29016/1‑100 (MQ=60)
                                                       CAAGGTAACAGAAAAATTTCCAACCCAATTTTTCCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATGAGAGCATCATTAAAA    <  SRR3722092.129484/100‑1 (MQ=60)
                                                         AGGTAACAGAAAAATTTCCAACCCAATTTTTCCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATGAGAGCATCATTAAAATT  <  SRR3722092.378940/100‑1 (MQ=60)
                                                                                        |                                                                    
GGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATGAGAGCATCATTAAAATT  >  NZ_CP009273/2753081‑2753237

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ 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: