Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 181,321 C→T 10.4% G146G (GGC→GGT 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,3210CT10.4% 129.8 / 4.7 48G146G (GGC→GGTdegPserine endoprotease (protease Do), membrane‑associated
Reads supporting (aligned to +/- strand):  ref base C (18/25);  new base T (3/2);  total (21/27)
Fisher's exact test for biased strand distribution p-value = 6.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATC  >  W3110S.gb/181254‑181384
                                                                   |                                                               
gCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGt                                                              <  1:360033/71‑1 (MQ=255)
 cTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGt                                                              <  1:3342679/70‑1 (MQ=255)
  tATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGt                                                              >  1:925897/1‑69 (MQ=255)
    tGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAg                                                           <  1:2857566/70‑1 (MQ=255)
           aCCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCg                                                       <  1:1907888/67‑1 (MQ=255)
               acaacCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGAcgcg                                                  >  1:196394/1‑68 (MQ=255)
                 aacCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGCCGCGAAGAt                                             >  1:90166/1‑71 (MQ=255)
                 aacCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGAt                                             >  1:1472912/1‑71 (MQ=255)
                 aacCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGAt                                             >  1:2290993/1‑71 (MQ=255)
                  acCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATg                                            <  1:3471866/71‑1 (MQ=255)
                    cACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGt                                          <  1:317192/71‑1 (MQ=255)
                     aCGTTGTTGATAACGCGACTGTCATTAAAGTTCAACTGAGCGATGGCCGTAAg                                                           <  1:1874611/53‑1 (MQ=255)
                          gttgATAACGCGACGGTGATTAAAGTCCAACTGAGCGATGGTCGTAAGTTCGACGCGAAGATGGTTGGCa                                     <  1:3791359/70‑1 (MQ=255)
                          gttgATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCa                                     >  1:324085/1‑70 (MQ=255)
                           ttgATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTa                                                             <  1:2990654/45‑1 (MQ=255)
                           ttgATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTa                                                             <  1:616415/45‑1 (MQ=255)
                           ttgATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGaa                                                <  1:3683389/58‑1 (MQ=255)
                           ttgATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGaa                                                <  1:2241351/58‑1 (MQ=255)
                           ttgATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGt                                          <  1:1769991/64‑1 (MQ=255)
                           ttgATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGt                                          <  1:480964/64‑1 (MQ=255)
                               tAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAGGATGGTTGGCAAAGATc                               >  1:3183268/1‑71 (MQ=255)
                                aaCGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTaa                                                            >  1:1651739/1‑41 (MQ=255)
                                   gcgACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGa                                              >  1:3298536/1‑52 (MQ=255)
                                   gcgACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGa                                              >  1:36847/1‑52 (MQ=255)
                                   gcgACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGa                                              >  1:1784485/1‑52 (MQ=255)
                                   gcgACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCcgcg                           <  1:2201299/71‑1 (MQ=255)
                                   gcgACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCcgcg                           <  1:2945699/71‑1 (MQ=255)
                                    cgACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGtt                                         >  1:1646363/1‑56 (MQ=255)
                                       cGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATc                               >  1:2834757/1‑63 (MQ=255)
                                       cGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATc                               >  1:1595330/1‑63 (MQ=255)
                                         gtgaTTAAAGTCCAACTGAGCGATGGTCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGa                     <  1:3010838/68‑1 (MQ=255)
                                          tCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGTTTGGCAAAg                                  <  1:2493645/57‑1 (MQ=255)
                                          tCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAg                                  <  1:2812566/57‑1 (MQ=255)
                                          tCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGat                    <  1:2195007/71‑1 (MQ=255)
                                           taTTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCaaa                                   >  1:1315576/2‑55 (MQ=255)
                                           cATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCaaa                                   >  1:1355468/1‑55 (MQ=255)
                                            aTTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGaa                                                <  1:3454965/41‑1 (MQ=255)
                                            aTTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGaa                                                <  1:1983799/41‑1 (MQ=255)
                                               aaaGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATcgc               >  1:462021/1‑71 (MQ=255)
                                                   tccAACTGAGCGATGGTCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATcgcg              >  1:2623947/3‑68 (MQ=255)
                                                   tccAACTGAGCGATGGTCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATcgcg              >  1:2791227/3‑68 (MQ=255)
                                                   tccAACTGAGCGATGGTCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATcgc               >  1:3851740/3‑67 (MQ=255)
                                                     cAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGAt         >  1:2178164/1‑71 (MQ=255)
                                                            gCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATc  <  1:338034/71‑1 (MQ=255)
                                                            gCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATc  <  1:2207718/71‑1 (MQ=255)
                                                            gCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATc  <  1:3242074/71‑1 (MQ=255)
                                                                tGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGa                                 <  1:1023260/36‑1 (MQ=255)
                                                                tGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGAt                                <  1:1056457/37‑1 (MQ=255)
                                                                   |                                                               
GCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATC  >  W3110S.gb/181254‑181384

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: