Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,955,291 G→A 9.5% A750A (GCC→GCT ptr ← protease III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,955,2910GA 9.5% 115.5 / 2.8 42A750A (GCC→GCTptrprotease III
Reads supporting (aligned to +/- strand):  ref base G (13/25);  new base A (2/2);  total (15/27)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

GGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAA  >  W3110S.gb/2955228‑2955357
                                                               |                                                                  
ggTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCttttt                                                             >  1:2634083/1‑71 (MQ=255)
    tATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCaaa                                                         <  1:3589634/71‑1 (MQ=255)
      ttCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCaaa                                                         <  1:3496313/69‑1 (MQ=255)
        cATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGAt                                                      <  1:1054435/70‑1 (MQ=255)
        cATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGAt                                                      <  1:2812270/70‑1 (MQ=255)
        cATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGAt                                                      <  1:729317/70‑1 (MQ=255)
        cATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGAt                                                      <  1:3342925/70‑1 (MQ=255)
        cATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGAt                                                      <  1:585794/70‑1 (MQ=255)
         aTCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATg                                                     <  1:265535/70‑1 (MQ=255)
         aTCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGa                                                    >  1:286390/1‑71 (MQ=255)
           cGTAGCCAGTCGGTACAAATATCGCTGCCAGTGCGGAGTCGGTGCTGTTACCAGCTTTTTCAAAGATg                                                     <  1:3731385/68‑1 (MQ=255)
           cGTAGCCAGTCGGTACAAATATCGCTGCCAGTGCGGAGTCGGTGCTGTTACCAGCTTTTTCAAAGATg                                                     <  1:3169414/68‑1 (MQ=255)
            gTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGa                                                    >  1:2407925/1‑68 (MQ=255)
                 cAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGt                                            >  1:1116246/1‑71 (MQ=255)
                  aGTCGGTACAAATATCGCTGCCAGTGCGGAGTCGGTGCTGTTACCAGCTTTTTCAAAGATGACGGATTGtt                                           >  1:505570/1‑71 (MQ=255)
                  aGTCGGTACAAATATCGCTGCCAGTGCGGAGTCGGTGCTGTTACCAGCTTTTTCAAAGATGACGGATTGtt                                           >  1:3371704/1‑71 (MQ=255)
                    tCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGtttt                                         >  1:3444058/1‑71 (MQ=255)
                          cAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATc                                    <  1:1020292/70‑1 (MQ=255)
                          cAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATc                                    >  1:3149756/1‑70 (MQ=255)
                          cAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCg                                   <  1:717820/71‑1 (MQ=255)
                                ccGCTGCCAGTGCGGAGTCTGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCAc                              <  1:865504/70‑1 (MQ=255)
                                ccGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCAc                              >  1:2167982/1‑70 (MQ=255)
                                ccGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCAc                              <  1:285899/70‑1 (MQ=255)
                                    tGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTa                            <  1:3199137/68‑1 (MQ=255)
                                    tGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTa                            <  1:930608/68‑1 (MQ=255)
                                    tGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTa                            <  1:2803858/68‑1 (MQ=255)
                                      ccAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGtttttt                                       <  1:1280645/55‑1 (MQ=255)
                                       cAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGAt                                                      <  1:2510762/39‑1 (MQ=255)
                                       cAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCtt                      <  1:2032910/71‑1 (MQ=255)
                                             ggAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTc                <  1:2225264/71‑1 (MQ=255)
                                             ggAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTc                <  1:1777816/71‑1 (MQ=255)
                                                gTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACt                             <  1:1929187/55‑1 (MQ=255)
                                                 tCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGAc                                 <  1:1356620/50‑1 (MQ=255)
                                                   ggTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGacac           >  1:1569699/1‑70 (MQ=255)
                                                      gCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATc                        >  1:287742/1‑54 (MQ=255)
                                                       cTGTTACCGGCTTTTTCAAAGATGACGGATTGttttt                                        >  1:2130675/1‑37 (MQ=255)
                                                         gTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATc                                    >  1:1788958/1‑39 (MQ=255)
                                                          ttACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGtt                  <  1:3727523/56‑1 (MQ=255)
                                                          ttACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTg    <  1:2589661/70‑1 (MQ=255)
                                                           tACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGaa  >  1:2300334/1‑71 (MQ=255)
                                                            aCCGGCTTTTTCAAAGATGACGGATTGTTTTTTATc                                    >  1:2804625/1‑36 (MQ=255)
                                                              cGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTc                <  1:1403911/54‑1 (MQ=255)
                                                               ggCTTTTTCAAAGATGACGGATTGTTTTTTATCGAc                                 <  1:3441592/36‑1 (MQ=255)
                                                               |                                                                  
GGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAA  >  W3110S.gb/2955228‑2955357

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: