Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,850 T→G 25.0% F373V (TTC→GTC)  thrC → threonine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,8500TG25.0% 41.8 / 27.6 40F373V (TTC→GTC) thrCthreonine synthase
Reads supporting (aligned to +/- strand):  ref base T (15/15);  new base G (6/4);  total (21/19)
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.83e-01

CGGAGCCGCACGCTGCCGTAGCTTATCGTGCGCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCT  >  W3110S.gb/4782‑4920
                                                                    |                                                                      
ctgAGCCGCACGCTGCCGTAGCTTATCGTGCTCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGt                                                                        <  1:2363989/67‑1 (MQ=255)
 ggAGCCGCACGCTGCCGTAGCTTATCGTGCGCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTc                                                                      <  1:1543263/70‑1 (MQ=255)
               ccGTAGCTTATCGTGCGCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCAt                                                       <  1:1102245/71‑1 (MQ=255)
                    gCTTATCGTGCGCTTCGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGc                                                  >  1:2389867/1‑71 (MQ=255)
                    gCTTATCGTGCGCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGc                                                  >  1:3489525/1‑71 (MQ=255)
                    gCTTATCGTGCGCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGGTCCTCgg                                                                 >  1:3404354/1‑56 (MQ=255)
                    gCTTATCGTGCGCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGGTCCTCgg                                                                 >  1:4572402/1‑56 (MQ=255)
                       tATCGTGCGCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATc                                                      >  1:4086254/1‑64 (MQ=255)
                        aTCGTGCGCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGaa                                               <  1:2154964/70‑1 (MQ=255)
                         tCGTGCGCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGGTCCTCGGCACCGCGCATCCgg                                                   >  1:2256919/1‑65 (MQ=255)
                         tCGTGCGCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGGTCCTCGGCACCGCGCATCCgg                                                   >  1:3034520/1‑65 (MQ=255)
                           gTGCGCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGc                                                  >  1:2701437/1‑64 (MQ=255)
                                cTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTa                                          <  1:3246719/67‑1 (MQ=255)
                                   cGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATc                                                      <  1:1327608/52‑1 (MQ=255)
                                   cGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAagagag                                   <  1:742495/71‑1 (MQ=255)
                                       aTCAGTTGAATCCAGGCGAATATGGCTTGGTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTg                               <  1:3826775/71‑1 (MQ=255)
                                        tCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAAtt                                            >  1:4135130/1‑57 (MQ=255)
                                            ttGAATCCAGGCGAATATGGCTTGGTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGc                          <  1:2513836/71‑1 (MQ=255)
                                              gAATCCAGGCGAATATGGCTTGGTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAg                           <  1:1116254/68‑1 (MQ=255)
                                                aTCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGAt                       >  1:2251480/1‑70 (MQ=255)
                                                   cAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATtct                    <  1:3322467/70‑1 (MQ=255)
                                                   cAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATtct                    <  1:2073582/70‑1 (MQ=255)
                                                    aGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCg                  <  1:2952984/71‑1 (MQ=255)
                                                      gCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTa                                          <  1:2814279/45‑1 (MQ=255)
                                                       cGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGa                        >  1:3187360/1‑62 (MQ=255)
                                                        gAATATGGCTTGGTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTg               >  1:902775/1‑70 (MQ=255)
                                                          atatGGCTTGGTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGaaa            <  1:75809/71‑1 (MQ=255)
                                                           tatGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTaa                                         <  1:3906844/41‑1 (MQ=255)
                                                           tatGGCTTGGTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAAc           >  1:2750164/1‑71 (MQ=255)
                                                           tatGGCTTATTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAAc           >  1:187923/1‑71 (MQ=255)
                                                             tGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTa                                          >  1:2255422/1‑38 (MQ=255)
                                                             tGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGAtt                      >  1:1303070/1‑58 (MQ=255)
                                                             tGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCgg                 >  1:285401/1‑63 (MQ=255)
                                                             tGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGt         >  1:728769/1‑71 (MQ=255)
                                                              ggCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGtt        >  1:180637/1‑71 (MQ=255)
                                                              ggCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGtt        >  1:4318316/1‑71 (MQ=255)
                                                                cTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTg                               <  1:1454482/46‑1 (MQ=255)
                                                                cTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTgg      <  1:3260754/71‑1 (MQ=255)
                                                                cTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGAGAAACGTTgg      <  1:932127/71‑1 (MQ=255)
                                                                 ttGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTg                               <  1:1608035/45‑1 (MQ=255)
                                                                  tGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGa     >  1:3817065/1‑70 (MQ=255)
                                                                    ttCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCt  <  1:3819219/71‑1 (MQ=255)
                                                                    |                                                                      
CGGAGCCGCACGCTGCCGTAGCTTATCGTGCGCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCT  >  W3110S.gb/4782‑4920

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: