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 I229 R1 214 17.4 943020 97.5% 919444 87.1

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% 87.8 / 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)

TTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGATCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACA  >  NZ_CP009273/4073085‑4073245
                                                                              |                                                                                  
ttATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGc                                                                         >  1:132515/1‑90 (MQ=255)
  aTGACCTCTCTAGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGCGGTGCAg                                                                       <  2:330420/90‑1 (MQ=255)
         ctctGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATc                                                                >  2:174777/1‑90 (MQ=255)
                   tataGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGc                                                      >  1:16116/1‑90 (MQ=255)
                     taGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCgg                                                    <  2:433622/90‑1 (MQ=255)
                       gATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGAc                                                  <  2:346611/90‑1 (MQ=255)
                        atatGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGAcc                                                 <  1:216723/90‑1 (MQ=255)
                            gggCAGAATTTTACTCGGTTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGc                                                      >  1:336096/1‑81 (MQ=255)
                            gggCAGAATTTTACTCGGTTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGc                                                      <  2:336096/81‑1 (MQ=255)
                            gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                             >  1:15396/1‑90 (MQ=255)
                            gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                             >  2:92723/1‑90 (MQ=255)
                            gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                             >  1:205524/1‑90 (MQ=255)
                            gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                             >  1:181171/1‑90 (MQ=255)
                                 ggaTTTTACTCGATTTATCGAATGAGGTGATTGAGCAACTTGATATTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAggg                                        >  2:44946/3‑90 (MQ=255)
                                 gAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAggg                                        >  2:283366/1‑90 (MQ=255)
                                  aaTTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGa                                       <  2:437665/90‑1 (MQ=255)
                                  aaTTTTACTCGATTTATCGAATGAGGTGATAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGaa                                      <  2:320068/90‑1 (MQ=255)
                                         cTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTa                                <  2:171217/90‑1 (MQ=255)
                                             aTTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGAc                                                  >  2:15747/1‑68 (MQ=255)
                                             aTTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGAc                                                  <  1:15747/68‑1 (MQ=255)
                                                   cGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATAcc                      <  2:132515/90‑1 (MQ=255)
                                                               ttAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAaatcaat          >  1:382769/1‑90 (MQ=255)
                                                                    cAACTTGATGTTCTCGAGGAGCAGCGTAATCTTCCGCGGGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGGTAATTCACTCGCaa     >  2:445226/1‑90 (MQ=255)
                                                                       cTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACa  >  1:70997/1‑90 (MQ=255)
                                                                              |                                                                                  
TTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGATCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACA  >  NZ_CP009273/4073085‑4073245

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ 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
AAATGATATTGGCGAGATATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGATCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACAGCAAGAACCAGTGTTCC  >  NZ_CP009273/4073065‑4073262
                                                                                                  |                                                                                                   
AAATGATATTGGCGAGATATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTT                                                                                                    >  SRR3722113.27889/1‑100 (MQ=60)
       ATTGGCGAGATATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGG                                                                                             <  SRR3722113.235768/100‑1 (MQ=60)
          tgggctcggagatgtgtataagagaCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGC                                                                                          <  SRR3722113.457839/75‑1 (MQ=60)
          GGCGAGATATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGC                                                                                          >  SRR3722113.133630/1‑100 (MQ=60)
                             CTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGC                                                                       >  SRR3722113.16262/1‑100 (MQ=60)
                                      GTATAGATATGGGCAGAATTTTACTCGGTTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATT                                                              >  SRR3722113.339505/1‑100 (MQ=60)
                                      GTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATT                                                              >  SRR3722113.15537/1‑100 (MQ=60)
                                      GTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATT                                                              >  SRR3722113.182603/1‑100 (MQ=60)
                                      GTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATT                                                              >  SRR3722113.207107/1‑100 (MQ=60)
                                           tataagagacagATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGG                                                         <  SRR3722113.15890/88‑1 (MQ=60)
                                            ATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGA                                                        <  SRR3722113.218440/100‑1 (MQ=60)
                                                      tagtctcgtgggctcggagatgtgtataagagacagACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGAC                                              <  SRR3722113.208164/64‑1 (MQ=60)
                                                                         AATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAAT                           >  SRR3722113.386769/1‑100 (MQ=60)
                                                                                 GATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACA                   >  SRR3722113.71616/1‑100 (MQ=60)
                                                                                                 GTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACAGCAAGAACCAGTGTTC   >  SRR3722113.261013/1‑100 (MQ=60)
                                                                                                  TTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACAGCAAGAACCAGTGTTCC  >  SRR3722113.134764/1‑100 (MQ=60)
                                                                                                  |                                                                                                   
AAATGATATTGGCGAGATATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGATCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACAGCAAGAACCAGTGTTCC  >  NZ_CP009273/4073065‑4073262

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