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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb9,4420TC53.1% 2.4 / 41.8 32L46S (TTA→TCA) mogpredicted molybdochelatase
Reads supporting (aligned to +/- strand):  ref base T (2/13);  new base C (4/13);  total (6/26)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGATGAAATGAGTTGCC  >  W3110S.gb/9381‑9507
                                                             |                                                                 
ccTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTCAATc                                                               >  1:1274984/1‑66 (MQ=255)
ccTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTCAATc                                                               >  1:1848901/1‑66 (MQ=255)
                 gCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTCAATCCCCGATGAGCAGGCGatcatc                                          >  1:752397/1‑70 (MQ=255)
                 gCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTCAATCCCCGATGAGCAGGCGatcatc                                          >  1:1281953/1‑70 (MQ=255)
                     aCATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCa                                     <  1:1005571/71‑1 (MQ=255)
                     aCATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCa                                     <  1:714408/71‑1 (MQ=255)
                     aCATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCa                                     <  1:694805/71‑1 (MQ=255)
                     aCATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCa                                     <  1:288218/71‑1 (MQ=255)
                     aCATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCa                                     <  1:94240/71‑1 (MQ=255)
                     aCATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTCAATCCCCGATGAGCAGGCGATCATCGAGCa                                     <  1:1571948/71‑1 (MQ=255)
                     aCATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTCAATCCCCGATGAGCAGGCGATCATCGAGCa                                     <  1:1152781/71‑1 (MQ=255)
                     aCATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTCAATCCCCGATGAGCAGGCGATCATCGAGCa                                     <  1:1851932/71‑1 (MQ=255)
                     aCATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTCAATCCCCGATGAGCAGGCGATCATCGAGCa                                     <  1:1287168/71‑1 (MQ=255)
                     aCATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTCAATCCCCGATGAGCAGGCGATCATCGAGCa                                     <  1:1425643/71‑1 (MQ=255)
                     aCATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTCAATCCCCGATGAGCAGGCGATCATCGAGCa                                     <  1:1752227/71‑1 (MQ=255)
                          ggCGCTAACCACGCCGTTTGAACTGGAAACCCGCTCAATCCCCGATGAGCAGGCGATCATCGAGCAAACGt                                <  1:1719673/71‑1 (MQ=255)
                          ggCGCTAACCACGCCGTTTGAACTGGAAACCCGCTCAATCCCCGATGAGCAGGCGATCATCGAGCAAACGt                                <  1:1518547/71‑1 (MQ=255)
                           gcgcTAACCACGCCGTTTGAACTGGAAACCCGCTCAATCCCCGATGAGCAGGCGATCATCGAGCa                                     <  1:94944/65‑1 (MQ=255)
                                        cGTTTGAACTGGAAACCCGCTCAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCtggtg                  <  1:1708325/71‑1 (MQ=255)
                                        cGTTTGAACTGGAAACCCGCTCAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCtggtg                  <  1:238671/71‑1 (MQ=255)
                                        cGTTTGAACTGGAAACCCGCTCAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCtggtg                  <  1:1072701/71‑1 (MQ=255)
                                            tGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGAt               >  1:722908/1‑70 (MQ=255)
                                            tGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGAt               >  1:1019015/1‑70 (MQ=255)
                                              aaCTGGAAACCCGCTCAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGATGa             <  1:362567/70‑1 (MQ=255)
                                                        ccGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGATGAAATGAGTTGcc  <  1:901682/71‑1 (MQ=255)
                                                        ccGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGATGAAATGAGTTGcc  <  1:980055/71‑1 (MQ=255)
                                                        ccGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGATGAAATGAGTTGcc  <  1:488957/71‑1 (MQ=255)
                                                        ccGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGATGAAATGAGTTGcc  <  1:382704/71‑1 (MQ=255)
                                                        ccGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGATGAAATGAGTTGcc  <  1:177247/71‑1 (MQ=255)
                                                        ccGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGATGAAATGAGTTGcc  <  1:170965/71‑1 (MQ=255)
                                                        ccGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGATGAAATGAGTTGcc  <  1:1013020/71‑1 (MQ=255)
                                                         cGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGATGAAATGAGTTGcc  <  1:230649/70‑1 (MQ=255)
                                                             |                                                                 
CCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGATGAAATGAGTTGCC  >  W3110S.gb/9381‑9507

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: