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,661,508 T→G T250P (ACC→CCC)  hcaR ← DNA‑binding transcriptional regulator HcaR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,661,5080TG100.0% 79.5 / NA 25T250P (ACC→CCC) hcaRDNA‑binding transcriptional regulator HcaR
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (13/12);  total (13/12)

CATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGAATATTCGTTGC  >  NZ_CP009273/2661428‑2661592
                                                                                |                                                                                    
cATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAatt                                                                             >  1:186298/1‑90 (MQ=255)
 aTCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAatta                                                                            <  2:240908/90‑1 (MQ=255)
          ttcttcCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATaa                                                                   >  1:294465/1‑90 (MQ=255)
             ttcCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAattatt                                                                          >  2:487465/1‑80 (MQ=255)
             ttcCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAattatt                                                                          <  1:487465/80‑1 (MQ=255)
             ttcCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCg                                                                <  2:553485/90‑1 (MQ=255)
             ttcCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCg                                                                <  2:535778/90‑1 (MQ=255)
                 aCGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTa                                                            >  1:286102/1‑90 (MQ=255)
                     cATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCaa                                                        <  2:299178/90‑1 (MQ=255)
                      aTCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCaaa                                                       >  1:448542/1‑90 (MQ=255)
                          gTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGa                                                   <  2:485256/90‑1 (MQ=255)
                               gCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGccc                                              <  2:416075/90‑1 (MQ=255)
                                  aTGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGccc                                              >  2:566152/1‑87 (MQ=255)
                                  aTGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGccc                                              <  1:566152/87‑1 (MQ=255)
                                  aTGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGc                                           >  1:40008/1‑90 (MQ=255)
                                       aGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCAt                                      >  2:337430/1‑90 (MQ=255)
                                       aGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCAt                                      >  2:181043/1‑90 (MQ=255)
                                         gCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGc                                    >  2:465616/1‑90 (MQ=255)
                                         gCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGc                                    >  1:271127/1‑90 (MQ=255)
                                              tttCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCccacc                               <  2:462359/90‑1 (MQ=255)
                                              tttCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCccacc                               <  1:382510/90‑1 (MQ=255)
                                              tttCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCccacc                               <  2:536913/90‑1 (MQ=255)
                                                             cGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCAcc                <  2:310740/90‑1 (MQ=255)
                                                                          gTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGAATATTCGTTg   >  2:482415/1‑90 (MQ=255)
                                                                           tCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGAATATTCGTTGc  >  1:302152/1‑90 (MQ=255)
                                                                                |                                                                                    
CATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGAATATTCGTTGC  >  NZ_CP009273/2661428‑2661592

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ 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
ATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGAATATTCGTTGCCACCTGGAC  >  NZ_CP009273/2661418‑2661601
                                                                                          |                                                                                             
ATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATT                                                                                      >  SRR3722111.188624/1‑100 (MQ=60)
          CATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAA                                                                            >  SRR3722111.298242/1‑100 (MQ=60)
             CTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCG                                                                         <  SRR3722111.494950/100‑1 (MQ=60)
                 CCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTA                                                                     >  SRR3722111.289734/1‑100 (MQ=60)
                      CTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAA                                                                >  SRR3722111.455329/1‑100 (MQ=60)
                                  CAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGC                                                    >  SRR3722111.40495/1‑100 (MQ=60)
                                         GCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGC                                             >  SRR3722111.274445/1‑100 (MQ=60)
                                         GCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGC                                             <  SRR3722111.575053/100‑1 (MQ=60)
                                                        TTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGG                              <  SRR3722111.388004/100‑1 (MQ=60)
                                                                           AAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGAATATTCGTTGC           >  SRR3722111.306066/1‑100 (MQ=60)
                                                                                    GTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGAATATTCGTTGCCACCTGGAC  >  SRR3722111.326675/1‑100 (MQ=60)
                                                                                          |                                                                                             
ATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGAATATTCGTTGCCACCTGGAC  >  NZ_CP009273/2661418‑2661601

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: