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 F5 I102 R1 278 39.1 2251254 95.8% 2156701 87.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 871,743 A→G E318G (GAA→GGA)  dgcI → GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273871,7430AG100.0% 100.2 / NA 28E318G (GAA→GGA) dgcIGGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/9);  total (19/9)

ATGTATAACTCGGCGTTCCTGGACGACCTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGAAAAAGGCGGTTATTTGTGCCTGTTTGACGTCGATAAATTCAAAAACATTAACGACACCTTCGGTCATTTGCTGGGCGATGAAGTGT  >  NZ_CP009273/871661‑871829
                                                                                  |                                                                                      
atgtatAACTCGGCGTTCCTGGACGACCTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCg                                                                                >  2:362079/1‑91 (MQ=255)
   tataACTCGGCGTTCCTGGACGACCTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCGGtt                                                                             >  1:224019/1‑91 (MQ=255)
     taACTCGGCGTTCCTGGACGACCTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCGGTTAt                                                                           >  1:586360/1‑91 (MQ=255)
       aCTCGGCGTTCCTGGACGACCTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCGGTTAttt                                                                         <  1:274880/91‑1 (MQ=255)
       aCTCGGCGTTCCTGGACGACCTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCGGTTAttt                                                                         >  1:717511/1‑91 (MQ=255)
       aCTCGGCGTTCCTGGACGACCTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCGGTTAttt                                                                         >  1:627834/1‑91 (MQ=255)
        cTCGGCGTTCCTGGACGACCTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCGGTTATTtg                                                                        >  1:1007784/1‑91 (MQ=255)
        cTCGGCGTTCCTGGACGACCTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCGGTTATTtg                                                                        >  2:580574/1‑91 (MQ=255)
          cGGCGTTCCTGGACGACCTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCGGTTATTtgtg                                                                      >  1:377448/1‑91 (MQ=255)
          cGGCGTTCCTGGACGACCTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCGGTTATTtgtg                                                                      >  1:713042/1‑91 (MQ=255)
          cGGCGTTCCTGGACGACCTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCGGTTATTtgtg                                                                      >  1:371160/1‑91 (MQ=255)
          cGGCGTTCCTGGACGACCTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCGGTTATTtgtg                                                                      >  2:767774/1‑91 (MQ=255)
          cGGCGTTCCTGGACGACCTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCGGTTATTtgtg                                                                      >  2:20226/1‑91 (MQ=255)
            gCGTTCCTGGACGACCTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCGGTTATTTGTGcc                                                                    >  2:287852/1‑91 (MQ=255)
               ttCCTGGACGACCTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCGGTTATTTGTGCCTGt                                                                 >  2:395376/1‑91 (MQ=255)
                         acCTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCGGTTATTTGTGCCTGTTTGACGTCGa                                                       >  2:681370/1‑91 (MQ=255)
                         acCTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCGGTTATTTGTGCCTGTTTGACGTCGa                                                       >  1:1067372/1‑91 (MQ=255)
                           cTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCGGTTATTTGTGCCTGTTTGACGTCGATa                                                     <  2:503824/91‑1 (MQ=255)
                            tGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCGGTTATTTGTGCCTGTTTGACGTCGATaa                                                    <  1:92332/91‑1 (MQ=255)
                                     cAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCGGTTATTTGTGCCTGTTTGACGTCGATAAATTCaaaaa                                           <  1:312483/91‑1 (MQ=255)
                                          ggCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCGGTTATTTGTGCCTGTTTGACGTCGATAAATTCAAAAACATTa                                      <  2:422587/91‑1 (MQ=255)
                                                     cTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCGGTTATTTGTGCCTGTTTGACGTCGATAAATTCAAAAACATTAACGACACCTTc                           <  2:100105/91‑1 (MQ=255)
                                                             aaGATTTAAAAGCGCTGCCGGGAAAAGGCGGTTATTTGTGCCTGTTTGACGTCGATAAATTCAAAAACATTAACGACACCTTCGGTCAttt                   >  2:881486/1‑91 (MQ=255)
                                                                         cgcTGCCGGGAAAAGGCGGTTATTTGTGCCTGTTTGACGTCGATAAATTCAAAAACATTAACGACACCTTCGGTCATTTGCTGGGCGATGa       <  1:654070/91‑1 (MQ=255)
                                                                         cgcTGCCGGGAAAAGGCGGTTATTTGTGCCTGTTTGACGTCGATAAATTCAAAAACATTAACGACACCTTCGGTCATTTGCTGGGCGATGa       >  2:699060/1‑91 (MQ=255)
                                                                           cTGCCGGGAAAAGGCGGTTATTTGTGCCTGTTTGACGTCGATAAATTCAAAAACATTAACGACACCTTCGGTCATTTGCTGGGCGATGAAg     >  2:220631/1‑91 (MQ=255)
                                                                           cTGCCGGGAAAAGGCGGTTATTTGTGCCTGTTTGACGTCGATAAATTCAAAAACATTAACGACACCTTCGGTCATTTGCTGGGCGATGAAg     <  2:742033/91‑1 (MQ=255)
                                                                              ccGGGAAAAGGCGGTTATTTGTGCCTGTTTGACGTCGATAAATTCAAAAACATTAACGACACCTTCGGTCATTTGCTGGGCGATGAAgtgt  <  2:205667/91‑1 (MQ=255)
                                                                                  |                                                                                      
