Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,301,561 G→A 24.5% D429D (GAC→GAT fdhF ← formate dehydrogenase‑H, selenopolypeptide subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,301,5610GA24.5% 93.1 / 28.7 49D429D (GAC→GATfdhFformate dehydrogenase‑H, selenopolypeptide subunit
Reads supporting (aligned to +/- strand):  ref base G (12/25);  new base A (2/10);  total (14/35)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

CATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTG  >  W3110S.gb/4301494‑4301629
                                                                   |                                                                    
cATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTcc                                                                   >  1:1970494/1‑71 (MQ=255)
cATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTcc                                                                   >  1:227414/1‑71 (MQ=255)
cATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTcc                                                                   >  1:1327687/1‑71 (MQ=255)
cATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTcc                                                                   >  1:1337641/1‑71 (MQ=255)
cATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTcc                                                                   >  1:571442/1‑71 (MQ=255)
  tGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTg                                                                 <  1:412802/71‑1 (MQ=255)
  tGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTg                                                                 <  1:1688316/71‑1 (MQ=255)
  tGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTg                                                                 <  1:1185168/71‑1 (MQ=255)
  tGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTg                                                                 <  1:1246938/71‑1 (MQ=255)
  tGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTg                                                                 <  1:825698/71‑1 (MQ=255)
  tGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTg                                                                 <  1:1368213/71‑1 (MQ=255)
       gCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAAc                                                              <  1:1843014/69‑1 (MQ=255)
                  ttGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAAcc                                                       <  1:1366027/65‑1 (MQ=255)
                  ttGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAAcc                                                       <  1:1812299/65‑1 (MQ=255)
                  ttGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAAcc                                                       <  1:1961421/65‑1 (MQ=255)
                  ttGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAAcc                                                       <  1:58832/65‑1 (MQ=255)
                    gACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTc                                                  >  1:1919853/1‑68 (MQ=255)
                    gACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTc                                                  >  1:1387296/1‑68 (MQ=255)
                    gACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTc                                                  >  1:1175911/1‑68 (MQ=255)
                    gACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTc                                                  >  1:116864/1‑68 (MQ=255)
                     aCGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAg                                               <  1:653263/70‑1 (MQ=255)
                     aCGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAg                                               <  1:849737/70‑1 (MQ=255)
                     aCGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATATCCTGAACGATAACCAGTTCCAg                                               <  1:390978/70‑1 (MQ=255)
                     aCGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATATCCTGAACGATAACCAGTTCCAg                                               <  1:1899036/70‑1 (MQ=255)
                     aCGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATATCCTGAACGATAACCAGTTCCAg                                               <  1:484194/70‑1 (MQ=255)
                     aCGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATATCCTGAACGATAACCAGTTCCAg                                               <  1:168768/70‑1 (MQ=255)
                     aCGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATATCCTGAACGATAACCAGTTCCAg                                               <  1:532422/70‑1 (MQ=255)
                     aCGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATATCCTGAACGATAACCAGTTCCAg                                               <  1:1532908/70‑1 (MQ=255)
                     aCGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATATCCTGAACGATAACCAGTTCCAg                                               <  1:1530413/70‑1 (MQ=255)
                     aCGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATATCCTGAACGATAACCAGTTCCAg                                               <  1:1489510/70‑1 (MQ=255)
                     aCGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATATCCTGAACGATAACCAGTTCCAg                                               <  1:1137450/70‑1 (MQ=255)
                     aCGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATATCCTGAACGATAACCAGTACCAg                                               <  1:1124943/70‑1 (MQ=255)
                            aaTAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAAc                                                        <  1:1601519/54‑1 (MQ=255)
                                     cgccgccgACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTa                               <  1:851219/70‑1 (MQ=255)
                                     cgccgccgACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTa                               <  1:396958/70‑1 (MQ=255)
                                     cgccgccgACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTa                               <  1:1017355/70‑1 (MQ=255)
                                     cgccgccgACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTa                               <  1:262435/70‑1 (MQ=255)
                                     cgccgccgACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTa                               <  1:1857017/70‑1 (MQ=255)
                                     cgccgccgACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTa                               <  1:1295052/70‑1 (MQ=255)
                                     cgccgccgACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTa                               <  1:734631/70‑1 (MQ=255)
                                     cgccgccgACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTa                               <  1:1624141/70‑1 (MQ=255)
                                     cgccgccgACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTa                               <  1:586034/70‑1 (MQ=255)
                                      gccgccgACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACg                             <  1:76551/71‑1 (MQ=255)
                                      gccgccgACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACg                             <  1:547938/71‑1 (MQ=255)
                                               gcgGTTTTGGTCATAAAGATATCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGa                    >  1:980315/1‑71 (MQ=255)
                                               gcgGTTTTGGTCATAAAGATATCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGa                    >  1:524719/1‑71 (MQ=255)
                                                             aaaGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCa       >  1:566071/1‑70 (MQ=255)
                                                             aaaGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCa       >  1:412557/1‑70 (MQ=255)
                                                             aaaGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAg      >  1:357872/1‑71 (MQ=255)
                                                                   gTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTg  >  1:648759/1‑69 (MQ=255)
                                                                   gTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTg  >  1:1394097/1‑69 (MQ=255)
                                                                   gTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTg  >  1:1651206/1‑69 (MQ=255)
                                                                   gTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTg  >  1:32208/1‑69 (MQ=255)
                                                                   gTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTg  >  1:971224/1‑69 (MQ=255)
                                                                   |                                                                    
CATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTG  >  W3110S.gb/4301494‑4301629

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: