Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,609,568 T→C 17.5% V69V (GTT→GTC hyfI → hydrogenase 4, Fe‑S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,609,5680TC17.5% 145.5 / 25.1 63V69V (GTT→GTChyfIhydrogenase 4, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base T (21/31);  new base C (4/7);  total (25/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

AAATCGAAATTTTTGCTGCCATTACACCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTG  >  W3110S.gb/2609504‑2609635
                                                                |                                                                   
aaaTCGAAATTTTTGCTGCCATTACACCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTCtcat                                                                 <  1:2557622/69‑1 (MQ=255)
aaaTCGAAATTTTTGCTGCCATTACACCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTCtcat                                                                 <  1:2455890/69‑1 (MQ=255)
  aTCGAAATTTTTGCTGCCATTACACCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTTTCtcat                                                                 <  1:476721/67‑1 (MQ=255)
     gAAATTTTTGCTGCCATTACACCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTtcatc                                                                >  1:885453/1‑65 (MQ=255)
     gAAATTTTTGCTGCCATTACACCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTtcatc                                                                >  1:3454325/1‑65 (MQ=255)
     gAAATTTTTGCTGCCATTACACCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTtcatc                                                                >  1:289200/1‑65 (MQ=255)
               ctgcCATTACACCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGcc                                                   <  1:1327972/68‑1 (MQ=255)
                tgcCATTACACCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGcc                                                   <  1:142531/67‑1 (MQ=255)
                     ttACACCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATAttt                                            <  1:1022977/69‑1 (MQ=255)
                     ttACACCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATAttt                                            <  1:1450623/69‑1 (MQ=255)
                     ttACACCAGTATTCGACGCAGAACGTTTTGGCATTAAAGTTGTTTCATCACCGCGTCACGCCGATAttt                                            <  1:3392070/69‑1 (MQ=255)
                       acacCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGttt                                                                    <  1:851677/43‑1 (MQ=255)
                       acacCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATAtttt                                           >  1:3221815/1‑68 (MQ=255)
                       acacCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATAtttt                                           >  1:2339983/1‑68 (MQ=255)
                       acacCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTCtcat                                                                 <  1:2995340/46‑1 (MQ=255)
                       acacCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTCtcat                                                                 <  1:2697548/46‑1 (MQ=255)
                            aGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTa                                       <  1:3468880/67‑1 (MQ=255)
                               aTTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTAc                                  >  1:784552/1‑69 (MQ=255)
                                ttCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACt                                 >  1:2645469/1‑69 (MQ=255)
                                ttCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACt                                 >  1:3604487/1‑69 (MQ=255)
                                ttCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACt                                 >  1:2082024/1‑69 (MQ=255)
                                ttCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACt                                 >  1:2162647/1‑69 (MQ=255)
                                ttCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACt                                 >  1:2225702/1‑69 (MQ=255)
                                ttCGACGCAGAACGTTTTGGCATTAAGGTTGTCtcatc                                                                >  1:1392368/1‑38 (MQ=255)
                                ttCGACGCAGAACGTTTTGGCATTAAGGTTGTCtcatc                                                                >  1:819569/1‑38 (MQ=255)
                                ttCGACGCAGAACGTTTTGGCATTAAGGTTGTCtcat                                                                 >  1:1905536/1‑37 (MQ=255)
                                ttCGACGCAGAACGTTTAGGCATTAAGGTTGTCtcatc                                                                >  1:3217948/1‑38 (MQ=38)
                                  cGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTgg                               <  1:3615044/69‑1 (MQ=255)
                                  cGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTgg                               <  1:3285248/69‑1 (MQ=255)
                                  cGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTgg                               <  1:1137784/69‑1 (MQ=255)
                                  cGACGCAGAACGTTTTGGCATTAAGGTTGTCTCATCACCGCGTCACGCCGATATTTTGt                                         <  1:361642/59‑1 (MQ=255)
                                   gACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTgg                               >  1:433859/1‑68 (MQ=255)
                                     cGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGgcgc                            <  1:67951/69‑1 (MQ=255)
                                     cGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGgcgc                            <  1:3517225/69‑1 (MQ=255)
                                     cGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGgcgc                            <  1:2586128/69‑1 (MQ=255)
                                     cGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGgcgc                            <  1:3096621/69‑1 (MQ=255)
                                     cGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGgcgc                            <  1:365426/69‑1 (MQ=255)
                                      gCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGgcgc                            <  1:719876/68‑1 (MQ=255)
                                      gCAGAACGTTTTGGCATTAAGGTTGTCTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCa                           <  1:1008378/69‑1 (MQ=255)
                                        aGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGt                         <  1:1304423/69‑1 (MQ=255)
                                           aCGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTa                                       <  1:1428552/52‑1 (MQ=255)
                                           aCGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCAc                      <  1:3260357/69‑1 (MQ=255)
                                           aCGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCAc                      >  1:3258331/1‑69 (MQ=255)
                                             gTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCa                       <  1:1131858/66‑1 (MQ=255)
                                             gTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCa                       <  1:89926/66‑1 (MQ=255)
                                                  ggCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGc                <  1:1277820/68‑1 (MQ=255)
                                                  ggCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGc                <  1:1629957/68‑1 (MQ=255)
                                                  ggCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGc                <  1:2324891/68‑1 (MQ=255)
                                                  ggCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGc                <  1:1745151/68‑1 (MQ=255)
                                                     aTTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATg            >  1:1686581/1‑69 (MQ=255)
                                                       tAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCg          <  1:3470180/69‑1 (MQ=255)
                                                       tAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCg          <  1:3716177/69‑1 (MQ=255)
                                                        aaGGTTGTTTCATCACCGCGTCACGCCGATATTTTGtt                                        <  1:13260/38‑1 (MQ=255)
                                                        aaGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCa                           >  1:3405955/1‑51 (MQ=255)
                                                        aaGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCa                           <  1:281001/51‑1 (MQ=255)
                                                        aaGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGc                >  1:331882/1‑62 (MQ=255)
                                                        aaGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGc                >  1:1113470/1‑62 (MQ=255)
                                                        aaGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCg          >  1:1669459/1‑68 (MQ=255)
                                                        aaGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCg          <  1:1142441/68‑1 (MQ=255)
                                                        aaGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGt         >  1:1300519/1‑69 (MQ=255)
                                                        aaGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGt         >  1:818596/1‑69 (MQ=255)
                                                        aaGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTAATTACTGGCGCAGTCACCCGGGCGATGCGt         >  1:2258063/1‑69 (MQ=255)
                                                          ggTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTa        <  1:2886980/68‑1 (MQ=255)
                                                          ggTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTa        <  1:3671398/68‑1 (MQ=255)
                                                              gttTCATCACCGCGTCACGCCGATATTTTGTTAttt                                    >  1:174435/1‑36 (MQ=255)
                                                                ttCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTg  <  1:981471/68‑1 (MQ=255)
                                                                |                                                                   
AAATCGAAATTTTTGCTGCCATTACACCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTG  >  W3110S.gb/2609504‑2609635

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: