Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,615,349 A→G 23.8% P369P (CCT→CCC ubiB ← 2‑octaprenylphenol hydroxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,615,3490AG23.8% 34.0 / 8.5 21P369P (CCT→CCCubiB2‑octaprenylphenol hydroxylase
Reads supporting (aligned to +/- strand):  ref base A (10/6);  new base G (2/3);  total (12/9)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

ATTAAACAGATTTAACAGTACATGTCCAAACGAAATTTCGGCCAGCGGTTTCTCAAAGATAGGTTCACAGACCGTACGAATGGCAAATTCGAACTCTTCAACGTTGGTATCTGGTGGCACCCAGC  >  W3110S.gb/3615289‑3615413
                                                            |                                                                
aTTAAACAGATTTAACAGTACATGTCCAAACGAAATTTCGGCCAGCGGTTTCTCAAAGATAGGTTCACAGa                                                        >  1:653029/1‑71 (MQ=255)
 ttAAACAGATTTAACAGTACATGTCCAAACGAAATTTCGGCCAGCGGTTTCTCAAAGATAGGTTCACAGAc                                                       >  1:1271875/1‑71 (MQ=255)
  tAAACAGATTTAACAGTACATGTCCAAACGAAATTTCGGCCAGCGGTTTCTCAAAGATGGGTTCACAGa                                                        >  1:940421/1‑69 (MQ=255)
    aaCAGATTTAACAGTACATGTCCAAACGAAATTTCGGCCAGCGGTTTCTCAAAGATAGGTTCACAGACCg                                                     <  1:489581/70‑1 (MQ=255)
      cAGATTTAACAGTACATGTCCAAACGAAATTTCGGCCAGCGGTTTCTCAAAGATAGGTTCACAGACCGTAc                                                  >  1:150339/1‑71 (MQ=255)
         aTTTAACAGTACATGTCCAAACGAAATTTCGGCCAGCGGTTTCTCAAAGATAGGTTCACAGACCGTACGaa                                               >  1:380096/1‑71 (MQ=255)
              aCAGTACATGTCCAAACGAAATTTCGGCCAGCGGTTTCTCAAAGATAGGTTca                                                            >  1:1694138/1‑53 (MQ=255)
                    cATGTCCAAACGAAATTTCGGCCAGCGGTTTCTCAAAGATAGGTTCACAGACCGTACGAATGGc                                           >  1:1007571/1‑64 (MQ=255)
                        tCCAAACGAAATTTCGGCCAGCGGTTTCTCAAAGATAGGTTCACAGACCGTACGAATGGCAAATTCGa                                   <  1:1369423/68‑1 (MQ=255)
                        tCCAAACGAAATTTCGGCCAGCGGTTTCTCAAAGATAGGTTCACAGACCGTACGAATGGCAAATTCGa                                   <  1:747787/68‑1 (MQ=255)
                              cGAAATTTCGGCCAGCGGTTTCTCAAAGATAGGTTCACAGACCGTACGAATGGCAAATTCGAACTCTTCaa                          <  1:432951/71‑1 (MQ=255)
                                  aTTTCGGCCAGCGGTTTCTCAAAGATGGGTTCACAGACCGTACGAATGGCAAATTCGAACTCTTCAACGtt                      <  1:990294/71‑1 (MQ=255)
                                   tttCGGCCAGCGGTTTCTCAAAGATAGGTTCACAGACCGTACGAATGGCAAATTCGAACTCTTCAACGTTg                     <  1:181324/71‑1 (MQ=255)
                                      cGGCCAGCGGTTTCTCAAAGATGGGTTCACAGACCGTACGAATGGCAAATTCGAACTCTTCAACGTTGGt                   <  1:1250504/70‑1 (MQ=255)
                                          cAGCGGTTTCTCAAAGATGGGTTCACAGACCGTACGAATgg                                            <  1:655892/41‑1 (MQ=255)
                                          cAGCGGTTTCTCAAAGATAGGTTCACAGACCGTACGAATGGCAAATTCGAACTCTTCAACGTTGGTATCtg              <  1:1093103/71‑1 (MQ=255)
                                             cGGTTTCTCAAAGATAGGTTCACAGACCGTACGAATGGCAAATTCGAACTCTTCAACGTTGGTATCtggtg           >  1:578975/1‑71 (MQ=255)
                                                    tcAAAGATGGGTTCACAGACCGTACGAATGGCAAATTCGAACTCTTCAACGTTGGTATCTGGTGGCAcc      >  1:580532/1‑69 (MQ=255)
                                                      aaaGATAGGTTCACAGACCGTACGAATGGCAAATTCGAACTCTTCAACGTTGGTATCTGGTGGCACccagc  >  1:630858/1‑71 (MQ=255)
                                                      aaaGATAGGTTCACAGACCGTACGAATGGCAAATTCGAACTCTTCAACGTTGGTATCTGGTGGCACccagc  >  1:150493/1‑71 (MQ=255)
                                                          aTAGGTTCACAGACCGTACGAATGGCAAATTCGAACTCTTCAACGTTGGTAt                 >  1:620865/1‑52 (MQ=255)
                                                            |                                                                
ATTAAACAGATTTAACAGTACATGTCCAAACGAAATTTCGGCCAGCGGTTTCTCAAAGATAGGTTCACAGACCGTACGAATGGCAAATTCGAACTCTTCAACGTTGGTATCTGGTGGCACCCAGC  >  W3110S.gb/3615289‑3615413

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: