Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,313,823 T→C 62.2% 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,8230TC62.2% 23.0 / 30.9 32E349G (GAA→GGA) yojLpredicted thiamine biosynthesis lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (6/6);  new base C (8/11);  total (14/18)
Fisher's exact test for biased strand distribution p-value = 7.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAGGATAGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGTTTTAAAGCT  >  W3110S.gb/2313755‑2313872
                                                                    |                                                 
tAAGGATAGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAAttttttt                                                   <  1:1026958/69‑1 (MQ=255)
tAAGGATAGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAAttttttt                                                   <  1:1652297/69‑1 (MQ=255)
      tAGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGCTg                                              >  1:3388837/1‑68 (MQ=255)
      tAGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGa                                             <  1:88635/69‑1 (MQ=255)
       aGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGCTGAc                                            <  1:1157512/69‑1 (MQ=255)
       aGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGCTGAc                                            <  1:2031049/69‑1 (MQ=255)
        gCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGa                                             <  1:34098/67‑1 (MQ=255)
        gCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGAc                                            <  1:135180/68‑1 (MQ=255)
             tAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAGGaa                                      >  1:1195032/1‑69 (MQ=255)
                      cGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTaaa                             >  1:1690802/1‑69 (MQ=255)
                      cGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTa                               <  1:2140165/67‑1 (MQ=255)
                          cACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACtgt                         >  1:1710816/1‑69 (MQ=255)
                           aCCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACtgt                         <  1:2725188/68‑1 (MQ=255)
                              aTCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTgtggtg                     <  1:333032/69‑1 (MQ=255)
                              aTCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTgtggt                      <  1:1526194/68‑1 (MQ=255)
                                cACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGAc                   >  1:2120448/1‑69 (MQ=255)
                                cACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGAc                   >  1:824916/1‑69 (MQ=255)
                                 aCGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAACtg                          >  1:3542648/1‑61 (MQ=255)
                                 aCGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAACtg                          >  1:54836/1‑61 (MQ=255)
                                 aCGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAACtg                          >  1:2257654/1‑61 (MQ=255)
                                 aCGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAACtg                          >  1:2790667/1‑61 (MQ=255)
                                 aCGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAACtg                          <  1:1475493/61‑1 (MQ=255)
                                 aCGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACa                  >  1:378522/1‑69 (MQ=255)
                                   gCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACtgt                         <  1:3535462/60‑1 (MQ=255)
                                   gCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATc                <  1:3557303/69‑1 (MQ=255)
                                    cAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTgtggt                      >  1:911937/1‑62 (MQ=255)
                                    cAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTgtgg                       >  1:1807624/1‑61 (MQ=255)
                                    cAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGATGACAAGGAAGCTTTTAAACTgtggt                      >  1:2206673/1‑62 (MQ=255)
                                                 tCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGTTTTAAAGCt  <  1:961880/69‑1 (MQ=255)
                                                 tCTTGCGCTTTAATTTTTTACGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGTTTTAAAGCt  <  1:1304488/69‑1 (MQ=255)
                                                  cTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAACTgtggt                      <  1:2045641/48‑1 (MQ=255)
                                                  cTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAACTgtggt                      <  1:3658078/48‑1 (MQ=255)
                                                  cTTGCGCTTTAATTTTTTCCGGTGACAATGAAGCTTTTAAACTGTGGTGACATCCAg             <  1:3467983/57‑1 (MQ=255)
                                                   ttGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATcc               <  1:212496/54‑1 (MQ=255)
                                                                    |                                                 
TAAGGATAGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGTTTTAAAGCT  >  W3110S.gb/2313755‑2313872

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: