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 I228 R1 214 26.5 1448192 97.1% 1406194 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,5550TC92.6% 71.7 / ‑0.4 27F232S (TTC→TCC) yqeItranscriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (12/13);  total (14/13)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.67e-01

AATTGATGGGTTGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTTCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCC  >  NZ_CP009273/2982469‑2982641
                                                                                      |                                                                                      
aaTTGATGGGTTGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTcccc                                                                                     <  1:286949/90‑1 (MQ=255)
      tGGGTTGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGc                                                                               >  2:12987/1‑90 (MQ=255)
         gTTGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCaaa                                                                            >  2:252357/1‑90 (MQ=255)
          ttGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCaaaa                                                                           <  2:155796/90‑1 (MQ=255)
                  gTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCAcc                                                                   >  1:49249/1‑90 (MQ=255)
                         aaaTATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGgatg                                                            <  1:678998/90‑1 (MQ=255)
                         aaaTATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGgatg                                                            <  2:337857/90‑1 (MQ=255)
                         aaaTATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGgatg                                                            <  1:95466/90‑1 (MQ=255)
                         aaaTATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGgatg                                                            <  2:60694/90‑1 (MQ=255)
                         aaaTATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGgatg                                                            <  1:23264/90‑1 (MQ=255)
                         aaaTATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGgatg                                                            <  2:679280/90‑1 (MQ=255)
                           atatAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGgatgat                                                          <  2:379814/90‑1 (MQ=255)
                                aaaCACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTa                                                     >  2:612674/1‑90 (MQ=255)
                                 aaCACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTAt                                                    >  2:156820/1‑90 (MQ=255)
                                     cGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCa                                                >  2:5399/1‑90 (MQ=255)
                                     cGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCa                                                >  2:42177/1‑90 (MQ=255)
                                      ggATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTTCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCAt                                               >  2:298635/1‑90 (MQ=255)
                                         tACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCa                                            <  1:406895/90‑1 (MQ=255)
                                                  cGGGGATTAACTGTAAAAAACATCCGTGGTTATACTTCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACt                                   >  1:52193/1‑90 (MQ=255)
                                                   ggggATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACta                                  >  1:641213/1‑90 (MQ=255)
                                                              gTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATg                       >  1:423889/1‑90 (MQ=255)
                                                                     aCATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGc                >  1:173720/1‑90 (MQ=255)
                                                                      cATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCa               >  2:216217/1‑90 (MQ=255)
                                                                      cATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCa               >  1:153805/1‑90 (MQ=255)
                                                                         ccGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAAc            <  1:612674/90‑1 (MQ=255)
                                                                                 ataCTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAAc    <  2:131290/90‑1 (MQ=255)
                                                                                   aCTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAAccc  <  2:153805/90‑1 (MQ=255)
                                                                                   aCTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAAccc  <  1:42177/90‑1 (MQ=255)
                                                                                      |                                                                                      
AATTGATGGGTTGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTTCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCC  >  NZ_CP009273/2982469‑2982641

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
TAACAGAAGATACAATTGATGGGTTGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTTCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCG  >  NZ_CP009273/2982456‑2982651
                                                                                                   |                                                                                                
TAACAGAAGATACAATTGATGGGTTGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTC                                                                                                  <  SRR3722112.189513/100‑1 (MQ=60)
TAACAGAAGATACAATTGATGGGTTGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTC                                                                                                  <  SRR3722112.200280/100‑1 (MQ=60)
TAACAGAAGATACAATTGATGGGTTGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTC                                                                                                  <  SRR3722112.366649/100‑1 (MQ=60)
TAACAGAAGATACAATTGATGGGTTGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTC                                                                                                  >  SRR3722112.427214/1‑100 (MQ=60)
             AATTGATGGGTTGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAA                                                                                     <  SRR3722112.290282/100‑1 (MQ=60)
                     GGTTGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACC                                                                             >  SRR3722112.49854/1‑100 (MQ=60)
                                      AAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGG                                                            <  SRR3722112.23558/100‑1 (MQ=60)
                                      AAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGG                                                            <  SRR3722112.96656/100‑1 (MQ=60)
                                      AAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGG                                                            <  SRR3722112.689615/100‑1 (MQ=60)
                                                CACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCtgtctcttatacacatctgacgctgccgacgaatgtcaa                                                  >  SRR3722112.98798/1‑61 (MQ=60)
                                                     ATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTTCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACT                                             >  SRR3722112.52842/1‑100 (MQ=60)
                                                      TACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTA                                            <  SRR3722112.412277/100‑1 (MQ=60)
                                                      TACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTA                                            >  SRR3722112.651183/1‑100 (MQ=60)
                                                                 GGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATG                                 >  SRR3722112.429543/1‑100 (MQ=60)
                                                                        ACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGC                          >  SRR3722112.175848/1‑100 (MQ=60)
                                                                         CTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCA                         >  SRR3722112.155694/1‑100 (MQ=60)
                                                                                      CCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCC            <  SRR3722112.622087/100‑1 (MQ=60)
                                                                                         TGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGT         >  SRR3722112.311057/1‑100 (MQ=60)
                                                                                                ACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCG  <  SRR3722112.42681/100‑1 (MQ=60)
                                                                                                   |                                                                                                
TAACAGAAGATACAATTGATGGGTTGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTTCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCG  >  NZ_CP009273/2982456‑2982651

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