Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,094,702 T→G 10.0% intergenic (‑430/+21) rhaD ← / ← rhaA rhamnulose‑1‑phosphate aldolase/L‑rhamnose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,094,7020TG10.0% 162.6 / 9.3 50intergenic (‑430/+21)rhaD/rhaArhamnulose‑1‑phosphate aldolase/L‑rhamnose isomerase
Reads supporting (aligned to +/- strand):  ref base T (27/18);  new base G (3/2);  total (30/20)
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

TCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTGCCACGGCAACGATTTCTGC  >  NC_000913/4094562‑4094849
                                                                                                                                            |                                                                                                                                                   
tCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATGCGGCGCGAGcgcctt                                                                                                                                             >  4:543644/1‑147 (MQ=255)
     tCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATGCGGCGCGAGCGCCTTATCCg                                                                                                                                        <  3:543644/149‑1 (MQ=11)
            tACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATTCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAgtgt                                                                                                                                 >  3:416112/1‑149 (MQ=255)
                       ccGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCgcggc                                                                                                                      >  1:558143/1‑149 (MQ=255)
                        cGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCgcggcg                                                                                                                     >  1:71973/1‑149 (MQ=255)
                            cTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTc                                                                                                                 >  3:294369/1‑149 (MQ=255)
                                    cGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAAtttc                                                                                                         >  4:642424/1‑149 (MQ=255)
                                    cGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAAtttc                                                                                                         <  4:294369/149‑1 (MQ=255)
                                      gATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAAtttctt                                                                                                       >  1:56517/1‑149 (MQ=255)
                                           ggCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTCCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCa                                                                                                  >  3:84102/1‑149 (MQ=255)
                                           ggCACTTGTGCCGCATCCGGCAATCAATGCCTGATGAGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCa                                                                                                  >  4:653222/1‑149 (MQ=255)
                                              aCTTGTGCTGCATCTGGCAATAAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCTGTCCCGTAGGCCGGACAAGGCACTCGCGCCGCATGCGGGAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATaa                                                                                               <  3:427296/149‑1 (MQ=18)
                                              aCTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATaa                                                                                               <  2:71973/149‑1 (MQ=255)
                                                  gtgCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGccc                                                                                           >  3:551576/1‑149 (MQ=255)
                                                     ccGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCa                                                                                        >  2:248114/1‑149 (MQ=255)
                                                           ccGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGctct                                                                                  >  1:262382/1‑149 (MQ=255)
                                                            cGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCtctc                                                                                 <  4:551576/149‑1 (MQ=255)
                                                                 atcaatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGcc                                                                            >  4:38023/1‑149 (MQ=255)
                                                                     atGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTc                                                                        <  1:248114/149‑1 (MQ=255)
                                                                       gCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGc                                                                      <  1:652798/149‑1 (MQ=255)
                                                                          tGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTAc                                                                   <  3:642424/149‑1 (MQ=255)
                                                                           gATGCGACGCTGTCGCGCCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGCCTCAAAATTTCTTTCTCAGAAGCCCGTACGCACTCCAGCCATTCGCTAcc                                                                  <  4:416112/149‑1 (MQ=25)
                                                                                gACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACctgctg                                                             >  3:125500/1‑149 (MQ=255)
                                                                                gACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACctgctg                                                             >  3:323617/1‑149 (MQ=255)
                                                                                  cGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGc                                                           <  2:56517/149‑1 (MQ=255)
                                                                                       tCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATc                                                      <  2:262382/149‑1 (MQ=255)
                                                                                       tCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATc                                                      <  4:323617/149‑1 (MQ=255)
                                                                                         gcgTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATGCGGCGCGAGCGCCTTATCCg                                                                                                                                        <  3:347787/65‑1 (MQ=14)
                                                                                         gcgTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATGCGGCGCGAGCGCCTTATCCg                                                                                                                                        >  4:347787/1‑65 (MQ=14)
                                                                                          cgTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTg                                                   >  4:391454/1‑149 (MQ=255)
                                                                                           gTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCa                                                                                        >  3:14640/1‑111 (MQ=255)
                                                                                           gTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCa                                                                                        <  4:14640/111‑1 (MQ=255)
                                                                                             cTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACg                                                >  2:408733/1‑149 (MQ=255)
                                                                                               tATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGtt                                              <  4:125500/149‑1 (MQ=255)
                                                                                                aTCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGc                                                                                     >  4:416249/1‑109 (MQ=255)
                                                                                                aTCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGc                                                                                     <  3:416249/109‑1 (MQ=255)
                                                                                                aTCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGa                                                  >  1:542468/1‑144 (MQ=255)
                                                                                                aTCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGa                                                  <  2:542468/144‑1 (MQ=255)
                                                                                                 tCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTgg                                            <  1:467193/149‑1 (MQ=255)
                                                                                                 tCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTgg                                            <  3:653222/149‑1 (MQ=255)
                                                                                                       cTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATAc                                      >  2:648335/1‑149 (MQ=255)
                                                                                                        tACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACa                                     >  1:8294/1‑149 (MQ=255)
                                                                                                           aaCTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACAttt                                  >  3:102806/1‑149 (MQ=255)
                                                                                                                     gACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGc                        >  4:563937/1‑149 (MQ=255)
                                                                                                                       cccGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCACACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCt                      <  2:558143/149‑1 (MQ=255)
                                                                                                                        ccGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTg                     >  2:639590/1‑149 (MQ=255)
                                                                                                                                cGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTGCCACGGCa             >  4:620612/1‑149 (MQ=255)
                                                                                                                                 ggATAAGCCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTGCCACGCCaa            >  4:284961/1‑149 (MQ=255)
                                                                                                                                  gatgcggCGCGAGCGCCTTATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTGCCACGGCAAc           >  4:259854/6‑149 (MQ=255)
                                                                                                                                           cTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTGCCACGGCAACGATTTCTGc  >  1:494047/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
TCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTGCCACGGCAACGATTTCTGC  >  NC_000913/4094562‑4094849

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: