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 802,549 T→G T63T (ACA→ACC ybhC ← putative acyl‑CoA thioester hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273802,5490TG100.0% 87.0 / NA 26T63T (ACA→ACCybhCputative acyl‑CoA thioester hydrolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (12/14);  total (12/14)

TACACCTTGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTAC  >  NZ_CP009273/802462‑802625
                                                                                       |                                                                            
tACACCTTGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGgggt                                                                            >  2:1152198/1‑90 (MQ=255)
 acacCTTGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTc                                                                           >  1:1073893/1‑90 (MQ=255)
  cacCTTGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCa                                                                          >  2:383353/1‑90 (MQ=255)
    ccTTGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAgg                                                                        <  2:884107/90‑1 (MQ=255)
    ccTTGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAgg                                                                        <  2:1129763/90‑1 (MQ=255)
      ttGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGa                                                                      >  2:882731/1‑90 (MQ=255)
      ttGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGa                                                                      <  2:1021025/90‑1 (MQ=255)
         aGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGc                                                                   <  1:1093299/90‑1 (MQ=255)
         aGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGc                                                                   <  1:263777/90‑1 (MQ=255)
            gCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCaaa                                                                >  2:877693/1‑90 (MQ=255)
              ccgccgGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAAt                                                              >  2:401845/1‑90 (MQ=255)
              ccgccgGGACGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGc                                                                   >  1:79974/1‑85 (MQ=255)
              ccgccgGGACGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGc                                                                   <  2:79974/85‑1 (MQ=255)
                gccgGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAg                                                            <  2:563100/90‑1 (MQ=255)
                   gggCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAgtgt                                                         <  1:96847/90‑1 (MQ=255)
                        gACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAGTGTTGAGc                                                    <  2:1159813/90‑1 (MQ=255)
                                    aGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTg                                        >  1:55283/1‑90 (MQ=255)
                                                          ggggAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATc                  <  2:1058419/90‑1 (MQ=255)
                                                            ggAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCgg                <  1:783738/90‑1 (MQ=255)
                                                            ggAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCgg                >  2:967540/1‑90 (MQ=255)
                                                                ggATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCga            <  1:741745/90‑1 (MQ=255)
                                                                  aTTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCgaga          >  2:1157547/1‑90 (MQ=255)
                                                                  aTTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCgaga          <  1:561179/90‑1 (MQ=255)
                                                                     ccACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGagaaga       >  2:1228688/1‑90 (MQ=255)
                                                                     ccACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGagaaga       >  2:977682/1‑90 (MQ=255)
                                                                          cTGCAGCACCTGGGGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTAc  >  2:370124/1‑90 (MQ=255)
                                                                                    tGGGGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCgaga          <  2:359061/72‑1 (MQ=255)
                                                                                       |                                                                            
TACACCTTGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTAC  >  NZ_CP009273/802462‑802625

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
GGTATGCGTTACACCTTGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGT  >  NZ_CP009273/802453‑802635
                                                                                                |                                                                                      
GGTATGCGTTACACCTTGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTC                                                                                     >  SRR3722089.1090033/1‑100 (MQ=60)
     GCGTTACACCTTGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGA                                                                                <  SRR3722089.189055/100‑1 (MQ=60)
             CCTTGAGTGCCCGCCGGGACGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAAT                                                                        >  SRR3722089.80950/1‑100 (MQ=60)
                  AGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAGTGT                                                                   <  SRR3722089.1109803/100‑1 (MQ=60)
                  AGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAGTGT                                                                   <  SRR3722089.267156/100‑1 (MQ=60)
                            GGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAGTGTTGAGCATCGA                                                         <  SRR3722089.98050/100‑1 (MQ=60)
                               cccccaaaaccgtctggtggagcaggtagggcaataggggaaggatgagACGCTGCAGCACCGGGGGTCAGGGATGCAAAATAGTGTAGAGCAGCGAAAT                                                      <  SRR3722089.1106293/51‑1 (MQ=30)
                                   CAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTG                                                  >  SRR3722089.55966/1‑100 (MQ=60)
                                                   ggtcgtcggcagcgtcagatgtgtataagagacaGCAGCACCTGGGGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGAC                                  <  SRR3722089.1022016/66‑1 (MQ=60)
                                                                     GGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAG                <  SRR3722089.794435/100‑1 (MQ=60)
                                                                         GGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTAC            <  SRR3722089.751491/100‑1 (MQ=60)
                                                                           ATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCA          <  SRR3722089.567731/100‑1 (MQ=60)
                                                                                   CTGCAGCACCTGGGGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCActgtctct  >  SRR3722089.363491/1‑92 (MQ=60)
                                                                                   CTGCAGCACCTGGGGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGT  >  SRR3722089.470065/1‑100 (MQ=60)
                                                                                   CTGCAGCACCTGGGGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGT  >  SRR3722089.697182/1‑100 (MQ=60)
                                                                                                |                                                                                      
GGTATGCGTTACACCTTGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGT  >  NZ_CP009273/802453‑802635

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