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 F2 I200 R1 216 30.8 1683622 97.2% 1636480 87.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 951,809 T→G T94P (ACC→CCC)  ycaO ← 30S ribosomal protein S12 methylthiotransferase accessory factor YcaO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273951,8090TG96.3% 83.4 / ‑3.5 27T94P (ACC→CCC) ycaO30S ribosomal protein S12 methylthiotransferase accessory factor YcaO
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (16/10);  total (17/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.18e-01

CAGCCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGTTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCGAGTGCAGAAGCCAGCGCCGCTTT  >  NZ_CP009273/951726‑951893
                                                                                   |                                                                                    
cagcCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCgtc                                                                                <  1:542370/90‑3 (MQ=255)
cagcCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGcc                                                                                <  1:239834/90‑1 (MQ=255)
 agcCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGccc                                                                               >  2:826360/1‑90 (MQ=255)
 agcCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGccc                                                                               >  1:659164/1‑90 (MQ=255)
 agcCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGccc                                                                               >  1:439956/1‑90 (MQ=255)
       tCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCccagccag                                                                         >  2:400179/1‑90 (MQ=255)
            cacaTCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTc                                                                    >  2:207790/1‑90 (MQ=255)
               aTCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGc                                                                 <  1:826360/90‑1 (MQ=255)
               aTCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTAGGCGATGGGTTCGCCCAGCCAGAAGTCCGc                                                                 >  1:370928/1‑90 (MQ=255)
                tCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCa                                                                >  2:228768/1‑90 (MQ=255)
                           ggTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAgtt                                                     <  1:603038/90‑1 (MQ=255)
                           ggTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAgtt                                                     <  1:794053/90‑1 (MQ=255)
                           ggTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAgtt                                                     <  1:400179/90‑1 (MQ=255)
                                gTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAgttggttg                                                >  2:804054/1‑90 (MQ=255)
                                     aaCCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTgagaga                                           >  2:834252/1‑90 (MQ=255)
                                                 ttGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAatat                               >  2:798133/1‑90 (MQ=255)
                                                       tAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCg                         <  1:146472/90‑1 (MQ=255)
                                                       tAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCg                         <  2:68050/90‑1 (MQ=255)
                                                             aCGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCGAGTGCa                   >  2:318040/1‑90 (MQ=255)
                                                             aCGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCGAGTGCa                   >  2:547439/1‑90 (MQ=255)
                                                                    gACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCGAGTGCAGAAGCCa            >  2:10019/1‑90 (MQ=255)
                                                                    gACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCGAGTGCAGAAGCCa            >  2:453554/1‑90 (MQ=255)
                                                                      ccGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCGAGTGCAGAAGCCAgc          >  2:228245/1‑90 (MQ=255)
                                                                        gTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCGAGTGCAGAAGCCAgcgc        >  1:530953/1‑90 (MQ=255)
                                                                          tGGCGATGGTTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCGAGTGCAGAAGCCAGcgccg      >  1:311985/1‑90 (MQ=255)
                                                                           ggCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCGAGTGCAGAAGCCAGcgccgc     <  2:401926/90‑1 (MQ=255)
                                                                              gATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCGAGTGCAGAAGCCAGCGCCGCttt  <  2:439956/90‑1 (MQ=255)
                                                                                   |                                                                                    
CAGCCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGTTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCGAGTGCAGAAGCCAGCGCCGCTTT  >  NZ_CP009273/951726‑951893

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGTTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCGAGTGCAGAAGCCAGCGCCGCTTTCTTGGTTGCGCCT  >  NZ_CP009273/951717‑951906
                                                                                            |                                                                                                 
GTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCC                                                                                            >  SRR3722079.443778/1‑100 (MQ=60)
GTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCC                                                                                            >  SRR3722079.666006/1‑100 (MQ=60)
         CAGCCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGTCCAGCCAGAAG                                                                                   <  SRR3722079.547512/100‑1 (MQ=60)
         CAGCCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAG                                                                                   <  SRR3722079.241856/100‑1 (MQ=60)
              CTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTAGGCGATGGGTTCGCCCAGCCAGAAGTCCGC                                                                              >  SRR3722079.373905/1‑100 (MQ=60)
                        ATCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAG                                                                    <  SRR3722079.835207/100‑1 (MQ=60)
                                    GGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGA                                                        <  SRR3722079.403549/100‑1 (MQ=60)
                                    GGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGA                                                        <  SRR3722079.609120/100‑1 (MQ=60)
                                    GGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGA                                                        <  SRR3722079.802574/100‑1 (MQ=60)
                                                                TAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCGAGTGCAGAAG                            <  SRR3722079.147738/100‑1 (MQ=60)
                                                                       CGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCGAGTGCAGAAGCCAGCGC                     >  SRR3722079.535905/1‑100 (MQ=60)
                                                                         AACGGACCGTTGGCGATGGTTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCGAGTGCAGAAGCCAGCGCCG                   >  SRR3722079.314514/1‑100 (MQ=60)
                                                                                 GTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCGAGGGCAGAAGCCAGCGCCGCTTTCTTG           >  SRR3722079.501276/1‑100 (MQ=60)
                                                                                          GGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCGctgtctcttatacacatctccgagcccacgagacca  >  SRR3722079.654701/1‑64 (MQ=60)
                                                                                            |                                                                                                 
GTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGTTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCGAGTGCAGAAGCCAGCGCCGCTTTCTTGGTTGCGCCT  >  NZ_CP009273/951717‑951906

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: