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,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% 147.9 / NA 43H147P (CAC→CCC) glpTglycerol‑3‑phosphate transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (22/21);  total (22/21)

AACAGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTC  >  NZ_CP009273/2345327‑2345500
                                                                                       |                                                                                      
aaCAGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGc                                                                                      <  2:1052337/90‑1 (MQ=255)
aaCAGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGc                                                                                      <  1:387420/90‑1 (MQ=255)
     cagcGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCAt                                                                                 <  1:949413/90‑1 (MQ=255)
        cggcggAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGt                                                                              <  2:1141971/90‑1 (MQ=255)
          gcggAATACCACCCCCGACGTTGTGCGCACAGTTCCACACTGACCCAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTAc                                                                            <  1:57247/90‑1 (MQ=255)
                       accGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACAcggcgg                                                               >  1:255815/1‑90 (MQ=255)
                          gACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCa                                                            >  2:998138/1‑90 (MQ=255)
                              ttGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCAcccc                                                        >  1:381569/1‑90 (MQ=255)
                                gtgCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCAt                                                      <  1:995946/90‑1 (MQ=255)
                                    gcACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATcccc                                                  >  2:229343/1‑90 (MQ=255)
                                    gcACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATcccc                                                  <  2:101577/90‑1 (MQ=255)
                                        aGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGa                                              >  1:183172/1‑90 (MQ=255)
                                         gTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGaa                                             >  2:529881/1‑90 (MQ=255)
                                           tCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGaacc                                           <  2:475073/90‑1 (MQ=255)
                                           tCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGaacc                                           <  2:1094308/90‑1 (MQ=255)
                                                  tGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCa                                    <  2:1054276/90‑1 (MQ=255)
                                                   gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGTACCAACCGCag                                   >  2:560010/1‑90 (MQ=255)
                                                   gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGc                                     <  2:1086632/88‑1 (MQ=255)
                                                   gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGc                                     >  1:1086632/1‑88 (MQ=255)
                                                   gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCag                                   >  1:336910/1‑90 (MQ=255)
                                                   gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCag                                   >  1:979393/1‑90 (MQ=255)
                                                   gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCag                                   >  2:838678/1‑90 (MQ=255)
                                                   gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCag                                   >  2:841362/1‑90 (MQ=255)
                                                   gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCag                                   >  2:397623/1‑90 (MQ=255)
                                                   gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCag                                   >  2:366141/1‑90 (MQ=255)
                                                       caATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGa                               <  2:1150652/90‑1 (MQ=255)
                                                         aTGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCaga                                  <  1:581300/85‑1 (MQ=255)
                                                         aTGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCaga                                  >  2:581300/1‑85 (MQ=255)
                                                         aTGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGaac                             <  2:212134/90‑1 (MQ=255)
                                                         aTGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGaac                             >  2:140798/1‑90 (MQ=255)
                                                         aTGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGaac                             >  1:830113/1‑90 (MQ=255)
                                                         aTGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGaac                             >  2:811854/1‑90 (MQ=255)
                                                         aTGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGaac                             <  1:658608/90‑1 (MQ=255)
                                                           gccgccACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGaacaa                           >  1:492709/1‑90 (MQ=255)
                                                            ccgccACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGaacaac                          <  2:1232831/90‑1 (MQ=255)
                                                              gccACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAg                        >  2:22302/1‑90 (MQ=255)
                                                                cACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTa                      <  2:958537/90‑1 (MQ=255)
                                                                   gTTCTTTCTGCGCCCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACaa                   <  1:658056/90‑1 (MQ=255)
                                                                   gTTCTTTCTGCGACCCCCAGGGCCCCATAGTACGACCACACGGCGGCCCCCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACaa                   <  2:734634/90‑1 (MQ=255)
                                                                        ttctGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATc              <  2:979393/90‑1 (MQ=255)
                                                                            gCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCg          >  2:517893/1‑90 (MQ=255)
                                                                              gACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCa        >  2:741632/1‑90 (MQ=255)
                                                                                    caGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTc  <  2:1082268/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
CCAGCAGGAACAGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCAT  >  NZ_CP009273/2345319‑2345512
                                                                                               |                                                                                                  
CCAGCAGGAACAGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCAC                                                                                                <  SRR3722089.656121/100‑1 (MQ=60)
     AGGAACAGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAG                                                                                           <  SRR3722089.22400/100‑1 (MQ=60)
     AGGAACAGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAG                                                                                           <  SRR3722089.236131/100‑1 (MQ=60)
        AACAGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTAC                                                                                        <  SRR3722089.392140/100‑1 (MQ=60)
           cggcgaccaccgagatctacactaagatggtcgtcggcagcgtcagatgtgtataagagacagCGTTCTTTCTGCGACCACCAGGGCACCATAGTACGAC                                                                                     <  SRR3722089.358784/37‑1 (MQ=60)
             CAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCA                                                                                   <  SRR3722089.963365/100‑1 (MQ=60)
                  GCGGAATACCACCCCCGACGTTGTGCGCACAGTTCCACACTGACCCAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGG                                                                              <  SRR3722089.57951/100‑1 (MQ=60)
                     GAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGG                                                                           >  SRR3722089.259099/1‑100 (MQ=60)
                            ACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCC                                                                    >  SRR3722089.386224/1‑100 (MQ=60)
                                      TTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGA                                                          >  SRR3722089.185461/1‑100 (MQ=60)
                                        GTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAAC                                                        <  SRR3722089.1010736/100‑1 (MQ=60)
                                                 GTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAG                                               >  SRR3722089.993885/1‑100 (MQ=60)
                                                 GTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAG                                               >  SRR3722089.1103004/1‑100 (MQ=60)
                                                 GTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAG                                               >  SRR3722089.341112/1‑100 (MQ=60)
                                                       CACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAAC                                         >  SRR3722089.841841/1‑100 (MQ=60)
                                                         CTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAA                                       >  SRR3722089.498544/1‑100 (MQ=60)
                                                            ACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAG                                    <  SRR3722089.588267/100‑1 (MQ=60)
                                                                 ATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAA                               <  SRR3722089.667095/100‑1 (MQ=60)
                                                                           GTTCTTTCTGCGCCCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGC                     <  SRR3722089.666536/100‑1 (MQ=60)
                                                                                           CCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCC     >  SRR3722089.592400/1‑100 (MQ=60)
                                                                                           CCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCC     >  SRR3722089.93508/1‑100 (MQ=60)
                                                                                              GGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACctgtctcttatacacatctgacgctgccgacgacca  >  SRR3722089.1064689/1‑64 (MQ=60)
                                                                                               |                                                                                                  
CCAGCAGGAACAGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCAT  >  NZ_CP009273/2345319‑2345512

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ 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: