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 I227 R1 223 21.3 1167702 97.0% 1132670 86.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,117,071 T→C F44S (TTC→TCC)  rpmE → 50S ribosomal protein L31

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,117,0710TC100.0% 63.4 / NA 20F44S (TTC→TCC) rpmE50S ribosomal protein L31
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (11/9);  total (11/9)

TTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAG  >  NZ_CP009273/4116991‑4117154
                                                                                |                                                                                   
ttGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCaa                                                                            <  1:502332/90‑1 (MQ=255)
      tAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCg                                                                      >  2:398083/1‑90 (MQ=255)
        aCGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTg                                                                    <  2:454760/90‑1 (MQ=255)
         cGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGa                                                                   <  1:251440/90‑1 (MQ=255)
               gAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGc                                                             >  2:250150/1‑90 (MQ=255)
                        cTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTgg                                                    <  2:51281/90‑1 (MQ=255)
                                 tGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGa                                           >  1:448174/1‑90 (MQ=255)
                                      aTGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCt                                      >  1:231254/1‑90 (MQ=255)
                                             gAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTcaacaa                               <  1:356299/90‑1 (MQ=255)
                                             gAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTcaacaa                               >  1:181943/1‑90 (MQ=255)
                                                 cTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGGCCGCTTCAACAAGCGt                           <  1:250150/90‑1 (MQ=255)
                                                  tCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGAc                                          >  1:489357/1‑74 (MQ=255)
                                                  tCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGAc                                          <  2:489357/74‑1 (MQ=255)
                                                  tCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGtt                          >  1:542482/1‑90 (MQ=255)
                                                      cGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCaa                      <  2:349583/90‑1 (MQ=255)
                                                         gtgCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACAt                   >  1:513285/1‑90 (MQ=255)
                                                                      cACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCaaa      >  1:419576/1‑90 (MQ=255)
                                                                      cACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCaaa      >  1:360097/1‑90 (MQ=255)
                                                                       aCCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAAt     >  2:123020/1‑90 (MQ=255)
                                                                          cGTTCTCCACTGGCAAACAGCGCGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAg  <  1:503442/90‑1 (MQ=255)
                                                                                |                                                                                   
TTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAG  >  NZ_CP009273/4116991‑4117154

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ 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
TTACTGCTAGCTGCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAA  >  NZ_CP009273/4116975‑4117170
                                                                                                |                                                                                                   
cggcagcgtcagatgtgtataagagacagATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCAC                                                                                                  <  SRR3722111.14068/71‑1 (MQ=60)
                TTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGAT                                                                                  <  SRR3722111.510082/100‑1 (MQ=60)
                         CGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACC                                                                         <  SRR3722111.254395/100‑1 (MQ=60)
                                       GCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGA                                                           >  SRR3722111.454954/1‑100 (MQ=60)
                                            ACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCT                                                      >  SRR3722111.234037/1‑100 (MQ=60)
                                                   GTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAA                                               >  SRR3722111.184217/1‑100 (MQ=60)
                                                        GACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACActgtct                                          >  SRR3722111.496879/1‑94 (MQ=60)
                                                        GACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTT                                          >  SRR3722111.551013/1‑100 (MQ=60)
                                                             GAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAAC                                     <  SRR3722111.361197/100‑1 (MQ=60)
                                                               ACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACAT                                   >  SRR3722111.521262/1‑100 (MQ=60)
                                                                 CTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGGCCGCTTCAACAAGCGTTTCAACATCC                                 <  SRR3722111.253088/100‑1 (MQ=60)
                                                                            CAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAA                      >  SRR3722111.365083/1‑100 (MQ=60)
                                                                            CAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAA                      >  SRR3722111.425773/1‑100 (MQ=60)
                                                                                        CCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATT          >  SRR3722111.172932/1‑100 (MQ=60)
                                                                                          CGTTCTCCACTGGCAAACAGCGCGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTC        <  SRR3722111.511210/100‑1 (MQ=60)
                                                                                                CCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACctgtctcttatacacatctgacgctgccgacgaat  >  SRR3722111.583347/1‑65 (MQ=60)
                                                                                                |                                                                                                   
TTACTGCTAGCTGCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAA  >  NZ_CP009273/4116975‑4117170

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