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 F14 I0 R1 248 61.6 2601574 98.7% 2567753 146.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947_CJ‑RC 4,946,923 G→A 33.3% intergenic (+580/+218) PP_5673 → / ← flhB D‑alanine‑‑D‑alanine ligase/flagellin export apparatus substrate specificity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947_CJ‑RC4,946,9230GA33.3% 44.7 / 29.9 39intergenic (+580/+218)PP_5673/flhBD‑alanine‑‑D‑alanine ligase/flagellin export apparatus substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base G (7/19);  new base A (3/10);  total (10/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

GCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGC  >  NC_002947_CJ‑RC/4946777‑4947024
                                                                                                                                                  |                                                                                                     
gCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAgct                                                                                                     >  1:56501/1‑149 (MQ=255)
 cAAAATTTTCACGCCCAGGGGGCGCCTGCGGCGGGCATTGCTAGAAGCTTTCCAACTTTGTTATCTGCTGTGACGCGGGGGGTGCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAgctc                                                                                                    <  1:265079/149‑1 (MQ=37)
 cAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGGGAATCGGGGTGTTCGCCGTGAGTTTTTTGTCGCCTGTGACATCGCGCGCCGCCCGCGCGGCGCATCGCGAGCAAgct                                                                                                     <  1:14636/148‑1 (MQ=255)
    aaCTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAgctcgct                                                                                                 >  5:207129/1‑149 (MQ=255)
      cTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTcc                                                                                               <  5:148718/149‑1 (MQ=255)
             gCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACAttt                                                                                        >  3:290394/1‑149 (MQ=255)
                 aGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGttt                                                                                    <  2:38960/149‑1 (MQ=255)
                  gggggCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTc                                                                                   <  2:106383/149‑1 (MQ=255)
                      gcgcCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCGCCTACAATTTTTTggggg                                                                               >  3:38783/1‑148 (MQ=255)
                             cGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTa                                                                         <  1:67886/148‑1 (MQ=255)
                               gCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGAGCCGCCCGCGCGGCGCATCGAGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTAtta                                                                      >  7:186537/1‑148 (MQ=255)
                                   gCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGc                                                                  <  5:14639/149‑1 (MQ=255)
                                      ttGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCgg                                                                >  1:30202/1‑148 (MQ=255)
                                            gAAGCTTTCCAACTTTGTGGTCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTa                                                          <  8:305250/148‑1 (MQ=255)
                                              aggtgttcgccgggggattttcagcgcctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccgcgcggcgcttcgcggGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTgcgc                                                       <  4:295723‑M2/53‑1 (MQ=255)
                                                gtgttcgccgggggattttcagcgcctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccgcgcggcgcttcgcggGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTgcgcgc                                                     <  1:237146‑M2/55‑1 (MQ=255)
                                                  gttcgccgggggattttcagcgcctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccgcgcggcgcttcgcggGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGa                                                   >  2:246236‑M2/93‑149 (MQ=255)
                                                   ttcgccgggggattttcagcgcctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccgcgcggcgcttcgcggGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGAc                                                  >  2:13009‑M2/92‑149 (MQ=255)
                                                    ccAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGAcc                                                 <  1:232281/149‑1 (MQ=255)
                                                          ggggattttcagcgcctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccgcgcggcgcttcgcggGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCcttg                                           <  2:85289‑M2/65‑2 (MQ=255)
                                                          ggggattttcagcgcctgggaggtgtgcgccgggagagtttcagcgcctgggagatcgagcgccgcccgcgcggcgcttcgcggGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCcttg                                           <  5:137023‑M2/65‑2 (MQ=255)
                                                           gggattttcagcgcctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccgcgcggcgcttcgcggGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGccttgg                                          >  5:100091‑M2/84‑147 (MQ=255)
                                                            ggattttcagcgcctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccgcgcggcgcttcgcggGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCgccttggg                                         <  5:223532‑M2/67‑4 (MQ=255)
                                                             tgATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTcgc                                         <  5:117085/148‑1 (MQ=255)
                                                             tgATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTcgc                                         <  8:277661/148‑1 (MQ=255)
                                                             tgATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTcgc                                         <  5:121981/148‑1 (MQ=255)
                                                                 ttcagcgcctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccgcgcggcgcttcgcggGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCgaccgccttggggcat                                     <  6:100091‑M2/71‑8 (MQ=255)
                                                                 ttcagcgcctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccgcgcggcgcttcgcggGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCgaccgccttggggcat                                     <  1:13009‑M2/71‑8 (MQ=255)
                                                                 ttcagcgcctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccgcgcggcgcttcgcggGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCgaccgccttggggcat                                     <  2:271339‑M2/71‑8 (MQ=255)
                                                                   gctgTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACAcc                                  >  2:180185/1‑149 (MQ=255)
                                                                    ctgTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTa                                      <  1:87682/144‑1 (MQ=255)
                                                                         aaTCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGata                            <  7:79175/149‑1 (MQ=255)
                                                                              gggTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTgcg                       <  6:207129/149‑1 (MQ=255)
                                                                                   ttCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCg                                                                 <  5:119904/102‑1 (MQ=255)
                                                                                   ttCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCg                                                                 >  6:119904/1‑102 (MQ=255)
                                                                                     cgccgggagagtttcagcgcctgggagatcgagcgccgcccgcgcggcgcttcgcggGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTac                 <  4:146378‑M2/91‑1 (MQ=255)
                                                                                     cgccgggagagtttcagcgcctgggagatcgagcgccgcccgcgcggcgcttcgcggGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTac                 <  2:153807‑M2/91‑1 (MQ=255)
                                                                                          tttgttttttGTCGCCTGTGATATGGTGCGCCGCCGGCGCGGCGCTTCGCGAGCAAGCTCGCTCCAACTTTTGTTTCGGGCCAGTTATTGCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTacaaacaa           <  4:6779/144‑1 (MQ=25)
                                                                                             ggTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAgg         <  2:235394/148‑1 (MQ=255)
                                                                                              gTTTTTGTCGCCTGTGAGATCGAGCGCGGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGTGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCg       <  7:310290/149‑1 (MQ=255)
                                                                                               tttttGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCgg      <  1:180185/149‑1 (MQ=255)
                                                                                                 tttGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTc    >  1:175576/1‑149 (MQ=255)
                                                                                                    gTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTcgc  <  4:290394/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                     
GCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGC  >  NC_002947_CJ‑RC/4946777‑4947024

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

GATK/CNVnator alignment

N/A