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,950 G→A 30.9% intergenic (+607/+191) 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,9500GA30.9% 56.5 / 29.3 42intergenic (+607/+191)PP_5673/flhBD‑alanine‑‑D‑alanine ligase/flagellin export apparatus substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base G (12/17);  new base A (3/10);  total (15/27)
Fisher's exact test for biased strand distribution p-value = 3.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

CGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCC  >  NC_002947_CJ‑RC/4946806‑4947089
                                                                                                                                                |                                                                                                                                           
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                                                                                                          <  8:277661/148‑1 (MQ=255)
                                tgATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTcgc                                                                                                          <  5:121981/148‑1 (MQ=255)
                                tgATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTcgc                                                                                                          <  5:117085/148‑1 (MQ=255)
                                    ttcagcgcctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccgcgcggcgcttcgcggGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCgaccgccttggggcat                                                                                                      <  1:13009‑M2/71‑8 (MQ=255)
                                    ttcagcgcctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccgcgcggcgcttcgcggGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCgaccgccttggggcat                                                                                                      <  6:100091‑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)
                                                                                                                       tcgctcCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGgagcgagc                  >  2:61004/1‑149 (MQ=255)
                                                                                                                          ctcCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTg               >  1:195018/1‑149 (MQ=255)
                                                                                                                            ccTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCt             >  6:162826/1‑149 (MQ=255)
                                                                                                                             cTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTacaaacaa                                                                            <  7:163650/85‑1 (MQ=255)
                                                                                                                             cTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTacaaacaa                                                                            >  8:163650/1‑85 (MQ=255)
                                                                                                                             cTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTc            >  7:60836/1‑149 (MQ=255)
                                                                                                                                 aTTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGa        >  8:329096/1‑149 (MQ=255)
                                                                                                                                       ttCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGcgcc  >  8:250331/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                           
CGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCC  >  NC_002947_CJ‑RC/4946806‑4947089

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