Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 9,442 T→C 28.3% L46S (TTA→TCA)  mog → predicted molybdochelatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb9,4420TC28.3% 67.2 / 33.9 46L46S (TTA→TCA) mogpredicted molybdochelatase
Reads supporting (aligned to +/- strand):  ref base T (18/15);  new base C (7/6);  total (25/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGATGAAATG  >  W3110S.gb/9375‑9500
                                                                   |                                                          
ggCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTa                                                           <  1:3350183/69‑1 (MQ=255)
   aTCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATc                                                        <  1:778355/69‑1 (MQ=255)
      ccTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATcccc                                                     >  1:1153999/1‑69 (MQ=255)
      ccTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATcccc                                                     >  1:3689850/1‑69 (MQ=255)
      ccTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATcccc                                                     >  1:564334/1‑69 (MQ=255)
      ccTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTCAATcccc                                                     >  1:3191207/1‑69 (MQ=255)
       cTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATcccc                                                     <  1:619688/68‑1 (MQ=255)
                     tGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATcccc                                                     >  1:218400/1‑54 (MQ=255)
                     tGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCg                                                    >  1:1650940/1‑55 (MQ=255)
                     tGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGatc                                      >  1:3661755/1‑69 (MQ=255)
                     tGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGatc                                      <  1:2071971/69‑1 (MQ=255)
                     tGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGatc                                      <  1:1922425/69‑1 (MQ=255)
                     tGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGatc                                      <  1:1117584/69‑1 (MQ=255)
                      ggCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGatc                                      <  1:340607/68‑1 (MQ=255)
                      ggCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGatc                                      <  1:641997/68‑1 (MQ=255)
                      ggCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGatc                                      <  1:483625/68‑1 (MQ=255)
                             aTCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTCAATCCCCGATGAGCAGGCGATCATCGAGCa                              <  1:1995228/69‑1 (MQ=255)
                                ggCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAg                                               >  1:129188/1‑49 (MQ=255)
                                ggCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGat                                       >  1:1322087/1‑57 (MQ=255)
                                ggCGCTAACCACGCCGTTTGAACTGGAAACCCGCTCAATCCCCGATGAGCAGGCGATCATCGAGCAAAc                           >  1:2702480/1‑69 (MQ=255)
                                ggCGCTAACCACGCCGTTTGAACTGGAAACCCGCTCAATCCCCGATGAGCAGGCGATCATCGAGCAAAc                           >  1:2519611/1‑69 (MQ=255)
                                ggCGCTAACCACGCCGTTTGAACTGGAAACCCGCTCAATCCCCGATGAGCAGGCGATCATCGAGCAAAc                           >  1:1797207/1‑69 (MQ=255)
                                ggCGCTAACCACGCCGTTTGAACTGGAAACCCGCTCAATACCCGATGAGCAGGCGATCATCGa                                 >  1:3714215/1‑63 (MQ=255)
                                  cgcTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGt                         >  1:2835431/1‑69 (MQ=255)
                                  cgcTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGt                         >  1:125308/1‑69 (MQ=255)
                                    cTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTtg                       <  1:953899/69‑1 (MQ=255)
                                       aCCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTtgt                      >  1:1863987/1‑67 (MQ=255)
                                           cGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGc                >  1:372677/1‑69 (MQ=255)
                                           cGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGc                >  1:3715155/1‑69 (MQ=255)
                                                  tGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCtggtgg          >  1:204958/1‑68 (MQ=255)
                                                  tGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCtggtg           >  1:118564/1‑67 (MQ=255)
                                                  tGAACTGGAAACCCGCTCAATCCCCGATGAGCAGGCGat                                       <  1:984997/39‑1 (MQ=255)
                                                  tGAACTGGAAACCCGCTCAATCCCCGATGAGCAGGCGat                                       >  1:2142874/1‑39 (MQ=255)
                                                  tGAACTGGAAACACGCTCAATCCCCGATGAGCAGGCGat                                       >  1:3237301/1‑39 (MQ=38)
                                                      cTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGATGa      <  1:653240/68‑1 (MQ=255)
                                                      cTGGAAACCCGCTCAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGATGa      <  1:3279411/68‑1 (MQ=255)
                                                      cTGGAAACCCGCTCAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGATGa      <  1:2737283/68‑1 (MQ=255)
                                                      cTGGAAACCCGCTCAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGATGa      <  1:2253874/68‑1 (MQ=255)
                                                      cTGGAAACCCGCTCAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGATGa      <  1:1385853/68‑1 (MQ=255)
                                                       tGGAAACCCGCTTAATCCCTGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGATGaa     >  1:666607/1‑68 (MQ=255)
                                                        ggAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGATGaa     >  1:1568749/1‑67 (MQ=255)
                                                         gAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGATGAAATg  >  1:2216239/1‑69 (MQ=255)
                                                           aaCCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGATGAAATg  <  1:3299244/67‑1 (MQ=255)
                                                           aaCCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGATGAAATg  <  1:2596363/67‑1 (MQ=255)
                                                           aaCCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGATGAAATg  <  1:504222/67‑1 (MQ=255)
                                                           aaCCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGATGAAATg  <  1:1849291/67‑1 (MQ=255)
                                                                   |                                                          
GGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGATGAAATG  >  W3110S.gb/9375‑9500

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: