Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,008,716 T→C 5.8% K102K (AAA→AAG hybO ← hydrogenase 2, small subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,008,7160TC 5.8% 194.2 / 2.4 69K102K (AAA→AAGhybOhydrogenase 2, small subunit
Reads supporting (aligned to +/- strand):  ref base T (23/42);  new base C (2/2);  total (25/44)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

ATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACTTCGTGATACTCCAGAGAGATAGTCTCC  >  minE/2008647‑2008785
                                                                     |                                                                     
aTCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTAtgt                                                                       >  1:231648/1‑70 (MQ=255)
aTCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTAtgt                                                                       <  1:1350584/70‑1 (MQ=255)
 tCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTAtgtt                                                                      >  1:601024/1‑70 (MQ=255)
 tCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTAtgtt                                                                      <  1:406831/70‑1 (MQ=255)
    tttAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTAtgtttgt                                                                   >  1:1835034/1‑70 (MQ=255)
     ttAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTAtgttt                                                                     >  1:803067/1‑67 (MQ=255)
      tAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTtc                                                                 <  1:1473791/70‑1 (MQ=255)
        aTGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTtct                                                                <  1:1162753/69‑1 (MQ=255)
           ggATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTc                                                            <  1:35525/70‑1 (MQ=255)
           ggATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCg                                                           <  1:2198104/71‑1 (MQ=255)
            gATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGa                                                          <  1:903197/71‑1 (MQ=255)
             aTGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTAtgtt                                                                      <  1:2115316/58‑1 (MQ=255)
              tGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGAcc                                                        <  1:3083264/71‑1 (MQ=255)
                gAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCtg                                                      <  1:1187964/71‑1 (MQ=255)
                 aaCCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCtg                                                      <  1:2925256/70‑1 (MQ=255)
                 aaCCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGCTTGTTCTCTTCGACCtga                                                     >  1:3815521/1‑71 (MQ=255)
                  aCCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCtgat                                                    <  1:1659327/71‑1 (MQ=255)
                  aCCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCtgat                                                    <  1:944460/71‑1 (MQ=255)
                   ccATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTAtgt                                                                       >  1:2591741/1‑51 (MQ=255)
                   ccATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTAtgt                                                                       >  1:2234938/1‑51 (MQ=255)
                   ccATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCtgatg                                                   <  1:112991/71‑1 (MQ=255)
                    cATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCtga                                                     <  1:1509735/68‑1 (MQ=255)
                      tCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTAtgtt                                                                      >  1:2047519/1‑49 (MQ=255)
                        ccacacTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCg                                               >  1:1136389/4‑70 (MQ=255)
                        caccacTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGa                                              <  1:1421039/71‑1 (MQ=255)
                         accacTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTAtgtttg                                                                    <  1:2191009/48‑1 (MQ=255)
                            acTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTAtgtttg                                                                    >  1:79502/1‑45 (MQ=255)
                            acTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAggc                                          >  1:2438029/1‑71 (MQ=255)
                             cTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTGTGTTTGTTCTCTTCGACCtg                                                      <  1:1388431/58‑1 (MQ=255)
                             cTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTtctc                                                               <  1:210122/49‑1 (MQ=255)
                                acacATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTtct                                                                <  1:971273/45‑1 (MQ=255)
                                acacATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCtt                                                             <  1:3784477/48‑1 (MQ=255)
                                acacATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCg                                                           <  1:3771554/50‑1 (MQ=255)
                                acacATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCt                                                       <  1:2096855/54‑1 (MQ=255)
                                acacATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGgcggcg                                      >  1:2700897/1‑71 (MQ=255)
                                 cacaTACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAggcgg                                        <  1:3169603/68‑1 (MQ=255)
                                 cacaTACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGcggcgg                                     <  1:2256705/71‑1 (MQ=255)
                                  acaTACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGa                                    <  1:1228325/71‑1 (MQ=255)
                                   cacaCTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGaa                                   <  1:1366125/68‑1 (MQ=255)
                                   caTACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAggcggc                                       >  1:166145/1‑67 (MQ=255)
                                   caTACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGa                                    <  1:1287523/70‑1 (MQ=255)
                                   caTACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGa                                    >  1:3629257/1‑70 (MQ=255)
                                    aTACTGCCCTTTGTACTTCTCGAGAGCGTTAtgtttgt                                                                   <  1:1903489/38‑1 (MQ=255)
                                    aTACTGCCCTTTGTACTTCTCGAGAGCGTTAtgtttgt                                                                   <  1:790964/38‑1 (MQ=255)
                                    aTACTGCCCTTTGTACTTCTCGAGAGCGTTAtgtttgt                                                                   <  1:1898941/38‑1 (MQ=255)
                                    aTACTGCCCTTTGTACTTCTCGAGAGCGTTATGCTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGa                                    >  1:2035965/1‑69 (MQ=255)
                                      aCTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTtc                                                                 >  1:242463/1‑38 (MQ=255)
                                       cTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCg                                                           <  1:1282077/43‑1 (MQ=255)
                                        tGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCAc                              >  1:303198/1‑71 (MQ=255)
                                         gCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACt                             >  1:303455/1‑71 (MQ=255)
                                          cccTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACtt                            <  1:3662327/71‑1 (MQ=255)
                                          cccTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACtt                            <  1:3190074/71‑1 (MQ=255)
                                          cccTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACt                             <  1:2462791/70‑1 (MQ=255)
                                            cTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGc                                >  1:1492222/1‑65 (MQ=255)
                                            cTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACTTCg                          >  1:1448836/1‑71 (MQ=255)
                                             tttGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGAcc                                                        <  1:2096882/40‑1 (MQ=255)
                                             tttGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCtgat                                                    <  1:2672840/44‑1 (MQ=255)
                                             tttGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGaa                                             <  1:1736319/51‑1 (MQ=255)
                                             tttGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACTTCGt                         >  1:2762159/1‑71 (MQ=255)
                                             tttGTACTTCTCGAGAGCGTTATGCTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCa                               <  1:3936076/65‑1 (MQ=255)
                                               tGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGAc                                                         >  1:1985273/1‑37 (MQ=255)
                                               tGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAgg                                           >  1:3933897/1‑51 (MQ=255)
                                               tGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCAc                              >  1:2006111/1‑64 (MQ=255)
                                               tGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACTTCGt                         <  1:2450982/69‑1 (MQ=255)
                                               tGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACTTCGTGa                       >  1:3481965/1‑71 (MQ=255)
                                                gTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACTTCGTGa                       <  1:3522687/70‑1 (MQ=255)
                                                 tACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGa                                              >  1:2466745/1‑46 (MQ=255)
                                                 tACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGa                                    <  1:3288367/56‑1 (MQ=255)
                                                   cTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACTTCGTGATAc                    <  1:501348/70‑1 (MQ=255)
                                                    ttCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGaa                                   <  1:933286/54‑1 (MQ=255)
                                                    ttCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACTTCGTGATACTc                  >  1:1872417/1‑71 (MQ=255)
                                                     tctcGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACTTCGTGATACTc                  >  1:2071974/1‑70 (MQ=255)
                                                      ctcGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAg                                 >  1:2299579/1‑54 (MQ=255)
                                                      ctcGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACTTCGTGATACTc                  >  1:2244729/1‑69 (MQ=255)
                                                      ctcGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACTTCGTGATACTc                  <  1:3356018/69‑1 (MQ=255)
                                                      ctcGAGAGCGTTATGCTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACTTCGTGATACTc                  <  1:3250036/69‑1 (MQ=255)
                                                         gagaGCGTTATGTTTGTTCTCTTCGACCTGATGACCGa                                              <  1:1889662/38‑1 (MQ=255)
                                                         gagaGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACTTCGTGATACTCCa                >  1:1791988/1‑68 (MQ=255)
                                                            agCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACTTCGTGATACTCCAgagaga          <  1:2740351/71‑1 (MQ=255)
                                                              cGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACTTCGTGATACTCCAGAGAGATa        <  1:2246740/71‑1 (MQ=255)
                                                               gTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACTTCGTGATACTCCagaga            <  1:525385/66‑1 (MQ=255)
                                                               gTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACTTCGTGATACTCCagaga            <  1:608650/66‑1 (MQ=255)
                                                                  aTGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACTTCGTGATACTCCAGAGAGATAGtct    >  1:2170018/1‑71 (MQ=255)
                                                                  aTGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACTTCGTGATACTCCAGAGAGATAGtct    >  1:1385984/1‑71 (MQ=255)
                                                                   tgtttgttCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACTTCGTGATACTCCAGAGAGATAGtctc   >  1:3009472/1‑71 (MQ=255)
                                                                     tttgttCTCTTCGACCTGATGACCGAAGGCGGCGGa                                    <  1:1078218/36‑1 (MQ=255)
                                                                     tttgttCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACTTc                           <  1:3400284/45‑1 (MQ=255)
                                                                     tttgttCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACTTc                           <  1:2509737/45‑1 (MQ=255)
                                                                     tttgttCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACTTCGTGATACTCCAGAGAGATAGtctc   <  1:1658724/69‑1 (MQ=255)
                                                                     tttgttCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACTTCGTGATACTCCAGAGAGATAGTCTcc  <  1:2060024/70‑1 (MQ=255)
                                                                     |                                                                     
ATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACTTCGTGATACTCCAGAGAGATAGTCTCC  >  minE/2008647‑2008785

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: