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
A4 F60 I1 R1 17 33.8 1212966 98.8% 1198410 132.7

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,599,160 C→T L94L (CTG→CTA opgB ← phosphoglycerol transferase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,599,1600CT96.7% 94.0 / NA 30L94L (CTG→CTAopgBphosphoglycerol transferase I
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (13/16);  total (13/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGATTTCACCAGTTCCGTTATCTGACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAG  >  NC_000913/4599031‑4599287
                                                                                                                                 |                                                                                                                               
ggATTTCACCAGTTCCGTTATCTGACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATccagcctagc                                                                                                                              >  1:84988/1‑133 (MQ=255)
ggATTTCACCAGTTCCGTTATCTGACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATccagcctagc                                                                                                                              <  1:364260/133‑1 (MQ=255)
           gTTCCGTTATCTGACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCTAGCGCACCGAacac                                                                                                                   >  1:21173/1‑133 (MQ=255)
            ttCCGTTATCTGACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCTAGCGCACCGAACAcc                                                                                                                  >  2:114343/1‑133 (MQ=255)
                  tATCTGACTAAATGCAGGGCTGGCGTCCGCAGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCTAGCGCACCGAACACCGCTGTc                                                                                                            <  1:273144/133‑1 (MQ=255)
                  tATCTGACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCTAGCGCACCGAACACCGCTGTc                                                                                                            <  1:386394/133‑1 (MQ=255)
                     cTGACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCTAGCGCACCGAACACCGCTGTCAGc                                                                                                         >  2:460129/1‑133 (MQ=255)
                      tGACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACCGCGCAGGAGCCGACCTGGCGCACCGAACACTGCTGTCAGcc                                                                                                        >  2:195955/1‑133 (MQ=255)
                       gACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCTAGCGCACCGAACACCGCTGTCAGccc                                                                                                       <  2:84988/133‑1 (MQ=255)
                              tGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCTGGCGGACCAAACACCGCAGTCAACCCCAGcaca                                                                                                >  2:202953/1‑133 (MQ=255)
                                             cACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCTAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGc                                                                                 <  1:114343/133‑1 (MQ=255)
                                               ccGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCTAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCaa                                                                               <  2:21173/133‑1 (MQ=255)
                                                  aaCCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCTAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAAt                                                                            >  1:29226/1‑133 (MQ=255)
                                                      cAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCTAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTAt                                                                        <  2:14031/133‑1 (MQ=255)
                                                       aGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCTAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCaa                                                                               <  1:347362/125‑1 (MQ=255)
                                                       aGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCTAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCaa                                                                               >  2:347362/1‑125 (MQ=255)
                                                           ccAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCTAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCt                                                                   <  2:29226/133‑1 (MQ=255)
                                                                   agCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCTAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGGTACCCGGAAAAATGTATTTGCTGACGCCAg                                                           >  1:290984/1‑133 (MQ=255)
                                                                       ccAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCTAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCAcc                                                       <  1:116598/133‑1 (MQ=255)
                                                                               cTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCTAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAGTGTATTTGCTGGCGGCAGCACAGGTCAGGc                                               >  1:172824/1‑133 (MQ=255)
                                                                                 gTAACCAAAATGGTGGGGGTGAGGGCGACGGCGGCGCAGGATCCAGCCTAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGTCAGGCTg                                             <  1:202953/133‑1 (MQ=255)
                                                                                 gTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCTAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGTCAGGCTg                                             <  2:122350/133‑1 (MQ=255)
                                                                                     ccAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCTAGCGCACCGCACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGACAAAATGTATTTGGTGGCGCCAGCAGCGGTTAGGCTGTTgg                                         >  2:164660/1‑133 (MQ=255)
                                                                                     ccAAAATGGTGCGGACGATGGCCACGACGGGGCAGGATCCAGCCTAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGTCAGGCTGTTgg                                         <  1:195955/133‑1 (MQ=255)
                                                                                      cAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCTAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGTCAGGCTGTTGGt                                        <  1:507107/133‑1 (MQ=255)
                                                                                       aaaaTGGTGCGGATGCTGGCGAAGACGGCGCAGGATCCAGCCTAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGTCAGGCTGTTGGtt                                       <  1:229719/133‑1 (MQ=255)
                                                                                       aaaaTGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCTAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGTCAGGCTGTTGGtt                                       <  2:469292/133‑1 (MQ=255)
                                                                                        aaaTGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCTAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGTCAGGCTGTTGGTTa                                      >  2:566750/1‑133 (MQ=255)
                                                                                        aaaTGGTGCGGATGAGGGCGACGACGGCTCAGTATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGTCAGGCTGTTGGTTa                                      <  2:242164/133‑1 (MQ=255)
                                                                                             gTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCTAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGTCAGGCTGTTGGTTAAGGta                                 >  1:298969/1‑133 (MQ=255)
                                                                                              tGCGGATGATGGCGACGACGGCGCAGGATCCAGCCTAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGTCAGGCTGTTGGTTAAGGtat                                >  1:87172/1‑133 (MQ=255)
                                                                                                         gCGACGACAGCGCAGCATGCAGTCTAGCGCTCCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGTCAGGCTGTTGGTTAAGGTATAGAGAACcgcg                     <  1:154589/133‑1 (MQ=255)
                                                                                                                            cagcctagcGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAg  <  2:87172/133‑1 (MQ=255)
                                                                                                                                 |                                                                                                                               
GGATTTCACCAGTTCCGTTATCTGACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAG  >  NC_000913/4599031‑4599287

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

GATK/CNVnator alignment

N/A