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 4,034,572 T→G Q74H (CAA→CAC mobB ← molybdopterin‑guanine dinucleotide biosynthesis protein MobB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,034,5720TG100.0% 101.4 / NA 31Q74H (CAA→CACmobBmolybdopterin‑guanine dinucleotide biosynthesis protein MobB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (14/17);  total (14/17)

CCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCA  >  NZ_CP009273/4034487‑4034651
                                                                                     |                                                                               
ccAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTgct                                                                             <  1:163486/90‑1 (MQ=255)
 cAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTgctg                                                                            >  2:506941/1‑90 (MQ=255)
 cAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTgctg                                                                            >  2:363205/1‑90 (MQ=255)
      ttGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGctggctgg                                                                       <  1:51714/90‑1 (MQ=255)
            tATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGa                                                                 <  1:188558/90‑1 (MQ=255)
            tATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGa                                                                 <  1:499232/90‑1 (MQ=255)
            tATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGa                                                                 <  1:115625/90‑1 (MQ=255)
             aTCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGAt                                                                >  1:194314/1‑90 (MQ=255)
             aTCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGAt                                                                >  1:378706/1‑90 (MQ=255)
             aTCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGAt                                                                >  2:440142/1‑90 (MQ=255)
                 aTTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGtt                                                            <  1:16179/90‑1 (MQ=255)
                  ttCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGttt                                                           >  2:38921/1‑90 (MQ=255)
                       cTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGcc                                                      >  2:251663/1‑90 (MQ=255)
                       cTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGcc                                                      <  2:230832/90‑1 (MQ=255)
                            gagaAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGcc                                                 <  1:363205/90‑1 (MQ=255)
                            gagaAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGcc                                                 <  1:76562/90‑1 (MQ=255)
                               aaaCTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGc                                              <  2:487879/90‑1 (MQ=255)
                                      aGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTgcgc                                       <  2:214694/90‑1 (MQ=255)
                                         tCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGcc                                                 >  1:170412/1‑77 (MQ=255)
                                         tCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGcc                                                 <  2:170412/77‑1 (MQ=255)
                                           cAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTc                                  <  1:506941/90‑1 (MQ=255)
                                           cAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTc                                  <  1:378684/90‑1 (MQ=255)
                                            aGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCa                                 >  1:431501/1‑90 (MQ=255)
                                                  tcttcGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTa                           >  1:523633/1‑90 (MQ=255)
                                                  tcttcGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTa                           >  1:13669/1‑90 (MQ=255)
                                                      cGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCtt                       >  2:137824/1‑90 (MQ=255)
                                                                      catcaAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAAc       >  2:106763/1‑90 (MQ=255)
                                                                       atcaAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACa      <  1:512443/90‑1 (MQ=255)
                                                                       atcaAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACa      <  2:378706/90‑1 (MQ=255)
                                                                         caAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATc    <  1:159328/90‑1 (MQ=255)
                                                                           aGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCa  >  2:277707/1‑90 (MQ=255)
                                                                                     |                                                                               
CCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCA  >  NZ_CP009273/4034487‑4034651

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ 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
CCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTA  >  NZ_CP009273/4034487‑4034671
                                                                                     |                                                                                                   
CCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAAC                                                                                       <  SRR3722114.165110/100‑1 (MQ=60)
   GCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGAT                                                                                    >  SRR3722114.196213/1‑100 (MQ=60)
   GCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGAT                                                                                    >  SRR3722114.383264/1‑100 (MQ=60)
      TTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGT                                                                                 <  SRR3722114.52232/100‑1 (MQ=60)
            TATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGC                                                                           <  SRR3722114.116846/100‑1 (MQ=60)
            TATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGC                                                                           <  SRR3722114.190401/100‑1 (MQ=60)
            TATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGC                                                                           <  SRR3722114.505493/100‑1 (MQ=60)
             ATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGctgtctcttatacacatc                                                                          >  SRR3722114.416687/1‑82 (MQ=60)
                 ATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGC                                                                      <  SRR3722114.16333/100‑1 (MQ=60)
                            GAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGC                                                           <  SRR3722114.77357/100‑1 (MQ=60)
                            GAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGC                                                           <  SRR3722114.367486/100‑1 (MQ=60)
                               AAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCctg                                                        >  SRR3722114.172106/1‑97 (MQ=60)
                                  CTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCA                                                     >  SRR3722114.436817/1‑100 (MQ=60)
                                        ATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTA                                               >  SRR3722114.13797/1‑100 (MQ=60)
                                        ATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTA                                               >  SRR3722114.530186/1‑100 (MQ=60)
                                           acggtctcgtgggctcggagatgtgtataagagacagCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCT                                            <  SRR3722114.285720/63‑1 (MQ=60)
                                           CAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCT                                            <  SRR3722114.513300/100‑1 (MQ=60)
                                           CAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCT                                            <  SRR3722114.383242/100‑1 (MQ=60)
                                                                       ATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCAT                <  SRR3722114.518869/100‑1 (MQ=60)
                                                                         CAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGG              <  SRR3722114.160920/100‑1 (MQ=60)
                                                                                     GTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTA  >  SRR3722114.440330/1‑100 (MQ=60)
                                                                                     |                                                                                                   
CCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTA  >  NZ_CP009273/4034487‑4034671

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