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 4,153,816 T→C S84P (TCA→CCA)  btuB → TonB‑dependent vitamin B12 receptor BtuB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,153,8160TC96.2% 74.9 / ‑4.0 26S84P (TCA→CCA) btuBTonB‑dependent vitamin B12 receptor BtuB
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (10/15);  total (10/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.47e-01

GACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGG  >  NZ_CP009273/4153729‑4153891
                                                                                       |                                                                           
gACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCa                                                                           <  1:395725/90‑1 (MQ=255)
            tCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTAtt                                                               <  2:198445/90‑1 (MQ=255)
               aCCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTcgc                                                            <  1:170672/90‑1 (MQ=255)
               aCCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTcgc                                                            <  2:171238/90‑1 (MQ=255)
               aCCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATCcgc                                                            <  1:298169/90‑1 (MQ=255)
                      tCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACaa                                                     <  2:46654/90‑1 (MQ=255)
                              gtgCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCCGGGCAGCTCCCCTCTATTTTTTTTTGCGGGACAAATGCCAGt                                             >  1:368334/1‑90 (MQ=255)
                                gctgcGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCa                                           <  1:173589/90‑1 (MQ=255)
                                gctgcGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCa                                           <  2:139470/90‑1 (MQ=255)
                                   gcgcCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTTAGGTCAGCTCCCCTCTATTTTTTTTTGCGGGACAAATGCCAGTCAtgt                                        >  2:417292/1‑90 (MQ=255)
                                   gcgcCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCCTCTATTTTTTTTTGCGGTACAAATGCCCGTCAtgt                                        >  2:26477/1‑90 (MQ=255)
                                   gcgcCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGGCAGCTCCCCTCTATTTTTTTTTGCGGGACAAAAGCCAGTCAtgt                                        >  2:461273/1‑90 (MQ=255)
                                   gcgcCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGGCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCAtgt                                        >  1:13820/1‑90 (MQ=255)
                                   gcgcCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGGCAGCTCCCATCTATTTTTATTCGCGGGACAAATGCCAGTCAtgt                                        >  1:262923/1‑90 (MQ=255)
                                     gcCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCCTCTTTTTTTTTTTGCGGGACAAAAGCCAGTCATgtgt                                      >  2:217130/1‑90 (MQ=255)
                                     gcCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATgtgt                                      <  1:417292/90‑1 (MQ=255)
                                        gTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTgg                                   >  1:128397/1‑90 (MQ=255)
                                             ccGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTTTTCGCGGGACAAAAGCCAGTCATGTGTTGGGGTTa                              >  2:195361/1‑90 (MQ=255)
                                                   gTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGAt                        <  2:190162/90‑1 (MQ=255)
                                                   gTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGAt                        <  1:195361/90‑1 (MQ=255)
                                                        tatCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGt                   <  2:522686/90‑1 (MQ=255)
                                                        tatCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGt                   <  1:286382/90‑1 (MQ=255)
                                                                       cggTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGc    >  1:225789/1‑90 (MQ=255)
                                                                        ggTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCg   <  1:398971/90‑1 (MQ=255)
                                                                         gTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCgg  <  2:128397/90‑1 (MQ=255)
                                                                         gTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCgg  <  1:63322/90‑1 (MQ=255)
                                                                                       |                                                                           
GACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGG  >  NZ_CP009273/4153729‑4153891

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 19 ≤ 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
CGTCAGGATATCGACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTG  >  NZ_CP009273/4153717‑4153906
                                                                                                   |                                                                                          
CGTCAGGATATCGACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATACCACCCAAAACGGCGGTTCAGGGCAGCTCC                                                                                            >  SRR3722088.591586/1‑100 (MQ=60)
    AGGATATCGACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATC                                                                                        <  SRR3722088.134581/100‑1 (MQ=60)
       ATATCGACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTAT                                                                                     <  SRR3722088.437094/100‑1 (MQ=60)
            GACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTA                                                                                <  SRR3722088.399992/100‑1 (MQ=60)
                           ACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATG                                                                 <  SRR3722088.172394/100‑1 (MQ=60)
                           ACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATCCGCGGTACAAATG                                                                 <  SRR3722088.301072/100‑1 (MQ=60)
                                GGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCCGGGCAGCTCCCCTCTATTTTTTTTTGCGGGACAAATGCCAGT                                                            >  SRR3722088.372248/1‑100 (MQ=60)
                                     ATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGGCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGT                                                       >  SRR3722088.13973/1‑100 (MQ=60)
                                     ATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGGCAGCTCCCATCTATTTTTATTCGCGGGACAAATGCCAGTCATGT                                                       >  SRR3722088.265473/1‑100 (MQ=60)
                                          GTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCCGGGCAGCTCCCATCTATTTTTTTTTGCGGGACAAAAGCCAGTCATGTGTTGG                                                  >  SRR3722088.664491/1‑100 (MQ=60)
                                          GTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGG                                                  >  SRR3722088.129730/1‑100 (MQ=60)
                                          GTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGGCAGCTCCCATCTATTTTTTTTTGCGGTAAAAATGCCAGgccttggttgg                                                  >  SRR3722088.332420/1‑89 (MQ=60)
                                            GCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTG                                                <  SRR3722088.175333/100‑1 (MQ=60)
                                                 GCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAAT                                           <  SRR3722088.421922/100‑1 (MQ=60)
                                                               GTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCC                             <  SRR3722088.197305/100‑1 (MQ=60)
                                                                    TATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAAT                        <  SRR3722088.289139/100‑1 (MQ=60)
                                                                         CCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGC                   >  SRR3722088.228044/1‑100 (MQ=60)
                                                                                    GGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTG        <  SRR3722088.403300/100‑1 (MQ=60)
                                                                                     GTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGG       <  SRR3722088.63984/100‑1 (MQ=60)
                                                                                          GGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTG  >  SRR3722088.65059/1‑100 (MQ=60)
                                                                                                   |                                                                                          
CGTCAGGATATCGACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTG  >  NZ_CP009273/4153717‑4153906

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