Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,063,843 A→C 13.2% intergenic (‑82/+32) yihR ← / ← yihS mutarotase superfamily protein, YihR family/aldose‑ketose isomerase; D‑mannose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,063,8430AC13.2% 961.7 / 23.1 350intergenic (‑82/+32)yihR/yihSmutarotase superfamily protein, YihR family/aldose‑ketose isomerase; D‑mannose isomerase
Reads supporting (aligned to +/- strand):  ref base A (252/31);  new base C (2/43);  total (254/96)
Fisher's exact test for biased strand distribution p-value = 3.72e-31
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 415 total aligned reads displayed.

GGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCA  >  CP009273/4063776‑4063908
                                                                   |                                                                 
 gCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGagag                                                                  >  1:3737764/1‑68 (MQ=255)
   tGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAggg                                                                >  1:33802981/1‑68 (MQ=255)
       gTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGaaa                                                            >  1:29646028/1‑68 (MQ=255)
       gTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGaaa                                                            >  1:27014571/1‑68 (MQ=255)
        tCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGaaa                                                            >  1:17220883/1‑67 (MQ=255)
             cccAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGc                                                      >  1:2050552/1‑68 (MQ=255)
              ccAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGc                                                      >  1:28468162/1‑67 (MQ=255)
                aaaacCTCCCCCGGCCGGCCCCCTCGCCCCTTTGGGGAGCGGGTTCGGCCGCGGGGAAAACCCGCAcc                                                   <  1:14676561/68‑1 (MQ=12)
                aaaaaCTGCACCGGCCGGTCCCCTCGCCGCTTTGGGTAGAGGGTTTGGGAGTGGGGAAAACCCGCAcc                                                   <  1:17697720/68‑1 (MQ=21)
                 aaaaCTCCCCCGCACGGCCCCCCCCCCCCTTTGGGGAGCGGGTTCGGGAGCGGGGAAAACCCGCAccc                                                  <  1:19198507/68‑1 (MQ=16)
                 aaaaCTCCACCGGACGGTCCCCTCGCCCCTTTGGGGAGTGGGTTTGGGAGAGGGGAAAACCCGCAccc                                                  <  1:6800746/68‑1 (MQ=21)
                    aCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAt                                               >  1:10130744/1‑68 (MQ=255)
                    aCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAt                                               >  1:21713237/1‑68 (MQ=255)
                    aCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAt                                               >  1:31587862/1‑68 (MQ=255)
                    aCTCCACCGGCGGGCCCCCTCCCCCCTTTGGGGTGTGGTTTTGGGCGTGGGGAAAACCCGCACCCTAt                                               <  1:11349691/68‑1 (MQ=12)
                     cTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtt                                              >  1:28963924/1‑68 (MQ=255)
                         aCCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCa                                          >  1:32825846/1‑68 (MQ=255)
                         aCCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCa                                          >  1:7089735/1‑68 (MQ=255)
                          ccGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTATCAg                                         >  1:29997083/1‑68 (MQ=255)
                              aCGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCCCCCTATTTTCAGATCa                                     >  1:23310508/1‑68 (MQ=255)
                              aCGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCa                                     >  1:36376002/1‑68 (MQ=255)
                              aCGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCa                                     >  1:29232762/1‑68 (MQ=255)
                              aCGGTCCCCTCGCCCCGTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCCCCCTTTTTTAAGATCa                                     >  1:23511327/1‑68 (MQ=18)
                               cGGGCCCCGCGCCCCTTGGGTTTGGGGTTTTGGTTGCGGGCAAAACCCGCACCCTATTTTCAGATCAg                                    <  1:24804726/68‑1 (MQ=255)
                                  cccccTCCCCCGTTTGGGGATTGGGTTAGGTTGTGGGGAAAACCCGCACCCTATTTTCAGATCAGtt                                  <  1:5364922/66‑1 (MQ=21)
                                   ccccTCGCCCCTTTGGGGAGAGGGTTAGGGAGCGGGGAACACCCGCACCCTATTTTCAGATCAGTTAt                                >  1:2203233/1‑68 (MQ=255)
                                   ccccTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTAt                                >  1:18726468/1‑68 (MQ=255)
                                      cTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTc                             >  1:34541001/1‑68 (MQ=255)
                                        cGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGc                           >  1:27748810/1‑68 (MQ=255)
                                        cGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGc                           >  1:31384723/1‑68 (MQ=255)
                                        cGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGc                           >  1:30448816/1‑68 (MQ=255)
                                        cGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGc                           >  1:34869452/1‑68 (MQ=255)
                                        cGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGc                           >  1:35696404/1‑68 (MQ=255)
                                        cGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGc                           >  1:3888930/1‑68 (MQ=255)
                                         cccccTTGGGGGAGTGGGTTCGGGAGCGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCa                          <  1:24707023/67‑1 (MQ=255)
                                          ccccTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCAt                         >  1:7779598/1‑68 (MQ=255)
                                           cccTTTGGGGAGAGGGTTAGGGACAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCAtt                        <  1:6453390/68‑1 (MQ=255)
                                               ttGGGGTGTGGGTTTGGGCGCGGGCAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAat                     <  1:25574459/67‑1 (MQ=255)
                                               ttGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAata                    >  1:36451628/1‑68 (MQ=255)
                                                tGGGGTGTGGTTTGGGGCGCGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAatat                   <  1:20808893/68‑1 (MQ=255)
                                                tGGGGAGAGGGTTAGGGAGCGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAatat                   <  1:23781719/68‑1 (MQ=255)
                                                 ggggAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAatat                   >  1:9622732/1‑67 (MQ=255)
                                                 ggggAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAatat                   >  1:25228845/1‑67 (MQ=255)
                                                 ggggAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAatat                   >  1:13084600/1‑67 (MQ=255)
                                                   ggAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCca                >  1:24171304/1‑68 (MQ=255)
                                                   ggAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCca                >  1:12006855/1‑68 (MQ=255)
                                                   ggAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCca                >  1:12426130/1‑68 (MQ=255)
                                                   ggAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCca                >  1:5892427/1‑68 (MQ=255)
                                                   ggAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCca                <  1:12883438/68‑1 (MQ=255)
                                                   ggAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCca                >  1:36725448/1‑68 (MQ=255)
                                                   ggAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCca                >  1:20993939/1‑68 (MQ=255)
                                                   ggAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCca                >  1:14920480/1‑68 (MQ=255)
                                                   cgagaGGGTTACGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCca                >  1:20083865/2‑68 (MQ=255)
                                                    gagtGGGTTTGGGAGCGGGCAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCca                <  1:35992392/67‑1 (MQ=255)
                                                    gagaGGGTTAGGGAGCGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCca                <  1:1550461/67‑1 (MQ=255)
                                                    gagaGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCca                >  1:30218977/1‑67 (MQ=255)
                                                     agagGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCcagc              >  1:7264133/1‑68 (MQ=255)
                                                     agagGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCcagc              >  1:32673978/1‑68 (MQ=255)
                                                     agagGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCcagc              >  1:16060978/1‑68 (MQ=255)
                                                     agagGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCcagc              >  1:4538876/1‑68 (MQ=255)
                                                        gggTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCcagcag            >  1:14068064/1‑67 (MQ=255)
                                                             aGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGAccc        >  1:34273713/1‑66 (MQ=255)
                                                              gggAGAGGGGAAAACCCGCACCCTATTTTCAGATCa                                     <  1:10711052/36‑1 (MQ=255)
                                                              gggAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGAccc        >  1:4735043/1‑65 (MQ=255)
                                                               ggggagGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCt     >  1:18044425/4‑67 (MQ=255)
                                                               ggAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCCGCAGACCctct     >  1:8752868/1‑67 (MQ=255)
                                                               ggAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCt     >  1:24962414/1‑67 (MQ=255)
                                                               ggAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCt     >  1:26055493/1‑67 (MQ=255)
                                                               ggAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCt     >  1:26670469/1‑67 (MQ=255)
                                                               ggAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCt     >  1:23121105/1‑67 (MQ=255)
                                                               ggAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCt     >  1:2272262/1‑67 (MQ=255)
                                                               ggAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCt     >  1:2744146/1‑67 (MQ=255)
                                                               ggAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCt     >  1:13242000/1‑67 (MQ=255)
                                                               ggAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCt     >  1:32177642/1‑67 (MQ=255)
                                                               ggAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCt     >  1:22394388/1‑67 (MQ=255)
                                                               ggAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCt     >  1:21456126/1‑67 (MQ=255)
                                                               ggAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCt     >  1:14424953/1‑67 (MQ=255)
                                                               ggAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCt     >  1:11088272/1‑67 (MQ=255)
                                                               ggAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCt     >  1:19497762/1‑67 (MQ=255)
                                                               ggAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCt     >  1:7617057/1‑67 (MQ=255)
                                                               ggAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCt     >  1:18444416/1‑67 (MQ=255)
                                                               ggAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCt     >  1:10901793/1‑67 (MQ=255)
                                                               ggAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCt     >  1:8494707/1‑67 (MQ=255)
                                                               ggAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCt     >  1:34051240/1‑67 (MQ=255)
                                                               ggAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCt     >  1:17480300/1‑67 (MQ=255)
                                                               ggAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCt     >  1:24519563/1‑67 (MQ=255)
                                                                gagaGGGGAAAACCCGCCCCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGCCCCGCTGc   >  1:33682356/1‑68 (MQ=255)
                                                                gagaGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGc   >  1:36831096/1‑68 (MQ=255)
                                                                gagaGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGc   >  1:9267820/1‑68 (MQ=255)
                                                                gagaGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGc   >  1:26359712/1‑68 (MQ=255)
                                                                gagaGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGc   >  1:5605321/1‑68 (MQ=255)
                                                                gagaGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGc   >  1:35166534/1‑68 (MQ=255)
                                                                gagaGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGc   >  1:15867951/1‑68 (MQ=255)
                                                                gagaGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGc   >  1:31021357/1‑68 (MQ=255)
                                                                gagaGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGc   >  1:32405216/1‑68 (MQ=255)
                                                                 agagGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCa  >  1:5217950/1‑68 (MQ=255)
                                                                 agagGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCa  >  1:36197234/1‑68 (MQ=255)
                                                                 agagGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCa  >  1:1965753/1‑68 (MQ=255)
                                                                 agagGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCa  >  1:33211065/1‑68 (MQ=255)
                                                                 agagGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCa  >  1:16818670/1‑68 (MQ=255)
                                                                   |                                                                 
GGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCA  >  CP009273/4063776‑4063908

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: