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
A6 F31 I1 R1 2 18.6 702272 97.0% 681203 135.2

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 4,173,025 A→T Q618L (CAG→CTG)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,173,0250AT87.0% 71.8 / 1.2 23Q618L (CAG→CTG) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base T (15/5);  total (16/7)
Fisher's exact test for biased strand distribution p-value = 2.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATC  >  CP009273/4172912‑4173139
                                                                                                                 |                                                                                                                  
gACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGaagaag                                                                                          >  1:142302/1‑140 (MQ=255)
 aCTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCGAACTTGGATGAAGAAgg                                                                                         >  2:307326/1‑140 (MQ=255)
 aCTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAgg                                                                                         >  2:46266/1‑140 (MQ=255)
 aCTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAgg                                                                                         >  2:33022/1‑140 (MQ=255)
  cTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGc                                                                                        <  2:48814/140‑1 (MQ=255)
           aCGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGt                                                                                <  1:156467/139‑1 (MQ=255)
           aCGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGt                                                                                >  2:156467/1‑139 (MQ=255)
                         aCTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTgg                                                                                                                  <  1:12500/91‑1 (MQ=255)
                         aCTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTgg                                                                                                                  >  2:12500/1‑91 (MQ=255)
                         aCTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAAc                                                                 >  1:98581/1‑140 (MQ=255)
                         aCTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAAc                                                                 >  2:145894/1‑140 (MQ=255)
                         aCTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAAc                                                                 >  2:181088/1‑140 (MQ=255)
                         aCTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAAc                                                                 >  2:21307/1‑140 (MQ=255)
                         aCTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAAc                                                                 >  2:238350/1‑140 (MQ=255)
                         aCTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAAc                                                                 <  2:307613/140‑1 (MQ=255)
                             cGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGc                                                             >  1:114363/1‑140 (MQ=255)
                                                   tgttgtAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCt                                       >  1:147665/1‑140 (MQ=255)
                                                   tgttgtAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCt                                       >  1:147609/1‑140 (MQ=255)
                                                      tgtAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGt                                    <  2:142302/140‑1 (MQ=255)
                                                             gACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCcg                             >  1:275437/1‑140 (MQ=255)
                                                                gAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGa                          <  1:181088/140‑1 (MQ=255)
                                                                          aCCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGAc                >  2:190912/1‑140 (MQ=255)
                                                                                        gaagaaGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATc  <  1:190912/140‑1 (MQ=255)
                                                                                                                 |                                                                                                                  
GACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATC  >  CP009273/4172912‑4173139

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

GATK/CNVnator alignment

N/A