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 2,849,315 A→C Y540S (TAC→TCC)  flhA → formate hydrogenlyase transcriptional activator FlhA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,849,3150AC100.0% 72.1 / NA 22Y540S (TAC→TCC) flhAformate hydrogenlyase transcriptional activator FlhA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (12/10);  total (12/10)

GACGGACGTGCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGCTGGCGAAAGCCTTTACC  >  NZ_CP009273/2849232‑2849397
                                                                                   |                                                                                  
gacggacgTGCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCt                                                                              <  1:378336/90‑1 (MQ=255)
        tGCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACGTAt                                                                      >  1:795857/1‑90 (MQ=255)
            tCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACGTATTccc                                                                  <  1:306972/90‑1 (MQ=255)
              tAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGa                                                                >  2:1033997/1‑90 (MQ=255)
                   gccgcGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCAc                                                           <  1:215376/90‑1 (MQ=255)
                    ccgcGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCAcc                                                          >  1:735250/1‑90 (MQ=255)
                    ccgcGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCAcc                                                          >  2:667056/1‑90 (MQ=255)
                    ccgcGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCAcc                                                          >  1:257162/1‑90 (MQ=255)
                     cgcgACAACCCGCGATCTAAAAAAAATGGCCGCCGACCGTGAGTCCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCt                                                         <  1:357295/90‑1 (MQ=255)
                       ggACTAACCGCGATCTGAAAAAAATGGCCCCCGCCCGTGAGTCCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTgc                                                       <  1:1214631/89‑1 (MQ=255)
                       cgACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTgc                                                       >  1:475322/1‑90 (MQ=255)
                           tacccGCGATCTGAAAAAAATGGTCGCCGCCCGTGAGTCCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTgccgcc                                                   <  1:891220/87‑1 (MQ=255)
                                                 gcccccGCCCGGGATTCCCGTAGCGATCTCTTTTCCCGCCTGAACTTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGat                             <  1:146255/86‑1 (MQ=255)
                                                      cGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTcc                        <  2:437037/90‑1 (MQ=255)
                                                        cccGTGATTCCCGTAGCGATCTCTTTCCCCGCCTGAACGTATTCCCGATTCCCCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCgc                      <  1:471878/89‑1 (MQ=37)
                                                                  ccGTGGCGATCTCTCTCCCCCCCTGAACGTTTTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGCTGGCGaa            <  1:931372/90‑1 (MQ=38)
                                                                    gTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGCTGGCGAAAg          >  1:1026386/1‑90 (MQ=255)
                                                                         gATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGCTGGCGAAAGCCttt     >  2:1067357/1‑90 (MQ=255)
                                                                         gATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGCTGGCGAAAGCCttt     >  2:341726/1‑90 (MQ=255)
                                                                         gATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGCTGGCGAAAGCCttt     >  2:508923/1‑90 (MQ=255)
                                                                         gATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGCTGGCGAAAGCCttt     >  2:584416/1‑90 (MQ=255)
                                                                            ctctATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGCTGGCGAAAGCCTTTAcc  >  1:599300/1‑90 (MQ=255)
                                                                                   |                                                                                  
GACGGACGTGCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGCTGGCGAAAGCCTTTACC  >  NZ_CP009273/2849232‑2849397

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGCTGGCGAAAGCCTTTACCTTCAAAATTGCC  >  NZ_CP009273/2849217‑2849409
                                                                                                  |                                                                                              
CAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCC                                                                                               <  SRR3722089.1083911/100‑1 (MQ=60)
CAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCC                                                                                               <  SRR3722089.502552/100‑1 (MQ=60)
CAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCC                                                                                               <  SRR3722089.878253/100‑1 (MQ=60)
CAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCC                                                                                               <  SRR3722089.29752/100‑1 (MQ=60)
    AAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCC                                                                                           <  SRR3722089.689135/100‑1 (MQ=60)
             CAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACGTAT                                                                                  >  SRR3722089.806829/1‑100 (MQ=60)
               GACGGACGTGCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACGTATTC                                                                                <  SRR3722089.382956/100‑1 (MQ=60)
                         CGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCACC                                                                      >  SRR3722089.744849/1‑100 (MQ=60)
                         CGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCACC                                                                      >  SRR3722089.260464/1‑100 (MQ=60)
                           TCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTG                                                                    <  SRR3722089.310866/100‑1 (MQ=60)
                            CTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGC                                                                   >  SRR3722089.480983/1‑100 (MQ=60)
                                  GCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCAC                                                             <  SRR3722089.218101/100‑1 (MQ=60)
                                    CGCGACAACCCGCGATCTAAAAAAAATGGCCGCCGACCGTGAGTCCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTA                                                           <  SRR3722089.361703/100‑1 (MQ=60)
                                      GGACTAACCGCGATCTGAAAAAAATGGCCCCCGCCCGTGAGTCCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACG                                                         <  SRR3722089.1233111/100‑1 (MQ=60)
                                          TACCCGCGATCTGAAAAAAATGGTCGCCGCCCGTGAGTCCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAG                                                     <  SRR3722089.904100/100‑1 (MQ=60)
                                                                gcccccgcccgggattcCCGTAGCGATCTCTTTTCCCGCCTGAACTTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGC                               <  SRR3722089.148067/83‑1 (MQ=60)
                                                                       cccgtgattcCCGTAGCGATCTCTTTCCCCGCCTGAACGTATTCCCGATTCCCCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGCTGGCGAA                        <  SRR3722089.477514/90‑1 (MQ=60)
                                                                         CGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGCTGGCGAAAG                      >  SRR3722089.1041697/1‑100 (MQ=60)
                                                                                 CCGTGGCGATCTCTCTCCCCCCCTGAACGTTTTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGCTGGCGAAAGCCTTTACC              <  SRR3722089.944976/100‑1 (MQ=60)
                                                                                 CCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGCTGGCGAAAGCCTTTACC              >  SRR3722089.606660/1‑100 (MQ=60)
                                                                                          TCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGCTGGCGAAAGCCTTTACCTTCAAAATT     >  SRR3722089.1082349/1‑100 (MQ=60)
                                                                                           CTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGCTGGCGAAAGCCTTTACCTTCAAAATTG    >  SRR3722089.952419/1‑100 (MQ=60)
                                                                                           CTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGCTGGCGAAAGCCTTTACCTTCAAAATTG    >  SRR3722089.807951/1‑100 (MQ=60)
                                                                                           CTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGCTGGCGAAAGCCTTTACCTTCAAAATTG    >  SRR3722089.735570/1‑100 (MQ=60)
                                                                                           CTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGCTGGCGAAAGCCTTTACCTTCAAAATTG    >  SRR3722089.298618/1‑100 (MQ=60)
                                                                                           CTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGCTGGCGAAAGCCTTTACCTTCAAAATTG    >  SRR3722089.128666/1‑100 (MQ=60)
                                                                                           CTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGCTGGCGAAAGCCTTTACCTTCAAAATTG    >  SRR3722089.1090289/1‑100 (MQ=60)
                                                                                             CTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGCTGGCGAAAGCCTTTACCTTCAAAATTGCC  >  SRR3722089.440627/1‑100 (MQ=60)
                                                                                                  |                                                                                              
CAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGCTGGCGAAAGCCTTTACCTTCAAAATTGCC  >  NZ_CP009273/2849217‑2849409

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