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 1,026,705 A→C V57G (GTG→GGG)  yccA ← FtsH protease modulator YccA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,026,7050AC100.0% 118.1 / NA 32V57G (GTG→GGG) yccAFtsH protease modulator YccA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (18/14);  total (18/14)

CAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTC  >  NZ_CP009273/1026617‑1026791
                                                                                        |                                                                                      
cAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATAcccccc                                                                                      >  1:135045/1‑91 (MQ=255)
cAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATAcccccc                                                                                      >  1:736902/1‑91 (MQ=255)
     aaCCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGt                                                                                 >  2:885888/1‑91 (MQ=255)
       ccGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCa                                                                               <  2:3349/91‑1 (MQ=255)
                   aaTGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCg                                                                   >  1:586145/1‑91 (MQ=255)
                      gCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGgag                                                                >  1:938447/1‑91 (MQ=255)
                            aTAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGgca                                                          <  2:206614/91‑1 (MQ=255)
                            aTAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGgca                                                          <  2:583218/91‑1 (MQ=255)
                               aTCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGgcagca                                                       >  2:931481/1‑91 (MQ=255)
                               aTCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGgcagca                                                       >  1:245338/1‑91 (MQ=255)
                                    ggTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGc                                                  >  1:458893/1‑91 (MQ=255)
                                     gTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCa                                                 >  2:23523/1‑91 (MQ=255)
                                      tCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCAc                                                <  2:586145/91‑1 (MQ=255)
                                      tCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCAc                                                <  2:461966/91‑1 (MQ=255)
                                           ttATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTAc                                           <  1:135555/91‑1 (MQ=255)
                                                  ccGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGt                                    <  1:475231/91‑1 (MQ=255)
                                                         ataGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTa                             <  1:460807/91‑1 (MQ=255)
                                                                 ggAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGgca                                                          >  1:310997/1‑54 (MQ=255)
                                                                 ggAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGgca                                                          <  2:310997/54‑1 (MQ=255)
                                                                    aCATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGaaaa                  >  2:145825/1‑91 (MQ=255)
                                                                      aTCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATgg                                                            <  2:945930/47‑1 (MQ=255)
                                                                      aTCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATgg                                                            >  1:945930/1‑47 (MQ=255)
                                                                      aTCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAgg                >  1:606946/1‑91 (MQ=255)
                                                                      aTCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAgg                >  1:1071562/1‑91 (MQ=255)
                                                                       tCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGc               >  1:393203/1‑91 (MQ=255)
                                                                        cAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGcc              <  2:413839/91‑1 (MQ=255)
                                                                               tacatacCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCa       <  2:155252/91‑1 (MQ=255)
                                                                               tacatacCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCa       >  2:142388/1‑91 (MQ=255)
                                                                               tacatacCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCa       >  1:477072/1‑91 (MQ=255)
                                                                                   tacCCCCCAGCGTCAGAATCAGACCTGGAGATGGCAGCATCAGCAc                                                >  2:337860/1‑46 (MQ=255)
                                                                                   tacCCCCCAGCGTCAGAATCAGACCTGGAGATGGCAGCATCAGCAc                                                <  1:337860/46‑1 (MQ=255)
                                                                                   tacCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCt   >  2:632825/1‑91 (MQ=255)
                                                                                   tacCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCt   >  1:289271/1‑91 (MQ=255)
                                                                                    acCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTc  <  2:785212/91‑1 (MQ=255)
                                                                                        |                                                                                      
CAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTC  >  NZ_CP009273/1026617‑1026791

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
GGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAA  >  NZ_CP009273/1026607‑1026800
                                                                                                  |                                                                                               
GGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCC                                                                                               <  SRR3722125.1021380/101‑1 (MQ=60)
GGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCC                                                                                               >  SRR3722125.137610/1‑101 (MQ=60)
GGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCC                                                                                               <  SRR3722125.158351/101‑1 (MQ=60)
GGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCC                                                                                               >  SRR3722125.752796/1‑101 (MQ=60)
GGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCC                                                                                               <  SRR3722125.563504/101‑1 (MQ=60)
                   GGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCG                                                                            >  SRR3722125.598736/1‑101 (MQ=60)
                      AAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAG                                                                         >  SRR3722125.959285/1‑101 (MQ=60)
                               TGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCA                                                                >  SRR3722125.250267/1‑101 (MQ=60)
                                    AGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGC                                                           >  SRR3722125.468577/1‑101 (MQ=60)
                                                agtctcgtgggctcggagatgtgtataagagacagAAACCATACATACCCCCCAGCGTCAGAATCAGACCTGGAGATGGCAGCATCAGCACCGTACTGGCA                                               <  SRR3722125.344887/66‑1 (MQ=60)
                                                     TTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGC                                          <  SRR3722125.138132/101‑1 (MQ=60)
                                                            CCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATC                                   <  SRR3722125.485250/101‑1 (MQ=60)
                                                                 TTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCtgtctcttatacacatctccgagccc                              >  SRR3722125.317394/1‑75 (MQ=60)
                                                                   ATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAA                            <  SRR3722125.470521/101‑1 (MQ=60)
                                                                      GGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCctgtctcttatacacatctccgagcccacgagac                         >  SRR3722125.966964/1‑67 (MQ=60)
                                                                      GGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGG                         >  SRR3722125.620016/1‑101 (MQ=60)
                                                                      GGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGG                         >  SRR3722125.1095703/1‑101 (MQ=60)
                                                                       GTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGC                        >  SRR3722125.401438/1‑101 (MQ=60)
                                                                               CATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCA                >  SRR3722125.487138/1‑101 (MQ=60)
                                                                                   AAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCT            >  SRR3722125.295169/1‑101 (MQ=60)
                                                                                      CCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAG         >  SRR3722125.762398/1‑101 (MQ=60)
                                                                                        ATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCA       >  SRR3722125.26798/1‑101 (MQ=60)
                                                                                         TACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCctgtctcttatacacat      >  SRR3722125.840474/1‑84 (MQ=60)
                                                                                         TACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAG      >  SRR3722125.613294/1‑101 (MQ=60)
                                                                                          ACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGA     >  SRR3722125.85841/1‑101 (MQ=60)
                                                                                             TACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAA  >  SRR3722125.860017/1‑101 (MQ=60)
                                                                                                  |                                                                                               
GGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAA  >  NZ_CP009273/1026607‑1026800

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