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 I206 R1 206 23.7 1314362 96.8% 1272302 86.2

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,527,151 T→C Q87Q (CAA→CAG nanS ← 9‑O‑acetyl‑N‑acetylneuraminic acid deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,527,1510TC100.0% 112.2 / NA 34Q87Q (CAA→CAGnanS9‑O‑acetyl‑N‑acetylneuraminic acid deacetylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/16);  total (18/16)

GAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAA  >  NZ_CP009273/4527072‑4527236
                                                                               |                                                                                     
gAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGTCAACGGTGCCGTACTGTGTCTGATGATTCg                                                                             <  1:407898/90‑1 (MQ=255)
 aaTCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGt                                                                            >  2:303527/1‑90 (MQ=255)
      cTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCa                                                                       >  1:132912/1‑90 (MQ=255)
      cTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCa                                                                       >  1:497242/1‑90 (MQ=255)
       tGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCag                                                                      <  2:578090/90‑1 (MQ=255)
           ttATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGa                                                                  <  1:511095/90‑1 (MQ=255)
             aTCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATg                                                                <  1:189414/90‑1 (MQ=255)
             aTCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATg                                                                <  1:515413/90‑1 (MQ=255)
                   aaTGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATa                                                          >  2:272786/1‑90 (MQ=255)
                      gAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAAcc                                                       <  2:610304/90‑1 (MQ=255)
                           gCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCa                                                  <  2:132912/90‑1 (MQ=255)
                           gCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCa                                                  <  2:32449/90‑1 (MQ=255)
                           gCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCa                                                  <  1:454946/90‑1 (MQ=255)
                           gCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCa                                                  <  1:209300/90‑1 (MQ=255)
                                   ttCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGa                                          >  2:441395/1‑90 (MQ=255)
                                      cGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACa                                       >  1:594621/1‑90 (MQ=255)
                                        tGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATc                                     <  1:166352/90‑1 (MQ=255)
                                        tGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATc                                     >  2:227714/1‑90 (MQ=255)
                                                   gTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGt                          >  1:169604/1‑90 (MQ=255)
                                                    tGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTg                         >  2:191797/1‑90 (MQ=255)
                                                    tGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTg                         >  1:176341/1‑90 (MQ=255)
                                                        tGGCCAACTGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGAt                                                           <  1:566200/52‑1 (MQ=255)
                                                        tGGCCAACTGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGAt                                                           >  2:566200/1‑52 (MQ=255)
                                                        tGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTc                     >  2:347758/1‑90 (MQ=255)
                                                        tGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTc                     >  2:515778/1‑90 (MQ=255)
                                                           ccAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGt                  <  2:648003/90‑1 (MQ=255)
                                                           ccAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGt                  <  2:378725/90‑1 (MQ=255)
                                                           ccAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGt                  <  1:303527/90‑1 (MQ=255)
                                                             aaCGGTGCCGTACTGTGTCTGATGATTCGTTCCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTgg                >  2:439345/1‑90 (MQ=255)
                                                               cGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGaa              >  1:101788/1‑90 (MQ=255)
                                                                           gtgtCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCa                                                  >  2:520419/1‑42 (MQ=255)
                                                                           gtgtCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCa                                                  <  1:520419/42‑1 (MQ=255)
                                                                           gtgtCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTaa  >  1:280366/1‑90 (MQ=255)
                                                                           gtgtCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTaa  >  2:576278/1‑90 (MQ=255)
                                                                               |                                                                                     
GAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAA  >  NZ_CP009273/4527072‑4527236

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGG  >  NZ_CP009273/4527056‑4527246
                                                                                               |                                                                                               
aagacggcatacgagatctagtacggtctcgtgggctcggagatgtgtataagagacagAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGAT                                                                                             <  SRR3722088.421745/41‑1 (MQ=60)
            ATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCA                                                                                 >  SRR3722088.134297/1‑100 (MQ=60)
            ATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCA                                                                                 >  SRR3722088.503196/1‑100 (MQ=60)
                GAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGTCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGG                                                                             <  SRR3722088.412386/100‑1 (MQ=60)
                           TTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAAC                                                                  <  SRR3722088.517252/100‑1 (MQ=60)
                             ATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCC                                                                <  SRR3722088.521630/100‑1 (MQ=60)
                             ATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCC                                                                <  SRR3722088.191310/100‑1 (MQ=60)
                                  tgggctcggagatgtgtataagagacagATGCAGTGCCTGGCCAACTGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCAT                                                           <  SRR3722088.573143/72‑1 (MQ=60)
                                           tacggtctcgtgggctcggagatgtgtataagagacaGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAAC                                                  <  SRR3722088.526700/63‑1 (MQ=60)
                                           GCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAAC                                                  <  SRR3722088.211393/100‑1 (MQ=60)
                                           GCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAAC                                                  <  SRR3722088.460257/100‑1 (MQ=60)
                                            CAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACA                                                 >  SRR3722088.601938/1‑100 (MQ=60)
                                                        TGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAG                                     <  SRR3722088.168036/100‑1 (MQ=60)
                                                         GCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGT                                    >  SRR3722088.171311/1‑100 (MQ=60)
                                                          CAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTG                                   >  SRR3722088.178106/1‑100 (MQ=60)
                                                                     GCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAA                        >  SRR3722088.102838/1‑100 (MQ=60)
                                                                           CCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGT                  <  SRR3722088.306531/100‑1 (MQ=60)
                                                                                 GTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAA            >  SRR3722088.283061/1‑100 (MQ=60)
                                                                                      GTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGA       >  SRR3722088.326756/1‑100 (MQ=60)
                                                                                      GTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGA       >  SRR3722088.188484/1‑100 (MQ=60)
                                                                                           GTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTctgtctc  >  SRR3722088.171514/1‑93 (MQ=60)
                                                                                           GTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGG  >  SRR3722088.440180/1‑100 (MQ=60)
                                                                                           GTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGG  >  SRR3722088.476666/1‑100 (MQ=60)
                                                                                               |                                                                                               
CAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGG  >  NZ_CP009273/4527056‑4527246

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