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 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,2370CT95.6% 159.3 / ‑2.2 45L14F (CTT→TTT) yliIaldose sugar dehydrogenase YliI
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (24/19);  total (25/20)
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.30e-01

ATCTCCTCCACTGACTACGCTTTAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTCTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACCATCCCTGGGCACTGGCCT  >  NZ_CP009273/874150‑874321
                                                                                       |                                                                                    
aTCTCCTCCACTGACTACGCTTTAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGtttt                                                                                    <  2:869268/90‑1 (MQ=255)
            gACTACGCTTTAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGctctc                                                                        >  1:767841/1‑90 (MQ=255)
            gACTACGCTTTAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGctctc                                                                        >  2:1059790/1‑90 (MQ=255)
            gACTACGCTTTAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGctctc                                                                        >  2:767678/1‑90 (MQ=255)
            gACTACGCTTTAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTCTTTCTTCCGctctc                                                                        >  2:512721/1‑90 (MQ=255)
                   cTTTAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGggcgg                                                                 >  1:139618/1‑90 (MQ=255)
                      tAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTc                                                              >  1:996151/1‑90 (MQ=255)
                      tAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTc                                                              >  2:1088102/1‑90 (MQ=255)
                      tAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTc                                                              >  1:148842/1‑90 (MQ=255)
                      tAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTc                                                              >  2:405824/1‑90 (MQ=255)
                           cAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCa                                                         <  1:1424055/90‑1 (MQ=255)
                            agagTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCaa                                                        <  1:1000195/90‑1 (MQ=255)
                             gagTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAAc                                                       >  2:217741/1‑90 (MQ=255)
                             gagTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAAc                                                       >  2:1581779/1‑90 (MQ=255)
                             gagTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAAc                                                       >  1:653855/1‑90 (MQ=255)
                                  aaTCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTCTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTaa                                                  <  2:567059/90‑1 (MQ=255)
                                    tCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAAt                                                <  2:541726/90‑1 (MQ=255)
                                       ggAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTc                                             >  2:1007149/1‑90 (MQ=255)
                                            cGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGt                                        <  1:1373596/90‑1 (MQ=255)
                                                    aTCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAg                                >  2:31028/1‑90 (MQ=255)
                                                     tCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGa                               <  2:94560/90‑1 (MQ=255)
                                                      cGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGAc                              >  1:1429565/1‑90 (MQ=255)
                                                       gACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGc                                                                            <  2:955524/43‑1 (MQ=255)
                                                       gACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGc                                                                            >  1:955524/1‑43 (MQ=255)
                                                       gACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACa                             >  1:81343/1‑90 (MQ=255)
                                                       gACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACa                             >  1:1013332/1‑90 (MQ=255)
                                                       gACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACa                             >  1:1127507/1‑90 (MQ=255)
                                                       gACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACa                             >  2:759267/1‑90 (MQ=255)
                                                       gACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACa                             >  1:15347/1‑90 (MQ=255)
                                                       gACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACa                             >  1:338643/1‑90 (MQ=255)
                                                                ttttCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGTCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGAcc                    <  1:767678/90‑1 (MQ=255)
                                                                ttttCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCTAGACAAACTCGAcc                    <  1:1426027/90‑1 (MQ=255)
                                                                ttttCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGAcc                    <  1:326920/90‑1 (MQ=255)
                                                                ttttCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGAcc                    <  1:217741/90‑1 (MQ=255)
                                                                ttttCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGAcc                    <  2:338643/90‑1 (MQ=255)
                                                                ttttCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGAcc                    <  1:209907/90‑1 (MQ=255)
                                                                ttttCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGAcc                    <  1:1581779/90‑1 (MQ=255)
                                                                ttttCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGAcc                    <  2:653855/90‑1 (MQ=255)
                                                                ttttCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGAcc                    <  2:1483000/90‑1 (MQ=255)
                                                                ttttCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGAcc                    >  1:409964/1‑90 (MQ=255)
                                                                ttttCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGAcc                    <  1:738707/90‑1 (MQ=255)
                                                                 tttCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACCa                   <  2:953761/90‑1 (MQ=255)
                                                                       tgtgCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGAc                              >  2:1509029/1‑73 (MQ=255)
                                                                       tgtgCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGAc                              <  1:1509029/73‑1 (MQ=255)
                                                                            cccTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACCATCCCTGGGCAc        <  1:1059790/90‑1 (MQ=255)
                                                                                  tttGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACCATCCCTGGGCACTGGCCt  >  2:916759/1‑90 (MQ=255)
                                                                                       |                                                                                    
ATCTCCTCCACTGACTACGCTTTAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTCTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACCATCCCTGGGCACTGGCCT  >  NZ_CP009273/874150‑874321

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TTGTTATCTCCTCCACTGACTACGCTTTAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTCTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACCATCCCTGGGCACTGGCCTTTTTACCCGATAA  >  NZ_CP009273/874145‑874334
                                                                                            |                                                                                                 
TTGTTATCTCCTCCACTGACTACGCTTTAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTC                                                                                            <  SRR3722078.1480561/100‑1 (MQ=60)
       CTCCTCCACTGACTACGCTTTAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTC                                                                                     >  SRR3722078.777744/1‑100 (MQ=60)
              ACTGACTACGCTTTAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGG                                                                              >  SRR3722078.141249/1‑100 (MQ=60)
                 GACTACGCTTTAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTC                                                                           >  SRR3722078.1010650/1‑100 (MQ=60)
                 GACTACGCTTTAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTC                                                                           >  SRR3722078.150573/1‑100 (MQ=60)
                        CTTTAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAAC                                                                    >  SRR3722078.661347/1‑100 (MQ=60)
                                CAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATG                                                            <  SRR3722078.1447188/100‑1 (MQ=60)
                                 AGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGT                                                           <  SRR3722078.1014802/100‑1 (MQ=60)
                                                 CGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGAC                                           <  SRR3722078.1395736/100‑1 (MQ=60)
                                                 CGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGAC                                           >  SRR3722078.1452806/1‑100 (MQ=60)
                                                  GTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGctgtctcttatacacatctgacgctgccgacgatggc                                          >  SRR3722078.969013/1‑63 (MQ=60)
                                                  GTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACA                                          >  SRR3722078.1144752/1‑100 (MQ=60)
                                                  GTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACA                                          >  SRR3722078.1028258/1‑100 (MQ=60)
                                                  GTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACA                                          >  SRR3722078.82295/1‑100 (MQ=60)
                                                  GTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACA                                          >  SRR3722078.342787/1‑100 (MQ=60)
                                                  GTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACA                                          >  SRR3722078.15535/1‑100 (MQ=60)
                                                           gagacagCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACC                                 <  SRR3722078.1533637/93‑1 (MQ=60)
                                                           CGACAATCCTTTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACC                                 >  SRR3722078.414920/1‑100 (MQ=60)
                                                                     tataagagacagCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACCATCCCTGGGC                       <  SRR3722078.1604433/88‑1 (MQ=60)
                                                                     TTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGTCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACCATCCCTGGGC                       <  SRR3722078.777579/100‑1 (MQ=60)
                                                                     TTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCTAGACAAACTCGACCATCCCTGGGC                       <  SRR3722078.1449203/100‑1 (MQ=60)
                                                                     TTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACCATCCCTGGGC                       <  SRR3722078.1607527/100‑1 (MQ=60)
                                                                     TTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACCATCCCTGGGC                       <  SRR3722078.212445/100‑1 (MQ=60)
                                                                     TTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACCATCCCTGGGC                       <  SRR3722078.220376/100‑1 (MQ=60)
                                                                     TTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACCATCCCTGGGC                       <  SRR3722078.330931/100‑1 (MQ=60)
                                                                     TTTTCCTTGTGCCCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACCATCCCTGGGC                       <  SRR3722078.748032/100‑1 (MQ=60)
                                                                                 CCCTTATTTGTTTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACCATCCCTGGGCACTGGCCTTTTT           <  SRR3722078.1075734/100‑1 (MQ=60)
                                                                                          GTTTTTCTTCCGCTTTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACCATCCCTGGGCACTGGCCTTTTTACCCGATAA  <  SRR3722078.1068085/100‑1 (MQ=60)
                                                                                            |                                                                                                 
TTGTTATCTCCTCCACTGACTACGCTTTAAGCCAGAGTCAATCCGGAGGCGTTATGCATCGACAATCCTTTTTCCTTGTGCCCCTTATTTGTCTTTCTTCCGCTCTCTGGGCGGCTCCTGCAACGGTAAATGTCGAAGTACTGCAAGACAAACTCGACCATCCCTGGGCACTGGCCTTTTTACCCGATAA  >  NZ_CP009273/874145‑874334

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