Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 181,477 T→C 9.7% G198G (GGT→GGC degP → serine endoprotease (protease Do), membrane‑associated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb181,4770TC 9.7% 100.6 / 3.3 41G198G (GGT→GGCdegPserine endoprotease (protease Do), membrane‑associated
Reads supporting (aligned to +/- strand):  ref base T (17/20);  new base C (2/2);  total (19/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

TAAGATGGCGGATTCTGATGCACTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGCCGAAAAC  >  W3110S.gb/181411‑181543
                                                                  |                                                                  
tAAGATGGCGGATTCTGAGGCACTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTgg                                                                >  1:1831906/1‑71 (MQ=255)
    aTGGCGGATTCTGATGCACTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCga                                                            <  1:378721/71‑1 (MQ=255)
     tGGCGGATTCTGATGCACTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCga                                                            <  1:191559/70‑1 (MQ=255)
      ggCGGATTCTGATGCACTGCGCGTGTGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCga                                                            <  1:47185/69‑1 (MQ=255)
          gATTCTGATGCACTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGt                                                      <  1:1092572/71‑1 (MQ=255)
          gATTCTGATGCACTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGt                                                      <  1:122811/71‑1 (MQ=255)
              cTGATGCACTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGAc                                                         >  1:2713829/1‑64 (MQ=255)
                      cTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTc                                                <  1:1323700/65‑1 (MQ=255)
                        gcgcGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTg                                        >  1:170139/1‑71 (MQ=255)
                            gTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTggg                                                               >  1:3831391/1‑44 (MQ=255)
                            gTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGtctc                                    >  1:577142/1‑71 (MQ=255)
                               ggTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTgc                                 <  1:1059852/71‑1 (MQ=255)
                                 tGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTgcgc                               <  1:3752716/71‑1 (MQ=255)
                                   aTTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTcc                                               <  1:3357908/53‑1 (MQ=255)
                                     tACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTggg                           >  1:1796357/1‑71 (MQ=255)
                                         ccGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGt                       <  1:889632/71‑1 (MQ=255)
                                            tAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTgc                                 <  1:3173072/58‑1 (MQ=255)
                                            tAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTgc                                 <  1:2719639/58‑1 (MQ=255)
                                            tAGCAATTGGTAACCCGTTCGGCCTGGGCGAGACAGTAACTTCCgg                                             <  1:2417259/46‑1 (MQ=25)
                                            tAGCAATTGGTAACCCGTTCGGCCTGGGCGAGACAGTAACTTCCgg                                             <  1:842966/46‑1 (MQ=25)
                                               cGATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGc                 <  1:3656703/71‑1 (MQ=255)
                                                  ttGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTgg                            >  1:3849395/1‑57 (MQ=255)
                                                  ttGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTgg                            >  1:695244/1‑57 (MQ=255)
                                                  ttGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTgg                            >  1:228397/1‑57 (MQ=255)
                                                  ttGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTgg                            >  1:2288504/1‑57 (MQ=255)
                                                  ttGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCt               >  1:2184672/1‑70 (MQ=255)
                                                  ttGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTg              <  1:2667050/71‑1 (MQ=255)
                                                     gTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGaa            <  1:2558634/70‑1 (MQ=255)
                                                      tAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCgg                                             <  1:1206905/36‑1 (MQ=255)
                                                       aaCCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGc         >  1:2711869/1‑71 (MQ=255)
                                                       aaCCCGTTCGGCCTGGGCGAGACAGTAACTTCCGGGATTGTCTCTGCGCTGGGACGTAGCGGCCTGAACGc         >  1:2472989/1‑71 (MQ=255)
                                                       aaCCCGTTCGGCCTGGGCGAGACAGTAACTTCCGGGATTGTCTCTGCGCTGGGACGTAGCGGCCTGAACGc         >  1:1337091/1‑71 (MQ=255)
                                                        aCCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGc                         >  1:2315034/1‑54 (MQ=255)
                                                        aCCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTACGCTGGGGc                         >  1:794684/1‑54 (MQ=255)
                                                           cGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGCCGaa     <  1:2133750/71‑1 (MQ=255)
                                                           cGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGCCGa      >  1:426476/1‑70 (MQ=255)
                                                             tttGGTCTGGGCGAGACGGTAACTTCCGGGATTGTctct                                   <  1:1692414/39‑1 (MQ=255)
                                                             tttGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTg              <  1:197667/60‑1 (MQ=255)
                                                             tttGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTg              <  1:1025532/60‑1 (MQ=255)
                                                               tGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGc         >  1:3545710/1‑63 (MQ=255)
                                                               tGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGCCGAAAAc  >  1:421294/1‑70 (MQ=255)
                                                                  |                                                                  
TAAGATGGCGGATTCTGATGCACTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGCCGAAAAC  >  W3110S.gb/181411‑181543

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: