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 5,196,590 T→G 20.3% intergenic (+29/+24) PP_4574 → / ← PP_4575 hypothetical protein/urea amidolyase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947_CJ‑RC5,196,5900TG20.3% 33.7 / 7.8 29intergenic (+29/+24)PP_4574/PP_4575hypothetical protein/urea amidolyase‑like protein
Reads supporting (aligned to +/- strand):  ref base T (12/11);  new base G (4/2);  total (16/13)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

CCAAGCAGGCCTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCGCCAGCGCCAGCGGCGTCAGCG  >  NC_002947_CJ‑RC/5196443‑5196738
                                                                                                                                                   |                                                                                                                                                    
ccAAGCAGGCCTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTc                                                                                                                                                     <  8:19922/149‑1 (MQ=255)
          cTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGaaa                                                                                                                                           <  8:174398/149‑1 (MQ=255)
          cTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGaaa                                                                                                                                           <  4:40501/149‑1 (MQ=255)
          cTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCCTTGGCCGGCGGGCCcccgcca                                                                                                                                           >  4:159117/1‑146 (MQ=255)
             ggAGTACACCGCCTCGTCCCATCTGGTGCGCGTCCTGATCGAGGACAATATGTATCCGCTGACCCTGCAAGGGTTTGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAggg                                                                                                                                        <  6:192667/149‑1 (MQ=255)
             ggAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTGCCGGCCCTTTCGCGGGCGAGCCCGCGAAAg                                                                                                                                          <  1:33159/147‑1 (MQ=255)
             ggAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAgg                                                                                                                                         <  2:74069/148‑1 (MQ=255)
             ggAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAg                                                                                                                                          <  1:281738/147‑1 (MQ=255)
                    acCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCGTTAGCGGGCTCGCCCGCGCAAGGGCccgccc                                                                                                                                 >  1:276287/1‑147 (MQ=255)
                               ccATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTACCTGTACCGGCCCATTGGCCGGCTCGCCCGCGCAAGGGCCAGCACTCTTTCGAcc                                                                                                                       >  5:111185/1‑148 (MQ=255)
                                       tGCGCGACCTGATCGAGTACAAGATGTATCCGCTGACCCTGCATGGGTTTGAGGACGGGATGCGCTACCTTTCCAAGTGAGGTTGCCTGTACGGGCCCTTTCGCGGGCTCGCCCGCGAATGGGCCAGCAGTCTTACGACCAACGCCCCa                                                                                                              <  6:159317/149‑1 (MQ=255)
                                                      aGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAgctgctgct                                                                                               >  8:82798/1‑149 (MQ=255)
                                                       cgtcatgATGTCTCCGCTGCCCCTGCACGGGTTGGCGCACGGGTGGCTCTACGTGTCCATGTGCGGTTGCCTGTACCGGCCCTTTCGCGGGGTCGCCCGCGAACGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGctc                                                                                              <  2:316196/143‑1 (MQ=255)
                                                             gATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCATTCGCGGGCTCGCCCGCCAAAGGGCCAGCAGTATTACGACCAACGCCCCAGGTAAAGCTGCAGCTCTAGCCa                                                                                        >  7:277561/1‑149 (MQ=255)
                                                                         gACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCCGGCGCGCCCGCGAAAGGGGCAGCAGTGGTACGACCCACGCCCCAGGTTAAGCTGCTGCTCTCGCCACGCCTCGTccag                                                                            >  5:60598/1‑147 (MQ=255)
                                                                            ccTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGcac                                                                         >  2:244373/1‑149 (MQ=255)
                                                                                 aaGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTAGCGGGCTCGGCCGCGAAAGGGCCAGCCGTCGTAAGACCAACCCCCCAGGTTAAGCAGCTGCTCTCGCCCGGCCGCGGGCGGCcccaccc                                                                    >  6:256879/1‑148 (MQ=255)
                                                                                          tgcaCGGGACCCGCTACCTGCCAATGTGCGGTGGCCGGAACCGGCCCTTCCGCGGCCTCGCCCGCTAATGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAAcc                                                           <  2:273522/146‑1 (MQ=255)
                                                                                          aggaCGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAAcc                                                           <  1:13472/149‑1 (MQ=255)
                                                                                                  aTGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCATTACCGGGCTAGCCCCCGAAAGGGCCAGCAATCTTAACACCAACCCCCCAGGTAAAACAGCAGCTAACCCCAGGCCAACTACAGCACCCCCCCCCAGAACCCcaagata                                                   >  4:250927/1‑144 (MQ=255)
                                                                                                      gcTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGGAAGGGCCAGCACTCGTACGCCCCACGCCCCAGGGAAAGCGGCTGCGCGCGCCAGGCCGCGTCCGGCACCCCCGCCCGGAACCCCCCCCCCGGcc                                               >  6:238022/1‑149 (MQ=255)
                                                                                                          ccTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCgg                                           >  2:18616/1‑149 (MQ=255)
                                                                                                               ccAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGc                                      <  5:79293/149‑1 (MQ=255)
                                                                                                                       ggTTGCCTGTACCGGCCCTTTAGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCGTACGGCCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCCGGCCGCGTCCTGCCCCCCCCCCCTGCACCCCACCACCCTCCCCCGCAAGCACAGCGcc                              >  6:61570/1‑149 (MQ=255)
                                                                                                                                      cccTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTgcgccagcgcc               >  7:80212/1‑149 (MQ=255)
                                                                                                                                         tttCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCgccagcgccagc            >  1:284793/1‑149 (MQ=255)
                                                                                                                                         tttCGCGGGCGAGCCCGCGAAAGGGCCAGCAGTCTTACGACCAATGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCgccagcgccagc            >  6:315965/1‑149 (MQ=255)
                                                                                                                                          ttCGCGGGCGAGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGCACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCGccagcgccagcg           >  1:257325/1‑149 (MQ=255)
                                                                                                                                            cgcgGGCGCGCCCGCTAATGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCGCCAGCGCCAgcggc         <  4:16893/149‑1 (MQ=255)
                                                                                                                                                   tCGCCCGCGAAAGGGCCAGCAGTCTTACGCCCAACCCCCCAGGTAAAGCTGCAGCACTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGACCCCGGCATGCACTGCGCCACTTGCGCCAGCGCCAGCGGCGTCAgcg  >  1:37239/1‑149 (MQ=255)
                                                                                                                                                   tCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCCCCGTCCCCGGCATGCACTGCGCCAGTTGCGCCAGCGCCAGCGGCGTCAgcg  >  6:196041/1‑149 (MQ=255)
                                                                                                                                                   tCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCGCCAGCGCCAGCGGCGTCAgcg  >  7:232024/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCAAGCAGGCCTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCGCCAGCGCCAGCGGCGTCAGCG  >  NC_002947_CJ‑RC/5196443‑5196738

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

GATK/CNVnator alignment

N/A