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 874,237 C→T L14F (CTT→TTT)  yliI → aldose sugar dehydrogenase YliI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273874,2370CT91.7% 74.2 / ‑1.3 24L14F (CTT→TTT) yliIaldose sugar dehydrogenase YliI
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (12/10);  total (13/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

ATCTCCTCCACTGACTACGCTTTAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTCTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACCATCCCTGGGCACTGGC  >  NZ_CP009273/874150‑874319
                                                                                       |                                                                                  
aTCTCCTCCACTGACTACGCTTTAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGtttt                                                                                  >  2:318791/1‑90 (MQ=255)
  ctcctcCACTGACTACGCTTTAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGtttttc                                                                                <  1:250444/90‑1 (MQ=255)
       ccACTGACTACGCTTTAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCg                                                                           >  2:477673/1‑90 (MQ=255)
        cACTGACTACGCTTTAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGc                                                                          <  1:140946/90‑1 (MQ=255)
          cTGACTACGCTTTAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGctc                                                                        <  1:477673/90‑1 (MQ=255)
            gACTACGCTTTAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGctctc                                                                      >  1:265083/1‑90 (MQ=255)
            gACTACGCTTTAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGctctc                                                                      >  2:104332/1‑90 (MQ=255)
              cTACGCTTTAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTg                                                                    >  2:417114/1‑90 (MQ=255)
                     ttAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCt                                                             >  2:441282/1‑90 (MQ=255)
                                       ggAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCaa                                                      <  1:228374/79‑1 (MQ=255)
                                       ggAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCaa                                                      >  2:228374/1‑79 (MQ=255)
                                       ggAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTc                                           >  2:433451/1‑90 (MQ=255)
                                             gTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTCTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTa                                     >  1:1421/1‑90 (MQ=255)
                                                    aTCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAg                              >  2:711761/1‑90 (MQ=255)
                                                     tCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTCTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGa                             <  2:1421/90‑1 (MQ=255)
                                                      cGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGAc                            <  2:385416/90‑1 (MQ=255)
                                                       gACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACa                           >  1:146250/1‑90 (MQ=255)
                                                           aTCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACt                       >  1:255289/1‑90 (MQ=255)
                                                                ttttCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGAcc                  <  1:643102/90‑1 (MQ=255)
                                                                  ttCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACCAt                <  2:265083/90‑1 (MQ=255)
                                                                       tgtgCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACCATCCCTg           <  1:318791/90‑1 (MQ=255)
                                                                       tgtgCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACCATCCCTg           >  1:17436/1‑90 (MQ=255)
                                                                       tgtgCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACCATCCCTg           >  2:459497/1‑90 (MQ=255)
                                                                               ttATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACCATCCCTGGGCACTgg   <  1:319686/90‑1 (MQ=255)
                                                                               ttATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACCATCCCTGGGCACTgg   <  1:20026/90‑1 (MQ=255)
                                                                                tATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACCATCCCTGGGCACTGGc  <  1:136543/90‑1 (MQ=255)
                                                                                       |                                                                                  
ATCTCCTCCACTGACTACGCTTTAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTCTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACCATCCCTGGGCACTGGC  >  NZ_CP009273/874150‑874319

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
CTCCTCCACTGACTACGCTTTAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTCTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACCATCCCTGGGCACTGGCCTTTTTACCC  >  NZ_CP009273/874152‑874329
                                                                                     |                                                                                            
CTCCTCCACTGACTACGCTTTAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTC                                                                                <  SRR3722112.253403/100‑1 (MQ=60)
CTCCTCCACTGACTACGCTTTAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTC                                                                                >  SRR3722112.268198/1‑100 (MQ=60)
      CACTGACTACGCTTTAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCG                                                                          <  SRR3722112.142710/100‑1 (MQ=60)
        CTGACTACGCTTTAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGC                                                                        <  SRR3722112.484515/100‑1 (MQ=60)
                          GGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGT                                                      <  SRR3722112.231119/100‑1 (MQ=60)
                                 ATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTCTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTA                                               >  SRR3722112.1437/1‑100 (MQ=60)
                                           GTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACA                                     >  SRR3722112.148065/1‑100 (MQ=60)
                                               TGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACT                                 >  SRR3722112.258304/1‑100 (MQ=60)
                                                           CCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACCATCCCTG                     >  SRR3722112.17660/1‑100 (MQ=60)
                                                              TTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACCATCCCTGGGC                  <  SRR3722112.653109/100‑1 (MQ=60)
                                                                     TGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACCATCCCTGGGCACTGGCC           <  SRR3722112.322515/100‑1 (MQ=60)
                                                                             TTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACCATCCCTGGGCACTGGCCTTTTTACC   <  SRR3722112.20275/100‑1 (MQ=60)
                                                                             TTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACCATCCCTGGGCACTGGCCTTTTTACC   <  SRR3722112.323424/100‑1 (MQ=60)
                                                                              TATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACCATCCCTGGGCACTGGCCTTTTTACCC  <  SRR3722112.138243/100‑1 (MQ=60)
                                                                                     |                                                                                            
CTCCTCCACTGACTACGCTTTAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTCTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACCATCCCTGGGCACTGGCCTTTTTACCC  >  NZ_CP009273/874152‑874329

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: