Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,609,011 T→C 46.7% N64N (AAT→AAC hyfH → hydrogenase 4, Fe‑S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,609,0110TC46.7% ‑0.9 / 39.9 30N64N (AAT→AAChyfHhydrogenase 4, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base T (4/12);  new base C (5/9);  total (9/21)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

CGGTGCCTGCGCCTGTGCTTGTCCGGCAAATGCGCTGACTATCCAGACCGACGACCAGCAAAATTCGCGCACCTGGCAGCTCTATCTGGGGCGTTGTATTTACTGCGGACGTTGTGAAGAAGTG  >  W3110S.gb/2608948‑2609071
                                                               |                                                            
cGGTGCCTGCGCCTGTGCTTGTCCGGCAAATGCGCTGACTATCCAGACCGACGACCAGCAAAATTCgcgc                                                        <  1:5099575/70‑1 (MQ=255)
cGGTGCCTGCGCCTGTGCTTGTCCGGCAAATGCGCTGACTATCCAGACCGACGACCAGCAAAATTCgcgc                                                        <  1:4369588/70‑1 (MQ=255)
  gTGCCTGCGCCTGTGCTTGTCCGGCAAATGCGCTGACTATCCAGACCGACGACCAGCAAAATTCGCGCAcc                                                     <  1:1171408/71‑1 (MQ=255)
         cgcCTGTGCTTGTCCGGCAAATGCGCTGACTATCCAGACCGACGACCAGCAAAATTCGCGCACCTGGCAGc                                              <  1:4078122/71‑1 (MQ=255)
           ccTGTGCTTGTCCGGCAAATGCGCTGACTATCCAGACCGACGACCAGCAAAACTCGCGCACCTGGCAGctc                                            <  1:4200124/71‑1 (MQ=255)
            cTGTGCTTGTCCGGCAAATGCGCTGACTATCCAGACCGACGACCAGCAAAACTCGCGCACCTGGCAGctct                                           <  1:3853418/71‑1 (MQ=255)
               tgCTTGTCCGGCAAATGCGCTGACTATCCAGACCGACGACCAGCAAAATTCGCGCACCTGGCAGCTCTATc                                        <  1:3666731/71‑1 (MQ=255)
               tgCTTGTCCGGCAAATGCGCTGACTATCCAGACCGACGACCAGCAAAATTCGCGCACCTGGCAGCTCTATc                                        <  1:1359492/71‑1 (MQ=255)
                 cTTGTCCGGCAAATGCGCTGACTATCCAGACCGACGACCAGCAAAACTCGCGCACCTGGCAGCTCTATCTg                                      <  1:2571886/71‑1 (MQ=255)
                   tGTCCGGCAAATGCGCTGACTATCCAGACCGACGACCAGCAAAACTCGCGCACCTGGCAGCTCTATCTggg                                    >  1:408625/1‑71 (MQ=255)
                     tCCGGCAAATGCGCTGACTATCCAGACCGACGACCAGCAAAATTCGCGCACCTGGCAGCTCTATCTgggg                                   <  1:1245347/70‑1 (MQ=255)
                      ccGGCAAATGCGCTGACTATCCAGACCGACGACCAGCAAAAtt                                                             >  1:3334767/1‑43 (MQ=255)
                      ccGGCAAATGCGCTGACTATCCAGACCGACGACCAGCAAAAtt                                                             >  1:1081316/1‑43 (MQ=255)
                       cGGCAAATGCGCTGACTATCCAGACCGACGACCAGCAAAATTCGCGCACCTGGCAGCTCTATCTGGGGCGt                                >  1:2359170/1‑71 (MQ=255)
                             aTGCTCTGACTATCCAGACCGACGACCAGCAAAATTCGCGCACCTGGCAGCTCTATCTGGGGCGTTGTAtt                          <  1:4685428/71‑1 (MQ=255)
                             aTGCGCTGACTATCCAGACCGACGACCAGCAAAATTCGCGCACCTGGCAGCTCTATc                                        <  1:2845254/57‑1 (MQ=255)
                             aTGCGCTGACTATCCAGACCGACGACCAGCAAAATTCGCGCACCTGGCAGCTCTATc                                        <  1:5102218/57‑1 (MQ=255)
                             aTGCGCTGACTATCCAGACCGACGACCAGCAAAATTCGCGCACCTGGCAGCTCTATCTGGGGCGTTGTAtt                          <  1:4408309/71‑1 (MQ=255)
                             aTGCGCTGACTATCCAGACCGACGACCAGCAAAACTCGCGCACCTTGCAGCTCTATCTGGGGCGTTGTAtt                          <  1:1902449/71‑1 (MQ=255)
                             aTGCGCTGACTATCCAGACCGACGACCAGCAAAACTCGCGCACCTGGCAGCTCTATCTGGGGCGTTGTAtt                          <  1:2976707/71‑1 (MQ=255)
                                           cAGCCCGACGACCAGCAAAACTCGCGCACCTGGCAGCTCTATCTGGGGCGTTGTATTTACTGCGGACGtt             >  1:4408708/1‑70 (MQ=255)
                                           cAGACCGACGACCAGCAAAACTCGCGCACCTGGCAGCTCTATCTGGGGCGTTGTATTTACTGCGGACGTtg            >  1:4273712/1‑71 (MQ=255)
                                             gACCGACGACCAGCAAAATTCGCGCACCTGGCAGCTCTATCTGGGGCGTTGTATTTACTGCGGACGTtgtg          <  1:2225819/71‑1 (MQ=255)
                                              cccGACGACCAGCAAAATTCGCGCACCTGGCAGCTCTATCTGGGGCGTTGTATTTACTGCGGACGTtgt           >  1:2221222/2‑69 (MQ=255)
                                               ccGACGGCCAGCAAAACTCGCGCACCTGGCAGCTCTa                                          >  1:661011/1‑37 (MQ=37)
                                               ccGACGACCAGCAAAACTCGCGCACCTGGCAGCTCTa                                          >  1:3811840/1‑37 (MQ=255)
                                                 gacgacCAGCAAAACTCGCGCACCTGGCAGCTCTATCTGGGGCGTTGTATTTACTGCGGACGTTGTgaaga      <  1:1701974/71‑1 (MQ=255)
                                                 gacgacCAGCAAAACTCGCGCACCTGGCAGCTCTATCTGGGGCGTTGTATTTACTGCGGACGTTGTgaaga      <  1:2005605/71‑1 (MQ=255)
                                                  acgacCAGCAAAATTCGCGCACCTGGCAGCTCTATCTGGGGc                                  >  1:2135303/1‑42 (MQ=255)
                                                      ccAGCAAAATTCGCGCACCTGGCAGCTCTATCTGGGGCGTTGTATTTACTg                     >  1:3566639/1‑51 (MQ=255)
                                                      ccAGCAAAACTCGCGCACCTGGCAGCTCTATCTGGGGCGTTGTATTTACTGCGGACGTTGTGAAGAAgtg  <  1:4303963/70‑1 (MQ=255)
                                                         gCAAAACTCGCGCACCTGGCAGCTCTATCTGGGGCGTTGTa                            <  1:983556/41‑1 (MQ=255)
                                                               |                                                            
CGGTGCCTGCGCCTGTGCTTGTCCGGCAAATGCGCTGACTATCCAGACCGACGACCAGCAAAATTCGCGCACCTGGCAGCTCTATCTGGGGCGTTGTATTTACTGCGGACGTTGTGAAGAAGTG  >  W3110S.gb/2608948‑2609071

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: