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 I228 R1 214 26.5 1448192 97.1% 1406194 86.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,439,101 T→C K255E (AAG→GAG)  mepA ← penicillin‑insensitive murein endopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,439,1010TC100.0% 70.5 / NA 22K255E (AAG→GAG) mepApenicillin‑insensitive murein endopeptidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/8);  total (14/8)

AAACAGGCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTTCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGG  >  NZ_CP009273/2439016‑2439176
                                                                                     |                                                                           
aaaCAGGCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCa                                                                         <  2:429417/90‑1 (MQ=255)
      acTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCttt                                                                   >  1:302385/2‑90 (MQ=255)
               cGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTccc                                                          <  2:302385/90‑1 (MQ=255)
                 tttCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCgg                                                        >  2:479690/1‑90 (MQ=255)
                              aCGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCa                                           >  1:166795/1‑90 (MQ=255)
                              aCGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCa                                           >  1:227680/1‑90 (MQ=255)
                              aCGTGCTCATCCAGTAGCGCCCGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCa                                           >  2:545130/1‑90 (MQ=255)
                                 tGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAAc                                        >  1:163010/1‑90 (MQ=255)
                                  gCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAAcc                                       >  2:466153/1‑90 (MQ=255)
                                   cTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCa                                                                         <  1:515230/55‑1 (MQ=255)
                                   cTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCa                                                                         >  2:515230/1‑55 (MQ=255)
                                   cTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCa                                      >  2:99936/1‑90 (MQ=255)
                                     cATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGc                                    <  2:558188/90‑1 (MQ=255)
                                          aGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGc                               <  2:405115/90‑1 (MQ=255)
                                          aGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGc                               <  2:169141/90‑1 (MQ=255)
                                             aGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGt                            >  2:391219/1‑90 (MQ=255)
                                               cgcCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAAcc                                       >  1:545273/1‑77 (MQ=255)
                                               cgcCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAAcc                                       <  2:545273/77‑1 (MQ=255)
                                               cgcCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTc                          >  1:562834/1‑90 (MQ=255)
                                                gcCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCt                         >  2:428460/1‑90 (MQ=255)
                                                gcCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCt                         >  1:287422/1‑90 (MQ=255)
                                                                       aGGCGGTGTCCTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATgg  <  1:207095/90‑1 (MQ=255)
                                                                                     |                                                                           
AAACAGGCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTTCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGG  >  NZ_CP009273/2439016‑2439176

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ 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
CCATAAACAGGCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTTCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACAC  >  NZ_CP009273/2439012‑2439198
                                                                                         |                                                                                                 
CCATAAACAGACTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTT                                                                                         >  SRR3722112.305865/1‑100 (MQ=60)
    agcgtcagatgtgtataagagacagAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTG                                                                                     <  SRR3722112.522775/75‑1 (MQ=60)
                        CCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCA                                                                 >  SRR3722112.168845/1‑100 (MQ=60)
                        CCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCA                                                                 >  SRR3722112.230420/1‑100 (MQ=60)
                           TTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAAC                                                              >  SRR3722112.165014/1‑100 (MQ=60)
                                      tacacattgacattcgtcggcagcgtcagatgtgtataagagacagCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAG                                                   <  SRR3722112.660348/54‑1 (MQ=60)
                                         CATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGctg                                                >  SRR3722112.553465/1‑97 (MQ=60)
                                         CATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTC                                                >  SRR3722112.571376/1‑100 (MQ=60)
                                          ATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCT                                               >  SRR3722112.290757/1‑100 (MQ=60)
                                                                           AGGCGGTGTCCTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGT              <  SRR3722112.209603/100‑1 (MQ=60)
                                                                                    CTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCA     >  SRR3722112.612575/1‑100 (MQ=60)
                                                                                       CTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAGCCATCGCCTGATGGCGGTAAAGGTTGATCTTCACAC  <  SRR3722112.555656/100‑1 (MQ=60)
                                                                                       CTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACAC  <  SRR3722112.553319/100‑1 (MQ=60)
                                                                                         |                                                                                                 
CCATAAACAGGCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTTCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACAC  >  NZ_CP009273/2439012‑2439198

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