Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,141,782 G→T 16.6% P347T (CCG→ACG)  hybB ← predicted hydrogenase 2 cytochrome b type component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,141,7820GT16.6% 85.6 / 4.7 30P347T (CCG→ACG) hybBpredicted hydrogenase 2 cytochrome b type component
Reads supporting (aligned to +/- strand):  ref base G (16/9);  new base T (2/3);  total (18/12)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01

TGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTG  >  W3110S.gb/3141713‑3141844
                                                                     |                                                              
tGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCggg                                                               <  1:926930/71‑1 (MQ=255)
tGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCggg                                                               <  1:807802/71‑1 (MQ=255)
      cGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTa                                                          <  1:1721159/70‑1 (MQ=255)
         tGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAg                                                            <  1:639658/65‑1 (MQ=255)
         tGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAg                                                            <  1:1303888/65‑1 (MQ=255)
         tGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAg                                                            <  1:1452738/65‑1 (MQ=255)
             aGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTa                                                          >  1:1614497/1‑63 (MQ=255)
             aGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAAc                                                  <  1:1808988/71‑1 (MQ=255)
             aGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAAc                                                  <  1:158982/71‑1 (MQ=255)
                     tctcAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGccc                                          >  1:388239/1‑71 (MQ=255)
                     tctcAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGccc                                          >  1:897050/1‑71 (MQ=255)
                           tAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTg                                    <  1:558462/71‑1 (MQ=255)
                           tAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTg                                    <  1:539792/71‑1 (MQ=255)
                             gCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAAc                                                  >  1:1788800/1‑55 (MQ=255)
                             gCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGaa                                  >  1:1434737/1‑71 (MQ=255)
                                 caAAACCAATAGAAATCAACAGTTCTTCCCAGGTCggg                                                               <  1:934157/38‑1 (MQ=255)
                                   aaaCCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTa                                                          >  1:751769/1‑41 (MQ=255)
                                   aaaCCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTa                                                          >  1:25456/1‑41 (MQ=255)
                                   aaaCCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTa                                                          >  1:1509649/1‑41 (MQ=255)
                                   aaaCCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTa                                                          >  1:129771/1‑41 (MQ=255)
                                          tAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGaa                     >  1:43713/1‑71 (MQ=255)
                                          tAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGaa                     >  1:1048288/1‑71 (MQ=255)
                                          tAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGa                      >  1:369171/1‑70 (MQ=255)
                                          tAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGa                      >  1:873424/1‑70 (MQ=255)
                                                  aaCAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGt                >  1:237494/1‑68 (MQ=255)
                                                  aaCAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGt                >  1:235312/1‑68 (MQ=255)
                                                  aaCAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGt                >  1:811349/1‑68 (MQ=255)
                                                      gTTCTTCCCAGGTCGTGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGc         >  1:275309/1‑71 (MQ=255)
                                                      gTTCTTCCCAGGTCGTGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGc         >  1:1223536/1‑71 (MQ=255)
                                                          ttcCCAGGTCGGGATGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCa       <  1:257456/69‑1 (MQ=255)
                                                             ccAGGTCGTGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTg  <  1:1586247/71‑1 (MQ=255)
                                                             ccAGGTCGTGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTg  <  1:827287/71‑1 (MQ=255)
                                                             ccAGGTCGTGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTg  <  1:1251219/71‑1 (MQ=255)
                                                                     |                                                              
TGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTG  >  W3110S.gb/3141713‑3141844

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: