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 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,5720TG97.1% 110.2 / ‑4.3 34Q74H (CAA→CACmobBmolybdopterin‑guanine dinucleotide biosynthesis protein MobB
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (11/22);  total (12/22)
Fisher's exact test for biased strand distribution p-value = 3.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCA  >  NZ_CP009273/4034488‑4034651
                                                                                    |                                                                               
cAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTgctg                                                                            >  1:476813/1‑90 (MQ=255)
cAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTgctg                                                                            >  2:714286/1‑90 (MQ=255)
cAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTgctg                                                                            >  1:65564/1‑90 (MQ=255)
           tATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGa                                                                 <  2:464732/90‑1 (MQ=255)
           tATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGa                                                                 <  1:709037/90‑1 (MQ=255)
            aTCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGAt                                                                >  1:562531/1‑90 (MQ=255)
                 ttCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGttt                                                           >  2:4837/1‑90 (MQ=255)
                 ttCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGttt                                                           >  1:86849/1‑90 (MQ=255)
                           gagaAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGcc                                                 <  2:72440/90‑1 (MQ=255)
                           gagaAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGcc                                                 <  1:36458/90‑1 (MQ=255)
                           gagaAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGcc                                                 <  2:499041/90‑1 (MQ=255)
                                   gTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTgc                                         >  2:716838/1‑90 (MQ=255)
                                     aGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTgcgc                                       <  2:86849/90‑1 (MQ=255)
                                     aGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTgcgc                                       <  1:221407/90‑1 (MQ=255)
                                     aGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTgcgc                                       <  2:255017/90‑1 (MQ=255)
                                     aGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTgcgc                                       <  1:592405/90‑1 (MQ=255)
                                          cAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTc                                  <  2:430706/90‑1 (MQ=255)
                                           aGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCt                                            >  1:694721/1‑79 (MQ=255)
                                           aGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCt                                            <  2:694721/79‑1 (MQ=255)
                                           aGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCa                                 >  2:311101/1‑90 (MQ=255)
                                                         tGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGcc                                                 <  1:647515/60‑1 (MQ=255)
                                                         tGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGcc                                                 >  2:647515/1‑60 (MQ=255)
                                                           gtgtTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTg                 <  2:355268/90‑1 (MQ=255)
                                                             gtTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGc               >  1:459981/1‑90 (MQ=255)
                                                              tttCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCt              <  2:156396/90‑1 (MQ=255)
                                                                    tcatcaAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCaa        <  1:357590/90‑1 (MQ=255)
                                                                    tcatcaAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCaa        <  1:360993/90‑1 (MQ=255)
                                                                     catcaAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAAc       <  1:396841/90‑1 (MQ=255)
                                                                      atcaAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACa      <  2:257024/90‑1 (MQ=255)
                                                                      atcaAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACa      <  1:410148/90‑1 (MQ=255)
                                                                       tcaAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACAt     >  1:558591/1‑90 (MQ=255)
                                                                        caAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATc    <  2:65564/90‑1 (MQ=255)
                                                                        caAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATc    <  1:75458/90‑1 (MQ=255)
                                                                          aGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCa  <  1:435233/90‑1 (MQ=255)
                                                                                    |                                                                               
CAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCA  >  NZ_CP009273/4034488‑4034651

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTA  >  NZ_CP009273/4034478‑4034671
                                                                                              |                                                                                                   
GAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTG                                                                                                >  SRR3722112.483635/1‑100 (MQ=60)
GAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTG                                                                                                >  SRR3722112.66353/1‑100 (MQ=60)
            GCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGAT                                                                                    >  SRR3722112.571065/1‑100 (MQ=60)
                 GAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTT                                                                               >  SRR3722112.87926/1‑100 (MQ=60)
                     TATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGC                                                                           <  SRR3722112.720123/100‑1 (MQ=60)
                                     cagatgtgtataagagacagCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGC                                                           <  SRR3722112.657604/80‑1 (MQ=60)
                                     GAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGC                                                           <  SRR3722112.36906/100‑1 (MQ=60)
                                           CTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCC                                                     >  SRR3722112.705570/1‑100 (MQ=60)
                                               AGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGC                                                 <  SRR3722112.224069/100‑1 (MQ=60)
                                               AGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGC                                                 <  SRR3722112.601452/100‑1 (MQ=60)
                                                             TTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGC                                   >  SRR3722112.466479/1‑100 (MQ=60)
                                                                       GTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACAT                         >  SRR3722112.567057/1‑100 (MQ=60)
                                                                              TCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATC                  <  SRR3722112.365613/100‑1 (MQ=60)
                                                                              TCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATC                  <  SRR3722112.362153/100‑1 (MQ=60)
                                                                               CATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCA                 <  SRR3722112.402053/100‑1 (MQ=60)
                                                                               CATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGGGCCGCGCCAGCCGTGCGCAGCTCCTAGCTATCTTTGCCTGnnnnnnnnnnnnnnnnnnnn                 >  SRR3722112.557706/1‑80 (MQ=60)
                                                                                ATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCAT                <  SRR3722112.415575/100‑1 (MQ=60)
                                                                                  CAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGG              <  SRR3722112.76375/100‑1 (MQ=60)
                                                                                    AGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTG            <  SRR3722112.441205/100‑1 (MQ=60)
                                                                                              GTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTA  >  SRR3722112.294864/1‑100 (MQ=60)
                                                                                              GTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTA  >  SRR3722112.449069/1‑100 (MQ=60)
                                                                                              GTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTA  >  SRR3722112.201171/1‑100 (MQ=60)
                                                                                              |                                                                                                   
GAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTA  >  NZ_CP009273/4034478‑4034671

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