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 I227 R1 223 21.3 1167702 97.0% 1132670 86.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,778,934 (T)6→5 intergenic (‑564/+187) BW25113_RS25705 ← / ← BW25113_RS13825 DUF5507 domain‑containing protein/tRNA‑Ile

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,778,9290T.100.0% 107.5 / NA 26intergenic (‑559/+192)BW25113_RS25705/BW25113_RS13825DUF5507 domain‑containing protein/tRNA‑Ile
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (10/16);  total (10/16)

AGGGGTATAATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCA  >  NZ_CP009273/2778841‑2779002
                                                                                        |                                                                         
aGGGGTATAATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCAtt                                                                          >  1:545324/1‑90 (MQ=255)
     tataATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCAtttttcc                                                                     >  2:452095/1‑88 (MQ=255)
            aaTGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACt                                                             >  2:21734/1‑90 (MQ=255)
                  aataCTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGaa                                                       <  1:203347/90‑1 (MQ=255)
                   ataCTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAAt                                                      >  2:544613/1‑90 (MQ=255)
                   ataCTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAAt                                                      >  1:55296/1‑90 (MQ=255)
                        ctGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTa                                                 >  1:416173/1‑90 (MQ=255)
                        ctGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTa                                                 >  2:194077/1‑90 (MQ=255)
                          ggCTTTATCGTTAATTACTTAATTCCACATGTAGCAATTTGCCCGCTTGGCAAAGCAGGCATTTTTCCAGGTACTTTTGAATGAGGACTg                                              >  1:552396/1‑90 (MQ=255)
                             tttATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGAt                                            <  2:322561/90‑1 (MQ=255)
                             tttATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGAt                                            <  1:452095/90‑1 (MQ=255)
                             tttATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGAt                                            <  1:111668/90‑1 (MQ=255)
                                  cGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATa                                       <  2:396584/90‑1 (MQ=255)
                                  cGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATa                                       <  1:21734/90‑1 (MQ=255)
                                  cGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATa                                       <  1:82688/90‑1 (MQ=255)
                                  cGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATa                                       <  2:177035/90‑1 (MQ=255)
                                    ttAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATaaa                                     <  2:211072/90‑1 (MQ=255)
                                       aTTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATAc                                  <  2:16173/90‑1 (MQ=255)
                                       aTTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATAc                                  <  2:20753/90‑1 (MQ=255)
                                       aTTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATAc                                  <  1:379792/90‑1 (MQ=255)
                                          aCTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACAtt                               >  1:252184/1‑90 (MQ=255)
                                                           aGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGt              <  2:545324/90‑1 (MQ=255)
                                                               aTTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGAt                                            >  2:361067/1‑56 (MQ=39)
                                                               aTTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGAt                                            <  1:361067/56‑1 (MQ=255)
                                                                 ttGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCaaa        >  1:419926/1‑90 (MQ=255)
                                                                 ttGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCaaa        <  2:55296/90‑1 (MQ=255)
                                                                      cGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACAcc   <  2:114161/90‑1 (MQ=255)
                                                                       gCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCa  >  2:376735/1‑90 (MQ=255)
                                                                                        |                                                                         
AGGGGTATAATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCA  >  NZ_CP009273/2778841‑2779002

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ 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
ATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAA  >  NZ_CP009273/2778850‑2778996
                                                                                    |                                                              
ATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAAT                                                >  SRR3722111.55959/1‑100 (MQ=60)
     TGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTA                                           >  SRR3722111.422301/1‑100 (MQ=60)
       ATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAGCAATTTGCCCGCTTGGCAAAGCAGGCATTTTTCCAGGTACTTTTGAATGAGGACTG                                        >  SRR3722111.561072/1‑100 (MQ=60)
         AATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGA                                       <  SRR3722111.205861/100‑1 (MQ=60)
                    gcgtcagatgtgtataagagacagCATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATAC                            <  SRR3722111.366081/76‑1 (MQ=60)
                    TTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATAC                            <  SRR3722111.113074/100‑1 (MQ=60)
                    TTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATAC                            <  SRR3722111.458953/100‑1 (MQ=60)
                       ATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATT                         >  SRR3722111.255150/1‑100 (MQ=60)
                         CGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGC                       <  SRR3722111.22004/100‑1 (MQ=60)
                         CGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGC                       <  SRR3722111.83706/100‑1 (MQ=60)
                              ATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGG                  <  SRR3722111.385230/100‑1 (MQ=60)
                                              TGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAA  >  SRR3722111.426132/1‑100 (MQ=60)
                                                                                    |                                                              
ATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAA  >  NZ_CP009273/2778850‑2778996

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