Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,141,409 C→T 12.6% A79A (GCG→GCA hybC ← hydrogenase 2, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,141,4090CT12.6% 266.7 / 24.1 95A79A (GCG→GCAhybChydrogenase 2, large subunit
Reads supporting (aligned to +/- strand):  ref base C (46/37);  new base T (6/6);  total (52/43)
Fisher's exact test for biased strand distribution p-value = 7.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

ATGCGTGGTGTGCGCAGCCAGAATGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGCAC  >  W3110S.gb/3141340‑3141477
                                                                     |                                                                    
aTGCGTGGTGTGCGCAGCCAGAATGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAgcg                                                                     <  1:3068077/71‑1 (MQ=255)
 tGCGTGGTGTGCGCAGCCAGAATGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAgtgc                                                                    <  1:1756499/71‑1 (MQ=255)
   cGTGGTGTGCGCAGCCAGAATGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCAc                                                                  >  1:2382887/1‑71 (MQ=255)
   cGTGGTGTGCGCAGCCAGAATGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCAc                                                                  >  1:69265/1‑71 (MQ=255)
    gtggtgTGCGCAGCCAGAATGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACt                                                                 <  1:3958026/71‑1 (MQ=255)
       gtgtGCGCAGCCAGAATGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTc                                                              <  1:4387133/71‑1 (MQ=255)
       gtgtGCGCAGCCAGAATGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTc                                                              <  1:4168451/71‑1 (MQ=255)
       gtgtGCGCAGCCAGAATGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTc                                                              <  1:3355649/71‑1 (MQ=255)
          tgCGCAGCCAGAATGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTg                                                            >  1:3671260/1‑70 (MQ=255)
          tgCGCAGCCAGAATGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGc                                                           >  1:2610042/1‑71 (MQ=255)
            cgcAGCCAGAATGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGTGCACTTTCTGCCg                                                         >  1:4600842/1‑71 (MQ=255)
            cgcAGCCAGAATGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGcc                                                          >  1:3225444/1‑70 (MQ=255)
            cgcAGCCAGAATGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCg                                                         >  1:1489630/1‑71 (MQ=255)
            cgcAGCCAGAATGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCg                                                         >  1:3855898/1‑71 (MQ=255)
                 ccAGAATGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTg                                                            <  1:2412401/63‑1 (MQ=255)
                  cAGAATGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGaa                                                   >  1:4087333/1‑71 (MQ=255)
                  cAGAATGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGa                                                    <  1:2233004/70‑1 (MQ=255)
                   aGAATGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGcc                                                          >  1:4325088/1‑63 (MQ=255)
                   aGAATGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAAc                                                  >  1:664196/1‑71 (MQ=255)
                     aaTGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACgg                                                <  1:3149561/71‑1 (MQ=255)
                     aaTGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACg                                                 >  1:196520/1‑70 (MQ=255)
                      atgatgTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGcc                                                          >  1:2953691/1‑60 (MQ=255)
                      atgatgTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGcc                                                          >  1:2921319/1‑60 (MQ=255)
                      atgatgTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACgg                                                >  1:422513/1‑70 (MQ=255)
                       tgatgTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCATCGCACTTTCTGCCGCACGAACgg                                                >  1:4642513/1‑69 (MQ=255)
                       tgatgTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACgg                                                >  1:2981939/1‑69 (MQ=255)
                       tgatgTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACgg                                                >  1:1000434/1‑69 (MQ=255)
                       tgatgTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGaa                                              <  1:3279507/71‑1 (MQ=255)
                       tgatgTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGaa                                              <  1:2375389/71‑1 (MQ=255)
                       tgatgTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGaa                                              <  1:176875/71‑1 (MQ=255)
                        gatgTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACgg                                                >  1:28657/1‑68 (MQ=255)
                         atgTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAgagc                                                                    >  1:1999933/1‑47 (MQ=255)
                         atgTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACttt                                                               >  1:3501417/1‑52 (MQ=255)
                         atgTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACttt                                                               >  1:907104/1‑52 (MQ=255)
                         atgTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGa                                            >  1:3806878/1‑71 (MQ=255)
                         atgTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGa                                            <  1:2446318/71‑1 (MQ=255)
                            ttACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGaa                                              <  1:3468759/66‑1 (MQ=255)
                              aCGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACt                                                                 >  1:3312786/1‑45 (MQ=255)
                               cGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAgcgc                                      >  1:4513776/1‑71 (MQ=255)
                                 gATGTATTGCGCGTTAACCGGAACGTCGATATTCAGTGCACTTTCTGCCGCACGAACGGAAGACAGCGCGt                                    <  1:2218711/71‑1 (MQ=255)
                                 gATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGt                                    <  1:2135509/71‑1 (MQ=255)
                                  aTGTATTGCGCGTTAACCGGAACGTCGATATTCAGTGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTg                                   <  1:3318218/71‑1 (MQ=255)
                                   tGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACg                                                     <  1:868928/52‑1 (MQ=255)
                                   tGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGt                                    <  1:2988867/69‑1 (MQ=255)
                                    gTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACt                                                                 >  1:3540421/1‑39 (MQ=255)
                                    gTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGa                                            <  1:802615/60‑1 (MQ=255)
                                     tATTGCGCGTTAACCGGAACGTCGATATTCAGTGCACTTTCTGCCGCACg                                                     >  1:3673716/1‑50 (MQ=255)
                                     tATTGCGCGTTAACCGGAACGTCGATATTCAGTGCACTTTCTGCCGCACg                                                     >  1:328522/1‑50 (MQ=255)
                                      aTTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTg                               >  1:1472606/1‑71 (MQ=255)
                                       ttGCGCGTTAACCGGAACGTCGATATTCAGTGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTgg                              <  1:3810600/71‑1 (MQ=255)
                                       ttGCGCGTTAACCGGAACGTCGATATTCAGTGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTgg                              <  1:4096009/71‑1 (MQ=255)
                                       ttGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCa                                                       >  1:3364863/1‑46 (MQ=255)
                                       ttGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGaa                                              <  1:2247739/55‑1 (MQ=255)
                                       ttGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGt                                    <  1:1051011/65‑1 (MQ=255)
                                       ttGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTgg                              >  1:158492/1‑71 (MQ=255)
                                       ttGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTgg                              <  1:3492300/71‑1 (MQ=255)
                                        tGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGcgcg                                     >  1:2646189/1‑63 (MQ=255)
                                        tGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGcgcg                                     >  1:699328/1‑63 (MQ=255)
                                         gcgcGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAg                                 >  1:828054/1‑66 (MQ=255)
                                         gcgcGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTgg                              >  1:4058378/1‑69 (MQ=255)
                                         gcgcGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTgg                              >  1:1353090/1‑69 (MQ=255)
                                         gcgcGTTAACCGGAACGTCGATATCCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAg                                 >  1:1749274/1‑66 (MQ=255)
                                          cgcgTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGgta                            <  1:3197250/70‑1 (MQ=255)
                                             gTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTAc                        <  1:523455/71‑1 (MQ=255)
                                              ttAACCGGAACGTCGATATTCAGTGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGgtagta                         <  1:4535210/69‑1 (MQ=255)
                                              ttAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACa                       <  1:2873408/71‑1 (MQ=255)
                                              ttAACCGGAACGTCGATATTCAGCGCACTTTCTGACGCACGAACGGAAGACAGCGCGTGAGTGgt                             <  1:3645211/65‑1 (MQ=255)
                                               tAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGaa                                                   >  1:148393/1‑42 (MQ=255)
                                                aaCCGGAACGTCGATATTCAGTGCACTTTCTGCCGCACGaa                                                   >  1:937665/1‑41 (MQ=255)
                                                aaCCGGAACGTCGATATTCAGTGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACAt                      <  1:4449104/70‑1 (MQ=255)
                                                aaCCGGAACGTCGATATTCAGTGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATa                     >  1:1031575/1‑71 (MQ=255)
                                                aaCCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACAt                      <  1:3690190/70‑1 (MQ=255)
                                                aaCCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATa                     >  1:3269936/1‑71 (MQ=255)
                                                aaCCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATa                     >  1:3264396/1‑71 (MQ=255)
                                                  ccGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGaa                                                   >  1:2527896/1‑39 (MQ=255)
                                                  ccGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGaa                                                   >  1:4120282/1‑39 (MQ=255)
                                                  ccGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAg                                 <  1:1680862/57‑1 (MQ=255)
                                                  ccGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGgt                             <  1:3193595/61‑1 (MQ=255)
                                                   cGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTg                                   <  1:3429517/54‑1 (MQ=255)
                                                   cGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGc                  <  1:1000796/71‑1 (MQ=255)
                                                    ggAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTg                                   <  1:1326966/53‑1 (MQ=255)
                                                    ggAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGc                  >  1:1073989/1‑70 (MQ=255)
                                                    ggAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGc                  >  1:3269958/1‑70 (MQ=255)
                                                       aCGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACa                                          >  1:999912/1‑43 (MQ=255)
                                                       aCGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGt                                    >  1:72254/1‑49 (MQ=255)
                                                         gTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGa            <  1:4035859/71‑1 (MQ=255)
                                                          tCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGgtagta                         <  1:1826006/57‑1 (MQ=255)
                                                            gATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGa                                  >  1:4070688/1‑46 (MQ=255)
                                                            gATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGgtag                           <  1:788438/53‑1 (MQ=255)
                                                             atatTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGgt                             <  1:868900/50‑1 (MQ=255)
                                                                ttCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAg                                 <  1:2120232/43‑1 (MQ=255)
                                                                ttCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGa            <  1:1925730/64‑1 (MQ=255)
                                                                 tCAGTGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTgg                              >  1:2760775/1‑45 (MQ=255)
                                                                  cAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGgta                            >  1:3644145/1‑46 (MQ=255)
                                                                  cAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACa                       >  1:3134581/1‑51 (MQ=255)
                                                                  cAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGca   >  1:2589706/1‑71 (MQ=255)
                                                                   aGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGcac  <  1:2571726/71‑1 (MQ=255)
                                                                   aGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGcac  <  1:1770267/71‑1 (MQ=255)
                                                                     cgcACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGgtagt                          <  1:3021071/45‑1 (MQ=255)
                                                                     cgcACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCcaca              <  1:2578560/57‑1 (MQ=255)
                                                                     |                                                                    
ATGCGTGGTGTGCGCAGCCAGAATGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGCAC  >  W3110S.gb/3141340‑3141477

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