Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,313,823 T→C 76.8% 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,8230TC76.8% 82.9 / 22.9 43E349G (GAA→GGA) yojLpredicted thiamine biosynthesis lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (6/4);  new base C (16/17);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

TTAAGGATAGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGTTTTAAAGCTATCGCC  >  W3110S.gb/2313754‑2313878
                                                                     |                                                       
ttAAGGATAGACTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACACTCTTGCGCTTTAAttttttt                                                         >  1:3448602/1‑70 (MQ=255)
  aaGGATAGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGc                                                      >  1:2234785/1‑71 (MQ=255)
  aaGGATAGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGc                                                      >  1:2468088/1‑71 (MQ=255)
        aGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGAc                                                  >  1:2889394/1‑69 (MQ=255)
              tAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGc                                          >  1:2162630/1‑71 (MQ=255)
              tAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGc                                          >  1:2331529/1‑71 (MQ=255)
              tAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGc                                          >  1:266785/1‑71 (MQ=255)
                     cccGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTaa                                    <  1:3182544/70‑1 (MQ=255)
                     cccGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTaa                                    <  1:10478/70‑1 (MQ=255)
                     cccGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTaa                                    <  1:195629/70‑1 (MQ=255)
                     cccGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTa                                     <  1:1561072/69‑1 (MQ=255)
                      ccGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAAc                                  <  1:3298038/71‑1 (MQ=255)
                        ggTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACtg                                <  1:2606158/71‑1 (MQ=255)
                            aCCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAACTgtggt                            <  1:990253/71‑1 (MQ=255)
                               aTCACGCAAAAACCAACAATCTTGCTCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTgtggt                            <  1:1567801/68‑1 (MQ=255)
                               aTCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTgtggt                            <  1:582597/68‑1 (MQ=255)
                                tCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAgg                                              <  1:3725306/49‑1 (MQ=255)
                                tCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAgg                                              <  1:139005/49‑1 (MQ=255)
                                 cACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACAt                       >  1:2700431/1‑71 (MQ=255)
                                 cACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACAt                       >  1:1580197/1‑71 (MQ=255)
                                  aCGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCt                                         >  1:3556041/1‑52 (MQ=255)
                                  aCGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATc                      >  1:756667/1‑71 (MQ=255)
                                  aCGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATc                      >  1:956518/1‑71 (MQ=255)
                                  aCGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGa                          >  1:853275/1‑67 (MQ=255)
                                       aaaaCCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAACTgtggtg                           <  1:2718590/61‑1 (MQ=255)
                                       aaaaCCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAACTgtggtg                           <  1:365739/61‑1 (MQ=255)
                                         aaCCAACAATCTTGCTCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAACTGTGGTGAc                         >  1:2885956/1‑61 (MQ=255)
                                         aaCCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAACTGTGGTGAc                         >  1:897608/1‑61 (MQ=255)
                                         aaCCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAACTGTGGTGAc                         >  1:938497/1‑61 (MQ=255)
                                         aaCCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGttt               >  1:59831/1‑71 (MQ=255)
                                         aaCCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGttt               >  1:265804/1‑71 (MQ=255)
                                                aaTCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGTTTTAAAGc         >  1:3001553/1‑70 (MQ=255)
                                                aaTCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGTTTTAAAGc         >  1:3261877/1‑70 (MQ=255)
                                                aaTCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGTTTTAAAGCt        <  1:371098/71‑1 (MQ=255)
                                                aaTCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGTTTTAAAGCt        <  1:66551/71‑1 (MQ=255)
                                                 aTCTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGTTTTAAAGCTa       <  1:504004/71‑1 (MQ=255)
                                                 aTCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGTTTTAAAGCTa       >  1:3092896/1‑71 (MQ=255)
                                                  tCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTa                                     <  1:2511026/40‑1 (MQ=255)
                                                   cTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTa                                     <  1:1787572/39‑1 (MQ=255)
                                                      gcgcTTTAATTTTTTCCGCTGACAAGGAAGCTTTTaaa                                   >  1:3262940/1‑38 (MQ=255)
                                                      gcgcTTTAATTTTTTCCGCTGACAAGGAAGCTTTTaaa                                   >  1:11664/1‑38 (MQ=255)
                                                       cgcTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGTTTTAAAGCTATCGcc  <  1:593434/70‑1 (MQ=255)
                                                       cgcTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGTTTTAAAGCTATCGcc  <  1:269172/70‑1 (MQ=255)
                                                       cgcTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGTTTTAAAGCTATCGcc  <  1:264869/70‑1 (MQ=255)
                                                                     |                                                       
TTAAGGATAGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGTTTTAAAGCTATCGCC  >  W3110S.gb/2313754‑2313878

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: