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,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% 209.6 / NA 50intergenic (‑559/+192)BW25113_RS25705/BW25113_RS13825DUF5507 domain‑containing protein/tRNA‑Ile
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (32/18);  total (32/18)

AGGGGTATAATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCAGCCCT  >  NZ_CP009273/2778841‑2779007
                                                                                        |                                                                              
aGGGGTATAATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCAtt                                                                               >  1:124136/1‑90 (MQ=255)
aGGGGTATAATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATCTGCCCGCTTGGCATAGCAGGCAtt                                                                               >  2:428861/1‑90 (MQ=255)
 ggggTATAATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCAttt                                                                              <  2:192053/90‑1 (MQ=255)
           aaaTGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTAc                                                                   >  2:228360/1‑90 (MQ=255)
            aacgATAATACTCTGGCTTAATCGTTAATTACTTAATGCCACATGTAAGCAATTTCCCCCCTTGACATAGCAGGCATTTTTCCATGTACt                                                                  >  2:33692/4‑90 (MQ=255)
              tGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACttt                                                                >  1:290551/1‑90 (MQ=255)
              tGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACttt                                                                >  2:82427/1‑90 (MQ=255)
                  aataCTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGaa                                                            <  2:443007/90‑1 (MQ=255)
                    taCTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATg                                                          <  1:60204/90‑1 (MQ=255)
                    taCTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATg                                                          <  2:838457/90‑1 (MQ=255)
                       tctGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGt                                                       >  2:495440/1‑90 (MQ=255)
                       tctGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGt                                                       >  2:785968/1‑90 (MQ=255)
                        ctGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGt                                                                     <  2:374049/75‑1 (MQ=255)
                        ctGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGt                                                                     >  1:374049/1‑75 (MQ=255)
                        ctGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGaa                                                            >  1:423969/1‑84 (MQ=255)
                        ctGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGaa                                                            <  2:423969/84‑1 (MQ=255)
                        ctGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTa                                                      >  2:669623/1‑90 (MQ=255)
                        ctGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTa                                                      >  2:23448/1‑90 (MQ=255)
                        ctGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTa                                                      >  1:103079/1‑90 (MQ=255)
                        ctGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTa                                                      >  1:555836/1‑90 (MQ=255)
                        ctGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTa                                                      >  1:456670/1‑90 (MQ=255)
                        ctGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTa                                                      >  1:56076/1‑90 (MQ=255)
                         tGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTAc                                                     <  1:869241/90‑1 (MQ=255)
                             tttATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGAt                                                 <  1:542449/90‑1 (MQ=255)
                             tttATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGAt                                                 <  2:124136/90‑1 (MQ=255)
                             tttATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGAt                                                 <  2:722818/90‑1 (MQ=255)
                             tttATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGAt                                                 <  2:578421/90‑1 (MQ=255)
                              ttATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATg                                                >  2:733006/1‑90 (MQ=255)
                               tATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATgg                                               >  1:647639/1‑90 (MQ=255)
                               tATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATgg                                               >  2:881202/1‑90 (MQ=255)
                                  cGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATa                                            <  1:754795/90‑1 (MQ=255)
                                    ttAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGt                                                                     <  2:213921/63‑1 (MQ=255)
                                    ttAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGt                                                                     >  1:213921/1‑63 (MQ=255)
                                    ttAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGAt                                                 <  1:855484/83‑1 (MQ=255)
                                    ttAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGAt                                                 >  2:855484/1‑83 (MQ=255)
                                    ttAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATaaa                                          >  1:76189/1‑90 (MQ=255)
                                    ttAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATaaa                                          >  1:507186/1‑90 (MQ=255)
                                    ttAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATaaa                                          >  1:124494/1‑90 (MQ=255)
                                       aTTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATAc                                       <  1:655798/90‑1 (MQ=255)
                                               aTTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGt                               >  2:126454/1‑90 (MQ=255)
                                                   cacaTGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACCTTTGAATGAGTACTAATGGATAAATACATTGCAGTGGCg                           >  2:392715/1‑90 (MQ=255)
                                                    acaTGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATaaa                                          >  1:595235/1‑74 (MQ=255)
                                                    acaTGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATaaa                                          <  2:595235/74‑1 (MQ=255)
                                                    acaTGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGt                          >  1:307551/1‑90 (MQ=255)
                                                    acaTGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGt                          >  1:781184/1‑90 (MQ=255)
                                                    acaTGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGt                          <  1:723234/90‑1 (MQ=255)
                                                    acaTGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGAACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGt                          >  2:915690/1‑90 (MQ=255)
                                                       tGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGcc                       >  1:365928/1‑90 (MQ=255)
                                                         tAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCAc                     >  2:625133/1‑90 (MQ=255)
                                                          aaGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACg                    >  2:47004/1‑90 (MQ=255)
                                                          aaGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACg                    >  2:437447/1‑90 (MQ=255)
                                                           aGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGt                   <  2:647639/90‑1 (MQ=255)
                                                                            gCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCAGCCCt  <  2:781184/90‑1 (MQ=255)
                                                                                        |                                                                              
AGGGGTATAATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCAGCCCT  >  NZ_CP009273/2778841‑2779007

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ 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
AATGAAGGGGTATAATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCAGCCCTCATTCGAAACCACCCACCGCACTTC  >  NZ_CP009273/2778836‑2779032
                                                                                                  |                                                                                                  
gaagacggcatacgagataatcggacgtctcgtgggctcggagatgtgtataagagacagTGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCC                                                                                                   <  SRR3722072.144757/40‑1 (MQ=60)
         GTATAATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTT                                                                                         >  SRR3722072.295571/1‑100 (MQ=60)
                   TGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATctgtc                                                                               >  SRR3722072.380714/1‑95 (MQ=60)
                   TGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTA                                                                               >  SRR3722072.104982/1‑100 (MQ=60)
                   TGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTA                                                                               >  SRR3722072.57100/1‑100 (MQ=60)
                   TGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTA                                                                               >  SRR3722072.431678/1‑100 (MQ=60)
                   TGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTA                                                                               >  SRR3722072.465031/1‑100 (MQ=60)
                   TGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTA                                                                               >  SRR3722072.565828/1‑100 (MQ=60)
                         TACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATG                                                                         <  SRR3722072.61302/100‑1 (MQ=60)
                          ACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGG                                                                        >  SRR3722072.659201/1‑100 (MQ=60)
                              TGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAA                                                                    <  SRR3722072.884212/100‑1 (MQ=60)
                               GGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATctgtctcttatacacat                                                                   >  SRR3722072.217551/1‑83 (MQ=60)
                               GGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAA                                                                   >  SRR3722072.516321/1‑100 (MQ=60)
                               GGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAA                                                                   >  SRR3722072.77605/1‑100 (MQ=60)
                               GGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAA                                                                   >  SRR3722072.126758/1‑100 (MQ=60)
                                  TTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATAC                                                                <  SRR3722072.552216/100‑1 (MQ=60)
                                  TTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATAC                                                                <  SRR3722072.870250/100‑1 (MQ=60)
                                       CGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGC                                                           <  SRR3722072.768119/100‑1 (MQ=60)
                                            ATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGG                                                      <  SRR3722072.667503/100‑1 (MQ=60)
                                               ACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGctgtct                                                   >  SRR3722072.605853/1‑94 (MQ=60)
                                               ACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGT                                                   >  SRR3722072.312886/1‑100 (MQ=60)
                                               ACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGT                                                   >  SRR3722072.794885/1‑100 (MQ=60)
                                                  TAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCC                                                >  SRR3722072.372416/1‑100 (MQ=60)
                                                         ACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACC                                         <  SRR3722072.736090/100‑1 (MQ=60)
                                                                               TGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCAGCCCTCATTCGAA                   >  SRR3722072.585254/1‑100 (MQ=60)
                                                                                            caTTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCAGCCCTCATTCGAAACCACCCACCGC       <  SRR3722072.799741/98‑1 (MQ=60)
                                                                                                 TTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCAGCCCTCATTCGAAACCACCCACCGCACTTC  >  SRR3722072.744375/1‑100 (MQ=60)
                                                                                                  |                                                                                                  
AATGAAGGGGTATAATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCAGCCCTCATTCGAAACCACCCACCGCACTTC  >  NZ_CP009273/2778836‑2779032

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