Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,233,407 T→C 12.9% G375G (GGA→GGG ygjO ← predicted methyltransferase small domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,233,4070TC12.9% 79.7 / 4.8 33G375G (GGA→GGGygjOpredicted methyltransferase small domain
Reads supporting (aligned to +/- strand):  ref base T (19/9);  new base C (2/2);  total (21/12)
Fisher's exact test for biased strand distribution p-value = 5.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAA  >  W3110S.gb/3233351‑3233472
                                                        |                                                                 
ccGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGt                                                                   <  1:1689259/57‑1 (MQ=255)
ccGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGGCCCAGCTTCACCGc                                                     <  1:172779/71‑1 (MQ=255)
  ggATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCtt                                                   >  1:15499/1‑71 (MQ=255)
   gATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTc                                                                  >  1:678396/1‑55 (MQ=255)
   gATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCttt                                                  <  1:1147051/71‑1 (MQ=255)
   gATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCttt                                                  <  1:1568638/71‑1 (MQ=255)
        ggCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAg                                               <  1:227046/69‑1 (MQ=255)
         gCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTc                                                 <  1:1686596/66‑1 (MQ=255)
          cGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAg                                               >  1:1173263/1‑67 (MQ=255)
           gTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACca                                          >  1:691863/1‑71 (MQ=255)
                  ccgcATCCGGCAGGTAACTCCCAATCTTACCGACGACGCCCCAGCTTCACCGCTTTCAg                                               >  1:1004500/1‑59 (MQ=255)
                   cgcATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCAc                                                        >  1:1183002/1‑49 (MQ=255)
                   cgcATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCAc                                                        >  1:923406/1‑49 (MQ=255)
                   cgcATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCAc                                                        >  1:40267/1‑49 (MQ=255)
                    ccATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCAc                                                        >  1:865759/2‑48 (MQ=255)
                       tCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGAACCACaa                                       >  1:973119/1‑62 (MQ=255)
                         cGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATtc                            >  1:87612/1‑71 (MQ=255)
                            cAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCttt                                                  <  1:1197404/46‑1 (MQ=255)
                                  aCTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCa                                                >  1:1346693/1‑42 (MQ=255)
                                     cccAATCTTACCGACGACGCCCCAGCTTCACCGCTTTc                                                 <  1:500582/38‑1 (MQ=255)
                                      ccAATCTTACCGACGACGCCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAgtgg               >  1:931164/1‑71 (MQ=255)
                                          tCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATa                   >  1:1344350/1‑63 (MQ=255)
                                           cTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAg          >  1:1186458/1‑71 (MQ=255)
                                            ttACCGACGACGTCCCAGCTTCACCGCTTTCAGTACcaca                                        >  1:1084070/1‑40 (MQ=255)
                                            ttACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGt         >  1:669530/1‑71 (MQ=255)
                                            ttACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGt         >  1:175961/1‑71 (MQ=255)
                                              aCCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAg          <  1:740531/68‑1 (MQ=255)
                                              aCCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGt         <  1:39170/69‑1 (MQ=255)
                                              aCCGACGACGCCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGTAg          <  1:662485/68‑1 (MQ=255)
                                                 gacgacGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGcc     >  1:946862/1‑70 (MQ=255)
                                                   cgacgTCCCAGCTTCACCGCTTTCAGTACCACAAAttt                                   >  1:211547/1‑38 (MQ=255)
                                                   cgacgTCCCAGCTTCACCGCTTTCAGTACCACAAAttt                                   >  1:896032/1‑38 (MQ=255)
                                                   cgacgTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGaa  <  1:1435638/71‑1 (MQ=255)
                                                     acgTCCCAGCTTCAGCGCTTTGAGTACCACAAATTTATTATTCGTCGCAATAgtgg               <  1:1696879/56‑1 (MQ=255)
                                                        |                                                                 
CCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAA  >  W3110S.gb/3233351‑3233472

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: