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 I194 R1 197 32.2 1845014 95.2% 1756453 85.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 575,614 A→G intergenic (+72/+93) BW25113_RS02885 → / ← ylcI hypothetical protein/YlcI/YnfO family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273575,6140AG100.0% 133.6 / NA 37intergenic (+72/+93)BW25113_RS02885/ylcIhypothetical protein/YlcI/YnfO family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/15);  total (22/15)

CTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAAC  >  NZ_CP009273/575528‑575697
                                                                                      |                                                                                   
cTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAg                                                                                  >  2:265892/1‑90 (MQ=255)
            tGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGc                                                                      >  2:270095/1‑90 (MQ=255)
             gATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCa                                                                     <  2:829939/90‑1 (MQ=255)
             gATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCa                                                                     <  2:357866/90‑1 (MQ=255)
                cACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACt                                                                  <  2:514684/90‑1 (MQ=255)
                 aCCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTc                                                                 >  1:57382/1‑90 (MQ=255)
                       aCTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTAc                                                           >  1:853160/1‑90 (MQ=255)
                           cTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCa                                                       >  1:833385/1‑90 (MQ=255)
                            tAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCaa                                                      <  1:464233/90‑1 (MQ=255)
                            tAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCaa                                                      >  1:196894/1‑90 (MQ=255)
                            tAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCaa                                                      >  1:672197/1‑90 (MQ=255)
                               tCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGc                                                   <  2:824621/90‑1 (MQ=255)
                                  ttCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTAc                                                           <  2:553727/79‑1 (MQ=255)
                                  ttCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTAc                                                           >  1:553727/1‑79 (MQ=255)
                                  ttCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCt                                                >  1:738708/1‑90 (MQ=255)
                                  ttCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCt                                                >  1:159944/1‑90 (MQ=255)
                                  ttCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCt                                                <  1:39994/90‑1 (MQ=255)
                                          actTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTa                                        >  2:294227/1‑90 (MQ=255)
                                           ctTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTaa                                       >  2:393726/1‑90 (MQ=255)
                                              aaCCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTg                                    <  1:44150/90‑1 (MQ=255)
                                              aaCCTGTGACAGAGCCAACACGCAGTCAGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTg                                    <  2:414841/90‑1 (MQ=255)
                                                  tgtgACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAAtt                                >  1:709514/1‑90 (MQ=255)
                                                    tgACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGt                                     >  2:822680/1‑83 (MQ=255)
                                                    tgACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGt                                     <  1:822680/83‑1 (MQ=255)
                                                    tgACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTa                              >  1:850263/1‑90 (MQ=255)
                                                    tgACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTa                              >  1:76476/1‑90 (MQ=255)
                                                    tgACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTa                              >  2:526338/1‑90 (MQ=255)
                                                     gACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTAc                             <  1:393726/90‑1 (MQ=255)
                                                     gACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTAc                             <  1:199586/90‑1 (MQ=255)
                                                         gagCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGt                                     <  2:182167/78‑1 (MQ=255)
                                                         gagCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGt                                     >  1:182167/1‑78 (MQ=255)
                                                         gagCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAata                         >  1:673027/1‑90 (MQ=255)
                                                             cAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGt                     >  1:417679/1‑90 (MQ=255)
                                                               acacGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTcc                   <  2:57382/90‑1 (MQ=255)
                                                                       tCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGa           >  1:861369/1‑90 (MQ=255)
                                                                        ctgtcactgtcaGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAg          <  1:714058/90‑1 (MQ=255)
                                                                                gtcaGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAAc  >  1:220012/1‑90 (MQ=255)
                                                                                      |                                                                                   
CTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAAC  >  NZ_CP009273/575528‑575697

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGG  >  NZ_CP009273/575518‑575701
                                                                                                |                                                                                       
GGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAG                                                                                      <  SRR3722072.638251/100‑1 (MQ=60)
                 TTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTC                                                                     >  SRR3722072.58429/1‑100 (MQ=60)
                       GATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTAC                                                               >  SRR3722072.867894/1‑100 (MQ=60)
                           ACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCA                                                           >  SRR3722072.847830/1‑100 (MQ=60)
                            CCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAA                                                          >  SRR3722072.200289/1‑100 (MQ=60)
                            CCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAA                                                          >  SRR3722072.684235/1‑100 (MQ=60)
                                  CTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCT                                                    >  SRR3722072.162730/1‑100 (MQ=60)
                                  CTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCT                                                    >  SRR3722072.751797/1‑100 (MQ=60)
                                  CTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCC                                                    >  SRR3722072.563677/1‑100 (MQ=60)
                                      TAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTA                                                <  SRR3722072.472723/100‑1 (MQ=60)
                                            TTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAG                                          <  SRR3722072.40724/100‑1 (MQ=60)
                                                  CTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATT                                    >  SRR3722072.722169/1‑100 (MQ=60)
                                                    ACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTA                                  >  SRR3722072.864957/1‑100 (MQ=60)
                                                    ACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTA                                  >  SRR3722072.77896/1‑100 (MQ=60)
                                                       TAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAA                               <  SRR3722072.836990/100‑1 (MQ=60)
                                                        AACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAAT                              <  SRR3722072.44966/100‑1 (MQ=60)
                                                         ACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAAct                             >  SRR3722072.185319/1‑98 (MQ=60)
                                                         ACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATA                             >  SRR3722072.896704/1‑100 (MQ=60)
                                                         ACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATA                             >  SRR3722072.685081/1‑100 (MQ=60)
                                                             GTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGT                         >  SRR3722072.425247/1‑100 (MQ=60)
                                                               GACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCC                       <  SRR3722072.203030/100‑1 (MQ=60)
                                                               GACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCC                       <  SRR3722072.400809/100‑1 (MQ=60)
                                                                       CAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGA               >  SRR3722072.876231/1‑100 (MQ=60)
                                                                                GTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAAC      >  SRR3722072.223753/1‑100 (MQ=60)
                                                                                  CTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGC    <  SRR3722072.726781/100‑1 (MQ=60)
                                                                                    GTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATctgtctcttatacacatctccgagccca  >  SRR3722072.107067/1‑72 (MQ=60)
                                                                                                |                                                                                       
GGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGG  >  NZ_CP009273/575518‑575701

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