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
A1 F113 I2 R1 11 53.1 1830242 95.2% 1742390 142.8

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,184,542 C→A P1100Q (CCG→CAG)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,184,5420CA95.0% 112.9 / ‑5.8 42P1100Q (CCG→CAG) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (21/17);  total (25/17)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

GGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAA  >  NC_000913/4184405‑4184672
                                                                                                                                         |                                                                                                                                  
ggcccttACACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGCTAACAAGGGTGTAATTTCTAAGACCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTgg                                                                                                                               <  2:217384/138‑1 (MQ=255)
ggCGGGTAAGCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTgg                                                                                                                               <  1:832475/143‑1 (MQ=255)
   ggTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCCGTAGACACCGTACAGAACCGACGGGGCg                                                                                                                            >  1:455132/1‑143 (MQ=255)
        aaCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAAGCGGCGGGGCGtcgcg                                                                                                                       >  2:584368/1‑139 (MQ=255)
        aaCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGCAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTAATGAACGAGCTGGGCGTGCCg                                                                                                                       >  2:106443/1‑143 (MQ=255)
          cGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGcgcggaacgcc                                                                                                                     >  1:516589/1‑136 (MQ=255)
                            gACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACAGCGGTCa                                                                                                   >  2:395035/1‑143 (MQ=255)
                             aCAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAg                                                                                                  >  2:178088/1‑143 (MQ=255)
                                 gATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTCTGTACATCGTTCTGCTcc                                                                                              >  2:14545/1‑143 (MQ=255)
                                  attgCAGGTCGTCGCGGTAACAAGGGTGTAAGTTCTAAGATCACCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCt                                                                                             <  2:346898/140‑1 (MQ=255)
                                       aGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGCAAACGGTACGCCGGGAGACATCGTACTGAACCCGCTGGGCGTCCCGTCTCGTCTGACCACCGGACCGCTCCCCGaca                                                                                        >  1:232009/1‑141 (MQ=255)
                                               aCGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGCTCCTCGAAAccccccgg                                                                                >  2:883406/1‑141 (MQ=255)
                                                 ggTAACATGGGTGTAATTTCGAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTggg                                                                              <  1:474648/143‑1 (MQ=255)
                                                 ggTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTggg                                                                              <  2:912410/143‑1 (MQ=255)
                                                  gcaaCAAGGTTGTAATTTCTACGATCAACCCGATCGAAGATATGCCTTACGATGAAGACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCTGTTAGATCCTCGAAACCCACCTGGGt                                                                             <  1:903209/141‑1 (MQ=255)
                                                       aaGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCCGGGTCTGGc                                                                        >  1:633025/1‑143 (MQ=255)
                                                              tAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                                 <  1:826333/143‑1 (MQ=255)
                                                                aTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGCCCCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTCCGACAgg                                                               >  1:684894/1‑143 (MQ=255)
                                                                    ccaagcacAACCCGATCGAAGATAGGCCTTACGATGAAAACGGTACGCCGGTAGTCTGCCTACTATACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATc                                                           <  2:666229/136‑1 (MQ=255)
                                                                          tCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGAc                                                     >  2:191099/1‑143 (MQ=255)
                                                                            aaccgcaTCGACGAGATGCCGTACGATGGAAACGGTACGGCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACaa                                                   <  1:294856/137‑1 (MQ=255)
                                                                            aaCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACAGCGGTCAGATCCTCGCAACCCACCTGCGTATGGCTGCGAAAGGTATCGGCGACaa                                                   >  1:479079/1‑143 (MQ=255)
                                                                              cccGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGa                                                 >  2:737210/1‑143 (MQ=255)
                                                                               ccGATCGAAGAGATGCCTTACGATGAAAGCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCGCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGAt                                                <  1:106904/143‑1 (MQ=255)
                                                                                   tCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCGCACTGGGTATGGCTGCGAAAGGTAGCGGCGCCAAGATCCAc                                            >  1:662548/1‑143 (MQ=255)
                                                                                    cGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGCCACCCCCCTGGGTATGGCTGCGAAAGGTATCGGCCACAAGATCAACg                                           >  2:271254/1‑143 (MQ=255)
                                                                                     gAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACACGATAGACGc                                          >  2:242278/1‑143 (MQ=255)
                                                                                        gATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGAcgcgctccaccccat                                       >  2:841754/1‑129 (MQ=255)
                                                                                             gCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGTTCAGATCCTCGAAACCCACCTGGGGATGGCTGCGCAAGGTCTCCGCGACCAGATCAAAGCCATGCTGa                                  >  1:518076/1‑143 (MQ=255)
                                                                                             gCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGa                                  >  2:474983/1‑143 (MQ=255)
                                                                                               gttACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGaaa                                <  2:326553/142‑1 (MQ=255)
                                                                                                ttACGATGAAACCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAc                               <  2:598173/143‑1 (MQ=255)
                                                                                                  aCGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGACACCCACCTGGGTATGGCTGCGGCAGGTATCGGCGCCCAGGTCACCGCCATGCTGAAAcag                             >  2:37662/1‑143 (MQ=255)
                                                                                                  aCGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcag                             >  2:321966/1‑143 (MQ=255)
                                                                                                    gATGAAAACGGTCCGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTATCGTATGAACATCGGTCAGATCCTCGAAACCCACCGGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagca                           <  2:908775/143‑1 (MQ=255)
                                                                                                    gATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCCCCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGCTGAACGCAATGCTGAAAcagca                           >  2:225848/1‑143 (MQ=255)
                                                                                                      tGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACagcagc                         <  2:490919/143‑1 (MQ=255)
                                                                                                        accaCGGTCCGCCGGTAGACATCGTACTGAACCAGCGGGGCTTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaa                       <  2:120606/140‑1 (MQ=255)
                                                                                                              gTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTc                 >  2:540931/1‑143 (MQ=255)
                                                                                                                  gCCGGTAGACATCGCACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGa             <  2:521128/143‑1 (MQ=255)
                                                                                                                            aTCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGCACATCGGTCAGATCCGCGAGACCCACCCGGCTATGCCTGCGACAGGTACCGGCGACCACATCAACCCCATGCCGCAGCAGCAGCCACCCGCCGCGAAAcaccccca   >  1:95938/1‑136 (MQ=255)
                                                                                                                             tCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGaa  <  2:893225/143‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                  
GGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAA  >  NC_000913/4184405‑4184672

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑

GATK/CNVnator alignment

N/A