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 I228 R1 214 26.5 1448192 97.1% 1406194 86.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,934,547 T→C intergenic (‑97/+23) lpxM ← / ← mepM lauroyl‑Kdo(2)‑lipid IV(A) myristoyltransferase/murein DD‑endopeptidase MepM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,934,5470TC100.0% 87.1 / NA 27intergenic (‑97/+23)lpxM/mepMlauroyl‑Kdo(2)‑lipid IV(A) myristoyltransferase/murein DD‑endopeptidase MepM
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/11);  total (16/11)

AGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGA  >  NZ_CP009273/1934470‑1934624
                                                                             |                                                                             
aGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGaa                                                                   >  1:514004/1‑90 (MQ=255)
aGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGaa                                                                   >  1:46863/1‑90 (MQ=255)
aGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGaa                                                                   >  1:375178/1‑90 (MQ=255)
 ggCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAAt                                                                  <  1:330233/90‑1 (MQ=255)
     aaaaTCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAATGGAt                                                              >  1:106372/1‑90 (MQ=255)
      aaaTCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAATGGATg                                                             <  2:676919/90‑1 (MQ=255)
           aatcTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAATGGATGttaat                                                        >  1:60199/1‑90 (MQ=255)
            atcTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAATGGATGttaatt                                                       <  2:388495/90‑1 (MQ=255)
                ggTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAATGGATGTTAATTAATc                                                   <  1:46917/90‑1 (MQ=255)
                 gTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAATGGATGTTAATTAATCa                                                  >  2:127463/1‑90 (MQ=255)
                     tAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAATGGATGTTAATTAATCAAAcc                                              <  1:264652/90‑1 (MQ=255)
                          tAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGc                                         >  1:491107/1‑90 (MQ=255)
                          tAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGc                                         >  1:177391/1‑90 (MQ=255)
                          tAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAATGGATGGTAATTAATCAAACCGTAGc                                         >  1:473950/1‑90 (MQ=255)
                           aGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCt                                        >  2:490342/1‑90 (MQ=255)
                              ggcgCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTg                                       >  2:507397/1‑88 (MQ=255)
                              ggcgCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTg                                       <  1:507397/88‑1 (MQ=255)
                                     cTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAAt                              <  2:491107/90‑1 (MQ=255)
                                     cTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAAt                              <  2:4441/90‑1 (MQ=255)
                                         ccccGCACCAAATAAAAAAAGCCGGTACTGACTGCGCAcc                                                                            >  1:10451/1‑40 (MQ=255)
                                         ccccGCACCAAATAAAAAAAGCCGGTACTGACTGCGCAcc                                                                            <  2:10451/40‑1 (MQ=39)
                                            cGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTg                       >  2:52717/1‑90 (MQ=255)
                                             gCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTgg                      >  1:256761/1‑90 (MQ=255)
                                             gCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTgg                      >  2:525132/1‑90 (MQ=255)
                                                  aaaTAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTgc                 >  1:298151/1‑90 (MQ=255)
                                                       aaaaaaGCCGGTACTGACTGCGCACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAg            <  2:375178/90‑1 (MQ=255)
                                                                 gTACTGACTGCGCACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGcga  <  2:527852/90‑1 (MQ=255)
                                                                          gcgcACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGcga  <  2:244011/81‑1 (MQ=255)
                                                                             |                                                                             
AGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGA  >  NZ_CP009273/1934470‑1934624

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAG  >  NZ_CP009273/1934450‑1934633
                                                                                                 |                                                                                      
TGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCAC                                                                                      <  SRR3722112.545041/100‑1 (MQ=60)
TGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCAC                                                                                      <  SRR3722112.526037/100‑1 (MQ=60)
          GTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAA                                                                            >  SRR3722112.521532/1‑100 (MQ=60)
          GTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAA                                                                            >  SRR3722112.47439/1‑100 (MQ=60)
          GTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAA                                                                            >  SRR3722112.380032/1‑100 (MQ=60)
               GGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAATGGAT                                                                       >  SRR3722112.107693/1‑100 (MQ=60)
                     GGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAATGGATGTTAAT                                                                 <  SRR3722112.334198/100‑1 (MQ=60)
                     GGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAATGGATGTTAAT                                                                 >  SRR3722112.60927/1‑100 (MQ=60)
                                    GGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGC                                                  <  SRR3722112.47493/100‑1 (MQ=60)
                                    GGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGC                                                  >  SRR3722112.498203/1‑100 (MQ=60)
                                    GGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGC                                                  >  SRR3722112.179565/1‑100 (MQ=60)
                                    GGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAATGGATGGTAATTAATCAAACCGTAGC                                                  >  SRR3722112.480723/1‑100 (MQ=60)
                                         TAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGG                                             <  SRR3722112.267761/100‑1 (MQ=60)
                                                GCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATC                                      <  SRR3722112.514785/100‑1 (MQ=60)
                                                   GCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAATctgtctcttatacacatctgacgctgccgacgaatgtcaa                                   >  SRR3722112.10589/1‑60 (MQ=60)
                                                       AACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGG                               >  SRR3722112.259793/1‑100 (MQ=60)
                                                            GCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGCACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGC                          >  SRR3722112.301586/1‑100 (MQ=60)
                                                                                    GGTACTGACTGCGCACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAG  >  SRR3722112.246905/1‑100 (MQ=60)
                                                                                                 |                                                                                      
TGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAG  >  NZ_CP009273/1934450‑1934633

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