Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 172,490 A→G 41.4% E343E (GAA→GAG fhuB → fused subunits of iron‑hydroxamate transporter and membrane components of ABC superfamily transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb172,4900AG41.4% 10.4 / 33.3 29E343E (GAA→GAGfhuBfused subunits of iron‑hydroxamate transporter and membrane components of ABC superfamily transporter
Reads supporting (aligned to +/- strand):  ref base A (5/12);  new base G (6/6);  total (11/18)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGT  >  W3110S.gb/172430‑172552
                                                            |                                                              
gCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCaa                                                          <  1:3517687/67‑1 (MQ=255)
gCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCaa                                                          <  1:2215455/67‑1 (MQ=255)
gCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACa                                                        >  1:1567737/1‑69 (MQ=255)
gCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACa                                                        >  1:189781/1‑69 (MQ=255)
            cagcaTTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGtt                                            <  1:3647261/69‑1 (MQ=255)
            cagcaTTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGtt                                            <  1:3202755/69‑1 (MQ=255)
                aTTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAGCGCCAACATGTGCTGGCGTTTGcc                                        >  1:1439426/1‑69 (MQ=255)
                       cgcCGGATATTAAGGTCAACGATCGTGTCGCGGCTGAGCGCCAACATGTGCTGGCGTTTGCCCTCgcgg                                 <  1:362235/69‑1 (MQ=255)
                       cgcCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAGCGCCAACATGTGCTGGCGTTTGCCCTCgcgg                                 <  1:1067822/69‑1 (MQ=255)
                       cgcCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAGCGCCAACATGTGCTGGCGTTTGCCCTCgcgg                                 <  1:524215/69‑1 (MQ=255)
                       cgcCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAGCGCCAACATGTGCTGGCGTTTGCCCTCgcgg                                 <  1:2926675/69‑1 (MQ=255)
                              tatGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGt                          <  1:789326/69‑1 (MQ=255)
                              tatGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGt                          <  1:388154/69‑1 (MQ=255)
                                 gAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGggcggc                            >  1:2169157/1‑64 (MQ=255)
                                       cAACGATCGTGTCGCGGCTGAGCGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGtt                 >  1:186055/1‑69 (MQ=255)
                                       cAACGATCGTGTCGCGGCTGAGCGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGtt                 <  1:3169235/69‑1 (MQ=255)
                                       cAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGtt                 <  1:1786806/69‑1 (MQ=255)
                                       cAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGtt                 <  1:1524793/69‑1 (MQ=255)
                                       cAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGtt                 <  1:1322732/69‑1 (MQ=255)
                                       cAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGtt                 <  1:1305606/69‑1 (MQ=255)
                                        aaCGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGt                                             <  1:2424863/40‑1 (MQ=255)
                                        aaCGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGt                                             <  1:3720028/40‑1 (MQ=255)
                                           gATCGTGTCGCGGCTGAGCGCCAACATGTGCTGGCGTTTGCCCTCGCGggcggc                            <  1:367441/54‑1 (MQ=255)
                                           gATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGtg                         >  1:2542230/1‑57 (MQ=255)
                                              cGTGTCGCGGCTGAGCGCCAACATGTGCTGGCGTTTg                                          >  1:2455856/1‑37 (MQ=255)
                                                    gcgGCTGAGCGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTgtggtg    >  1:2422622/1‑69 (MQ=255)
                                                    gcgGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTgtggtg    >  1:896382/1‑69 (MQ=255)
                                                       gCTGAGCGCCAACATGTGCTGGCGTTTGCCCTCgcgggc                               >  1:1199025/1‑39 (MQ=255)
                                                       gCTGAGCGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTgtggtggt  >  1:1372120/1‑68 (MQ=255)
                                                            |                                                              
GCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGT  >  W3110S.gb/172430‑172552

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