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 I230 R1 226 18.8 1048726 96.7% 1014118 86.4

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

AACAGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTC  >  NZ_CP009273/2345327‑2345500
                                                                                       |                                                                                      
aaCAGCAGCGGCGGAATACCACCGCCGACGTTGTGCGCACAGTTCCACACTGCCACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGc                                                                                      <  2:481232/90‑1 (MQ=255)
aaCAGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGc                                                                                      <  1:156908/90‑1 (MQ=255)
   agcagcGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCAcc                                                                                   <  2:412800/90‑1 (MQ=255)
   agcagcGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCAcc                                                                                   <  1:426170/90‑1 (MQ=255)
   agcagcGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCAcc                                                                                   <  2:232750/90‑1 (MQ=255)
          gcggAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTAc                                                                            <  2:384096/90‑1 (MQ=255)
          gcggAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTAc                                                                            <  2:222031/90‑1 (MQ=255)
                 accaccaccGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACcaca                                                                     >  2:252528/1‑90 (MQ=255)
                              ttGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCAcccc                                                        >  1:210086/1‑90 (MQ=255)
                                                 cTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCATCCCCTGGAACCAACCGCa                                    <  1:274436/90‑1 (MQ=255)
                                                   gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCag                                   >  2:204628/1‑90 (MQ=255)
                                                   gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCag                                   >  2:513373/1‑90 (MQ=255)
                                                   gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCag                                   >  1:245376/1‑90 (MQ=255)
                                                     cacaATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCagag                                 >  1:500467/1‑90 (MQ=255)
                                                         aTGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGaac                             >  1:186540/1‑90 (MQ=255)
                                                           gccgccACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGaacaa                           <  2:483537/90‑1 (MQ=255)
                                                             cgccACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAacaaca                         >  1:269351/1‑90 (MQ=255)
                                                                cACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTa                      <  2:245376/90‑1 (MQ=255)
                                                                   gTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACaa                   <  2:269351/90‑1 (MQ=255)
                                                                             cGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGc         <  1:296249/90‑1 (MQ=255)
                                                                                    caGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTc  <  1:480783/90‑1 (MQ=255)
                                                                                       |                                                                                      
AACAGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTC  >  NZ_CP009273/2345327‑2345500

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 27 ≤ 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
ATCCCCAGCAGGAACAGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCC  >  NZ_CP009273/2345315‑2345510
                                                                                                   |                                                                                                
ATCCCCAGCAGGAACAGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGG                                                                                                  >  SRR3722114.224144/1‑100 (MQ=60)
  CCCCAGCAGGAACAGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGC                                                                                                <  SRR3722114.163277/100‑1 (MQ=60)
       GCAGGAACAGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCAT                                                                                           <  SRR3722114.239560/100‑1 (MQ=60)
            AACAGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTAC                                                                                      <  SRR3722114.158482/100‑1 (MQ=60)
               AGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGAC                                                                                   <  SRR3722114.431416/100‑1 (MQ=60)
                      catacgagatctagtacggtctcgtgggctcggagatgtgtataagagacagGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGG                                                                            <  SRR3722114.450745/48‑1 (MQ=60)
                                ACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCC                                                                  >  SRR3722114.212119/1‑100 (MQ=60)
                                                     GTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGctgtctcttatacacatctccgagcccacgagaccgtactagat                                             >  SRR3722114.36210/1‑56 (MQ=60)
                                                     GTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAG                                             >  SRR3722114.247753/1‑100 (MQ=60)
                                                       TCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAG                                           >  SRR3722114.506741/1‑100 (MQ=60)
                                                           CACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAAC                                       >  SRR3722114.188371/1‑100 (MQ=60)
                                                             CTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCATCCCCTGGAACCAACCGCAGAGGAACAAC                                    <  SRR3722114.277294/100‑1 (MQ=60)
                                                               GACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACA                                   >  SRR3722114.272131/1‑100 (MQ=60)
                                                                                         CGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGAC         <  SRR3722114.299405/100‑1 (MQ=60)
                                                                                                CAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCC  <  SRR3722114.486805/100‑1 (MQ=60)
                                                                                                   |                                                                                                
ATCCCCAGCAGGAACAGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCC  >  NZ_CP009273/2345315‑2345510

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