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 I194 R1 197 32.2 1845014 95.2% 1756453 85.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,196,055 T→C V766V (GTT→GTC BW25113_RS11045 → DUF4132 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,196,0550TC100.0% 81.3 / NA 24V766V (GTT→GTCBW25113_RS11045DUF4132 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/10);  total (14/10)

CCGGCACGCATTCCAAAAGCCCCTGATTTCTGGCATCCGCAAATGTGGTCCAGGCCGCGCCTTATCACTAATAATCAAACTGTTACCAATGACGCTCTGGAAATTATCGGCGAAATGCTGCGCTTTACCCAGGGGGGACGTTTTTATAGCGGGCTGGAACAACTGAA  >  NZ_CP009273/2195972‑2196138
                                                                                   |                                                                                   
ccGGCACGCATTCCAAAAGCCCCTGATTTCTGGCATCCGCAAATGTGGTCCAGGCCGCGCCTTATCACTAATAATCAAACTGTCACCAAt                                                                               >  1:556613/1‑90 (MQ=255)
   gCACGCATTCCAAAAGCCCCTGATTTCTGGCATCCGCAAATGTGGTCCAGGCCGCGCCTTATCACTAATAATCAAACTGTCACCAATGAc                                                                            >  1:711515/1‑90 (MQ=255)
                                    ccGCAAATGTGGTCCAGGCCGCGCCTTATCACTAATAATCAAACTGTCACCAATGACGCTCTGGAAATTATCGGCGAAATGCTGCGCttt                                           <  2:596616/90‑1 (MQ=255)
                                              ggTCCAGGCCGCGCCTTATCACTAATAATCAAACTGTCACCAATGACGCTCTGGAAATTATCGGCGAAATGCTGCGCTTTACCCAggggg                                 >  2:370147/1‑90 (MQ=255)
                                               gTCCAGGCCGCGCCTTATCACTAATAATCAAACTGTCACCAATGACGCTCTGGAAATTATCGGCGAAATGCTGCGCTTTACCCAgggggg                                <  1:820060/90‑1 (MQ=255)
                                                 ccAGGCCGCGCCTTATCACTAATAATCAAACTGTCACCAATGACGCTCTGGAAATTATCGGCGAAATGCTGCGCTTTACCCAGGGGGGAc                              <  2:556613/90‑1 (MQ=255)
                                                 ccAGGCCGCGCCTTATCACTAATAATCAAACTGTCACCAATGACGCTCTGGAAATTATCGGCGAAATGCTGCGCTTTACCCAGGGGGGAc                              <  1:887785/90‑1 (MQ=255)
                                                 ccAGGCCGCGCCTTATCACTAATAATCAAACTGTCACCAATGACGCTCTGGAAATTATCGGCGAAATGCTGCGCTTTACCCAGGGGGGAc                              <  1:630009/90‑1 (MQ=255)
                                                    ggCCGCGCCTTATCACTAATAATCAAACTGTCACCAATGACGCTCTGGAAATTATCGGCGAAATGCTGCGCTTTACCCAGGGGGGACGtt                           >  2:721702/1‑90 (MQ=255)
                                                        gcgcCTTATCACTAATAATCAAACTGTCACCAATGACGCTCTGGAAATTATCGGCGAAAtgct                                                  >  2:203123/1‑63 (MQ=255)
                                                        gcgcCTTATCACTAATAATCAAACTGTCACCAATGACGCTCTGGAAATTATCGGCGAAAtgct                                                  <  1:203123/63‑1 (MQ=255)
                                                            cTTATCACTAATAATCAAACTGTCACCAATGACGCTCTGGAAATTATCGGCGAAATGCTGCGCTTTACCCAGGGGGGACGTTTTTATAGc                   <  2:397018/90‑1 (MQ=255)
                                                              tATCACTAATAATCAAACTGTCACCAATGACGCTCTGGAAATTATCgg                                                           <  2:103820/48‑1 (MQ=255)
                                                              tATCACTAATAATCAAACTGTCACCAATGACGCTCTGGAAATTATCgg                                                           >  1:103820/1‑48 (MQ=255)
                                                                     aataatCAAACTGTCACCAATGACGCTCTGGAAATTATCGGCGAAATGCTGCGCTTTAccc                                       >  1:570993/1‑61 (MQ=255)
                                                                     aataatCAAACTGTCACCAATGACGCTCTGGAAATTATCGGCGAAATGCTGCGCTTTAccc                                       <  2:570993/61‑1 (MQ=255)
                                                                     aataatCAAACTGTCACCAATGACGCTCTGGAAATTATCGGCGAAATGCTGCGCTTTACCCAGGGGGGAc                              >  2:236960/1‑70 (MQ=255)
                                                                     aataatCAAACTGTCACCAATGACGCTCTGGAAATTATCGGCGAAATGCTGCGCTTTACCCAGGGGGGAc                              <  1:236960/70‑1 (MQ=255)
                                                                     aataatCAAACTGTCACCAATGACGCTCTGGAAATTATCGGCGAAATGCTGCGCTTTACCCAGGGGGGACGTTTTTATAGCGGGCTGGaa          >  2:908827/1‑90 (MQ=255)
                                                                        aatCAAACTGTCACCAATGACGCTCTGGAAATTATCGGCGAAATGCTGCGCTTTACCCAGGGGGGACGTTTTTATAGCGGGCTGGaacaa       >  2:486985/1‑90 (MQ=255)
                                                                             aaCTGTCACCAATGACGCTCTGGAAATTATCGGCGAAATGCTGCGCTTTACCCAGGGGGGCCGTTTTTATAGCGGGCTGGAACAACTGaa  >  2:914468/1‑90 (MQ=255)
                                                                             aaCTGTCACCAATGACGCTCTGGAAATTATCGGCGAAATGCTGCGCTTTACCCAGGGGGGACGTTTTTATAGCGGGCTGGAACAACTGaa  >  2:28870/1‑90 (MQ=255)
                                                                             aaCTGTCACCAATGACGCTCTGGAAATTATCGGCGAAATGCTGCGCTTTACCCAGGGGGGACGTTTTTATAGCGGGCTGGAACAACTGaa  >  1:861166/1‑90 (MQ=255)
                                                                             aaCTGTCACCAATGACGCTCTGGAAATTATCGGCGAAATGCTGCGCTTTACCCAGGGGGGACGTTTTTATAGCGGGCTGGAACAACTGaa  >  2:95872/1‑90 (MQ=255)
                                                                                   |                                                                                   
CCGGCACGCATTCCAAAAGCCCCTGATTTCTGGCATCCGCAAATGTGGTCCAGGCCGCGCCTTATCACTAATAATCAAACTGTTACCAATGACGCTCTGGAAATTATCGGCGAAATGCTGCGCTTTACCCAGGGGGGACGTTTTTATAGCGGGCTGGAACAACTGAA  >  NZ_CP009273/2195972‑2196138

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ 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
AATGGACATTTACCCGGCACGCATTCCAAAAGCCCCTGATTTCTGGCATCCGCAAATGTGGTCCAGGCCGCGCCTTATCACTAATAATCAAACTGTTACCAATGACGCTCTGGAAATTATCGGCGAAATGCTGCGCTTTACCCAGGGGGGACGTTTTTATAGCGGGCTGGAACAACTGAAAACGTTCTGC  >  NZ_CP009273/2195959‑2196148
                                                                                                |                                                                                             
AATGGACATTTACCCGGCACGCATTCCAAAAGCCCCTGATTTCTGGCATCCGCAAATGTGGTCCAGGCCGCGCCTTATCACTAATAATCAAACTGTCACC                                                                                            <  SRR3722072.892306/100‑1 (MQ=60)
   GGACATTTACCCGGCACGCATTCCAAAAGCCCCTGATTTCTGGCATCCGCAAATGTGGTCCAGGCCGCGCCTTATCACTAATAATCAAACTGTCACCAAT                                                                                         >  SRR3722072.566618/1‑100 (MQ=60)
      CATTTACCCGGCACGCATTCCAAAAGCCCCTGATTTCTGGCATCCGCAAATGTGGTCCAGGCCGCGCCTTATCACTAATAATCAAACTGTCACCAATGAC                                                                                      >  SRR3722072.724195/1‑100 (MQ=60)
        TTTACCCGGCACGCATTCCAAAAGCCCCTGATTTCTGGCATCCGCAAATGTGGTCCAGGCCGCGCCTTATCACTAATAATCAAACTGTCACCAATGACGC                                                                                    <  SRR3722072.411004/100‑1 (MQ=60)
                  gagataatcggacgtctcgtgggctcggagatgtgtataagagacagGGCCGCGCCTTATCACTAATAATCAAACTGTCACCAATGACGCTCTGGAAATT                                                                          <  SRR3722072.682529/53‑1 (MQ=60)
                                          atgtgtataagagacagGGTCCAGGCCGCGCCTTATCACTAATAATCAAACTGTCACCAATGACGCTCTGGAAATTATCGGCGAAATGCTGCGCTTTACC                                                  <  SRR3722072.206618/83‑1 (MQ=60)
                                                            GTCCAGGCCGCGCCTTATCACTAATAATCAAACTGTCACCAATGACGCTCTGGAAATTATCGGCGAAATGCTGCGCTTTACCCAGGGGGGACGCTTTTAT                                <  SRR3722072.834339/100‑1 (MQ=60)
                                                              taagagacaGCCTTATCACTAATAATCAAACTGTCACCAATGACGCTCTGGAAATTATCGGCGAAATGCTGCGCTTTACCCAGGGGGGACGTTTTTATAG                              <  SRR3722072.240914/91‑1 (MQ=60)
                                                              CCAGGCCGCGCCTTATCACTAATAATCAAACTGTCACCAATGACGCTCTGGAAATTATCGGCGAAATGCTGCGCTTTACCCAGGGGGGACGTTTTTATAG                              <  SRR3722072.641265/100‑1 (MQ=60)
                                                              CCAGGCCGCGCCTTATCACTAATAATCAAACTGTCACCAATGACGCTCTGGAAATTATCGGCGAAATGCTGCGCTTTACCCAGGGGGGACGTTTTTATAG                              <  SRR3722072.903052/100‑1 (MQ=60)
                                                                 GGCCGCGCCTTATCACTAATAATCAAACTGTCACCAATGACGCTCTGGAAATTATCGGCGAAATGCTGCtgtctcttatacacatctccgagcccacgag                           >  SRR3722072.105734/1‑69 (MQ=60)
                                                                        CCTTATCACTAATAATCAAACTGTCACCAATGACGCTCTGGAAATTATCGGCGAAATGCTGCGCTTTACCCAGGGGGGACGctgtctcttatacacatct                    >  SRR3722072.581197/1‑81 (MQ=60)
                                                                                CTAATAATCAAACTGTCACCAATGACGCTCTGGAAATTATCGGCGAAATGCTGCGCTTTACCCAGGGGGGACGTTTTTATAGCGGGCTGGAACAACTGAA            >  SRR3722072.876025/1‑100 (MQ=60)
                                                                                          AACTGTCACCAATGACGCTCTGGAAATTATCGGCGAAATGCTGCGCTTTACCCAGGGGctgtctcttatacacatctccgagcccacgagacgtccgatt  >  SRR3722072.15458/1‑58 (MQ=60)
                                                                                                |                                                                                             
AATGGACATTTACCCGGCACGCATTCCAAAAGCCCCTGATTTCTGGCATCCGCAAATGTGGTCCAGGCCGCGCCTTATCACTAATAATCAAACTGTTACCAATGACGCTCTGGAAATTATCGGCGAAATGCTGCGCTTTACCCAGGGGGGACGTTTTTATAGCGGGCTGGAACAACTGAAAACGTTCTGC  >  NZ_CP009273/2195959‑2196148

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: