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
A4 F13 I0 R1 140 52.5 2164708 98.8% 2138731 146.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947_CJ‑RC 5,258,571 C→G 24.8% intergenic (+20/+25) PP_5703 → / ← PP_4634 type III effector protein HopJ/alpha/beta family hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947_CJ‑RC5,258,5710CG24.8% 47.3 / 2.7 27intergenic (+20/+25)PP_5703/PP_4634type III effector protein HopJ/alpha/beta family hydrolase
Reads supporting (aligned to +/- strand):  ref base C (5/14);  new base G (3/4);  total (9/18)
Fisher's exact test for biased strand distribution p-value = 6.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGCACTACCGTGATGTGTTGGCCACGCCTGAGGGTAGCGACCATGGCAACATCCGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCT  >  NC_002947_CJ‑RC/5258429‑5258717
                                                                                                                                               |                                                                                                                                                  
gAGCACTACCGTGATGTGTTGGCCACGCCTGAGGGTAGCGACCATGGCAACATCCGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGcc                                                                                                                                                <  4:192469/147‑1 (MQ=255)
 aGCACTACCGTGATGTGTTGGCCACGCCTGTGGGTAGCGACGATTGCAACATCCGTGCACTGATCTAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCTAGGCTTGAGATTGCCGGGGCCGCTTTGCGGccc                                                                                                                                               <  2:147868/147‑1 (MQ=255)
 aGCACTACCGTGATGTGTTGGCCACGCCTGAGGGTAGCGACCATGGCAACATCCGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGccc                                                                                                                                               <  8:190694/147‑1 (MQ=255)
  gCACTACCGTGATGTGTTGGCCACGCCTGAGGGTAGCGACCATGGCAACATCCGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGcccc                                                                                                                                              <  4:130257/147‑1 (MQ=255)
  gCACTACCGTGATGTGTCGGCCACGCCTGAGGGTAGCGACCATGGCAACATCCGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGcccc                                                                                                                                              <  1:222703/147‑1 (MQ=255)
       aCCGTGATGTGTTGGCCACGCCTGAGGGTAGCGACCATGGCAACATCCGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGGCGCGTTTGGGGcccaccc                                                                                                                                           >  6:17912/1‑142 (MQ=255)
                             tGAGGGTAGCGACCATGGCAACATCCGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGGCGCGTTTGGGGCCCCCGCTCCGGCGCCCCCCCACCAACCat                                                                                                                   >  6:43483/1‑147 (MQ=255)
                                      cGACCATGGCAACATCCGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGACTTGAGATTGCCGGGGCCGCTTTGAGGCCCCCCGAACCTGCGCCCACCCAGAATGCATGATGGTg                                                                                                            >  8:267733/1‑145 (MQ=255)
                                                  cATCCGTGCACTGATCAAACACGGTCTGGAGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTTGCGCCCCACGATCCTGCCCCCACCCGAAAAACATTATGGTGCCCGGCcccccccc                                                                                              >  6:262417/1‑147 (MQ=255)
                                                                          tCTGGCGGGCTTTAAGTTTGCTGCGCAGCCGATTGCGCGCATGGTTTGTTTTTGCGGGGGCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCa                                                                      <  8:49458/147‑1 (MQ=255)
                                                                               ggggcgTTATTTTTGCGGCGCCGCGGTTGGCGCGAATGGCTTTATATGGCCGGGGCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCa                                                                 <  1:150264/146‑1 (MQ=255)
                                                                                        tttttGCGGCCCCGCCGCTGGCGCTAATGGTTTTATATGGCGGGGGCCGTTTGGCGCCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCAcc                                                         <  5:231233/144‑1 (MQ=255)
                                                                                         gTTTGCTGCGCCGCCGCTGCCGCTAATGGCTTGAGATTGCCGGGGCCGGTTTGGGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCgg                                                       <  7:157506/147‑1 (MQ=255)
                                                                                                                 cAAGGCTTGAGATTGCCGGGGCCCCGATTCGGCCCCACCATACAGCGCCCCCCCGCAAAACATGATGGTGCCAGGCCCCCCCCCATACCGCCACCCCAGCCGCAACCACAGCAAAGCCCCCGGCCACACGCATTCGCGaaaaccgg                                >  2:162076/1‑141 (MQ=255)
                                                                                                                  aaGGCTTGAGATTGCCGGGGCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCa                              <  2:3009/147‑1 (MQ=255)
                                                                                                                    ggCTTGTGATTGCGGGGGCCGCTTTGGGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATg                            <  2:176823/147‑1 (MQ=255)
                                                                                                                           gaTGGCGGGGGCCGATTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGaaaaa                      <  7:72828/146‑1 (MQ=255)
                                                                                                                            atggCCGGGCCCGCTTGGGGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAAcg                    <  7:65828/144‑1 (MQ=255)
                                                                                                                              ggCCGGGGCCGCTTGGCGCCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAAcgc                   <  4:207259/145‑1 (MQ=255)
                                                                                                                                  ggggCCGGTTGGGGCCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAAt              <  5:132349/147‑1 (MQ=255)
                                                                                                                                  gggcCCCCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAAt              <  6:259295/147‑1 (MQ=255)
                                                                                                                                   gggCCGCTGTGCGGCCCCACGATCCTGCTCCCACCCTCAATCCATGCTGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATg             <  4:183444/147‑1 (MQ=255)
                                                                                                                                    ggCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCACCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATgg            >  1:164530/1‑147 (MQ=255)
                                                                                                                                    gcccGCTTGGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATgg            <  6:228836/145‑1 (MQ=255)
                                                                                                                                         cTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCCGCCCCCCCTGCCGCAACCACTGCAAGGCCACCGGCCTCAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAAt       >  3:96440/1‑147 (MQ=255)
                                                                                                                                         cTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAAt       >  6:236236/1‑147 (MQ=255)
                                                                                                                                             tGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTcc   >  2:33200/1‑147 (MQ=255)
                                                                                                                                               cGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCt  >  4:203789/1‑147 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GAGCACTACCGTGATGTGTTGGCCACGCCTGAGGGTAGCGACCATGGCAACATCCGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCT  >  NC_002947_CJ‑RC/5258429‑5258717

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