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 I195 R1 222 27.1 1551108 96.0% 1489063 86.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,038,183 (T)5→4 intergenic (+78/‑236) BW25113_RS10330 → / → yeeJ tRNA‑Asn/inverse autotransporter adhesin YeeJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,038,1790T.100.0% 110.4 / NA 27intergenic (+74/‑240)BW25113_RS10330/yeeJtRNA‑Asn/inverse autotransporter adhesin YeeJ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (12/15);  total (12/15)

CAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTC  >  NZ_CP009273/2038095‑2038254
                                                                                    |                                                                           
cAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGttttaa                                                                        <  2:699473/90‑3 (MQ=255)
 agagGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGttttaac                                                                       <  1:696015/90‑4 (MQ=255)
       gCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTa                                                                >  1:762863/1‑90 (MQ=255)
        ccAAATTCTAAAAATTTGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTaa                                                               >  2:63520/1‑90 (MQ=255)
                 aaaaaTTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCg                                                      <  2:722275/90‑1 (MQ=255)
                  aaaaTTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGa                                                     <  2:703222/90‑1 (MQ=255)
                  aaaaTTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGa                                                     <  1:509578/90‑1 (MQ=255)
                  aaaaTTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGa                                                     <  1:677844/90‑1 (MQ=255)
                       tCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCttt                                                >  2:94912/1‑90 (MQ=255)
                                               gACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGa                        >  2:273908/1‑90 (MQ=255)
                                               gACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGa                        <  1:94912/90‑1 (MQ=255)
                                                aCCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGaa                       <  2:699652/90‑1 (MQ=255)
                                                 ccTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAAc                      <  1:773439/90‑1 (MQ=255)
                                                   ttAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGc                    >  2:542901/1‑90 (MQ=255)
                                                   ttAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGc                    >  1:605050/1‑90 (MQ=255)
                                                        tGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCa               <  2:11159/90‑1 (MQ=255)
                                                         gAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCAc              >  2:210673/1‑90 (MQ=255)
                                                         gAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCAc              >  1:541256/1‑90 (MQ=255)
                                                           aCATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATc                >  2:414062/1‑86 (MQ=255)
                                                           aCATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATc                <  1:414062/86‑1 (MQ=255)
                                                              ttGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCAttt         >  1:745980/1‑90 (MQ=255)
                                                              ttGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCAttt         >  1:135772/1‑90 (MQ=255)
                                                                gTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGa                                                     >  2:281757/1‑44 (MQ=38)
                                                                gTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGa                                                     <  1:281757/44‑1 (MQ=38)
                                                                 ttttttAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTcc      <  1:750306/90‑1 (MQ=255)
                                                                   ttttAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTcccc    <  1:122197/90‑1 (MQ=255)
                                                                   ttttAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTcccc    <  1:67053/90‑1 (MQ=255)
                                                                   ttttAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTcccc    <  1:421828/90‑1 (MQ=255)
                                                                     ttAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTc  <  2:741914/90‑1 (MQ=255)
                                                                                    |                                                                           
CAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTC  >  NZ_CP009273/2038095‑2038254

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 33 ≤ 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
GAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGT  >  NZ_CP009273/2038086‑2038279
                                                                                                 |                                                                                                
GAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAAC                                                                                               <  SRR3722073.267590/100‑1 (MQ=60)
GAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAAC                                                                                               <  SRR3722073.109134/100‑1 (MQ=60)
 AGTCCACTCAGAGGAGCCAAATTCTAAAAAATCGTTTTTTTAGCGCAATGTCACTGACCCCAGTTGACCATTGATTCTTAACGGATAGCGGGTTTTTATC                                                                                               <  SRR3722073.188045/100‑1 (MQ=60)
      AGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTA                                                                                         >  SRR3722073.778263/1‑100 (MQ=60)
          AGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCG                                                                                     <  SRR3722073.710152/100‑1 (MQ=60)
                           atatcgtcggcagcgtcagatgtgtataagagacagGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGT                                                                    <  SRR3722073.287711/64‑1 (MQ=60)
                           AAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGT                                                                    <  SRR3722073.691697/100‑1 (MQ=60)
                           AAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGT                                                                    <  SRR3722073.520317/100‑1 (MQ=60)
                                                  GTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGC                                             >  SRR3722073.617647/1‑100 (MQ=60)
                                                        GACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCAC                                       >  SRR3722073.552567/1‑100 (MQ=60)
                                                        GACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCAC                                       <  SRR3722073.97012/100‑1 (MQ=60)
                                                          CCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCA                                     <  SRR3722073.789031/100‑1 (MQ=60)
                                                             TAGTTGAGCATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTT                                  >  SRR3722073.761098/1‑100 (MQ=60)
                                                             TAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTT                                  >  SRR3722073.138598/1‑100 (MQ=60)
                                                                TTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCC                               <  SRR3722073.422910/100‑1 (MQ=60)
                                                                          TTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTAT                     <  SRR3722073.765502/100‑1 (MQ=60)
                                                                            TTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAA                   <  SRR3722073.430860/100‑1 (MQ=60)
                                                                            TTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAA                   <  SRR3722073.68621/100‑1 (MQ=60)
                                                                            TTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAA                   <  SRR3722073.124787/100‑1 (MQ=60)
                                                                                              TTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGT  >  SRR3722073.296704/1‑100 (MQ=60)
                                                                                                 |                                                                                                
GAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGT  >  NZ_CP009273/2038086‑2038279

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