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 I225 R1 227 21.4 1179702 97.1% 1145490 86.5

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% 111.4 / NA 30D26V (GAT→GTT) BW25113_RS20185hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (17/13);  total (17/13)

GAGATATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGATCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAA  >  NZ_CP009273/4073078‑4073243
                                                                                     |                                                                                
gagaTATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTtctc                                                                              >  1:317526/1‑90 (MQ=255)
   atatTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAg                                                                           >  2:142583/1‑90 (MQ=255)
       ttATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGc                                                                       <  1:220907/90‑1 (MQ=255)
            aaCTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGt                                                                  <  1:201806/90‑1 (MQ=255)
               tctctGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTGTCGAATGAGGTGATTAAGCAACTTTATGTTCTCGAGGTGCAGCgcagt                                                               >  1:305114/1‑86 (MQ=255)
                ctctGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATc                                                              >  1:121121/1‑90 (MQ=255)
                       aGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCg                                                       >  2:452542/1‑90 (MQ=255)
                             aGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGa                                                 <  2:317526/90‑1 (MQ=255)
                               atatGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGAcc                                               <  1:430262/90‑1 (MQ=255)
                                  tGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGTTGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAt                                            >  2:440173/1‑90 (MQ=255)
                                  tGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAt                                            >  2:562197/1‑90 (MQ=255)
                                  tGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAt                                            <  2:468725/90‑1 (MQ=255)
                                  tGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAt                                            >  2:331125/1‑90 (MQ=255)
                                   gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                           >  2:569926/1‑90 (MQ=255)
                                   gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                           >  2:526296/1‑90 (MQ=255)
                                   gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                           >  1:114977/1‑90 (MQ=255)
                                   gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                           >  1:200003/1‑90 (MQ=255)
                                    ggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTa                                          <  1:544872/90‑1 (MQ=255)
                                        gAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGa                                                 >  2:94325/1‑79 (MQ=255)
                                        gAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGa                                                 <  1:94325/79‑1 (MQ=255)
                                          aTTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGaa                                    <  1:331125/90‑1 (MQ=255)
                                               aCTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGt                               <  1:578351/90‑1 (MQ=255)
                                                cTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTa                              <  2:470483/90‑1 (MQ=255)
                                                  cGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGa                            >  2:301868/1‑90 (MQ=255)
                                                                        aaGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCg      >  2:564036/1‑90 (MQ=255)
                                                                          gCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCa    >  1:20480/1‑90 (MQ=255)
                                                                            aaCTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCaaa  <  1:509321/90‑1 (MQ=255)
                                                                            aaCTTAATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCaaa  <  1:376434/90‑1 (MQ=255)
                                                                              cTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTa                              >  2:45874/1‑60 (MQ=255)
                                                                              cTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTa                              <  1:45874/60‑1 (MQ=255)
                                                                                     |                                                                                
GAGATATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGATCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAA  >  NZ_CP009273/4073078‑4073243

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ 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
TGATATTGGCGAGATATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGATCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACAGCAAGAAC  >  NZ_CP009273/4073068‑4073253
                                                                                               |                                                                                          
TGATATTGGCGAGATATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTC                                                                                        >  SRR3722109.321711/1‑100 (MQ=60)
     TTGGCGAGATATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGT                                                                                   <  SRR3722109.259695/100‑1 (MQ=60)
               ATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTGTCGAATGAGGTGATTAAGCAACTTTATGTTCTCGAGGTGCAGCGcagt                                                                         >  SRR3722109.309085/1‑96 (MQ=60)
                TTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATC                                                                        >  SRR3722109.122615/1‑100 (MQ=60)
                 TTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCT                                                                       <  SRR3722109.223547/100‑1 (MQ=60)
                      AACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTC                                                                  <  SRR3722109.204255/100‑1 (MQ=60)
                                   GTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATT                                                     >  SRR3722109.116400/1‑100 (MQ=60)
                                   GAATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATT                                                     >  SRR3722109.202431/1‑100 (MQ=60)
                                       GGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGG                                                 <  SRR3722109.95495/100‑1 (MQ=60)
                                         ATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGA                                               <  SRR3722109.436599/100‑1 (MQ=60)
                                              GGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGG                                          <  SRR3722109.553136/100‑1 (MQ=60)
                                                    ATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATC                                    <  SRR3722109.335505/100‑1 (MQ=60)
                                                         ACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATAC                               <  SRR3722109.587105/100‑1 (MQ=60)
                                                          cagatgtgtataagagacagGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACC                              <  SRR3722109.46425/80‑1 (MQ=60)
                                                                          AGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCA              >  SRR3722109.20749/1‑100 (MQ=60)
                                                                                      AACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACAGCAAGAAC  <  SRR3722109.517028/100‑1 (MQ=60)
                                                                                      AACTTAATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACAGCAAGAAC  <  SRR3722109.381762/100‑1 (MQ=60)
                                                                                               |                                                                                          
TGATATTGGCGAGATATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGATCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACAGCAAGAAC  >  NZ_CP009273/4073068‑4073253

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