Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 1,855,181 (A)7→8 coding (513/1359 nt) ydjE ← predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,855,1741.A97.1% 108.0 / ‑2.8 34coding (520/1359 nt)ydjEpredicted transporter
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base A (16/17);  total (16/18)
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

TGAGAGAAACCAGGCTAACAGTATGCCAATACCACCCAGCAGAAACATTATTCGCCAACTAAAAAAAGCGATAACCACCACGCCTATCGCCGCAGACAG  >  W3110S.gb/1855115‑1855213
                                                            |                                        
tGAGAGAAACCAGGCTAACAGTATGCCAATACCACCCAGCAGAAACATTATTCGCCAACTaaaa                                     >  1:2719668/1‑64 (MQ=255)
tGAGAGAAACCAGGCTAACAGTATGCCAATACCACCCAGCAGAAACATTATTCGCCAACTaaaa                                     >  1:2425335/1‑64 (MQ=255)
tGAGAGAAACCAGGCTAACAGTATGCCAATACCACCCAGCAGAAACATTATTCGCCAACTaaaa                                     >  1:882017/1‑64 (MQ=255)
tGAGAGAAACCAGGCTAACAGTATGCCAATACCACCCAGCAGAAACATTATTCGCCAACTaaaa                                     >  1:592018/1‑64 (MQ=255)
tGAGAGAAACCAGGCTAACAGTATGCCAATACCACCCAGCAGAAACATTATTCGCCAACTaaaa                                     >  1:1126980/1‑64 (MQ=255)
tGAGAGAAACCAGGCTAACAGTATGCCAATACCACCCAGCAGAAACATTATTCGCCAACTaaaa                                     >  1:180514/1‑64 (MQ=255)
tGAGAGAAACCAGGCTAACAGTATGCCAATACCACCCAGCAGAAACATTATTCGCCAACTaaaa                                     >  1:1624877/1‑64 (MQ=255)
tGAGAGAAACCAGGCTAACAGTATGCCAATACCACCCAGCAGAAACATTATTCGCCAACTaaaa                                     >  1:1535539/1‑64 (MQ=255)
          cAGGCTAACAGTATGCCAATACCACCCAGCAGAAAGATTATTCGCCAACTAAAAAAAAGCGATaa                           <  1:2353458/65‑1 (MQ=255)
          cAGGCTAACAGTATGCCAATACCACCCAGCAGAAACATTATTCGCCAACTAAAAAAAAGCGATaa                           <  1:2915885/65‑1 (MQ=255)
          cAGGCTAACAGTATGCCAATACCACCCAGCAGAAACATTATTCGCCAACTAAAAAAAAGCGATaa                           <  1:2298122/65‑1 (MQ=255)
          cAGGCTAACAGTATGCCAATACCACCCAGCAGAAACATTATTCGCCAACTAAAAAAAAGCGATaa                           <  1:2253195/65‑1 (MQ=255)
          cAGGCTAACAGTATGCCAATACCACCCAGCAGAAACATTATTCGCCAACTAAAAAAAAGCGATaa                           <  1:174282/65‑1 (MQ=255)
          cAGGCTAACAGTATGCCAATACCACCCAGCAGAAACATTATTCGCCAACTAAAAAAAAGCGATaa                           <  1:642100/65‑1 (MQ=255)
          cAGGCTAACAGTATGCCAATACCACCCAGCAGAAACATTATTCGCCAACTAAAAAAAAGCGATaa                           <  1:1235983/65‑1 (MQ=255)
                           aaTACCACCCAGCAGAAACATTATTCGCCAACTAAAAAAAAGCGATAACCACCAGGCCTATcgc           >  1:1231960/1‑64 (MQ=255)
                           aaTACCACCCAGCAGAAACATTATTCGCCAACTAAAAAAAAGCGATAACCACCACGCCTATcgc           >  1:796843/1‑64 (MQ=255)
                           aaTACCACCCAGCAGAAACATTATTCGCCAACTAAAAAAAAGCGATAACCACCACGCCTATcgc           >  1:2912944/1‑64 (MQ=255)
                           aaTACCACCCAGCAGAAACATTATTCGCCAACTAAAAAAAAGCGATAACCACCACGCCTATcgc           >  1:2660654/1‑64 (MQ=255)
                           aaTACCACCCAGCAGAAACATTATTCGCCAACTAAAAAAAAGCGATAACCACCACGCCTATcgc           >  1:2595306/1‑64 (MQ=255)
                           aaTACCACCCAGCAGAAACATTATTCGCCAACTAAAAAAAAGCGATAACCACCACGCCTATcgc           >  1:2451931/1‑64 (MQ=255)
                           aaTACCACCCAGCAGAAACATTATTCGCCAACTAAAAAAAAGCGATAACCACCACGCCTATcgc           >  1:916280/1‑64 (MQ=255)
                           aaTACCACCCAGCAGAAACATTATTCGCCAACTAAAAAAAAGCGATAACCACCACGCCTATcgc           >  1:1123674/1‑64 (MQ=255)
                           aaTACCACCCAGCAGAAACATTATTCGCCAACTAAAAAAAAGCGATAACCACCACGCCTATcgc           >  1:814711/1‑64 (MQ=255)
                           aaTACCACCCAGCAGAAACATTATTCGCCAACTAAAAAAAAGCGATAACCACCACGCCTATcgc           >  1:1301804/1‑64 (MQ=255)
                           aaTACCACCCAGCAGAAACATTATTCGCCAACTAAAAAAAAGCGATAACCACCACGCCTATcgc           >  1:804404/1‑64 (MQ=255)
                           aaTACCACCCAGCAGAAACATTATTCGCCAACTAAAAAAAAGCGATAACCACCACGCCTATcgc           >  1:1732199/1‑64 (MQ=255)
                               ccaccCAGCAGAAACATTATTCGCCAACTAAAAAAAAGCGATAACCACCACGCCTATCGCCGCAg      <  1:552769/65‑1 (MQ=255)
                               ccaccCAGCAGAAACATTATTCGCCAACTAAAAAAAAGCGATAACCACCACGCCTATCGCCGCAg      <  1:1019980/65‑1 (MQ=255)
                               ccaccCAGCAGAAACATTATTCGCCAACTAAAAAAAAGCGATAACCACCACGCCTATCGCCGCAg      <  1:721901/65‑1 (MQ=255)
                               ccaccCAGCAGAAACATTATTCGCCAACTAAAAAAAAGCGATAACCACCACGCCTATCGCCGCAg      <  1:936516/65‑1 (MQ=255)
                               ccaccCAGCAGAAACATTATTCGCCAACTAAAAAAAAGCGATAACCACCACGCCTATCGCCGCAg      <  1:539184/65‑1 (MQ=255)
                               ccaccCAGCAGAAACATTATTCGCCAACTAAAAAAAAGCGATAACCACCACGCCTATCGCCGCAg      <  1:2401952/65‑1 (MQ=255)
                               ccaccCAGCAGAAACATTATTCGCCAACTAAAAAAAAGCGATAACCACCACGCCTATCGCCGCAg      <  1:1560015/65‑1 (MQ=255)
                               ccaccCAGCAGAAACATTATTCGCCAACTAAAAAAAAGCGATAACCACCACGCCTATCGCCGCAg      <  1:153197/65‑1 (MQ=255)
                               ccaccCAGCAGAAACATTATTCGCCAACTAAAAAAAAGCGATAACCACCACGCCTATCGCCGCAg      <  1:1054404/65‑1 (MQ=255)
                               ccaccCAGCAGAAACATTATTCGCCAACTAAAAAAAAGCGATAACCACCACGCCTATCGCCGCAg      <  1:1034844/65‑1 (MQ=255)
                                   ccAGCAGAAACATTATTCGCCAACTAAAAAAAAGCGAt                             >  1:152981/1‑38 (MQ=38)
                                   ccAGCAGAAACATTATTCGCCAACTAAAAAAAAGCGATAACCACCACGCCTATCGCCGCAGACAg  >  1:2362451/1‑65 (MQ=255)
                                   ccAGCAGAAACATTATTCGCCAACTAAAAAAAAGCGATAACCACCACGCCTATCGCCGCAGACAg  >  1:145092/1‑65 (MQ=255)
                                   ccAGCAGAAACATTATTCGCCAACTAAAAAAAAGCGATAACCACCACGCCTATCGCCGCAGACAg  >  1:1107885/1‑65 (MQ=255)
                                                      ccAACTATAAAAAAGCGATAACCACCACGCCTATCGCCGCAg      <  1:1169472/42‑1 (MQ=38)
                                                            |                                        
TGAGAGAAACCAGGCTAACAGTATGCCAATACCACCCAGCAGAAACATTATTCGCCAACTAAAAAAAGCGATAACCACCACGCCTATCGCCGCAGACAG  >  W3110S.gb/1855115‑1855213

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: