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 I206 R1 206 23.7 1314362 96.8% 1272302 86.2

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,003,480 T→C F60S (TTC→TCC)  yedK → SOS response‑associated peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,003,4800TC96.0% 72.9 / NA 25F60S (TTC→TCC) yedKSOS response‑associated peptidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (8/16);  minor base A (1/0);  total (9/16)
Fisher's exact test for biased strand distribution p-value = 3.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACCCATTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTTCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATGTTTAAAC  >  NZ_CP009273/2003393‑2003566
                                                                                       |                                                                                      
aaCCCATTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCt                                                                                      <  2:133799/90‑1 (MQ=255)
    cATTGGCAGATACAACGTCGCGCCCGCGACTAAAGTCCTACTGCTCAATACACGTGATGAACACCTTCCTCTGACTCCGGTTTACTgctg                                                                                  >  2:335011/1‑87 (MQ=255)
      ttGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGat                                                                                >  1:468760/1‑90 (MQ=255)
       tGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGata                                                                               <  2:630281/90‑1 (MQ=255)
         gCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATg                                                                             <  2:165714/90‑1 (MQ=255)
              tACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTccc                                                                        <  2:115519/90‑1 (MQ=255)
                    gTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAtgg                                                                  <  2:272461/90‑1 (MQ=255)
                       gcgcCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGatat                                                                              >  2:75781/1‑75 (MQ=255)
                       gcgcCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGatat                                                                              <  1:75781/75‑1 (MQ=255)
                       gcgcCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAtggtgg                                                               >  1:250909/1‑90 (MQ=255)
                                    aaGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACcgccgc                                                  <  2:140275/90‑1 (MQ=255)
                                             tGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGACAAACCGCCGCTGATTAACg                                         >  2:61440/1‑90 (MQ=255)
                                              gCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGc                                        >  1:568653/1‑90 (MQ=255)
                                                tCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGatat                                                                              <  1:543693/50‑1 (MQ=255)
                                                tCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGatat                                                                              >  2:543693/1‑50 (MQ=255)
                                                tCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGccc                                      <  1:82083/90‑1 (MQ=255)
                                                                  acCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCaa                    >  2:365258/1‑90 (MQ=255)
                                                                    cTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAGACTGCGGCAAcc                  <  1:278566/90‑1 (MQ=255)
                                                                    cTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAAcc                  <  1:558745/90‑1 (MQ=255)
                                                                    cTTCATCTGGATCCGGTTTCCTGGGGATATGCCCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAAcc                  <  2:477458/90‑1 (MQ=255)
                                                                          cTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGt            <  1:515513/90‑1 (MQ=255)
                                                                           tGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTa           <  1:577745/90‑1 (MQ=255)
                                                                           tGGATCCGGTTTCCTGGGGATATCCCCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTa           <  2:609560/90‑1 (MQ=255)
                                                                             gATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGAt                                                            <  2:45945/39‑1 (MQ=255)
                                                                             gATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGAt                                                            >  1:45945/1‑39 (MQ=39)
                                                                                    tttcCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATGTTTAAAc  <  2:568653/90‑1 (MQ=255)
                                                                                       |                                                                                      
AACCCATTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTTCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATGTTTAAAC  >  NZ_CP009273/2003393‑2003566

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 24 ≤ 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
CTATGATCCCGAACCCATTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTTCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATGTTTAAACCGCTCTGGCAA  >  NZ_CP009273/2003382‑2003577
                                                                                                  |                                                                                                 
CTATGATCCCGAACCCATTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCC                                                                                                  >  SRR3722088.269352/1‑100 (MQ=60)
       CCCGAACCCATTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGAT                                                                                           >  SRR3722088.474269/1‑100 (MQ=60)
                   cgtgggctcggagatgtgtataagagacagGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAT                                                                               <  SRR3722088.550322/70‑1 (MQ=60)
                   GACAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAT                                                                               <  SRR3722088.76549/100‑1 (MQ=60)
                        ATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGG                                                                          >  SRR3722088.253357/1‑100 (MQ=60)
                                               AAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGC                                                   >  SRR3722088.575623/1‑100 (MQ=60)
                                                           TCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAAC                                       <  SRR3722088.82907/100‑1 (MQ=60)
                                                                              CCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCctgtctcttatacacatctccgagcccacgagaccgtacta                    >  SRR3722088.46415/1‑59 (MQ=60)
                                                                               CTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAGACTGCGGCAACCAGTCGTATGT                   <  SRR3722088.281251/100‑1 (MQ=60)
                                                                               CTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATGT                   <  SRR3722088.565606/100‑1 (MQ=60)
                                                                                     CTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATGTTTAAAC             <  SRR3722088.521733/100‑1 (MQ=60)
                                                                                      TGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATGTTTAAACC            <  SRR3722088.584829/100‑1 (MQ=60)
                                                                                                TTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATGTTTAAACCGCTCTGGCAA  >  SRR3722088.560476/1‑100 (MQ=60)
                                                                                                  |                                                                                                 
CTATGATCCCGAACCCATTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTTCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATGTTTAAACCGCTCTGGCAA  >  NZ_CP009273/2003382‑2003577

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