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 I210 R1 226 19.6 1130256 95.4% 1078264 84.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,202,171 G→C W242C (TGG→TGC yehP → VWA domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,202,1710GC100.0% 64.5 / NA 21W242C (TGG→TGCyehPVWA domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (7/14);  total (7/14)

AACTGGTCTTACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTA  >  NZ_CP009273/2202092‑2202252
                                                                               |                                                                                 
aaCTGGTCTTACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGCCAGTTACCCg                                                                         >  2:134160/1‑90 (MQ=255)
     gTCTTACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGCCAGTTACCCGGCAtt                                                                    <  1:443305/90‑1 (MQ=255)
       cTTACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGCCAGTTACCCGGCATTCg                                                                  >  2:91600/1‑90 (MQ=255)
         tACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGCCAGTTACCCGGCATTCGTa                                                                >  1:287142/1‑90 (MQ=255)
          aCTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGCCAGTTACCCGGCATTCGTAc                                                               <  2:92673/90‑1 (MQ=255)
           cTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGCCAGTTACCCGGCATTCGTAcc                                                              <  1:182833/90‑1 (MQ=255)
           cTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGCCAGTTACCCGGCATTCGTAcc                                                              <  2:328182/90‑1 (MQ=255)
           cTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGCCAGTTACCCGGCATTCGTAcc                                                              <  1:249904/90‑1 (MQ=255)
                           gATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGCCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTg                                              <  1:134160/90‑1 (MQ=255)
                               gATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGCCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGacac                                          <  2:195400/90‑1 (MQ=255)
                               gATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGCCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGacac                                          <  2:357755/90‑1 (MQ=255)
                                aTGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGCCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGAcaca                                         >  2:199892/1‑90 (MQ=255)
                                     cGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGCCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGcgt                                    <  1:61313/90‑1 (MQ=255)
                                     cGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGCCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGcgt                                    >  2:38305/1‑90 (MQ=255)
                                          cGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGCCAGATACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTg                               <  2:429650/90‑1 (MQ=255)
                                               aTCCACTCTGCGGTGATGGCGGCCTGTTTGTGCCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGAtctc                          <  1:38305/90‑1 (MQ=255)
                                                      ctGCGGTGATGGCGGCCTGTTTGTGCCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAg                   >  1:473744/1‑90 (MQ=255)
                                                              aTGGCGGCCTGTTTGTGCCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGcc           <  1:91600/90‑1 (MQ=255)
                                                                   ggcCTGTTTGTGCCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATcc      <  2:540475/90‑1 (MQ=255)
                                                                      cTGTTTGTGCCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGt   <  1:199892/90‑1 (MQ=255)
                                                                       tGTTTGTGCCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTa  >  1:356601/1‑90 (MQ=255)
                                                                               |                                                                                 
AACTGGTCTTACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTA  >  NZ_CP009273/2202092‑2202252

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 33 ≤ 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
CGCCAAAGCGAACAATGGCAACTGGTCTTACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAAT  >  NZ_CP009273/2202073‑2202263
                                                                                                  |                                                                                            
CGCCAAAGCGAACAATGGCAACTGGTCTTACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGCC                                                                                             <  SRR3722092.121157/100‑1 (MQ=60)
                  CAACTGGTCTTACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGCCAGTTACCCGGCATTCGTA                                                                           >  SRR3722092.292100/1‑100 (MQ=60)
                        GTCTTACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGCCAGTTACCCGGCATTCGTACCCATC                                                                     <  SRR3722092.452012/100‑1 (MQ=60)
                              CTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGCCAGTTACCCGGCATTCGTACCCATCTGGTGG                                                               <  SRR3722092.185876/100‑1 (MQ=60)
                              CTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGCCAGTTACCCGGCATTCGTACCCATCTGGTGG                                                               <  SRR3722092.253951/100‑1 (MQ=60)
                                              GATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGCCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGT                                               <  SRR3722092.136417/100‑1 (MQ=60)
                                                        CGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGCCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTC                                     <  SRR3722092.62311/100‑1 (MQ=60)
                                                             CGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGCCAGTTACCCGGCATctgtctcttatacacatctgacgctgccgacgaccatcttagtgtaga                                >  SRR3722092.31194/1‑52 (MQ=60)
                                                               GTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGCCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAG                              >  SRR3722092.483145/1‑100 (MQ=60)
                                                                  ATCCACTCTGCGGTGATGGCGGCCTGTTTGTGCCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACG                           <  SRR3722092.38976/100‑1 (MQ=60)
                                                                                GATGGCGGCCTGTTTGTGCCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTA             >  SRR3722092.363140/1‑100 (MQ=60)
                                                                                 ATGGCGGCCTGTTTGTGCCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAG            <  SRR3722092.93169/100‑1 (MQ=60)
                                                                                         CTGTTTGTGCCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTA    <  SRR3722092.203196/100‑1 (MQ=60)
                                                                                           GTTTGTGCCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCCCGGCAGACGTTGCCGATCCGGAAGAGTTATTAAT  >  SRR3722092.33238/1‑100 (MQ=60)
                                                                                                  |                                                                                            
CGCCAAAGCGAACAATGGCAACTGGTCTTACTGGTTGATCAAAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAAT  >  NZ_CP009273/2202073‑2202263

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