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,819,095 T→C intergenic (‑160/‑96) mltB ← / → srlA lytic murein transglycosylase B/PTS glucitol/sorbitol transporter subunit IIC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,819,0950TC100.0% 88.6 / NA 27intergenic (‑160/‑96)mltB/srlAlytic murein transglycosylase B/PTS glucitol/sorbitol transporter subunit IIC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/13);  total (14/13)

GACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAG  >  NZ_CP009273/2819024‑2819181
                                                                       |                                                                                      
gACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATa                                                                      >  2:1005499/1‑90 (MQ=255)
  cTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAac                                                                    <  2:135208/90‑1 (MQ=255)
  cTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAac                                                                    <  2:254799/90‑1 (MQ=255)
    ttACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAacac                                                                  >  1:12299/1‑90 (MQ=255)
     tACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACt                                                                 >  1:1210435/1‑90 (MQ=255)
          aactcaactcATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTaa                                                            >  1:626309/1‑90 (MQ=255)
          aactcaactcATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTaa                                                            >  1:104447/1‑90 (MQ=255)
              caactcaTTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCt                                                        >  2:1225591/1‑90 (MQ=255)
              caactcaTTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCt                                                        >  1:760198/1‑90 (MQ=255)
                actcaTTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCttt                                                      >  1:666416/1‑90 (MQ=255)
                actcaTTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCttt                                                      >  2:1040652/1‑90 (MQ=255)
                     ttCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATc                                                 <  2:203189/90‑1 (MQ=255)
                     ttCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATc                                                 <  2:12299/90‑1 (MQ=255)
                     ttCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATc                                                 <  1:703876/90‑1 (MQ=255)
                     ttCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATc                                                 >  2:635861/1‑90 (MQ=255)
                        cccTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTga                                              >  1:992207/1‑90 (MQ=255)
                         ccTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTgat                                             <  1:294834/90‑1 (MQ=255)
                         ccTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTgat                                             <  2:721727/90‑1 (MQ=255)
                         ccTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTCTGCGATCAAAATAACACTTTTAAATCTTTCAATCTgat                                             <  1:838069/90‑1 (MQ=255)
                                      ttATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGc                                >  2:96160/1‑90 (MQ=255)
                                      ttATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGc                                >  1:998560/1‑90 (MQ=255)
                                        aTTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGcct                              >  2:747610/1‑89 (MQ=255)
                                                     aaTTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATaaaa                 <  2:1093582/90‑1 (MQ=255)
                                                     aaTTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATaaaa                 <  2:544791/90‑1 (MQ=255)
                                                          aaaTACTTATCTTCCGGTTTGCCAACAAAATAACACAAATAAATCTTTCACCCTCACTAGATTAGGTTGCCGTTTGGTaacaaaacaaaa            <  1:1112974/90‑3 (MQ=255)
                                                             tatTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAAt         <  1:604034/90‑1 (MQ=255)
                                                                    cTTCCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAg  <  1:1102391/90‑1 (MQ=255)
                                                                       |                                                                                      
GACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAG  >  NZ_CP009273/2819024‑2819181

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAA  >  NZ_CP009273/2819018‑2819191
                                                                             |                                                                                                
GATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACAC                                                                            >  SRR3722089.12444/1‑100 (MQ=60)
 ATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACT                                                                           >  SRR3722089.1228839/1‑100 (MQ=60)
      GACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAA                                                                      >  SRR3722089.105738/1‑100 (MQ=60)
      GACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAA                                                                      >  SRR3722089.634234/1‑100 (MQ=60)
          TTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCT                                                                  >  SRR3722089.770348/1‑100 (MQ=60)
            ACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCTTT                                                                >  SRR3722089.675026/1‑100 (MQ=60)
                    CAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGA                                                        >  SRR3722089.1006929/1‑100 (MQ=60)
                           TTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATT                                                 <  SRR3722089.713046/100‑1 (MQ=60)
                               CCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGT                                             <  SRR3722089.298598/100‑1 (MQ=60)
                               CCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTCTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGT                                             <  SRR3722089.849932/100‑1 (MQ=60)
                                  CGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGC                                          >  SRR3722089.1013404/1‑100 (MQ=60)
                                                                AAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAG            >  SRR3722089.957200/1‑100 (MQ=60)
                                                                   TATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAG         <  SRR3722089.611496/100‑1 (MQ=60)
                                                                    ATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACctgtctcttatacacatc        >  SRR3722089.456572/1‑82 (MQ=60)
                                                                          CTTCCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAA  <  SRR3722089.1119058/100‑1 (MQ=60)
                                                                             |                                                                                                
GATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAA  >  NZ_CP009273/2819018‑2819191

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