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 3,179,271 G→A A167T (GCC→ACC)  ygiL → fimbrial‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,179,2710GA100.0% 63.4 / NA 19A167T (GCC→ACC) ygiLfimbrial‑like protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (10/9);  total (10/9)

CTCCGGTAGTAATTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATGCCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCC  >  NZ_CP009273/3179188‑3179350
                                                                                   |                                                                               
cTCCGGTAGTAATTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCAccc                                                                           >  2:543026/1‑90 (MQ=255)
cTCCGGTAGTAATTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCAccc                                                                           >  2:539720/1‑90 (MQ=255)
cTCCGGTAGTAATTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCAccc                                                                           >  2:511321/1‑90 (MQ=255)
            ttACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGttt                                                               <  2:275734/90‑1 (MQ=255)
            ttACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGttt                                                               <  1:539720/90‑1 (MQ=255)
                gTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGc                                                           <  2:241266/90‑1 (MQ=255)
                   caacaaTACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCa                                                        <  1:33990/90‑1 (MQ=255)
                             aaCAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGc                                                     >  2:199211/1‑83 (MQ=255)
                             aaCAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGc                                                     <  1:199211/83‑1 (MQ=255)
                             aaCAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTAc                                              >  1:414112/1‑90 (MQ=255)
                                   tATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGAtt                                        >  2:12907/1‑90 (MQ=255)
                                      cAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGc                                                     <  2:391605/74‑1 (MQ=255)
                                      cAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGc                                                     >  1:391605/1‑74 (MQ=255)
                                           aCTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGa                                >  2:269140/1‑90 (MQ=255)
                                                 tttCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATc                          <  1:558652/90‑1 (MQ=255)
                                                               tGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGt            <  2:530041/90‑1 (MQ=255)
                                                                   aTCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGc        >  1:488252/1‑90 (MQ=255)
                                                                       ccAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTGATCGCATTTGGCCGGGTTCGCCCgg    >  2:329277/1‑90 (MQ=255)
                                                                         aaGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGcc  <  2:421238/90‑1 (MQ=255)
                                                                                   |                                                                               
CTCCGGTAGTAATTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATGCCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCC  >  NZ_CP009273/3179188‑3179350

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ 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
GAAGTAACTTTAGGTACTCCGGTAGTAATTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATGCCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAA  >  NZ_CP009273/3179172‑3179368
                                                                                                   |                                                                                                 
GAAGTAACTTTAGGTACTCCGGTAGTAATTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATA                                                                                                   >  SRR3722092.209954/1‑100 (MQ=60)
           gcgtcagatgtgtataagagacagCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGC                                                                                        <  SRR3722092.210383/76‑1 (MQ=60)
                            TTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGC                                                                       <  SRR3722092.550560/100‑1 (MQ=60)
                                   CAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTAC                                                                <  SRR3722092.34575/100‑1 (MQ=60)
                                   CAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTAC                                                                >  SRR3722092.422102/1‑100 (MQ=60)
                                      CAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAG                                                             <  SRR3722092.202502/100‑1 (MQ=60)
                                            CAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGctgtct                                                       >  SRR3722092.399087/1‑94 (MQ=60)
                                                                 TTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCC                                  <  SRR3722092.569862/100‑1 (MQ=60)
                                                                         CTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGC                          >  SRR3722092.497957/1‑100 (MQ=60)
                                                                                           GCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGctgtctcttatacacatctgacgctgccgacgaccatctta        >  SRR3722092.499494/1‑59 (MQ=60)
                                                                                           GCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACG        <  SRR3722092.391881/100‑1 (MQ=60)
                                                                                                GATACCACCCCGGGCAACGTTTACGCTCATGCTGACTACAAGATctgtctcttatacacatctgacgctgccgccgaccatcttagtgtagatctctacg   >  SRR3722092.420736/1‑44 (MQ=60)
                                                                                                 ATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACctgtctc  >  SRR3722092.178196/1‑93 (MQ=60)
                                                                                                   |                                                                                                 
GAAGTAACTTTAGGTACTCCGGTAGTAATTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATGCCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAA  >  NZ_CP009273/3179172‑3179368

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