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
A18 F53 I0 R2 221 66.0 2680176 94.9% 2543487 126.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,725,232 A→C 30.4% intergenic (‑64/‑465) U712_19205 ← / → U712_19210 Uncharacterized protein ywcI/Minor extracellular protease vpr

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,725,2320AC30.4% 21.8 / 5.2 29intergenic (‑64/‑465)U712_19205/U712_19210Uncharacterized protein ywcI/Minor extracellular protease vpr
Reads supporting (aligned to +/- strand):  ref base A (15/5);  new base C (3/6);  total (18/11)
Fisher's exact test for biased strand distribution p-value = 4.78e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGAAG  >  CP006881/3725113‑3725366
                                                                                                                          |                                                                                                                                       
aCCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCCGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACgg                                                                                                                      >  6:74602/1‑139 (MQ=255)
 cccAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGc                                                                                                                     >  7:157968/1‑139 (MQ=255)
      ttCATCAAAAAGACCATCCACCCCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTGGCAGGCCGCCGTTTTAGCCGCCCGGTTTCTCAGCCGCCGCCCGCAAGCAAAGACGGCCGAtt                                                                                                                >  6:72144/1‑139 (MQ=255)
      ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCCGCCTTTTTATCAGCCCGTTTTCTCAGCCCCACCCCGCAAGAAAAGACGGCCGAtt                                                                                                                >  1:245536/1‑139 (MQ=255)
      ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCCGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtt                                                                                                                >  8:141976/1‑139 (MQ=255)
           caAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtttcttt                                                                                                           <  3:258373/139‑1 (MQ=255)
           caAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGACCGTTTTCTCAGCCGCCGCCCCCAAGAAAAGACGGCCGAtttcgtt                                                                                                           >  5:244875/1‑139 (MQ=255)
           caAAAAGAACATCCACCCCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTCGCAGGCAGCCGTTTTAGCAGCCCGTTTTCTCAGCCGCAGACCGCAAGAAAATACGACCGAtatattt                                                                                                           >  7:42767/1‑139 (MQ=255)
                           acCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCgg                                                                                           >  2:135/1‑139 (MQ=255)
                             ccGAAAAGATACACGCAAACTTTTAATACAAAAACAGCCCTCCGTTTTTTTTTGGCCGGCCTCCTTTTTACCCGCCCGTTTTCCCCGCCGCAGCCCCCAACAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGct                                                                                         <  6:244875/139‑1 (MQ=255)
                                tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCGTTTTAGCAGCCCGTTTTCTCATCCGCAGCCCGCAAGAAAAGACGGCCGATTTGTTTTTTAGCCTACTATCggggc                                                                                        >  3:8176/1‑134 (MQ=255)
                                   aGATAAAACCAAACTTTACATACAAAAACTTCCCTCGGTTTTTTTTTGGCTGGCTGCTTTTTTAGCAGCCCGTTTTCCCCGCCGCAGCCCGAAATAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa                                                                                   <  4:293952/139‑1 (MQ=255)
                                     aTACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCaa                                                                                  >  3:136817/1‑138 (MQ=255)
                                       aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTGTTTTTGGCAGGCCGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCCGCCCGCAAGAAAAGACGGCCGAGTTCGTTTGTCGCCAACTAGCGGCTCGAGGaaaa                                                                                >  2:134569/1‑138 (MQ=255)
                                       aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTATCAGCCCGTTTTCTCCGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTATCGGCTCTAGCAAAAc                                                                               >  1:42814/1‑139 (MQ=255)
                                       aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCaaaa                                                                                >  2:176775/1‑138 (MQ=255)
                                           gaaaaCTTTTAATAAAAAAACTCCCCACGGATTCTTTTTGGCCGGCCGCTTTTTTCCCCGCCGGTTTGCCCCGCCCCCGCCGCAAAGAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                                                           <  4:320349/137‑1 (MQ=255)
                                           ccAAACGTGCCATCCAAAACCATCCCTCCGCTGCTTTTTGGCCGGCAGCCTTTTTAGCCGCCGGTGTTCTCCGCCGCCCCCCCCAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                                                           <  8:188762/138‑1 (MQ=255)
                                            cAAACGTTTCAAACAAAAACACCCCACCGCGTCTTTTTGGCAGGCACCCTTTTTACCCCCCCGGTTTCCCCCCCCCCCCCCCCAACACAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                                                           <  7:79951/139‑1 (MQ=255)
                                             aaccttttCATATAAAAACCTCGCGCCGTTTTTTTTGGGCCGGCCTCATTTTTCCCAGCCCTTTTTCCCCCCCGCCCCCCGAAAAAAACCACGGCCGTTTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGa                                                                         <  8:136109/134‑1 (MQ=255)
                                               ccGTTTCATACAAAAACATCCCACCGTTTCTTTTTGGCTGGCAGCCTTTTTAGCAGCCCGTTTTCTCACCCGCCGCCCGCAATAAAAGACGGCCGATTTCTTTTCTCGAAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                       <  7:21023/138‑1 (MQ=255)
                                               aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                       >  1:279150/1‑139 (MQ=255)
                                               aCGTTTCATACAAAAACACCCCTCCGCTTCTTTTGGGCTGGCCCCATTTTACGCCGCCGTTTTTCTCCGCCGCCGCCCCCAAAAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                       <  7:141976/139‑1 (MQ=255)
                                                         gAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCa                                                                                                    <  2:223499/100‑1 (MQ=255)
                                                         gAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCa                                                                                                    >  1:223499/1‑100 (MQ=255)
                                                               caccccTCCGCTTCTTTTTGGCCGGCCGCCTTTTTAGCAGCCGGTTTTCTCCCCCCCAGCCCCCAATAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTgg                                                       <  1:135/136‑1 (MQ=255)
                                                                                       gcagCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAatat                               >  1:247157/1‑139 (MQ=255)
                                                                                                         cccGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACACATGACTATATGATAGACCGTGGAAAGAGAATTCATTGTCACAATAAGAGTGAAAAAACCaa                >  5:144876/1‑138 (MQ=255)
                                                                                                                      ccGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTATCGCAAACTAGCGGCTCTAGCAAAACTCCTGACTCTATGATAGACCGTGGAAAGAGAAATCATTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGaac  >  2:3393/1‑138 (MQ=255)
                                                                                                                          |                                                                                                                                       
ACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGAAG  >  CP006881/3725113‑3725366

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

GATK/CNVnator alignment

N/A