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 I199 R1 220 57.5 3163680 97.0% 3068769 86.1

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,5080TG97.1% 103.4 / ‑3.7 34T250P (ACC→CCC) hcaRDNA‑binding transcriptional regulator HcaR
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (15/18);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

GGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGAATATTCGTT  >  NZ_CP009273/2661423‑2661590
                                                                                     |                                                                                  
ggCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTa                                                                                >  2:822016/1‑90 (MQ=255)
 gCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTaa                                                                               <  1:455619/90‑1 (MQ=255)
   ttCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTaaaa                                                                             <  1:925026/90‑1 (MQ=255)
     cATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAatt                                                                           >  1:948775/1‑90 (MQ=255)
      aTCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAatta                                                                          <  2:987006/90‑1 (MQ=255)
       tctcCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAattat                                                                         <  1:1097578/90‑1 (MQ=255)
              tttCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCAtata                                                                  <  2:483946/90‑1 (MQ=255)
               ttcttcCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATaa                                                                 >  1:542806/1‑90 (MQ=255)
                tcttcCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAAc                                                                <  2:884117/90‑1 (MQ=255)
                tcttcCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAAc                                                                <  1:640412/90‑1 (MQ=255)
                  ttcCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCg                                                              <  1:623295/90‑1 (MQ=255)
                  ttcCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCg                                                              <  2:846036/90‑1 (MQ=255)
                  ttcCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCg                                                              <  1:110462/90‑1 (MQ=255)
                  ttcCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCg                                                              <  2:345244/90‑1 (MQ=255)
                    ccACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCggg                                                            >  1:1093101/1‑90 (MQ=255)
                         ccATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCa                                                       >  2:521741/1‑90 (MQ=255)
                           aTCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCaaa                                                     >  1:1578168/1‑90 (MQ=255)
                               gTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGa                                                 <  2:674011/90‑1 (MQ=255)
                               gTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGa                                                 <  2:1130019/90‑1 (MQ=255)
                                 aaaGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACg                                               >  1:1503148/1‑90 (MQ=255)
                                 aaaGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACg                                               >  2:17161/1‑90 (MQ=255)
                                    gCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACACCCTGTCCGGGATTAAAATTATTCATATAACTGGGTATCAAAGTGACGccc                                            <  1:1479447/90‑1 (MQ=255)
                                     cAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCa                                           >  1:1460233/1‑90 (MQ=255)
                                      aaTGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAg                                          >  2:1044315/1‑90 (MQ=255)
                                       aTGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGc                                         >  2:1552280/1‑90 (MQ=255)
                                       aTGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGc                                         >  2:1485974/1‑90 (MQ=255)
                                            aGGCACATTTCCGGCAATGGGACGAAAAACAACCTTTCCGGGATTAAAATTATTCATAGAACCGGGTATCAAAGTGACGCCCACCCCCAt                                    >  2:1479447/1‑90 (MQ=255)
                                                           aaTGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCAt                     <  2:962133/90‑1 (MQ=255)
                                                                  cGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCAcc              <  2:1460233/90‑1 (MQ=255)
                                                                  cGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCAcc              <  2:839980/90‑1 (MQ=255)
                                                                   gAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCa             <  2:246119/90‑1 (MQ=255)
                                                                     aaaaCAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGa           >  1:935478/1‑90 (MQ=255)
                                                                              tGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGAATATTCGtt  >  1:1202669/1‑90 (MQ=255)
                                                                              tGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGAATATTCGtt  >  1:1502261/1‑90 (MQ=255)
                                                                                     |                                                                                  
GGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGAATATTCGTT  >  NZ_CP009273/2661423‑2661590

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ 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
ATCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGAATATTCGTTGCCACCTGGAC  >  NZ_CP009273/2661410‑2661601
                                                                                                  |                                                                                             
ATCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGA                                                                                              >  SRR3722078.1537928/1‑100 (MQ=60)
 TCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGAT                                                                                             <  SRR3722078.299176/100‑1 (MQ=60)
      TAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAA                                                                                        <  SRR3722078.1425325/100‑1 (MQ=60)
        ATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATT                                                                                      >  SRR3722078.962072/1‑100 (MQ=60)
              GCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCAT                                                                                <  SRR3722078.461063/100‑1 (MQ=60)
                TTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATAT                                                                              <  SRR3722078.937715/100‑1 (MQ=60)
                  CATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAA                                                                            >  SRR3722078.549218/1‑100 (MQ=60)
                    TCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACC                                                                          <  SRR3722078.1114276/100‑1 (MQ=60)
                       CCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGG                                                                       >  SRR3722078.1109708/1‑100 (MQ=60)
                             TCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAA                                                                 <  SRR3722078.647781/100‑1 (MQ=60)
                              CTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAA                                                                >  SRR3722078.1603854/1‑100 (MQ=60)
                               TTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAG                                                               <  SRR3722078.111766/100‑1 (MQ=60)
                               TTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAG                                                               <  SRR3722078.630477/100‑1 (MQ=60)
                                    CGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACG                                                          >  SRR3722078.1527666/1‑100 (MQ=60)
                                        ATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCA                                                      >  SRR3722078.1484050/1‑100 (MQ=60)
                                                 nnnnnnnnnnnnnnnnttaaGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGC                                             <  SRR3722078.1452875/80‑1 (MQ=60)
                                                 GCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACACCCTGTCCGGGATTAAAATTATTCATATAACTGGGTATCAAAGTGACGCCCAGCCCCATGC                                             <  SRR3722078.1503573/100‑1 (MQ=60)
                                                                tcgtcggcagcgtcagatgtgtataagagacaGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGG                              <  SRR3722078.1522781/68‑1 (MQ=60)
                                                                        AATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGA                      >  SRR3722078.948422/1‑100 (MQ=60)
                                                                                 AAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGAATATTCGTT             >  SRR3722078.1526764/1‑100 (MQ=60)
                                                                                 AAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGAATATTCGTT             >  SRR3722078.1221531/1‑100 (MQ=60)
                                                                                            GTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGAATATTCGTTGCCACCTGGAC  >  SRR3722078.34327/1‑100 (MQ=60)
                                                                                                  |                                                                                             
ATCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGAATATTCGTTGCCACCTGGAC  >  NZ_CP009273/2661410‑2661601

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