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 I199 R1 220 57.5 3163680 97.0% 3068769 86.1

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% 143.4 / NA 34intergenic (+74/‑240)BW25113_RS10330/yeeJtRNA‑Asn/inverse autotransporter adhesin YeeJ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (16/18);  total (16/18)

CAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCG  >  NZ_CP009273/2038095‑2038255
                                                                                    |                                                                            
cAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGttttaa                                                                         <  1:23313/90‑3 (MQ=255)
 agagGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGttttaac                                                                        <  1:625252/90‑4 (MQ=255)
          aaaTTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAgc                                                              <  2:1093819/90‑1 (MQ=255)
          aaaTTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAgc                                                              <  1:873749/90‑1 (MQ=255)
           aaTTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAgcg                                                             <  2:495429/90‑1 (MQ=255)
                 aaaaaTTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCg                                                       <  1:895533/90‑1 (MQ=255)
                 aaaaaTTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCg                                                       <  2:473653/90‑1 (MQ=255)
                 aaaaaTTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCg                                                       <  2:698905/90‑1 (MQ=255)
                    aaTTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGAcc                                                    <  1:529566/90‑1 (MQ=255)
                    aaTTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGAcc                                                    <  2:1069583/90‑1 (MQ=255)
                                       aTGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGttt                                 <  1:568794/90‑1 (MQ=255)
                                       aTGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGttt                                 >  2:713250/1‑90 (MQ=255)
                                           cACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAgc                                                              >  2:700003/1‑57 (MQ=255)
                                           cACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAgc                                                              <  1:700003/57‑1 (MQ=255)
                                               gACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGa                         >  1:1512146/1‑90 (MQ=255)
                                               gACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGa                         <  1:1471299/90‑1 (MQ=255)
                                               gACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGa                         >  2:740940/1‑90 (MQ=255)
                                               gACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGa                         >  2:8818/1‑90 (MQ=255)
                                                 ccTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAAc                       >  2:1001000/1‑90 (MQ=255)
                                                 ccTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAAc                       >  2:826305/1‑90 (MQ=255)
                                                   ttAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGc                     >  1:706000/1‑90 (MQ=255)
                                                        tGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTGTGAACGCCATCa                <  1:458979/90‑1 (MQ=255)
                                                        tGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCa                <  1:8818/90‑1 (MQ=255)
                                                         gAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGAcc                                                    >  2:1007022/1‑53 (MQ=255)
                                                         gAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGAcc                                                    <  1:1007022/53‑1 (MQ=255)
                                                         gAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCAc               >  2:1538459/1‑90 (MQ=255)
                                                         gAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCAc               >  2:71638/1‑90 (MQ=255)
                                                              ttGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCAttt          >  2:213196/1‑90 (MQ=255)
                                                              ttGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCAttt          >  1:208280/1‑90 (MQ=255)
                                                              ttGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCAttt          >  1:1354077/1‑90 (MQ=255)
                                                                   ttttAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTcccc     >  1:1456634/1‑90 (MQ=255)
                                                                     ttAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTc   <  2:1456634/90‑1 (MQ=255)
                                                                      tAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCg  <  2:1258361/90‑1 (MQ=255)
                                                                      tAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCg  <  1:1313174/90‑1 (MQ=255)
                                                                               gCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTcccc     >  2:1411514/1‑78 (MQ=255)
                                                                               gCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTcccc     <  1:1411514/78‑1 (MQ=255)
                                                                                    |                                                                            
CAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCG  >  NZ_CP009273/2038095‑2038255

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
CAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAG  >  NZ_CP009273/2038095‑2038281
                                                                                        |                                                                                                  
CAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGC                                                                                        <  SRR3722078.23607/100‑1 (MQ=60)
 AGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCG                                                                                       <  SRR3722078.632448/100‑1 (MQ=60)
          cgtcagatgtgtataagagacagAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACC                                                                              <  SRR3722078.708434/77‑1 (MQ=60)
          AAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACC                                                                              <  SRR3722078.885470/100‑1 (MQ=60)
                 AAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGG                                                                       <  SRR3722078.907566/100‑1 (MQ=60)
                    gcagcgtcagatgtgtataagagacagGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTG                                                                    <  SRR3722078.1021807/73‑1 (MQ=60)
                    AATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTG                                                                    <  SRR3722078.535823/100‑1 (MQ=60)
                                     CAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGA                                                   >  SRR3722078.1536809/1‑100 (MQ=60)
                                       ATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAAC                                                 <  SRR3722078.575461/100‑1 (MQ=60)
                                         GTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGC                                               >  SRR3722078.714544/1‑100 (MQ=60)
                                               GACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCAC                                         <  SRR3722078.1495305/100‑1 (MQ=60)
                                                    TAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTT                                    >  SRR3722078.210797/1‑100 (MQ=60)
                                                    TAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTT                                    >  SRR3722078.1375811/1‑100 (MQ=60)
                                                        TGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTGTGAACGCCATCACCATTTTCCC                                <  SRR3722078.464461/100‑1 (MQ=60)
                                                        TGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCC                                <  SRR3722078.8910/100‑1 (MQ=60)
                                                         GAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCC                               >  SRR3722078.1480380/1‑100 (MQ=60)
                                                                   agTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAA                     <  SRR3722078.1434421/98‑1 (MQ=60)
                                                                      TAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAAC                  <  SRR3722078.1334154/100‑1 (MQ=60)
                                                                                   GGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCA      >  SRR3722078.140331/1‑100 (MQ=60)
                                                                                       TTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAG  <  SRR3722078.765039/100‑1 (MQ=60)
                                                                                       TTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCATCGATTATAAAACTTGAGTTATTCAGTAG  <  SRR3722078.1222864/100‑1 (MQ=60)
                                                                                        |                                                                                                  
CAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAG  >  NZ_CP009273/2038095‑2038281

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