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
A20 F53 I0 R4 206 63.7 2487030 95.4% 2372626 127.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,725,205:1 +T 100% intergenic (‑37/‑492) 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,2051.T85.0% 82.0 / 14.3 40intergenic (‑37/‑492)U712_19205/U712_19210Uncharacterized protein ywcI/Minor extracellular protease vpr
Reads supporting (aligned to +/- strand):  ref base . (2/4);  new base T (24/5);  total (31/9)
Fisher's exact test for biased strand distribution p-value = 2.73e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

TGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAA  >  CP006881/3725068‑3725315
                                                                                                                                           |                                                                                                                 
tGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCt                                                                                                                  >  2:21001/1‑139 (MQ=255)
 gAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCGTTTTGGCAGgccgccgt                                                                                                                  >  5:236256/1‑137 (MQ=255)
  aGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTATCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCttt                                                                                                               >  7:30215/1‑139 (MQ=255)
            accgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTTGCAGccc                                                                                                    >  3:307744/4‑139 (MQ=255)
            accgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTGTTTAGCAGccc                                                                                                    >  3:265042/4‑139 (MQ=255)
            accgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCGTTTTGGCAGGCAGCCTTTTTAGCATccc                                                                                                    >  8:246969/4‑139 (MQ=255)
              cgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCCGCCCGt                                                                                                  >  6:108247/2‑139 (MQ=255)
              cgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGt                                                                                                  >  1:10077/2‑139 (MQ=255)
              cgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGt                                                                                                  >  4:230315/2‑139 (MQ=255)
              cgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTGGCAGGCAGCCGTTTTAGCAGCCCGt                                                                                                  >  3:291795/2‑139 (MQ=255)
              cgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTAGGGCAGGCAGCCTTTTATGCAGCCCGt                                                                                                  >  4:297440/2‑139 (MQ=255)
                 ggCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGCCCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTGGCAGGCAGCCGTTTTAGCAGCCCGtttt                                                                                               >  1:6945/1‑139 (MQ=255)
                 ggCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGgttt                                                                                               >  5:262056/1‑139 (MQ=255)
                      ccGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTGTTTTTTGGAGGCAGCCGATTTAGCATCGCGGTTTATCAg                                                                                          >  3:4654/1‑139 (MQ=255)
                                    tCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGcccgc                                                                               >  8:178287/1‑134 (MQ=255)
                                          gTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGa                                                                      <  7:178287/139‑1 (MQ=255)
                                                   ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGcccg                                                                                <  2:36034/120‑2 (MQ=255)
                                                   ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGcccg                                                                                >  1:36034/1‑119 (MQ=255)
                                                   ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGGCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGGGTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtg                                                             >  3:289871/1‑138 (MQ=255)
                                                             aGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCCGGCAGCCTTTTTAGCCGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCTTTTCTCg                                                   >  7:201611/1‑139 (MQ=255)
                                                             aGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGCCGGCCGATGTCTTTTCTCg                                                   >  8:141885/1‑139 (MQ=255)
                                                                             tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGCCGGCCGATTTCTTTTATCGCAAACTAGCGGGTCTa                                   >  4:65951/1‑139 (MQ=255)
                                                                             tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCTTTTCTCGCAAACTAGCGGctc                                     >  1:165876/1‑137 (MQ=255)
                                                                                aGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGAGTCTTTTTGGCTGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGAAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa                                <  8:67399/139‑1 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCCGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGCGTTGATTTGTCGGCAACTAGCCGGGGTAGCaaaa                             >  1:132172/1‑138 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCCGCCGTTTTTGCAGCCCGTTTTCTCAGCCGCAGCCCGCCAGAAAAGACGGCCGATTTCGTTTGTCGCAAACTTGCGGCTCTTGCaaaa                             >  3:18788/1‑138 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTACCGGCTCGAGCAAAAc                            >  8:83991/1‑139 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTAGCAAACTAGCGGCTCTAGCaaaa                             >  8:242659/1‑138 (MQ=255)
                                                                                      aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGACAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCGCTAGCAAAACt                           >  7:134659/1‑138 (MQ=255)
                                                                                      aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACt                           >  1:76065/1‑138 (MQ=255)
                                                                                      aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTATCGCAAACTAGCGGCTCTAGCAAAACt                           >  7:173887/1‑138 (MQ=255)
                                                                                       aGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTTGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGGCGCACCCCGCAAGAAAAGACGGCCGATTTATTTTGTAGCAACCAAGCGGCTGGAACAAAAcaaa                         >  6:3270/1‑136 (MQ=255)
                                                                                       aGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGCCGGCCGCTTTGTTTTCTCGCAAACGATCGGAGCTAGCaaaagtga                         >  6:103842/1‑135 (MQ=255)
                                                                                        ccAAAGGTTTCATACAAAAACACCCCACCGCTTCTTTTTGCCCGCCCCCTTTTTACCCCGCCCGTTTACTCCGCCCCCCCCCCCAATAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                        <  3:230315/138‑1 (MQ=255)
                                                                                        caaaccGTTCCATACAAAAACAGCCCACCGCTTCTTTTGGGCAGCCCCCTTTTTTACCCCCCGGTTTCCTCCCCCGCACCCCCAAACAAACCACGGCCGATTTCTTTTCTCGAAAACTAGCGGCTCTAGCAAAACTCAt                        <  6:262056/134‑1 (MQ=255)
                                                                                          aaaCTTTTAATATAAAAAAATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCCGCCCTTTTTCTCAGCCGCAGCCCGAAATAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGa                      <  7:246969/139‑1 (MQ=255)
                                                                                          aaaCGTTTCATAAAAAAACATCCCTCCGCTTTTTTTTGGCTGGCAGCCTTTTTAGCCGCCGTTTTTGTCCGCCGCAGCCCGCAATAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGa                      <  3:145451/139‑1 (MQ=255)
                                                                                            ccGTTACAAACAAAAACACCCCTCCGCTTCTTTTTGCCCGCCCCCCTTTTAACCCGCCCGTTTTCCCCCCCCCCCCCCCAAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                    <  2:10077/138‑1 (MQ=255)
                                                                                            acttttAATACAAAAACATCCCACCGCTTCTTTTTGGCTCGCAGCGTTTTTAACAGCCGGTTTTCTCAGCCCCCGCCCGCAAGAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                    <  7:192021/136‑1 (MQ=255)
                                                                                            aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATGTCTTTTCTCGCAAACGTCCGGCTCAAGCAAAACTCCTGAct                    >  3:56788/1‑139 (MQ=255)
                                                                                            aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCATTTTAGCAGCCCGTTTTCTAAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTATTCTCGCAAACTAGCGGCTAAAGAAAAACACATGCct                    >  4:253740/1‑139 (MQ=255)
                                                                                                ttCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTATATTCTCGCAAACTAGCGGCTCTAGCAAAACTCAAGACTCAAt                >  7:200463/1‑139 (MQ=255)
                                                                                                               cccTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCATTTATCGCAAACAAGCCGCTCGAGCAAAACTCATTCATCTATTATAGACCGTggga  >  5:189382/1‑136 (MQ=255)
                                                                                                                                           |                                                                                                                 
TGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAA  >  CP006881/3725068‑3725315

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

GATK/CNVnator alignment

N/A