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
A13 F56 I0 R1 198 55.3 2323794 95.3% 2214575 126.9

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.4% 109.7 / 11.4 41intergenic (‑37/‑492)U712_19205/U712_19210Uncharacterized protein ywcI/Minor extracellular protease vpr
Reads supporting (aligned to +/- strand):  ref base . (1/5);  new base T (21/11);  total (25/16)
Fisher's exact test for biased strand distribution p-value = 6.47e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.45e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

GAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAA  >  CP006881/3725069‑3725321
                                                                                                                                          |                                                                                                                     
gAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCtt                                                                                                                    >  3:119013/1‑139 (MQ=255)
             cgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGt                                                                                                        >  7:21807/2‑139 (MQ=255)
             cgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGt                                                                                                        >  7:169870/2‑139 (MQ=255)
             cgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTGGCAGGCAGCCGTTTTAGCAGCCCgg                                                                                                        >  1:254570/2‑138 (MQ=255)
                ggCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCCCACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCGTTTTAGCAGCCCGgtt                                                                                                      >  1:71154/1‑138 (MQ=255)
                ggCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCCGCCTTTTTAGCAGCCCGttt                                                                                                      >  5:200318/1‑138 (MQ=255)
                ggCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGtttt                                                                                                     >  4:179971/1‑139 (MQ=255)
                ggCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCGTTTTAGCAGCCCGttt                                                                                                      >  5:196234/1‑138 (MQ=255)
                ggCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGAGACAAGCACACGTTTCATAGATAAACATCCAGCCGCGTCGGTTTGGCAGGCAGCCTGTTTTGCAGCCCGtttg                                                                                                     >  8:213519/1‑138 (MQ=255)
                  cAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTtct                                                                                                   <  3:179971/139‑1 (MQ=255)
                                cTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGcccg                                                                                      <  7:236836/138‑2 (MQ=255)
                                                  ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCGGAtt                                                                   >  8:78097/1‑139 (MQ=255)
                                                    catcaAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtttc                                                                 >  4:81683/1‑139 (MQ=255)
                                                                 aTCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCaaa                                                     >  8:187875/1‑138 (MQ=255)
                                                                  tCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCCAAc                                                    >  4:240591/1‑138 (MQ=255)
                                                                       acCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTtct                                                                                                   >  1:225813/1‑86 (MQ=255)
                                                                       acCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTtct                                                                                                   <  2:225813/86‑1 (MQ=255)
                                                                            tAAAGATACAAGCAAACGTTTCATAGAAAAAGATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCGGCCCGTTTTCTCAGCCGCAGCCCGCAACAAACGAAGGCCGATTTCTTTTATCGCAAACTACCGGctc                                           >  7:140166/1‑137 (MQ=255)
                                                                            tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTGGCAGCCCGGTTTCTCAGCCGCCGCCCGCAAGAAAAGACGGCCGAGTTCTTTTCTCGCAAACTAGCGGctc                                           >  8:118372/1‑137 (MQ=255)
                                                                            tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGctc                                           >  7:173841/1‑137 (MQ=255)
                                                                            tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCATTTATCGCAAACTAGCGGctg                                           >  3:31481/1‑136 (MQ=255)
                                                                            tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCCAGAAAAGACGGCCGAGTTCTTTTAGCGCCAACAATCGGcgcaa                                         >  4:14852/1‑135 (MQ=255)
                                                                              aaGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAg                                        >  3:108276/1‑138 (MQ=255)
                                                                               ctaTACAACCAACCTTTTAATACAAAAACAACCCTCCGCTTCTTTTTGGCAGGCAGCTTTTTTCGCAGCCCTTTTCCTCCGCCGCAGCCCCAAATAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa                                      <  8:173841/137‑1 (MQ=255)
                                                                               aGATACACGAAAACGTTTCATAGAAAAACACCCCTCCGCTTCTTTTTCGCCGCCAGCGTTTTTCGCAGCCCGTTTTCTCAGCCGCACCCCGCAACAAAACACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa                                      <  7:118372/139‑1 (MQ=255)
                                                                                   aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCaaaa                                   >  6:107549/1‑138 (MQ=255)
                                                                                   aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCaaaa                                   >  7:19065/1‑138 (MQ=255)
                                                                                   aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAAc                                  >  5:24653/1‑139 (MQ=255)
                                                                                   aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAAc                                  >  5:289979/1‑139 (MQ=255)
                                                                                   aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTATCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTACTTTTCTCGCAAACTACCGGCTCTAGCaaaa                                   >  8:187100/1‑138 (MQ=255)
                                                                                       gaaacggTTTCATCTAAAAACATCCCTCCGATTCTTTTTGGCTGGCAGCTTTTTTCGCAGCGCGTTTTCTCAGCCGCAGCCCGAAACAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                              <  7:187100/133‑1 (MQ=255)
                                                                                       gCAAACGTTTCATATAAAACCACCCCTCCGCTTCTTTTTGCCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCCGCCCCAGCCCGCAATAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                              <  3:81683/139‑1 (MQ=255)
                                                                                       ccAAACTTTTCATAAAAAAACATCCCACCGCTTCTTTTTGCCAGGAAGCCTTTTTCTCACCCGGTTTTCTCCGCCGCCGCCCCCAATAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                              <  3:240591/138‑1 (MQ=255)
                                                                                          aaCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGAAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAc                           <  7:187875/139‑1 (MQ=255)
                                                                                           acttttCATAAAAAAACATCCCTCCGCTTCTTTTTGGCCGGCAGCTTTTTTAGCCGCCCGTTTTCTCAGCCGCCGCCCGAAATAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                          <  4:108276/136‑1 (MQ=255)
                                                                                           aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                          <  5:107549/139‑1 (MQ=255)
                                                                                           aCGTTTAATATAAAAAAATCCATCGGCTTATTTTTGGATGGCAGCGTTTTTTGATCCCCTTTTTCTCAGCCCCAGCCCGAATTAAAAGACGGCCGATTTCTTTTCTCGAAAACTAGCGGCTCTAGCAAAACTCATGAct                          <  8:19065/139‑1 (MQ=255)
                                                                                           aCGTTACAAACAAAAACATCCCACCGCTTCTTTTTGGCAGGCACCCTTTTTAGCCCCCGGTTTGCTCAGCCCCCGCCCCCAACAAAACACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                          <  7:177750/139‑1 (MQ=255)
                                                                                                      aaaaaCATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGAc               <  6:289979/139‑1 (MQ=255)
                                                                                                                cTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa                                      <  7:72249/106‑1 (MQ=255)
                                                                                                                cTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa                                      >  8:72249/1‑106 (MQ=255)
                                                                                                                   cGCTTCTTTTTGGCCGGCCGCCTTTTTCGCAGCCCGTTTCCTCAGCCGCAGCCCCCAACAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGaa  <  6:120890/139‑1 (MQ=255)
                                                                                                                                          |                                                                                                                     
GAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAA  >  CP006881/3725069‑3725321

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