Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,924,649 T→C 7.0% S191S (TCT→TCC bisC → biotin sulfoxide reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,924,6490TC 7.0% 153.6 / 2.6 57S191S (TCT→TCCbisCbiotin sulfoxide reductase
Reads supporting (aligned to +/- strand):  ref base T (30/23);  new base C (2/2);  total (32/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CTAACCCACTCAATACGCTGAAAATTGCGTGGAATGCATCCGATGAGCAGGGGCTTTCTTACTTTTCTGCACTGCGTGACAGCGGGAAAAAGCTGATCTGCATTGATCCAATGCGATCGGAAACCGTCGATTTCTT  >  W3110S.gb/3924582‑3924717
                                                                   |                                                                    
cTAACCCACTCAATACGCTGAAAATTGCGTGGAATGCATCCGATGAGCAGGGGCTTTCTTACTTTTCTGCa                                                                   <  1:3743087/71‑1 (MQ=255)
cTAACCCACTCAATACGCTGAAAATTGCGTGGAATGCATCCGATGAGCAGGGGCTTTCTTACTTTTCTGCa                                                                   >  1:3620678/1‑71 (MQ=255)
       aCTCAATACGCTGAAAATTGCGTGGAATGCATCCGATGAGCAGGGGCTTTCTTACTTTTCTGCACTGCGt                                                             >  1:1286281/1‑70 (MQ=255)
         tCAATACGCTGAAAATTGCGTGGAATGCATCCGATGAGCAGGGGCTTTCTTACTTTTCTg                                                                     <  1:1172088/60‑1 (MQ=255)
          cAATACGCTGAAAATTGCGTGGAATGCATCCGATGAGCAGGGGCTTTCTTACTTTTCTGCACTGCGTGa                                                           >  1:125389/1‑69 (MQ=255)
          cAATACGCTGAAAATTGCGTGGAATGCATCCGATGAGCAGGGGCTTTCTTACTTTTCTGCACTGCGTGa                                                           >  1:637418/1‑69 (MQ=255)
          cAATACGCTGAAAATTGCGTGGAATGCATCCGATGAGCAGGGGCTTTCTTACTTTTCTGCACTGCGTGa                                                           >  1:875854/1‑69 (MQ=255)
           aaTACGCTGAAAATTGCGTGGAATGCATCCGATGAGCAGGGGCTTTCTTACTTTTCTGCAc                                                                  >  1:710315/1‑61 (MQ=255)
           aaTACGCTGAAAATTGCGTGGAATGCATCCGATGAGCAGGGGCTTTCTTACTTTTCTGCACTGCGTGACAg                                                        >  1:229161/1‑71 (MQ=255)
           aaTACGCTGAAAATTGCGTGGAATGCATCCGATGAGCAGGGGCTTTCTTACTTTTCTGCACTGCGTGACAg                                                        >  1:2091586/1‑71 (MQ=255)
              aCGCTGAAAATTGCGTGGAATGCATCCGATGAGCAGGGGCTTTCTTACTTTTCTGCACTGCGTGACAGCgg                                                     >  1:1054019/1‑71 (MQ=255)
              aCGCTGAAAATTGCGTGGAATACATCCGATGAGCAGGGGCTTTCTTACTTTTCTGCACTGCGTGACAGCgg                                                     >  1:3794054/1‑71 (MQ=255)
               cGCTGAAAATTGCGTGGAATGCATCCGATGAGCAGGGGCTTTCTTACTTTTCTGCACTGCGTGACAGCgg                                                     >  1:3020840/1‑70 (MQ=255)
                 cTGAAAATTGCGTGGAATGCATCCGATGAGCAGGGGCTTTCTTACTTTTCTGCACTGCGTGACAGCGGGaa                                                  <  1:603554/71‑1 (MQ=255)
                  tGAAAATTGCGTGGAATGCATCCGATGAGCAGGGGCTTTCTTACTTTTCCGCACTGCg                                                              <  1:268783/58‑1 (MQ=255)
                  tGAAAATTGCGTGGAATGCATCCGATGAGCAGGGGCTTTCTTACTTTTCCGCACTGCg                                                              <  1:3619763/58‑1 (MQ=255)
                    aaaaTTGCGTGGAATGCATCCGATGAGCAGGGGCTTTCTTACTTTTCTGCa                                                                   >  1:879699/1‑51 (MQ=255)
                    aaaaTTGCGTGGAATGCATCCGATGAGCAGGGGCTTTCTTACTTTTCTGCACTGCGTGACAGCGGGaaaa                                                <  1:18898/70‑1 (MQ=255)
                    aaaaTTGCGTGGAATGCATCCGATGAGCAGGGGCTTTCTTACTTTTCTGCACTGCGTGACAGCGGGaaa                                                 <  1:21378/69‑1 (MQ=255)
                    aaaaTTGCGTGGAATGCATCCGATGAGCAGGGGCTTTCTTACTTTTCTGCACTGCGTGACAGCGGGaaa                                                 <  1:2373050/69‑1 (MQ=255)
                      aaTTGCGTGGAATGCATCCGATGAGCAGGGGCTTTCTTACTTTTCTGCACTGCGTGACAGCGGGAAAAAGc                                             >  1:373902/1‑71 (MQ=255)
                       aTTGCGTGGAATGCATCCGATGAGCAGGGGCTTTCTTACTTTTCTGCACTGCGTGACAGCGGGAAAAAGCt                                            <  1:1291959/71‑1 (MQ=255)
                           cGTGGAATGCATCCGATGAGCAGGGGCTTTCTTACTTTTCCGCACTGCGTGACAGCGGGAAAAAGCTGATc                                        >  1:1599275/1‑71 (MQ=255)
                           cGTGGAATGCATCCGATGAGCAGGGGCTTTCTTACTTTTCCGCACTGCGTGACAGCGGGAAAAAGCTGATc                                        >  1:2722735/1‑71 (MQ=255)
                            gTGGAATGCATCCGATGAGCAGGGGCTTTCTTACTTTTCTGCACTGCGTGACAGCGGGaa                                                  <  1:3319756/60‑1 (MQ=255)
                                  tGCATCCGATGAGCAGGGGCTTTCTTACTTTTCTGCACTGCGTGACAGCGGGAAAAAGCTGATCTGCAtt                                  >  1:1590588/1‑70 (MQ=255)
                                  tGCATCCGATGAGCAGGGGCTTTCTTACTTTTCTGCACTGCGTGACAGCGGGAAAAAGCTGATCTGCAtt                                  >  1:1390619/1‑70 (MQ=255)
                                  tGCATCCGATGAGCAGGGGCTTTCTTACTTTTCTGCACTGCGTGACAGCGGGAAAAAGCTGATCTGCATTg                                 <  1:1434412/71‑1 (MQ=255)
                                       ccGATGAGCAGGGGCTTTCTTACTTTTCTGCACTGCGTGACAGCGGGAAAAAGCTGATCTGCATTGATcc                             <  1:3204386/70‑1 (MQ=255)
                                        cGATGAGCAGGGGCTTTCTTACTTTTCTGCACTGCGTGa                                                           >  1:3819227/1‑39 (MQ=255)
                                        cGATGAGCAGGGGCTTTCTTACTTTTCTGCACTGCGTGACAGCGGGAAAAAGCTGATCTGCa                                    >  1:910890/1‑62 (MQ=255)
                                              gCAGGGGCTTTCTTACTTTTCTGCACTGCGTGACAGCGGGAAAAAGCTGATCTGCATTGATCCAATGCGAt                     <  1:1286081/71‑1 (MQ=255)
                                                aGGGGCTTTCTTACTTTTCTGCACTGCGTGACAGCGGGAAAAAGCTGATCTGCATTGATCCAATGCGATCg                   <  1:40409/71‑1 (MQ=255)
                                                aGGGGCTTTCTTACTTTTCTGCACTGCGTGACAGCGGGAAAAAGCTGATCTGCATTGATCCAATGCGATCg                   <  1:1825373/71‑1 (MQ=255)
                                                  gggCTTTCTTACTTTTCTGCACTGCGTGACAGCGGGAAAAAGCTGATCTGCATTGATCCAATGCGATCgg                  >  1:1602765/1‑70 (MQ=255)
                                                  gggCTTTCTTACTTTTCTGCACTGCGTGACAGCGGGAAAAAGCTGATCTGCATTGATCCAATGCGATCGGa                 <  1:347204/71‑1 (MQ=255)
                                                    gCTTTCTTACTTTTCTGCACTGCGTGACAGCGGGaaa                                                 <  1:1866329/37‑1 (MQ=255)
                                                     cTTTCTTACTTTTCTGCACTGCGTGACAGCGGGAAAAAGCTGATCTGCATTGATCCAATGCGATCGGa                 >  1:2530089/1‑68 (MQ=255)
                                                     cTTTCTTACTTTTCTGCACTGCGTGACAGCGGGAAAAAGCTGATCTGCATTGATCCAATGCGATCGGAAAc              >  1:3595556/1‑71 (MQ=255)
                                                     cTTTCTTACTTTTCTGCACTGCGTGACAGCGGGAAAAAGCTGATCTGCATTGATCCAATGCGATCGGAAAc              >  1:976129/1‑71 (MQ=255)
                                                       ttCTTACTTTTCTGCACTGCGTGACAGCGGGAAAAAGCTg                                           <  1:3308797/40‑1 (MQ=255)
                                                        tCTTACTTTTCTGCACTGCGTGACAGCGGGAAAAAGCTGATc                                        >  1:1110704/1‑42 (MQ=255)
                                                         cTTACTTTTCTGCACTGCGTGACAGCGGGAAAAAGCTGAt                                         >  1:1198452/1‑40 (MQ=255)
                                                          ttACTTTTCTGCACTGCGTGACAGCGGGAAAAAGCTGATCTg                                      <  1:1658596/42‑1 (MQ=255)
                                                          ttACTTTTCTGCACTGCGTGACAGCGGGAAAAAGCTGATCTGCATTGATCCAATGCGATCGGAAACCGTc          >  1:663991/1‑70 (MQ=255)
                                                            aCTTTTCTGCACTGCGTGACAGCGGGAAAAAGCTGATCTGCa                                    >  1:911113/1‑42 (MQ=255)
                                                            aCTTTTCTGCACTGCGTGACAGCGGGAAAAAGCTGATCTGCATTGATTCAATGCGATCGGAAACCGTCGAt       >  1:598430/1‑71 (MQ=255)
                                                            aCTTTTCTGCACTGCGTGACAGCGGGAAAAAGCTGATCTGCATTGATCCAATGCGATCGGAAACCGTCg         <  1:2592041/69‑1 (MQ=255)
                                                            aCTTTTCTGCACTGCGTGACAGCGGGAAAAAGCTGATCTGCATTGATCCAATGCGATCGGAAACCGTCg         <  1:2363061/69‑1 (MQ=255)
                                                            aCTTTTCTGCACTGCGTGACAGCGGGAAAAAGCTGATCTGCATTGATCCAATGCGATCGGAAACCGTCGAt       >  1:2697367/1‑71 (MQ=255)
                                                            aCTTTTCTGCACTGCGTGACAGCGGGAAAAAGCTGATCTGCATTGATCCAATGCGATCGGAAACCGTCGAt       >  1:1658473/1‑71 (MQ=255)
                                                            aCTTTTCTGCACTGCGTGACAGCGGGAAAAAGCTGATCTGCATTGATCCAATGCGATCGGAAACCGTCGAt       >  1:977743/1‑71 (MQ=255)
                                                              ttttCTGCACTGCGTGACAGCGGGAAAAAGCTGATCTGCATTGATCCAATGCGATCGGAAACCGTCGAt       <  1:3064603/69‑1 (MQ=255)
                                                              ttttCTGCACTGCGTGACAGCGGGAAAAAGCTGATCTGCATTGATCCAATGCGATCGGAAACCGTCGAt       <  1:2850073/69‑1 (MQ=255)
                                                              ttttCTGCACTGCGTGACAGCGGGAAAAAGCTGATCTGCATTGATCCAATGCGATCGGAAACCGTCGAt       <  1:857409/69‑1 (MQ=255)
                                                                ttCTGCACTGCGTGACAGCGGGAAAAAGCTGATCTGCATTGATCCAATGCGATCGGAAACCGTCGATTTCt   >  1:3074618/1‑71 (MQ=255)
                                                                  cTGCACTGCGTGACAGCGGGAAAAAGCTGATCTGCATTGATCCAATGCGATCGGAAACCGTCGATTTCtt  <  1:1392883/70‑1 (MQ=255)
                                                                   |                                                                    
CTAACCCACTCAATACGCTGAAAATTGCGTGGAATGCATCCGATGAGCAGGGGCTTTCTTACTTTTCTGCACTGCGTGACAGCGGGAAAAAGCTGATCTGCATTGATCCAATGCGATCGGAAACCGTCGATTTCTT  >  W3110S.gb/3924582‑3924717

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: