Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 442,006 G→A D412N (GAT→AAT)  thiI → sulfurtransferase required for thiamine and 4‑thiouridine biosynthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb442,0060GA91.0% 217.9 / 11.4 78D412N (GAT→AAT) thiIsulfurtransferase required for thiamine and 4‑thiouridine biosynthesis
Reads supporting (aligned to +/- strand):  ref base G (7/0);  new base A (49/22);  total (56/22)
Fisher's exact test for biased strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

GCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGGCTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCGAAGGGATTGATGTGGTTTC  >  W3110S.gb/441939‑442070
                                                                   |                                                                
gcagACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGTCTTCGGCCCGAACGACGTGATCCTCAata                                                               <  1:645398/71‑1 (MQ=255)
gcagACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGGCTTCGGCCCGAACGACGTGATCCTCAata                                                               <  1:1053470/71‑1 (MQ=255)
gcagACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGGCTTCGGCCCGAACGACGTGATCCTCAata                                                               <  1:806762/71‑1 (MQ=255)
gcagACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGGCTTCGGCCCGAACGACGTGATCCTCAata                                                               <  1:1496683/71‑1 (MQ=255)
gcagACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGGCTTCGGCCCGAACGACGTGATCCTCAata                                                               <  1:1963283/71‑1 (MQ=255)
gcagACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGGCTTCGGCCCGAACGACGTGATCCTCAata                                                               <  1:3494350/71‑1 (MQ=255)
gcagACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGGCTTCGGCCCGAACGACGTGATCCTCAata                                                               <  1:3524139/71‑1 (MQ=255)
gcagACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGGCTTCGGCCCGAACGACGTGATCCTCAata                                                               <  1:1315848/71‑1 (MQ=255)
 cagACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGGCTTCGGCCCGAACGACGTGATCCTCAata                                                               <  1:3796526/70‑1 (MQ=255)
                        agtggaaACCGTCAATGGCTTCGGCCCGAACGACGTGATCCTCAata                                                               <  1:2675013/47‑1 (MQ=255)
                         gtggaaACCGTCAATGGCTTCGGCCCGAACGACGTGATCCTCAata                                                               <  1:3487255/46‑1 (MQ=255)
                              aaCCGTCAATGGCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTATCGATGAACAGGAAGATAAGc                                 >  1:3013065/1‑71 (MQ=255)
                              aaCCGTCAATGGCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTATCGATGAACAGGAAGATAAGc                                 >  1:787481/1‑71 (MQ=255)
                                   tCAATGGCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTc                                                        <  1:16274/43‑1 (MQ=255)
                                   tCAATGGCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTc                                                        <  1:2686941/43‑1 (MQ=255)
                                     aaTGGCTTCGGCCCGAACGACGTGATCCTCAATATCCGtt                                                         >  1:3463558/1‑40 (MQ=255)
                                     aaTGGCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCt                                                       >  1:2529127/1‑42 (MQ=255)
                                     aaTGGCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTa                                                      >  1:3013502/1‑43 (MQ=255)
                                     aaTGGCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTa                                                      >  1:1303415/1‑43 (MQ=255)
                                     aaTGGCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTa                                                      >  1:1841594/1‑43 (MQ=255)
                                     aaTGGCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTa                                                      >  1:184714/1‑43 (MQ=255)
                                     aaTGGCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTa                                                      >  1:2730510/1‑43 (MQ=255)
                                     aaTGGCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTa                                                      >  1:3812504/1‑43 (MQ=255)
                                     aaTGGCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTa                                                      >  1:2693029/1‑43 (MQ=255)
                                     aaTGGCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTa                                                      >  1:755216/1‑43 (MQ=255)
                                     aaTGGCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTa                                                      >  1:2076330/1‑43 (MQ=255)
                                     aaTGGCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTa                                                      >  1:2099435/1‑43 (MQ=255)
                                     aaTGGCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTa                                                      >  1:889960/1‑43 (MQ=255)
                                     aaTGGCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTa                                                      >  1:230541/1‑43 (MQ=255)
                                     aaTGGCTTCGGCCCGAACGACGTGATACTCAATATCCGt                                                          >  1:2954475/1‑39 (MQ=255)
                                       tGGCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTATCGATg                                                <  1:504243/47‑1 (MQ=255)
                                       tGGCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTATCGATg                                                <  1:1017968/47‑1 (MQ=255)
                                       tGGCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTATCGATg                                                <  1:2596721/47‑1 (MQ=255)
                                       tGGCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTATCGATg                                                <  1:3103090/47‑1 (MQ=255)
                                       tGGCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTATCGATg                                                <  1:1017862/47‑1 (MQ=255)
                                       tGGCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTATCGATg                                                <  1:1087978/47‑1 (MQ=255)
                                       tGGCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTATCGATg                                                <  1:1055034/47‑1 (MQ=255)
                                        ggCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTATCGATg                                                <  1:3173365/46‑1 (MQ=255)
                                         gCTTCGGCCCGAACGAGGTGATCCTCAATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTc                      >  1:193060/1‑71 (MQ=255)
                                         gCTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAg                                           >  1:1207269/1‑50 (MQ=255)
                                         gCTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGt                       >  1:2460793/1‑70 (MQ=255)
                                         gCTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTc                      >  1:824002/1‑71 (MQ=255)
                                         gCTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTc                      >  1:2964952/1‑71 (MQ=255)
                                         gCTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTc                      >  1:3780797/1‑71 (MQ=255)
                                         gCTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTc                      >  1:3585540/1‑71 (MQ=255)
                                         gCTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTc                      >  1:1547033/1‑71 (MQ=255)
                                         gCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTATCGATGaa                                              >  1:2065495/1‑47 (MQ=255)
                                         gCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTATCGATGaa                                              >  1:1893728/1‑47 (MQ=255)
                                         gCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTATCGATGAACAGGaa                                        >  1:1314161/1‑53 (MQ=255)
                                         gCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTATCGATGAACAGGaa                                        >  1:93232/1‑53 (MQ=255)
                                         gCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTATCGATGAACAGGaa                                        >  1:3703133/1‑53 (MQ=255)
                                         gCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTATCGATGAACAGGaa                                        >  1:3586314/1‑53 (MQ=255)
                                         gCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTATCGATGAACAGGaa                                        >  1:1424030/1‑53 (MQ=255)
                                         gCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGt                       >  1:4001122/1‑70 (MQ=255)
                                         gCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTc                      >  1:890910/1‑71 (MQ=255)
                                         gCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTc                      >  1:1347308/1‑71 (MQ=255)
                                         gCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTc                      >  1:2740814/1‑71 (MQ=255)
                                         gCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTc                      >  1:3197448/1‑71 (MQ=255)
                                         gCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTc                      >  1:1683042/1‑71 (MQ=255)
                                         gCTTCGGCCCGAACGACGTGATCCTCAATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTc                      >  1:2747452/1‑71 (MQ=255)
                                           ttCGGCCCGAACGACGTGATCCTCAATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGaa                          >  1:737728/1‑65 (MQ=255)
                                           ttCGGCCCGAACGACGTGATCCTCAATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCg                     >  1:3210905/1‑70 (MQ=255)
                                           ttCGGCCCGAACGACGTGATCCTCAATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCg                     >  1:1080898/1‑70 (MQ=255)
                                           ttCGGCCCGAACGACGTGATCCTCAATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCg                     >  1:384616/1‑70 (MQ=255)
                                           ttCGGCCCGAACGACGTGATCCTCAATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCg                     >  1:1679913/1‑70 (MQ=255)
                                           ttCGGCCCGAACGACGTGATCCTCAATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCg                     >  1:1277844/1‑70 (MQ=255)
                                                 ccGAACGACGTGATCCTCAATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCGAAGGGAt              >  1:1631578/1‑71 (MQ=255)
                                                  cGAACGACGTGATCCTCAATATCCGTTCTATCGATg                                                <  1:1341041/36‑1 (MQ=255)
                                                          gTGATCCTCAATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCGaa                   >  1:2015584/1‑57 (MQ=255)
                                                          gTGATCCTCAATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCGAAggg                >  1:3086983/1‑60 (MQ=255)
                                                          gTGATCCTCAATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCGAAGGGAt              >  1:2180671/1‑62 (MQ=255)
                                                          gTGATCCTCAATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCGAAGGGAt              >  1:24513/1‑62 (MQ=255)
                                                          gTGATCCTCAATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCGAAGGGAt              >  1:262002/1‑62 (MQ=255)
                                                          gTGATCCTCAATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCGAAGGGAt              >  1:3505692/1‑62 (MQ=255)
                                                          gTGATCCTCAATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCGAAGGGAt              >  1:76715/1‑62 (MQ=255)
                                                          gTGATCCTCAATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCGAAGGGAt              >  1:3623159/1‑62 (MQ=255)
                                                          gTGATCCTCAATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCGAAGGGAt              >  1:1610474/1‑62 (MQ=255)
                                                             aTCCTCAATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCGAAGGGATTGATGTGGTTtc  >  1:3983272/1‑71 (MQ=255)
                                                                   |                                                                
GCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGGCTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCGAAGGGATTGATGTGGTTTC  >  W3110S.gb/441939‑442070

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: