Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,352,814 G→A 75.0% intergenic (‑61/+307) fabI ← / ← ycjD enoyl‑[acyl‑carrier‑protein] reductase, NADH‑dependent/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,352,8140GA75.0% 52.4 / 18.2 32intergenic (‑61/+307)fabI/ycjDenoyl‑[acyl‑carrier‑protein] reductase, NADH‑dependent/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (4/4);  new base A (11/13);  total (15/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAAACCCATAGCTTTAATCCTTATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATA  >  W3110S.gb/1352744‑1352871
                                                                      |                                                         
gaaaACCCATAGCTTTAATCCTTATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAg                                                           <  1:67556/71‑1 (MQ=255)
  aaaCCCATAGCTTTAATCCTTATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAATa                                                         <  1:2365990/71‑1 (MQ=255)
  aaaCCCATAGCTTTAATCCTTATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAATa                                                         <  1:1895486/71‑1 (MQ=255)
       cATAGCTTTAATCCTTATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAATAAAc                                                      <  1:436598/69‑1 (MQ=255)
       cATAGCTTTAATCCTTATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAATAAAc                                                      <  1:2461086/69‑1 (MQ=255)
       cATAGCTTTAATCCTTATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAATAAACAg                                                    <  1:1010675/71‑1 (MQ=255)
         tAGCTTTAATCCTTATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAg                                                    <  1:253704/69‑1 (MQ=255)
         tAGCTTTAATCCTTATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAATAAACa                                                     <  1:1141355/68‑1 (MQ=255)
           gCTTTAATCCTTATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACg                                                <  1:2810558/71‑1 (MQ=255)
           gCTTTAATCCTTATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAATAAACAGTAc                                                 >  1:2036753/1‑70 (MQ=255)
                   ccTTATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAATAAACAGTACGAACAGAt                                         >  1:1992355/1‑70 (MQ=255)
                    cTTATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAATAAACAGTACGAACAGATa                                        >  1:888792/1‑70 (MQ=255)
                      tATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAATAAACAGTACGAACa                                            >  1:2073339/1‑64 (MQ=255)
                       aTTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAATAAACAGTACGAACAGATAAACg                                    >  1:810354/1‑71 (MQ=255)
                          gttgATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATaaa                                      >  1:177776/1‑66 (MQ=255)
                          gttgATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAATAAACAGTACGaa                                              <  1:964968/58‑1 (MQ=255)
                              aTGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAATAAACAGTACGAACAGATAAACGGt                                  <  1:708689/66‑1 (MQ=255)
                               tGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAATAAACAGTACGAACAGATAAACGGt                                  <  1:860697/65‑1 (MQ=255)
                                gCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAATAAACAGTACGAACAGATAAACGGTtattat                            >  1:2247000/1‑70 (MQ=255)
                                 cTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATaa                          >  1:2165417/1‑71 (MQ=255)
                                 cTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAATAAACAGTACGAACAGATAAACGGtta                                >  1:865057/1‑65 (MQ=255)
                                       gtgCCTGAAAATCAGGCGATTCGTCGTTTTAATAAACAGTACGAACAGATAAACGGttat                               >  1:360708/1‑60 (MQ=255)
                                              aaaaTCAGGCGATTCGTCGTTTTAATAAACAGTACGAACAGATAAACGGttatta                             >  1:110923/1‑55 (MQ=255)
                                              aaaaTCAGGCGATTCGTCGTTTTAATAAACAGTACGAACAGATAAACGGttatta                             >  1:1465805/1‑55 (MQ=255)
                                                   cAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAg          >  1:2057212/1‑69 (MQ=255)
                                                      gCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATc                        <  1:1601896/52‑1 (MQ=255)
                                                      gCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCt     <  1:1322057/71‑1 (MQ=255)
                                                       cGATTCGTCGTTTTAATAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCta    <  1:2735548/71‑1 (MQ=255)
                                                       cGATTCGTCGTTTTAATAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCta    <  1:190199/71‑1 (MQ=255)
                                                        gATTCGTCGTTTTAATAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAg          >  1:1955685/1‑64 (MQ=255)
                                                         aTTCGTCGTTTTAATAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCtata  <  1:1608272/71‑1 (MQ=255)
                                                              tcgtTTTAATAAACAGTACGAACAGATAAACGGttat                               <  1:983774/37‑1 (MQ=255)
                                                                gtTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCtg              >  1:2578733/1‑52 (MQ=255)
                                                                      |                                                         
GAAAACCCATAGCTTTAATCCTTATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATA  >  W3110S.gb/1352744‑1352871

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