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 I225 R1 227 21.4 1179702 97.1% 1145490 86.5

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,345,414 T→G H147P (CAC→CCC)  glpT ← glycerol‑3‑phosphate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,345,4140TG100.0% 87.2 / NA 27H147P (CAC→CCC) glpTglycerol‑3‑phosphate transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (10/17);  total (10/17)

AACAGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTC  >  NZ_CP009273/2345327‑2345500
                                                                                       |                                                                                      
aaCAGCAGCGGCGGAATACCCCCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGc                                                                                      <  1:346983/90‑1 (MQ=255)
aaCAGCAGCGGCGGAATACCACCCCCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGc                                                                                      <  1:330765/90‑1 (MQ=255)
aaCAGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGc                                                                                      <  2:72247/90‑1 (MQ=255)
aaCAGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGc                                                                                      <  2:502910/90‑1 (MQ=255)
aaCAGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGc                                                                                      <  2:49570/90‑1 (MQ=255)
     cagcGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCAt                                                                                 <  1:231908/90‑1 (MQ=255)
              aaTACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGAcc                                                                        <  2:405806/90‑1 (MQ=255)
                   caccaccGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACAcg                                                                   <  2:6443/90‑1 (MQ=255)
                          gACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCa                                                            >  1:501932/1‑90 (MQ=255)
                          gACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCa                                                            >  1:251519/1‑90 (MQ=255)
                                gtgCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCAt                                                      <  1:32555/90‑1 (MQ=255)
                                     cacaGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCt                                                 <  1:287478/90‑1 (MQ=255)
                                                  tGACACAATGCCGCCACGTTCTTTCTGCGACCCCCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCa                                    <  1:44619/90‑1 (MQ=255)
                                                   gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCag                                   >  1:208397/1‑90 (MQ=255)
                                                   gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCag                                   >  1:123523/1‑90 (MQ=255)
                                                   gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCag                                   >  1:520072/1‑90 (MQ=255)
                                                   gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCag                                   >  2:410113/1‑90 (MQ=255)
                                                   gACACAATACCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGaacca                                          <  1:18032/83‑1 (MQ=255)
                                                   gACACAATACCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGaacca                                          >  2:18032/1‑83 (MQ=255)
                                                       caATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGa                               <  2:316209/90‑1 (MQ=255)
                                                           gccgccACGTTCTTTCTGCGACCCCCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGaacaa                           <  1:349899/90‑1 (MQ=255)
                                                             cgccACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAacaaca                         >  2:200722/1‑90 (MQ=255)
                                                                   gTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACaa                   <  2:520072/90‑1 (MQ=255)
                                                                        ttctGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATc              >  2:286046/1‑90 (MQ=255)
                                                                              gACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCa        >  1:424851/1‑90 (MQ=255)
                                                                               accaccAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCaa       <  1:410113/90‑1 (MQ=255)
                                                                                    caGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTc  <  2:319731/90‑1 (MQ=255)
                                                                                       |                                                                                      
AACAGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTC  >  NZ_CP009273/2345327‑2345500

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ 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
AACAGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCAT  >  NZ_CP009273/2345327‑2345512
                                                                                       |                                                                                                  
AACAGCAGCGGCGGAATACCCCCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTAC                                                                                        <  SRR3722109.351638/100‑1 (MQ=60)
AACAGCAGCGGCGGAATACCACCCCCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTAC                                                                                        <  SRR3722109.335140/100‑1 (MQ=60)
     CAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCA                                                                                   <  SRR3722109.234654/100‑1 (MQ=60)
                TACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCA                                                                        >  SRR3722109.254490/1‑100 (MQ=60)
                TACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCA                                                                        >  SRR3722109.509513/1‑100 (MQ=60)
                                GTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAAC                                                        <  SRR3722109.32940/100‑1 (MQ=60)
                                     CACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACC                                                   <  SRR3722109.291147/100‑1 (MQ=60)
                                         GTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAG                                               >  SRR3722109.125043/1‑100 (MQ=60)
                                         GTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAG                                               >  SRR3722109.210912/1‑100 (MQ=60)
                                         GTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAG                                               >  SRR3722109.527992/1‑100 (MQ=60)
                                            CCACACTGACACAATACCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGG                                            <  SRR3722109.18266/100‑1 (MQ=60)
                                                  TGACACAATGCCGCCACGTTCTTTCTGCGACCCCCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAAC                                      <  SRR3722109.45149/100‑1 (MQ=60)
                                                           GCCGCCACGTTCTTTCTGCGACCCCCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAAC                             <  SRR3722109.354635/100‑1 (MQ=60)
                                                                    TTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCA                    >  SRR3722109.431099/1‑100 (MQ=60)
                                                                               ACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGT         <  SRR3722109.416087/100‑1 (MQ=60)
                                                                                   CCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACActgtctcttatacacatctgacg     >  SRR3722109.471342/1‑77 (MQ=60)
                                                                                      GGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCAT  >  SRR3722109.582024/1‑100 (MQ=60)
                                                                                       |                                                                                                  
AACAGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCAT  >  NZ_CP009273/2345327‑2345512

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