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 I207 R1 218 45.0 2482936 97.1% 2410930 86.1

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,8090TG100.0% 97.6 / NA 29T94P (ACC→CCC) ycaO30S ribosomal protein S12 methylthiotransferase accessory factor YcaO
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (17/12);  total (17/12)

GAGCAGCCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGTTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCGAGTGCAGAAGCCAGCGCCGCTTT  >  NZ_CP009273/951723‑951893
                                                                                      |                                                                                    
gAGCAGCCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTc                                                                                   >  2:1152526/1‑90 (MQ=255)
   cagcCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGcc                                                                                <  1:59430/90‑1 (MQ=255)
   cagcCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGcc                                                                                <  1:667104/90‑1 (MQ=255)
    agcCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGccc                                                                               >  1:1171850/1‑90 (MQ=255)
    agcCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGccc                                                                               >  1:562975/1‑90 (MQ=255)
    agcCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGccc                                                                               >  1:630660/1‑90 (MQ=255)
                 caTCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCg                                                                  >  2:1221246/1‑90 (MQ=255)
                  aTCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGc                                                                 >  2:335245/1‑90 (MQ=255)
                              ggTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAgtt                                                     <  2:1080982/90‑1 (MQ=255)
                              ggTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAgtt                                                     <  2:396011/90‑1 (MQ=255)
                              ggTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAgtt                                                     <  1:111558/90‑1 (MQ=255)
                              ggTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAgtt                                                     <  1:991678/90‑1 (MQ=255)
                              ggTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAgtt                                                     <  1:731487/90‑1 (MQ=255)
                              ggTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAgtt                                                     >  1:15317/1‑90 (MQ=255)
                               gTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAgttg                                                    >  2:637289/1‑90 (MQ=255)
                                   gTGGGAACCATTTTTCGTTGGGATAATGCACGAATGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAgttggttg                                                >  2:530940/1‑90 (MQ=255)
                                           cATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGc                                        <  2:671784/90‑1 (MQ=255)
                                             tttttCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTc                                      >  1:884173/1‑90 (MQ=255)
                                             tttttCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTc                                      >  1:18629/1‑90 (MQ=255)
                                                ttCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGaa                                   <  1:837615/90‑1 (MQ=255)
                                                          tAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCg                         <  1:95873/90‑1 (MQ=255)
                                                                      ggACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCGAGTGCAGAAGcc             >  1:704954/1‑90 (MQ=255)
                                                                       gACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCGAGTGCAGAAGCCa            >  2:596016/1‑90 (MQ=255)
                                                                           gTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCGAGTGCAGAAGCCAgcgc        >  2:495917/1‑90 (MQ=255)
                                                                           gTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCGAGTGCAGAAGCCAgcgc        <  2:51167/90‑1 (MQ=255)
                                                                             tGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCGAGTGCAGAAGCCAGcgccg      >  2:1121094/1‑90 (MQ=255)
                                                                                 gATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCGAGTGCAGAAGCCAGCGCCGCttt  >  2:278258/1‑90 (MQ=255)
                                                                                 gATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCGAGTGCAGAAGCCAGCGCCGCttt  >  2:261983/1‑90 (MQ=255)
                                                                                   tGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGaa                                   <  2:283921/55‑1 (MQ=255)
                                                                                      |                                                                                    
GAGCAGCCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGTTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCGAGTGCAGAAGCCAGCGCCGCTTT  >  NZ_CP009273/951723‑951893

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ 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                                                                                            >  SRR3722089.1189688/1‑100 (MQ=60)
GTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCC                                                                                            >  SRR3722089.638668/1‑100 (MQ=60)
GTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCC                                                                                            >  SRR3722089.569566/1‑100 (MQ=60)
         CAGCCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAG                                                                                   <  SRR3722089.675727/100‑1 (MQ=60)
         CAGCCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAG                                                                                   <  SRR3722089.60157/100‑1 (MQ=60)
                          CGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTT                                                                  >  SRR3722089.15510/1‑100 (MQ=60)
                                    GGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGA                                                        <  SRR3722089.112925/100‑1 (MQ=60)
                                    GGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGA                                                        <  SRR3722089.741004/100‑1 (MQ=60)
                                    GGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGA                                                        <  SRR3722089.1006392/100‑1 (MQ=60)
                                         GTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTC                                                   >  SRR3722089.18866/1‑100 (MQ=60)
                                         GTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTC                                                   >  SRR3722089.896913/1‑100 (MQ=60)
                                                      TTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCG                                      <  SRR3722089.849472/100‑1 (MQ=60)
                                                                TAATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCGAGTGCAGAAG                            <  SRR3722089.97064/100‑1 (MQ=60)
                                                                  ATGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCGAGTGCAGAAGCC                          >  SRR3722089.714139/1‑100 (MQ=60)
                                                                   TGCACGAACGGACCGTTGGCGATGGGTTCGCCCAGCCAGAAGctgtctcttatacacatctgacgctgccgacgaccatcttagtgtagatctcggtggt                         >  SRR3722089.411593/1‑42 (MQ=60)
                                                                               CCGTTGGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCGctgtctcttatacacatctgacgct             >  SRR3722089.287556/1‑75 (MQ=60)
                                                                                    GGCGATGGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCGAGTGCAGAAGCCAGCGCCGCTTTCTTGGTT        >  SRR3722089.319084/1‑100 (MQ=60)
                                                                                          GGGTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCGAGTGCAGAAGCtgtctcttatacacatctgacgctg  >  SRR3722089.535817/1‑75 (MQ=60)
                                                                                            |                                                                                                 
GTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTCGTTGGGATAATGCACGAACGGACCGTTGGCGATGGTTTCGCCCAGCCAGAAGTCCGCAAAAAAGTAGTTGGTTGAGAGACGCTCGAAATATTCACCGAGTGCAGAAGCCAGCGCCGCTTTCTTGGTTGCGCCT  >  NZ_CP009273/951717‑951906

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