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 I227 R1 223 21.3 1167702 97.0% 1132670 86.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,982,555 T→C F232S (TTC→TCC)  yqeI → transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,982,5550TC100.0% 86.0 / NA 27F232S (TTC→TCC) yqeItranscriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/15);  total (12/15)

GATGGGTTGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTTCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCC  >  NZ_CP009273/2982473‑2982641
                                                                                  |                                                                                      
gATGGGTTGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCtt                                                                                 >  2:224839/1‑90 (MQ=255)
    ggTTGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCaa                                                                             <  2:427321/90‑1 (MQ=255)
              gTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCAcc                                                                   >  2:572234/1‑90 (MQ=255)
              gTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCAcc                                                                   >  1:386154/1‑90 (MQ=255)
              gTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCAcc                                                                   >  1:392814/1‑90 (MQ=255)
              gTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCAcc                                                                   >  2:45279/1‑90 (MQ=255)
              gTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCAcc                                                                   >  2:353544/1‑90 (MQ=255)
               ttCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCg                                                                  <  2:442019/90‑1 (MQ=255)
               ttCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCg                                                                  >  2:372995/1‑90 (MQ=255)
                     aaaTATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGgatg                                                            <  2:460641/90‑1 (MQ=255)
                     aaaTATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGgatg                                                            <  2:458367/90‑1 (MQ=255)
                      aaTATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGgatga                                                           >  2:111074/1‑90 (MQ=255)
                       atatAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGgatgat                                                          <  2:59319/90‑1 (MQ=255)
                       atatAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGgatgat                                                          <  2:330598/90‑1 (MQ=255)
                                         ggATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATaa                                        >  1:519087/1‑90 (MQ=255)
                                         ggATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATaa                                        <  2:519087/90‑1 (MQ=255)
                                           aTTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCGATaaaa                                      <  1:299978/90‑1 (MQ=255)
                                                                  cATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTggta               >  2:403416/1‑88 (MQ=255)
                                                                  cATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCa               >  2:437465/1‑90 (MQ=255)
                                                                    tCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCaaa             >  2:554935/1‑90 (MQ=255)
                                                                     ccGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAAc            <  1:318087/90‑1 (MQ=255)
                                                                     ccGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAAc            <  1:201060/90‑1 (MQ=255)
                                                                     ccGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAAc            <  1:104575/90‑1 (MQ=255)
                                                                               aCTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAAccc  <  2:38272/90‑1 (MQ=255)
                                                                               aCTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAAccc  <  2:386154/90‑1 (MQ=255)
                                                                               aCTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAAccc  <  1:572234/90‑1 (MQ=255)
                                                                               aCTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAAccc  <  1:452211/90‑1 (MQ=255)
                                                                                  |                                                                                      
GATGGGTTGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTTCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCC  >  NZ_CP009273/2982473‑2982641

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ 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
TAACAGAAGATACAATTGATGGGTTGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTTCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCG  >  NZ_CP009273/2982456‑2982651
                                                                                                   |                                                                                                
TAACAGAAGATACAATTGATGGGTTGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTC                                                                                                  <  SRR3722111.507050/100‑1 (MQ=60)
   CAGAAGATACAATTGATGGGTTGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCC                                                                                               <  SRR3722111.112324/100‑1 (MQ=60)
                     GGTTGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACC                                                                             >  SRR3722111.391716/1‑100 (MQ=60)
                     GGTTGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACC                                                                             >  SRR3722111.398509/1‑100 (MQ=60)
                                                CACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAA                                                  >  SRR3722111.527196/1‑100 (MQ=60)
                                                            ATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCGATAAAAACTATAACCA                                      <  SRR3722111.303852/100‑1 (MQ=60)
                                                                              ataagagacaGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTG                    <  SRR3722111.303642/90‑1 (MQ=60)
                                                                                      CCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCC            <  SRR3722111.105874/100‑1 (MQ=60)
                                                                                      CCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCC            <  SRR3722111.203559/100‑1 (MQ=60)
                                                                                      CCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCC            <  SRR3722111.322291/100‑1 (MQ=60)
                                                                                                ACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCG  <  SRR3722111.459070/100‑1 (MQ=60)
                                                                                                ACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCG  <  SRR3722111.581222/100‑1 (MQ=60)
                                                                                                   |                                                                                                
TAACAGAAGATACAATTGATGGGTTGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTTCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCG  >  NZ_CP009273/2982456‑2982651

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