Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 2,313,823 T→C E349G (GAA→GGA)  yojL ← predicted thiamine biosynthesis lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,313,8230TC95.5% 232.2 / 0.7 67E349G (GAA→GGA) yojLpredicted thiamine biosynthesis lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (0/3);  new base C (38/26);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 7.63e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

ATAGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGTTTTAAAGCTATC  >  W3110S.gb/2313760‑2313875
                                                               |                                                    
aTAGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGa                                                >  1:4017693/1‑70 (MQ=255)
aTAGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGa                                                >  1:1856215/1‑70 (MQ=255)
 tAGGCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGa                                                >  1:3749376/1‑69 (MQ=255)
 tAGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGa                                                >  1:349748/1‑69 (MQ=255)
 tAGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGa                                                >  1:3899405/1‑69 (MQ=255)
 tAGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGa                                                >  1:1879617/1‑69 (MQ=255)
 tAGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGa                                                >  1:570683/1‑69 (MQ=255)
 tAGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGa                                                >  1:28531/1‑69 (MQ=255)
 tAGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGAc                                               >  1:1743476/1‑70 (MQ=255)
 tAGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGAc                                               >  1:321256/1‑70 (MQ=255)
 tAGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGAc                                               >  1:1754224/1‑70 (MQ=255)
 tAGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGAc                                               >  1:712750/1‑70 (MQ=255)
 tAGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGAc                                               >  1:1384873/1‑70 (MQ=255)
      tttAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAg                                            <  1:2838663/68‑1 (MQ=255)
      tttAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAg                                            <  1:225277/68‑1 (MQ=255)
      tttAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAg                                            <  1:1456856/68‑1 (MQ=255)
        tAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAg                                        <  1:438777/70‑1 (MQ=255)
        tAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAg                                        <  1:533845/70‑1 (MQ=255)
        tAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAg                                        <  1:276161/70‑1 (MQ=255)
             tGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCtttt                                   >  1:2467753/1‑70 (MQ=255)
             tGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCtttt                                   >  1:2132426/1‑70 (MQ=255)
             tGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCtttt                                   >  1:2020674/1‑70 (MQ=255)
             tGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCtttt                                   >  1:359598/1‑70 (MQ=255)
             tGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCttt                                    >  1:1317065/1‑69 (MQ=255)
                   gTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACtgt                            >  1:2648996/1‑71 (MQ=255)
                   gTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACtgt                            >  1:3744469/1‑71 (MQ=255)
                   gTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACtgt                            >  1:255575/1‑71 (MQ=255)
                     caccaccacGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACt                              >  1:1778574/1‑67 (MQ=255)
                     cACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACt                              >  1:2712625/1‑67 (MQ=255)
                     cACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACt                              >  1:3275936/1‑67 (MQ=255)
                     cACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACt                              >  1:3426351/1‑67 (MQ=255)
                     cACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACt                              >  1:2007420/1‑67 (MQ=255)
                     cACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTgtgg                          >  1:1297261/1‑71 (MQ=255)
                         aTCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGa                       <  1:3899835/70‑1 (MQ=255)
                         aTCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGa                       <  1:506332/70‑1 (MQ=255)
                         aTCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGa                       <  1:1670790/70‑1 (MQ=255)
                         aTCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGa                       <  1:1632763/70‑1 (MQ=255)
                           cACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACAt                    >  1:177243/1‑71 (MQ=255)
                           cACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACAt                    >  1:591726/1‑71 (MQ=255)
                            aCGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACa                     >  1:2541622/1‑69 (MQ=255)
                            aCGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACa                     >  1:2110415/1‑69 (MQ=255)
                            aCGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATc                   >  1:1327770/1‑71 (MQ=255)
                               cAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACtgt                            <  1:672841/59‑1 (MQ=255)
                               cAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACtgt                            <  1:845036/59‑1 (MQ=255)
                               cAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACtgt                            <  1:1200915/59‑1 (MQ=255)
                               cAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACtgt                            <  1:3264217/59‑1 (MQ=255)
                               cAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACtgt                            <  1:1008183/59‑1 (MQ=255)
                               cAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACtgt                            <  1:2385583/59‑1 (MQ=255)
                               cAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACtgt                            <  1:1825348/59‑1 (MQ=255)
                               cAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACtgt                            <  1:2180340/59‑1 (MQ=255)
                                aaaaaCCAACAATCTTGCGCTTTTTTTTTTTCCGCTGACAAGGAAGGTTTTAAACtgt                            <  1:659988/58‑1 (MQ=255)
                                 aaaaCCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGt              <  1:2768638/71‑1 (MQ=255)
                                    aCCAACAATCTTGCGCTTTAATTTTTTCCGCTGCCAAGGAAGCTTTTAAACTGTGGTGACATCCAGGttt            >  1:2144844/1‑70 (MQ=255)
                                    aCCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGttt            >  1:3598901/1‑70 (MQ=255)
                                    aCCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGttt            >  1:1675980/1‑70 (MQ=255)
                                    aCCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGttt            >  1:2210206/1‑70 (MQ=255)
                                     ccAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGTTTTa          <  1:2983026/71‑1 (MQ=255)
                                     ccAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGTTTTa          <  1:2789387/71‑1 (MQ=255)
                                     ccAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGTTTTa          <  1:4090692/71‑1 (MQ=255)
                                         caaTCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACtgt                            <  1:3099054/49‑1 (MQ=255)
                                          aaTCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGTTTTAAAGCt     <  1:2384441/71‑1 (MQ=255)
                                          aaTCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGTTTTAAAGCt     <  1:2769612/71‑1 (MQ=255)
                                          aaTCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGTTTTAAAGCt     <  1:3634693/71‑1 (MQ=255)
                                          aaTCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGTTTTAAAGCt     <  1:914645/71‑1 (MQ=255)
                                             cTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGTTTTAAAGCTATc  >  1:1703970/1‑71 (MQ=255)
                                             cTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGTTTTAAAGCTATc  >  1:1553019/1‑71 (MQ=255)
                                                    tttAATTTTTTCGGCTGACAAGGAAGCTTTTAAACtgt                            <  1:3930158/38‑1 (MQ=38)
                                                               |                                                    
ATAGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGTTTTAAAGCTATC  >  W3110S.gb/2313760‑2313875

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