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 F4 I88 R1 217 30.5 1742422 95.9% 1670982 87.3

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 781,531 T→C I148T (ATC→ACC)  aroG → 3‑deoxy‑7‑phosphoheptulonate synthase AroG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273781,5310TC100.0% 114.5 / NA 33I148T (ATC→ACC) aroG3‑deoxy‑7‑phosphoheptulonate synthase AroG
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/14);  total (19/14)

TCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTG  >  NZ_CP009273/781444‑781613
                                                                                       |                                                                                  
tCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCAc                                                                                 <  1:845665/91‑1 (MQ=255)
  aaCGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCAccc                                                                               >  1:861692/1‑91 (MQ=255)
     gacgGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCcac                                                                            <  2:366782/91‑1 (MQ=255)
       cgGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACaa                                                                          <  1:657713/91‑1 (MQ=255)
             gCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATAtctc                                                                    >  1:300529/1‑91 (MQ=255)
               gTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGc                                                                  <  1:538641/91‑1 (MQ=255)
                  taGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGa                                                               <  2:861692/91‑1 (MQ=255)
                       cGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCtga                                                          >  1:94270/1‑91 (MQ=255)
                       cGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGAGAGATTCCCGCTATGACCACCCCACAATATCTAGCTGACCtga                                                          >  2:217370/1‑91 (MQ=255)
                        gTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCtgat                                                         >  1:133345/1‑91 (MQ=255)
                        gTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCtgat                                                         >  2:36820/1‑91 (MQ=255)
                        gTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCtgat                                                         >  1:343776/1‑91 (MQ=255)
                              tgctgcTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTg                                                   <  1:770677/91‑1 (MQ=255)
                               gctgctTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCCCAATATCTCGCTGACCTGATGAGCTgg                                                  >  2:30321/1‑91 (MQ=255)
                                        tatTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAAtt                                         <  1:710733/91‑1 (MQ=255)
                                              cGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCa                                   <  1:36820/91‑1 (MQ=255)
                                                 cAGCGGTCTGCCCGCGGCAGGGGGGTTTCTCGATATGACGACCCCCCAATATCCCGCTGCCCTGACGGCCTGGGGCGCAAATTGCGCACGt                                >  2:380409/1‑91 (MQ=255)
                                                 cAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGGGCAATTGGCGCACGt                                >  2:164534/1‑91 (MQ=255)
                                                 cAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGt                                >  2:298371/1‑91 (MQ=255)
                                                 cAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGt                                <  1:132398/91‑1 (MQ=255)
                                                 cAGCGGTCTGCCAGCGGCAGGGGGGTTTCTCGATATGACCACCCCACAGTATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCCCGt                                >  1:665069/1‑91 (MQ=255)
                                                   gCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGAac                              >  1:862702/1‑91 (MQ=255)
                                                   gCGGTCTGCCAGCGGCAGGGGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTac                              >  1:766384/1‑91 (MQ=255)
                                                     ggTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGc                                                     <  1:831713/66‑1 (MQ=255)
                                                     ggTCTGCCAGCGGCAGGGGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGc                                                     >  2:831713/1‑66 (MQ=255)
                                                      gTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTaccac                           >  1:471932/1‑91 (MQ=255)
                                                        cTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCg                         >  1:716361/1‑91 (MQ=255)
                                                        cTGCCAGCGGCAGGGGAGGTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGGCCGTACCACCg                         >  2:193361/1‑91 (MQ=255)
                                                            cAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATc                     <  2:187946/91‑1 (MQ=255)
                                                                   aGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTg              >  2:391446/1‑91 (MQ=255)
                                                                        aGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACcg         >  1:570745/1‑91 (MQ=255)
                                                                          tttCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACcgcg       <  1:673157/91‑1 (MQ=255)
                                                                           ttCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGa      >  2:427490/1‑91 (MQ=255)
                                                                             ctcGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAAc    <  2:495154/91‑1 (MQ=255)
                                                                               cGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTg  <  2:266606/91‑1 (MQ=255)
                                                                                       |                                                                                  
TCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTG  >  NZ_CP009273/781444‑781613

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCAGGGC  >  NZ_CP009273/781436‑781623
                                                                                               |                                                                                            
CTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCC                                                                                         >  SRR3722118.880860/1‑101 (MQ=60)
 TTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCC                                                                                        <  SRR3722118.759179/101‑1 (MQ=60)
      GATCAACGACGGTCTGCGTATAGCCCGTANNNNNNNNNTTGATATTAACGACNGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAAT                                                                                   <  SRR3722118.690252/101‑1 (MQ=60)
      GATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAAT                                                                                   <  SRR3722118.150118/101‑1 (MQ=60)
        TCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATAT                                                                                 <  SRR3722118.864362/101‑1 (MQ=60)
           cgagataggagtccgtctcgtgggctcggagatgtgtataagagacagGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTC                                                                              <  SRR3722118.628517/53‑1 (MQ=60)
           ACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTC                                                                              >  SRR3722118.306701/1‑101 (MQ=60)
               CGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTG                                                                          <  SRR3722118.671948/101‑1 (MQ=60)
                     GCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGA                                                                    >  SRR3722118.95996/1‑101 (MQ=60)
                      CGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGAT                                                                   >  SRR3722118.135871/1‑101 (MQ=60)
                      CGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGAT                                                                   >  SRR3722118.350902/1‑101 (MQ=60)
                       GTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATG                                                                  <  SRR3722118.550069/101‑1 (MQ=60)
                                    gtgtataagagacagTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAA                                                     <  SRR3722118.850041/86‑1 (MQ=60)
                                      TGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATT                                                   <  SRR3722118.787609/101‑1 (MQ=60)
                                               ATATTAACGACAGCGGTCTGCCAGCGGCAGGGGGGTTTCTCGATATGACCACCCCACAGTATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCCCGT                                          >  SRR3722118.679455/1‑101 (MQ=60)
                                                TATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTA                                         <  SRR3722118.726308/101‑1 (MQ=60)
                                                 ATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGAAC                                        >  SRR3722118.881896/1‑101 (MQ=60)
                                                 ATTAACGACAGCGGTCTGCCAGCGGCAGGGGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTAC                                        >  SRR3722118.783208/1‑101 (MQ=60)
                                                    AACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCAC                                     >  SRR3722118.481901/1‑101 (MQ=60)
                                                      CGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCG                                   >  SRR3722118.732076/1‑101 (MQ=60)
                                                      CGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCG                                   <  SRR3722118.37491/101‑1 (MQ=60)
                                                         CAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAAT                                <  SRR3722118.134909/101‑1 (MQ=60)
                                                                      GCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCG                   >  SRR3722118.582853/1‑101 (MQ=60)
                                                                                  TTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATC       <  SRR3722118.687751/101‑1 (MQ=60)
                                                                                       CGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCAGGGC  >  SRR3722118.738157/1‑101 (MQ=60)
                                                                                               |                                                                                            
CTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCAGGGC  >  NZ_CP009273/781436‑781623

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