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 I210 R1 226 19.6 1130256 95.4% 1078264 84.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,410,868 A→T G103G (GGA→GGT yfcC → putative basic amino acid antiporter YfcC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,410,8680AT100.0% 110.3 / NA 30G103G (GGA→GGTyfcCputative basic amino acid antiporter YfcC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (18/12);  total (18/12)

ACCTGAGTATCACCGCGTACAGCTGTTCACGACGGGCGATGAACGCCCGGGCCTGATGAACTTCCCGTTTGAAGGATTAACCTCAGGATCGAAATACGGGACAGCCGTTGGCATCATCATGTTTATGCTGGTGATTGGCGGCGCGTTTGGCATTGTGATGCGTACAGGAACCA  >  NZ_CP009273/2410781‑2410953
                                                                                       |                                                                                     
aCCTGAGTATCACCGCGTACAGCTGTTCACGACGGGCGATGAACGCCCGGGCCTGATGAACTTCCCGTTTGAAGGATTAACCTCAGGTTc                                                                                     <  2:47296/90‑1 (MQ=255)
aCCTGAGTATCACCGCGTACAGCTGTTCACGACGGGCGATGAACGCCCGGGCCTGATGAACTTCCCGTTTGAAGGATTAACCTCAGGTTc                                                                                     >  2:378966/1‑90 (MQ=255)
           aCCGCGTACAGCTGTTCACGACGGGCGATGAACGCCCGGGCCTGATGAACTTCCCGTTTGAAGGATTAACCTCAGGTTCGa                                                                                   <  1:320686/81‑1 (MQ=255)
           aCCGCGTACAGCTGTTCACGACGGGCGATGAACGCCCGGGCCTGATGAACTTCCCGTTTGAAGGATTAACCTCAGGTTCGa                                                                                   >  2:320686/1‑81 (MQ=255)
              gcgTACAGCTGTTCACGACGGGCGATGAACGCCCGGGCCTGATGAACTTCCCGTTTGAAGGATTAGCCTCAGGTTCGAAATACGGGACAg                                                                       >  1:464502/1‑90 (MQ=255)
                gTACAGCTGTTCACGACGGGCGATGAACGCCCGGGCCTGATGAACTTCCCGTTTGAAGGATTAACCTCAGGTTCGAAATACGGGACAGcc                                                                     >  2:507881/1‑90 (MQ=255)
                        gTTCACGACGGGCGATGAACGCCCGGGCCTGATGAACTTCCCGTTTGAAGGATTAACCTCAGGTTCGAAATACGGGACAGCCGTTGGcat                                                             >  2:417348/1‑90 (MQ=255)
                        gTTCACGACGGGCGATGAACGCCCGGGCCTGATGAACTTCCCGTTTGAAGGATTAACCTCAGGTTCGAAATACGGGACAGCCGTTGGcat                                                             >  1:11689/1‑90 (MQ=255)
                         ttCACGACGGGCGATGAACGCCCGGGCCTGATGAACTTCTCGTTTGAAGGATTAACCTCAGGTTCGAAATACGGGCCAGCCCTTGGcatc                                                            >  1:333035/1‑90 (MQ=255)
                            acgacgGGCGATGAACGCCCGGGCCTGATGAACTTCCCGTTTGAAGGATTAACCTCAGGTTCGAAATACGGGACAGCCGTTGGcatcatc                                                         <  1:43670/90‑1 (MQ=255)
                            acgacgGGCGATGAACGCCCGGGCCTGATGAACTTCCCGTTTGAAGGATTAACCTCAGGTTCGAAATACGGGACAGCCGTTGGcatcatc                                                         <  2:246438/90‑1 (MQ=255)
                             cgacgGGCGATGAACGCCCGGGCCTGATGAACTTCCCGTTTGAAGGATTAACCTCAGGTTCGAAATACGGGACAGCCGTTGGcatcatca                                                        >  1:410136/1‑90 (MQ=255)
                               acgGGCGATGAACGCCCGGGCCTGATGAACTTCCCGTTTGAAGGATTAACCTCAGGTTCGAAATACGGGACAGCCGTTGGCATCATCATg                                                      >  2:208214/1‑90 (MQ=255)
                                  ggCGATGAACGCCCGGGCCTGATGAACTTCCCGTTTGAAGGATTAACCTCAGGTTCGAAATACGGGACAGCCGTTGGCATCATCATGttt                                                   >  2:240398/1‑90 (MQ=255)
                                  ggCGATGAACGCCCGGGCCTGATGAACTTCCCGTTTGAAGGATTAACCTCAGGTTCGAAATACGGGACAGCCGTTGGCATCATCATGttt                                                   >  1:428239/1‑90 (MQ=255)
                                                      gatgaACTTCCCGTTTGAAGGATTAACCTCAGGTTCGAAATACGGGAc                                                                         >  1:437640/1‑48 (MQ=255)
                                                      gatgaACTTCCCGTTTGAAGGATTAACCTCAGGTTCGAAATACGGGAc                                                                         <  2:437640/48‑1 (MQ=255)
                                                      gatgaACTTCCCGTTTGAAGGATTAACCTCAGGTTCGAAATACGGGACAGCCGTTGGcatcatc                                                         <  1:4312/64‑1 (MQ=255)
                                                      gatgaACTTCCCGTTTGAAGGATTAACCTCAGGTTCGAAATACGGGACAGCCGTTGGcatcatc                                                         >  2:4312/1‑64 (MQ=255)
                                                      gatgaACTTCCCGTTTGAAGGATTAACCTCAGGTTCGAAATACGGGACAGCCGTTGGCATCATCATGTTTATGCTGGTGATTGgcggcg                                <  2:386536/89‑1 (MQ=255)
                                                      gatgaACTTCCCGTTTGAAGGATTAACCTCAGGTTCGAAATACGGGACAGCCGTTGGCATCATCATGTTTATGCTGGTGATTGgcggcg                                >  1:386536/1‑89 (MQ=255)
                                                      gatgaACTTCCCGTTTGAAGGATTAACCTCAGGTTCGAAATACGGGACAGCCGTTGGCATCATCATGTTTATGCTGGTGATTGGCGgcgc                               >  2:217325/1‑90 (MQ=255)
                                                      gatgaACTTCCCGTTTGAAGGATTAACCTCAGGTTCGAAATACGGGACAGCCGTTGGCATCATCATGTTTATGCTGGTGATTGGCGgcgc                               >  1:91960/1‑90 (MQ=255)
                                                        tgaACTTCCCGTTTGAAGGATTAACCTCAGGTTCGAAATACGGGACAGCCGTTGGCATCATCATGTTTATGCTGGTGATTGGCGGCGCGt                             <  2:91960/90‑1 (MQ=255)
                                                          aaCTTCCCGTTTGAAGGATTAACCTCAGGTTCGAAATACGGGACAGCCGTTGGCATCATCATGTTTATGCTGGTGATTGGCGGCGCGttt                           <  2:11689/90‑1 (MQ=255)
                                                              tCCCGTTTGAAGGATTAACCTCAGGTTCGAAATACGGGACAGCCGTTGGCATCATCATGTTTATGCTGGTGATTGGCGGCGCGTTTGGCa                       >  2:540939/1‑90 (MQ=255)
                                                                            ttAACCTCAGGTTCGAAATACGGGACAGCCGTTGGCATCATCATGTTTATGCTGGTGATTGGCGGCGCGTTTGGCATTGTGATGCGTACa         >  2:253352/1‑90 (MQ=255)
                                                                                 cTCAGGTTCGAAATACGGGACAGCCGTTGGCATCATCATGTTTATGCTGGTGATTGGCGGCGCGTTTGGCATTGTGATGCGTACAGGAAc    <  2:342395/90‑1 (MQ=255)
                                                                                 cTCAGGTTCGAAATACGGGACAGCCGTTGGCATCATCATGTTTATGCTGGTGATTGGCGGCGCGTTTGGCATTGTGATGCGTACAGGAAc    <  2:443262/90‑1 (MQ=255)
                                                                                   cAGGTTCGAAATACGGGACAGCCGTTGGCATCATCATGTTTATGCTGGTGATTGGCGGCGCGTTTGGCATTGTGATGCGTACAGGAACCa  <  1:507881/90‑1 (MQ=255)
                                                                                       |                                                                                     
ACCTGAGTATCACCGCGTACAGCTGTTCACGACGGGCGATGAACGCCCGGGCCTGATGAACTTCCCGTTTGAAGGATTAACCTCAGGATCGAAATACGGGACAGCCGTTGGCATCATCATGTTTATGCTGGTGATTGGCGGCGCGTTTGGCATTGTGATGCGTACAGGAACCA  >  NZ_CP009273/2410781‑2410953

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
CTGAGTATCACCGCGTACAGCTGTTCACGACGGGCGATGAACGCCCGGGCCTGATGAACTTCCCGTTTGAAGGATTAACCTCAGGATCGAAATACGGGACAGCCGTTGGCATCATCATGTTTATGCTGGTGATTGGCGGCGCGTTTGGCATTGTGATGCGTACAGGAACCATTGATAACGG  >  NZ_CP009273/2410783‑2410963
                                                                                     |                                                                                               
CTGAGTATCACCGCGTACAGCTGTTCACGACGGGCGATGAACGCCCGGGCCTGATGAACTTCCCGTTTGAAGGATTAACCTCAGGTTCGAAATACGGGAC                                                                                   <  SRR3722092.326354/100‑1 (MQ=60)
  GAGTATCACCGCGTACAGCTGTTCACGACGGGCGATGAACGCCCGGGCCTGATGAACTTCCCGTTTGAAGGATTAGCCTCAGGTTCGAAATACGGGACAG                                                                                 >  SRR3722092.473701/1‑100 (MQ=60)
            GCGTACAGCTGTTCACGACGGGCGATGAACGCCCGGGCCTGATGAACTTCCCGTTTGAAGGATTAACCTCAGGTTCGAAATACGGGACAGCCGTTGGCAT                                                                       >  SRR3722092.11898/1‑100 (MQ=60)
             CGTACAGCTGTTCACGACGGGCGATGAACGCCCGGGCCTGATGAACTTCTCGTTTGAAGGATTAACCTCAGGTTCGAAATACGGGCCAGCCCTTGGCATC                                                                      >  SRR3722092.338922/1‑100 (MQ=60)
                 CAGCTGTTCACGACGGGCGATGAACGCCCGGGCCTGATGAACTTCCCGTTTGAAGGATTAACCTCAGGTTCGAAATACGGGACAGCCGTTGGCATCATCA                                                                  >  SRR3722092.418046/1‑100 (MQ=60)
                      GTTCACGACGGGCGATGAACGCCCGGGCCTGATGAACTTCCCGTTTGAAGGATTAACCTCAGGTTCGAAATACGGGACAGCCGTTGGCATCATCATGTTT                                                             >  SRR3722092.436549/1‑100 (MQ=60)
                          tgtgtataagagacagGCCCGGGCCTGATGAACTTCCCGTTTGAAGGATTAACCTCAGGTTCGAAATACGGGACAGCCGTTGGCATCATCATGTTTATGC                                                         <  SRR3722092.4394/84‑1 (MQ=60)
                          ACGACGGGCGATGAACGCCCGGGCCTGATGAACTTCCCGTTTGAAGGATTAACCTCAGGTTCGAAATACGGGACAGCCGTTGGCATCATCATGTTTATGC                                                         <  SRR3722092.44410/100‑1 (MQ=60)
                                          GCCCGGGCCTGATGAACTTCCCGTTTGAAGGATTAACCTCAGGTTCGAAATACGGGACAGCCGTTGGCctgtctcttatacacatctgacgctgccgacg                                         >  SRR3722092.446208/1‑68 (MQ=60)
                                          GCCCGGGCCTGATGAACTTCCCGTTTGAAGGATTAACCTCAGGTTCGAAATACGGGACAGCCGTTGGCATCATCATGTTTATGCTGGTGATTGGCGGCGC                                         >  SRR3722092.393897/1‑100 (MQ=60)
                                          GCCCGGGCCTGATGAACTTCCCGTTTGAAGGATTAACCTCAGGTTCGAAATACGGGACAGCCGTTGGCATCATCATGTTTATGCTGGTGATTGGCGGCGC                                         >  SRR3722092.93532/1‑100 (MQ=60)
                                                                                 CAGGTTCGAAATACGGGACAGCCGTTGGCATCATCATGTTTATGCTGGTGATTGGCGGCGCGTTTGGCATTGTGATGCGTACAGGAACCATTGATAACGG  <  SRR3722092.518043/100‑1 (MQ=60)
                                                                                     |                                                                                               
CTGAGTATCACCGCGTACAGCTGTTCACGACGGGCGATGAACGCCCGGGCCTGATGAACTTCCCGTTTGAAGGATTAACCTCAGGATCGAAATACGGGACAGCCGTTGGCATCATCATGTTTATGCTGGTGATTGGCGGCGCGTTTGGCATTGTGATGCGTACAGGAACCATTGATAACGG  >  NZ_CP009273/2410783‑2410963

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: