Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 181,489 G→A 8.9% T202T (ACG→ACA 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,4890GA 8.9% 126.6 / 2.6 45T202T (ACG→ACAdegPserine endoprotease (protease Do), membrane‑associated
Reads supporting (aligned to +/- strand):  ref base G (20/21);  new base A (2/2);  total (22/23)
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.97e-01

GATTCTGATGCACTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGCCGAAAACTACGAAAACTTCATC  >  W3110S.gb/181421‑181558
                                                                    |                                                                     
gATTCTGATGCACTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGt                                                                     <  1:1092572/71‑1 (MQ=255)
gATTCTGATGCACTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGt                                                                     <  1:122811/71‑1 (MQ=255)
            cTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTc                                                               <  1:1323700/65‑1 (MQ=255)
              gcgcGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTg                                                       >  1:170139/1‑71 (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:2288504/1‑57 (MQ=255)
                                        ttGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTgg                                           >  1:228397/1‑57 (MQ=255)
                                        ttGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTgg                                           >  1:695244/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:1337091/1‑71 (MQ=255)
                                             aaCCCGTTCGGCCTGGGCGAGACAGTAACTTCCGGGATTGTCTCTGCGCTGGGACGTAGCGGCCTGAACGc                        >  1:2472989/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:1025532/60‑1 (MQ=255)
                                                   tttGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTg                             <  1:197667/60‑1 (MQ=255)
                                                     tGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGc                        >  1:3545710/1‑63 (MQ=255)
                                                     tGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGCCGAAAAc                 >  1:421294/1‑70 (MQ=255)
                                                          tGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGt                                      <  1:3569028/44‑1 (MQ=255)
                                                               gagaCGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTg                             <  1:1237469/48‑1 (MQ=255)
                                                               gagaCGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGCCGaaa                   <  1:3055141/58‑1 (MQ=255)
                                                               gagaCGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGCCGAAAACTACGAAAACtt      <  1:929507/71‑1 (MQ=255)
                                                                agaCGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGCCGAAAACTACGAAAACtt      >  1:2556068/1‑70 (MQ=255)
                                                                 gaCGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGaa                           >  1:1967557/1‑48 (MQ=255)
                                                                 gaCGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGaa                           >  1:3660921/1‑48 (MQ=255)
                                                                 gaCGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGaa                           >  1:3226322/1‑48 (MQ=255)
                                                                 gaCGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGaa                           >  1:1580209/1‑48 (MQ=255)
                                                                 gaCGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGa                            >  1:3149899/1‑47 (MQ=255)
                                                                    ggTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCt                              <  1:3723187/42‑1 (MQ=255)
                                                                    ggTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCt                              <  1:2475878/42‑1 (MQ=255)
                                                                    ggTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGCCGaa                    <  1:3693462/52‑1 (MQ=255)
                                                                    ggTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGCCGAAAACTACGAAAACTTCATc  >  1:3465725/1‑70 (MQ=255)
                                                                    ggTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGCCGAAAACTACGAAAACTTCATc  >  1:3224364/1‑70 (MQ=255)
                                                                    |                                                                     
GATTCTGATGCACTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGCCGAAAACTACGAAAACTTCATC  >  W3110S.gb/181421‑181558

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: