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 I231 R1 218 12.8 703680 97.2% 683976 86.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,073,163 A→T D26V (GAT→GTT)  BW25113_RS20185 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,073,1630AT100.0% 92.3 / NA 24D26V (GAT→GTT) BW25113_RS20185hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (14/10);  total (14/10)

GATATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGATCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAAC  >  NZ_CP009273/4073080‑4073244
                                                                                   |                                                                                 
gaTATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGa                                                                             >  2:246184/1‑90 (MQ=255)
     ttATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGc                                                                        <  1:246184/90‑1 (MQ=255)
     ttATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGc                                                                        <  2:215194/90‑1 (MQ=255)
     ttATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGc                                                                        <  2:16150/90‑1 (MQ=255)
         gAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCg                                                                    >  2:270737/1‑90 (MQ=255)
           aCTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTa                                                                  >  2:261536/1‑90 (MQ=255)
              ctctGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATc                                                               >  1:42430/1‑90 (MQ=255)
              ctctGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATc                                                               >  2:17796/1‑90 (MQ=255)
                        tataGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGc                                                     <  1:104049/90‑1 (MQ=255)
                           aGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGa                                                  <  2:7698/90‑1 (MQ=255)
                                tGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGAc                                                 <  1:134923/86‑1 (MQ=255)
                                tGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGAc                                                 >  2:134923/1‑86 (MQ=255)
                                 gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                            >  2:222819/1‑90 (MQ=255)
                                 gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                            >  1:308313/1‑90 (MQ=255)
                                 gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                            >  1:168640/1‑90 (MQ=255)
                                 gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                            >  2:279019/1‑90 (MQ=255)
                                      gaactttACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAggg                                       >  2:292778/5‑90 (MQ=255)
                                             aCTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGt                                <  2:42430/90‑1 (MQ=255)
                                              cTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTa                               <  1:160713/90‑1 (MQ=255)
                                              cTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTa                               <  1:141825/90‑1 (MQ=255)
                                                  aTTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATc                           >  2:18613/1‑90 (MQ=255)
                                                                    ttAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCgg                                                   <  1:45705/48‑1 (MQ=255)
                                                                    ttAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCgg                                                   >  2:45705/1‑48 (MQ=255)
                                                                           aCTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAAc  >  2:209337/1‑90 (MQ=255)
                                                                                   |                                                                                 
GATATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGATCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAAC  >  NZ_CP009273/4073080‑4073244

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AAAATGATATTGGCGAGATATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGATCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACAGCAAGAACCAGTGTTC  >  NZ_CP009273/4073064‑4073261
                                                                                                   |                                                                                                  
ggcagcgtcagatgtgtataagagacagCTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGT                                                                                                    <  SRR3722116.34548/72‑1 (MQ=60)
        ATTGGCGAGATATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGG                                                                                            <  SRR3722116.89458/100‑1 (MQ=60)
                    TTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATC                                                                                >  SRR3722116.42932/1‑100 (MQ=60)
                     TTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCT                                                                               <  SRR3722116.249623/100‑1 (MQ=60)
                                       GTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATT                                                             >  SRR3722116.170623/1‑100 (MQ=60)
                                       GTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATT                                                             >  SRR3722116.312867/1‑100 (MQ=60)
                                        TATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTA                                                            <  SRR3722116.105324/100‑1 (MQ=60)
                                          cgtcggcagcgtcagatgtgtataagagacaGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAG                                                          <  SRR3722116.46253/69‑1 (MQ=60)
                                            GATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGG                                                        <  SRR3722116.136498/100‑1 (MQ=60)
                                                              CTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACC                                      <  SRR3722116.143469/100‑1 (MQ=60)
                                                              CTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACC                                      <  SRR3722116.162542/100‑1 (MQ=60)
                                                                                                ATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACAGCAAGAACCAGTGT    >  SRR3722116.337054/1‑100 (MQ=60)
                                                                                                  GTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACAGCAAGAACCAGTGTTC  <  SRR3722116.6502/100‑1 (MQ=60)
                                                                                                   |                                                                                                  
AAAATGATATTGGCGAGATATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGATCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACAGCAAGAACCAGTGTTC  >  NZ_CP009273/4073064‑4073261

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