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 4,593,447 T→C S52P (TCC→CCC)  yjjQ → response regulator transcription factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,593,4470TC95.5% 64.0 / ‑3.8 22S52P (TCC→CCC) yjjQresponse regulator transcription factor
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (12/9);  total (13/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.49e-01

AGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTTCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATC  >  NZ_CP009273/4593364‑4593534
                                                                                   |                                                                                       
agGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAt                                                                                   >  1:311508/1‑90 (MQ=255)
agGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAt                                                                                   >  1:664439/1‑90 (MQ=255)
 ggTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAtt                                                                                  >  2:486727/1‑90 (MQ=255)
  gTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTa                                                                                 >  2:518733/1‑90 (MQ=255)
                             gAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGt                                                      <  1:677694/90‑1 (MQ=255)
                                       ccAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCAc                                            <  2:166183/90‑1 (MQ=255)
                                          gcgcCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGtt                                         >  1:586761/1‑90 (MQ=255)
                                            gcCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCt                                       <  1:515016/90‑1 (MQ=255)
                                             ccTCTGCGGAGGACCTTACCTTATTACGATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTg                                      >  1:452803/1‑90 (MQ=255)
                                                      aggaCCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAAc                             <  1:37313/90‑1 (MQ=255)
                                                      aggaCCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAAc                             <  1:658988/90‑1 (MQ=255)
                                                      aggaCCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAAc                             <  2:586761/90‑1 (MQ=255)
                                                         accttaccttATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACatt                          <  2:227867/90‑1 (MQ=255)
                                                            ttaccttaTTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACattatt                       >  1:77069/1‑90 (MQ=255)
                                                            ttaccttaTTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACattatt                       >  2:53276/1‑90 (MQ=255)
                                                            ttaccttaTTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACattatt                       >  2:671237/1‑90 (MQ=255)
                                                                  tattaCAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCtt                 >  2:379666/1‑90 (MQ=255)
                                                                        aaTTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAAtct           >  2:715637/1‑90 (MQ=255)
                                                                         aTTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAAtctc          >  2:268938/1‑90 (MQ=255)
                                                                          ttACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTca         <  1:273736/90‑1 (MQ=255)
                                                                           tACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTcac        <  1:271242/90‑1 (MQ=255)
                                                                              gtcgtcCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAat     >  2:533053/1‑90 (MQ=255)
                                                                                 gtTCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATc  >  1:320543/1‑90 (MQ=255)
                                                                                   |                                                                                       
AGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTTCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATC  >  NZ_CP009273/4593364‑4593534

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GGGTTAAAAGAGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTTCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGA  >  NZ_CP009273/4593354‑4593538
                                                                                             |                                                                                           
GGGTTAAAAGAGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAT                                                                                       >  SRR3722112.315077/1‑100 (MQ=60)
GGGTTAAAAGAGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAT                                                                                       >  SRR3722112.674800/1‑100 (MQ=60)
                                       GAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCAC                                                <  SRR3722112.688287/100‑1 (MQ=60)
                                          ATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTT                                             >  SRR3722112.595722/1‑100 (MQ=60)
                                             ATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACGATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTG                                          >  SRR3722112.459161/1‑100 (MQ=60)
                                                      GCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAAC                                 <  SRR3722112.522558/100‑1 (MQ=60)
                                                            GCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATT                           >  SRR3722112.78001/1‑100 (MQ=60)
                                                                AGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTC                       <  SRR3722112.37772/100‑1 (MQ=60)
                                                                AGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTC                       <  SRR3722112.669274/100‑1 (MQ=60)
                                                                                 CAATTACGTCGTTCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATC      >  SRR3722112.324296/1‑100 (MQ=60)
                                                                                 CAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATC      >  SRR3722112.296232/1‑100 (MQ=60)
                                                                                 CAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATC      >  SRR3722112.324641/1‑100 (MQ=60)
                                                                                 CAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATC      >  SRR3722112.586532/1‑100 (MQ=60)
                                                                                 CAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATC      >  SRR3722112.669347/1‑100 (MQ=60)
                                                                                    TTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGG   <  SRR3722112.276937/100‑1 (MQ=60)
                                                                                     TACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGA  <  SRR3722112.274412/100‑1 (MQ=60)
                                                                                             |                                                                                           
GGGTTAAAAGAGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTTCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGA  >  NZ_CP009273/4593354‑4593538

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