ATGTATAACTCGGCGTTCCTGGACGACCTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGAAAAAGGCGGTTATTTGTGCCTGTTTGACGTCGATAAATTCAAAAACATTAACGACACCTTCGGTCATTTGCTGGGCGATGAAGTGT  >  NZ_CP009273/871661‑871829

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ 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
AAAAAACTCTACAGTATGTATAACTCGGCGTTCCTGGACGACCTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGAAAAAGGCGGTTATTTGTGCCTGTTTGACGTCGATAAATTCAAAAACATTAACGACACCTTCGGTCATTTGCTGGGCGATGAAGTGTTGATGAAAG  >  NZ_CP009273/871646‑871838
                                                                                                 |                                                                                               
AAAAAACTCTACAGTATGTATAACTCGGCGTTCCTGGACGACCTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAA                                                                                              <  SRR3722125.207142/101‑1 (MQ=60)
   NAACTCTACAGTATGTATAACTCGGCGTTCCTGGACGACCTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGG                                                                                           >  SRR3722125.657139/1‑101 (MQ=60)
        CTACAGTATGTATAACTCGGCGTTCCTGGACGACCTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCGGTT                                                                                      >  SRR3722125.228461/1‑101 (MQ=60)
          ACAGTATGTATAACTCGGCGTTCCTGGACGACCTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCGGTTAT                                                                                    >  SRR3722125.598953/1‑101 (MQ=60)
            AGTATGTATAACTCGGCGTTCCTGGACGACCTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCGGTTATTT                                                                                  >  SRR3722125.733001/1‑101 (MQ=60)
            AGTATGTATAACTCGGCGTTCCTGGACGACCTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCGGTTATTT                                                                                  >  SRR3722125.641363/1‑101 (MQ=60)
             ttttttaaatCTCGGCGTTCCTGGACGACCTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCGGTTATTTG                                                                                 <  SRR3722125.700373/91‑1 (MQ=60)
             GTATGTATAACTCGGCGTTCCTGGACGACCTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCGGTTATTTG                                                                                 >  SRR3722125.1030318/1‑101 (MQ=60)
               ATGTATAACTCGGCGTTCCTGGACGACCTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCGGTTATTTGTG                                                                               >  SRR3722125.378935/1‑101 (MQ=60)
               ATGTATAACTCGGCGTTCCTGGACGACCTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCGGTTATTTGTG                                                                               >  SRR3722125.385359/1‑101 (MQ=60)
               ATGTATAACTCGGCGTTCCTGGACGACCTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCGGTTATTTGTG                                                                               >  SRR3722125.728417/1‑101 (MQ=60)
                      ACTCGGCGTTCCTGGACGACCTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCGGTTATTTGTGCCTGTTT                                                                        <  SRR3722125.280437/101‑1 (MQ=60)
                              TTCCTGGACGACCTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCGGTTATTTGTGCCTGTTTGACGTCGA                                                                >  SRR3722125.1091415/1‑101 (MQ=60)
                                           TGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCGGTTATTTGTGCCTGTTTGACGTCGATAAATTCAAAAAC                                                   <  SRR3722125.94074/101‑1 (MQ=60)
                                                    CAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGGAAAAGGCGGTTATTTGTGCCTGTTTGACGTCGATAAATTCAAAAACATTAACGAC                                          <  SRR3722125.318913/101‑1 (MQ=60)
                                                                                        CGCTGCCGGGAAAAGGCGGTTATTTGTGCCTGTTTGACGTCGATAAATTCAAAAACATTAACGACACCTTCGGTCATTTGCTGGGCGATGAAGTGTTGATG      <  SRR3722125.668170/101‑1 (MQ=60)
                                                                                            GCCGGGAAAAGGCGGTTATTTGTGCCTGTTTGACGTCGATAAATTCAAAAACATTAACGACACCTTCGGTCATTTGCTGGGCGATGAAGTGTTGATGAAAG  >  SRR3722125.413095/1‑101 (MQ=60)
                                                                                            GCCGGGAAAAGGCGGTTATTTGTGCCTGTTTGACGTCGATAAATTCAAAAACATTAACGACACCTTCGGTCATTTGCTGGGCGATGAAGTGTTGATGAAAG  >  SRR3722125.785526/1‑101 (MQ=60)
                                                                                                 |                                                                                               
AAAAAACTCTACAGTATGTATAACTCGGCGTTCCTGGACGACCTGACCAAAGCAATGGGCAGAAAATCCTTTGACGAAGATTTAAAAGCGCTGCCGGAAAAAGGCGGTTATTTGTGCCTGTTTGACGTCGATAAATTCAAAAACATTAACGACACCTTCGGTCATTTGCTGGGCGATGAAGTGTTGATGAAAG  >  NZ_CP009273/871646‑871838

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