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 I211 R1 223 18.0 1007644 96.3% 970361 85.7

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% 94.7 / NA 29A167T (GCC→ACC) ygiLfimbrial‑like protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (10/19);  total (10/19)

CTCCGGTAGTAATTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATGCCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACAT  >  NZ_CP009273/3179188‑3179354
                                                                                   |                                                                                   
cTCCGGTAGTAATTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCAccc                                                                               >  2:473046/1‑90 (MQ=255)
cTCCGGTAGTAATTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCAccc                                                                               >  2:280899/1‑90 (MQ=255)
cTCCGGTAGTAATTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCAccc                                                                               >  1:80699/1‑90 (MQ=255)
      tagtaATTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCa                                                                         <  1:178265/90‑1 (MQ=255)
      tagtaATTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCa                                                                         <  1:490692/90‑1 (MQ=255)
      tagtaATTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCa                                                                         <  1:434899/90‑1 (MQ=255)
        gtaATTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAAc                                                                       <  2:429538/90‑1 (MQ=255)
        gtaATTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAAc                                                                       <  1:484540/90‑1 (MQ=255)
          aaTTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGt                                                                     >  2:262125/1‑90 (MQ=255)
                   caacaaTACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCa                                                            >  2:390124/1‑90 (MQ=255)
                      caaTACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGc                                                         <  1:160075/90‑1 (MQ=255)
                        aTACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTg                                                       <  2:48488/90‑1 (MQ=255)
                        aTACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTg                                                       <  1:104356/90‑1 (MQ=255)
                        aTACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTg                                                       <  2:393677/90‑1 (MQ=255)
                        aTACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTg                                                       <  1:482737/90‑1 (MQ=255)
                           ccAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACt                                                    <  1:214122/90‑1 (MQ=255)
                             aaCAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTAc                                                  >  2:216452/1‑90 (MQ=255)
                               cAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACaa                                                >  1:146023/1‑90 (MQ=255)
                                aGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAg                                               >  1:334760/1‑90 (MQ=255)
                                  cTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGAt                                             <  1:176347/90‑1 (MQ=255)
                                  cTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGAt                                             <  1:169750/90‑1 (MQ=255)
                                           aCTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGa                                    >  1:369861/1‑90 (MQ=255)
                                            cTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAg                                   <  1:5184/90‑1 (MQ=255)
                                                      aaGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCAtt                         <  2:221746/90‑1 (MQ=255)
                                                          tCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGc                     <  2:289197/90‑1 (MQ=255)
                                                               tGGAATCCCCAAGCAAAGATACCCCCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGt                <  1:113798/90‑1 (MQ=255)
                                                               tGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGt                <  1:178224/90‑1 (MQ=255)
                                                                          aGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCca     <  2:252035/90‑1 (MQ=255)
                                                                             aaaGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCAcat  >  1:246048/1‑90 (MQ=255)
                                                                                   |                                                                                   
CTCCGGTAGTAATTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATGCCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACAT  >  NZ_CP009273/3179188‑3179354

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ 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
GTAACTTTAGGTACTCCGGTAGTAATTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATGCCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACAT  >  NZ_CP009273/3179175‑3179354
                                                                                                |                                                                                   
GTAACTTTAGGTACTCCGGTAGTAATTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCA                                                                                  <  SRR3722094.173823/100‑1 (MQ=60)
   ACTTTAGGTACTCCGGTAGTAATTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCC                                                                               >  SRR3722094.81852/1‑100 (MQ=60)
   ACTTTAGGTACTCCGGTAGTAATTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCC                                                                               <  SRR3722094.416018/100‑1 (MQ=60)
                   TAGTAATTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGC                                                               <  SRR3722094.499343/100‑1 (MQ=60)
                   TAGTAATTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGC                                                               <  SRR3722094.180732/100‑1 (MQ=60)
                   TAGTAATTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGC                                                               <  SRR3722094.442495/100‑1 (MQ=60)
                     GTAATTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTC                                                             <  SRR3722094.493089/100‑1 (MQ=60)
                                  ACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAA                                                >  SRR3722094.148074/1‑100 (MQ=60)
                                   CAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAG                                               <  SRR3722094.162309/100‑1 (MQ=60)
                                   CAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAG                                               >  SRR3722094.340251/1‑100 (MQ=60)
                                     ATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGAT                                             <  SRR3722094.105853/100‑1 (MQ=60)
                                     ATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGAT                                             <  SRR3722094.491252/100‑1 (MQ=60)
                                        CCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACCAGATTGC                                          <  SRR3722094.217005/100‑1 (MQ=60)
                                              GCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGA                                    >  SRR3722094.376098/1‑100 (MQ=60)
                                               CTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAG                                   <  SRR3722094.178794/100‑1 (MQ=60)
                                               CTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAG                                   <  SRR3722094.172118/100‑1 (MQ=60)
                                                         CTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATT                         <  SRR3722094.5265/100‑1 (MQ=60)
                                                                            TGGAATCCCCAAGCAAAGATACCCCCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCC      <  SRR3722094.115443/100‑1 (MQ=60)
                                                                            TGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCC      <  SRR3722094.180690/100‑1 (MQ=60)
                                                                                ATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACAT  >  SRR3722094.249609/1‑100 (MQ=60)
                                                                                                |                                                                                   
GTAACTTTAGGTACTCCGGTAGTAATTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATGCCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACAT  >  NZ_CP009273/3179175‑3179354

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ 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: