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 F5 I102 R1 278 39.1 2251254 95.8% 2156701 87.1

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% 132.4 / NA 38I148T (ATC→ACC) aroG3‑deoxy‑7‑phosphoheptulonate synthase AroG
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/19);  total (19/19)

TCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATC  >  NZ_CP009273/781444‑781618
                                                                                       |                                                                                       
tCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGGCAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCAc                                                                                      >  1:859066/1‑91 (MQ=255)
tCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCAc                                                                                      <  1:631153/91‑1 (MQ=255)
tCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCAc                                                                                      <  1:45617/91‑1 (MQ=255)
  aaCGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCAccc                                                                                    <  1:204171/91‑1 (MQ=255)
  aaCGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCAccc                                                                                    <  1:508276/91‑1 (MQ=255)
   acgacgGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCAcccc                                                                                   >  2:731877/1‑91 (MQ=255)
     gacgGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCcac                                                                                 <  2:42616/91‑1 (MQ=255)
        ggTCTGCGTATAGCCCGTAACTTGCTGCTTGATATTAACGCCCGCGGTCTGCCAGCGGCAGGGGAGTTTCTCGATATGACCACCCCACAat                                                                              >  2:1070381/1‑91 (MQ=255)
        ggTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGGGAGTTTCTCGATATGACCACCCCACAat                                                                              >  1:314500/1‑91 (MQ=255)
         gTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGGGAGTTTCTCGATATGACCACCCCACAata                                                                             >  2:489956/1‑91 (MQ=255)
             gCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATAtctc                                                                         >  1:480901/1‑91 (MQ=255)
               gTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGc                                                                       <  1:275367/91‑1 (MQ=255)
                   aGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTgcc                                                                   >  1:926473/1‑89 (MQ=255)
                       cGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCtga                                                               >  1:505515/1‑91 (MQ=255)
                       cGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCtga                                                               >  2:1045236/1‑91 (MQ=255)
                        gTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCtgat                                                              >  1:753216/1‑91 (MQ=255)
                        gTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGGGAGTTTTTCGATATGACCACCCCCCAAAATCTCGCTGACCtgat                                                              >  2:687174/1‑91 (MQ=255)
                         tAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCtgatg                                                             <  2:265116/91‑1 (MQ=255)
                            aTTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGc                                                          <  1:128235/91‑1 (MQ=255)
                             ttGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCt                                                         <  1:772031/91‑1 (MQ=255)
                                  gctTGATATTAACGACAGCGGTCTGCCAGCGGCAGGGGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGgc                                                    >  2:562515/1‑91 (MQ=255)
                                    ttGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGgcgc                                                  <  1:562515/91‑1 (MQ=255)
                                        tatTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAAtt                                              <  1:852178/91‑1 (MQ=255)
                                              cGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCa                                        <  2:963736/91‑1 (MQ=255)
                                               gACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCAc                                       <  2:924337/91‑1 (MQ=255)
                                                 cAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGt                                     >  2:493385/1‑91 (MQ=255)
                                                 cAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGt                                     >  1:392519/1‑91 (MQ=255)
                                                 cAGCGGTCTGCCAGCGGCAGGGGAGTTTTTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGGGCAATTGGCGCACGt                                     >  2:6068/1‑91 (MQ=255)
                                                 cAGCGGTCTGCCAGCGGCAGGGGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGt                                     >  2:870434/1‑91 (MQ=255)
                                                     ggTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTacca                                 >  2:233535/1‑91 (MQ=255)
                                                     ggTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTacca                                 >  1:770183/1‑91 (MQ=255)
                                                      gTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTaccac                                >  1:944297/1‑91 (MQ=255)
                                                       tCTTCCAGCGGCAAGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTaccacc                               <  2:923182/91‑1 (MQ=255)
                                                       tCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTaccacc                               <  2:103066/91‑1 (MQ=255)
                                                        cTGCCAGCGGCAGGTGGGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCg                              >  2:214061/1‑91 (MQ=255)
                                                            cAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATc                          <  2:1021647/91‑1 (MQ=255)
                                                                           ttCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGa           >  1:511296/1‑91 (MQ=255)
                                                                               cGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTg       <  2:480901/91‑1 (MQ=255)
                                                                                   atGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCATGTGCACCGCGAACTGGCAt   <  1:1088395/91‑1 (MQ=255)
                                                                                    tGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATc  <  2:759665/91‑1 (MQ=255)
                                                                                       |                                                                                       
TCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATC  >  NZ_CP009273/781444‑781618

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCAGGGCTTTC  >  NZ_CP009273/781434‑781627
                                                                                                 |                                                                                                
AGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGGCAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCAC                                                                                               >  SRR3722125.878002/1‑101 (MQ=60)
   gggctcggagatgtgtataagagacagGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCC                                                                                            <  SRR3722125.439309/74‑1 (MQ=60)
   agatgtgtataagagacagTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCC                                                                                            <  SRR3722125.58212/82‑1 (MQ=60)
        ataagagacaGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAAT                                                                                       <  SRR3722125.73262/91‑1 (MQ=60)
        GATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAAT                                                                                       <  SRR3722125.221044/101‑1 (MQ=60)
        GATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGGGAGTTTCTCGATATGACCACCCCACAAT                                                                                       >  SRR3722125.320963/1‑101 (MQ=60)
          TCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATAT                                                                                     <  SRR3722125.46442/101‑1 (MQ=60)
          TCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATAT                                                                                     <  SRR3722125.644742/101‑1 (MQ=60)
            AACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCT                                                                                   <  SRR3722125.208201/101‑1 (MQ=60)
            AACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCT                                                                                   <  SRR3722125.519058/101‑1 (MQ=60)
             ACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTC                                                                                  >  SRR3722125.491042/1‑101 (MQ=60)
                   GTGTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGCC                                                                            >  SRR3722125.947040/1‑101 (MQ=60)
                       GCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGA                                                                        >  SRR3722125.516230/1‑101 (MQ=60)
                        CGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGAT                                                                       >  SRR3722125.769495/1‑101 (MQ=60)
                         GTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATG                                                                      <  SRR3722125.280931/101‑1 (MQ=60)
                                      ATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAA                                                         <  SRR3722125.130650/101‑1 (MQ=60)
                                       TTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAAT                                                        <  SRR3722125.788806/101‑1 (MQ=60)
                                              TTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCA                                                 <  SRR3722125.574552/101‑1 (MQ=60)
                                                 ATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGT                                              >  SRR3722125.400738/1‑101 (MQ=60)
                                                  TATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTA                                             <  SRR3722125.870950/101‑1 (MQ=60)
                                                     TAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCA                                          >  SRR3722125.786907/1‑101 (MQ=60)
                                                      AACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCAC                                         >  SRR3722125.965303/1‑101 (MQ=60)
                                                                           GCAGGTGAGTTTCTCGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGA                    >  SRR3722125.522146/1‑101 (MQ=60)
                                                                                         CGATATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCAGGGC      >  SRR3722125.725965/1‑101 (MQ=60)
                                                                                             ATGACCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCATGTGCACCGCGAACTGGCATCAGGGCTTTC  <  SRR3722125.1112891/101‑1 (MQ=60)
                                                                                                 |                                                                                                
AGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCAGGGCTTTC  >  NZ_CP009273/781434‑781627

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