Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,530,282 G→A 25.2% I105I (ATC→ATT pdxK ← pyridoxal‑pyridoxamine kinase/hydroxymethylpyrimidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,530,2820GA25.2% 1255.9 / inf 694I105I (ATC→ATTpdxKpyridoxal‑pyridoxamine kinase/hydroxymethylpyrimidine kinase
Reads supporting (aligned to +/- strand):  ref base G (321/198);  new base A (109/66);  total (430/264)
Fisher's exact test for biased strand distribution p-value = 9.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.09e-01

Only 100 of 702 total aligned reads displayed.

TCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCAC  >  CP009273/2530215‑2530349
                                                                   |                                                                   
tCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAgg                                                                     <  1:16290882/68‑1 (MQ=255)
   ccGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGAtt                                                                  >  1:26646984/1‑68 (MQ=255)
       aTCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGa                                                              >  1:17962493/1‑68 (MQ=255)
       aTCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGa                                                              >  1:27648874/1‑68 (MQ=255)
       aTCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGa                                                              >  1:24496789/1‑68 (MQ=255)
       aTCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGa                                                              >  1:19737113/1‑68 (MQ=255)
       aTCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGa                                                              >  1:3740793/1‑68 (MQ=255)
       aTCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGa                                                              >  1:15749305/1‑68 (MQ=255)
       aTCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGa                                                              >  1:14131204/1‑68 (MQ=255)
       aTCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGa                                                              >  1:10715219/1‑68 (MQ=255)
       aTCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGAATTTTGa                                                              >  1:10128902/1‑68 (MQ=255)
       aTCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGAATTTTGa                                                              >  1:8141199/1‑68 (MQ=255)
            ccATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGc                                                         >  1:11842651/1‑68 (MQ=255)
             cATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCg                                                        >  1:33751749/1‑68 (MQ=255)
             cATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCg                                                        >  1:5281306/1‑68 (MQ=255)
             cATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCg                                                        >  1:11445522/1‑68 (MQ=255)
               tGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGAATTTTGATTTGCGAt                                                      <  1:16924151/68‑1 (MQ=255)
                  tCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGcc                                                   >  1:23375178/1‑68 (MQ=255)
                  tCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGcc                                                   >  1:31211875/1‑68 (MQ=255)
                  tCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGcc                                                   >  1:30661723/1‑68 (MQ=255)
                  tCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGcc                                                   >  1:29936978/1‑68 (MQ=255)
                  tCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGcc                                                   >  1:18726923/1‑68 (MQ=255)
                  tCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGcc                                                   >  1:28206560/1‑68 (MQ=255)
                  tCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGcc                                                   >  1:4002895/1‑68 (MQ=255)
                  tCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGcc                                                   >  1:17521326/1‑68 (MQ=255)
                  tCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGcc                                                   >  1:6583033/1‑68 (MQ=255)
                  tCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGcc                                                   >  1:10464589/1‑68 (MQ=255)
                  tCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGAATTTTGATTTGCGATGcc                                                   >  1:33373759/1‑68 (MQ=255)
                    aaTAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGAATTTTGATTTGCGATGCCGt                                                 <  1:25628367/68‑1 (MQ=255)
                     aTAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGtt                                                <  1:29249578/68‑1 (MQ=255)
                     aTAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGtt                                                <  1:35248677/68‑1 (MQ=255)
                     aTAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGtt                                                <  1:13742349/68‑1 (MQ=255)
                     aTAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGt                                                 >  1:17104117/1‑67 (MQ=255)
                     aTAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGt                                                 >  1:201608/1‑67 (MQ=255)
                     aTAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGAATTTTGATTTGCGATGCCGtt                                                <  1:2480887/68‑1 (MQ=255)
                     aTAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGAATTTTGATTTGCGATGCCGtt                                                <  1:12486490/68‑1 (MQ=255)
                     aTAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGAATTTTGATTTGCGATGCCGt                                                 >  1:27974467/1‑67 (MQ=255)
                     aTAGGGCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGtt                                                <  1:22992182/68‑1 (MQ=255)
                      tAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGAt                                                                   <  1:7445749/48‑1 (MQ=255)
                      tAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGtt                                                >  1:898840/1‑67 (MQ=255)
                       atgTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTcc                                              <  1:35646905/66‑1 (MQ=255)
                       aGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTcc                                              <  1:10960775/68‑1 (MQ=255)
                       aGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTcc                                              <  1:7222744/68‑1 (MQ=255)
                       aGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTcc                                              <  1:935299/68‑1 (MQ=255)
                       aGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTc                                               >  1:4555297/1‑67 (MQ=255)
                       aGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTc                                               >  1:26305030/1‑67 (MQ=255)
                       aGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGAATTTTGATTTGCGATGCCGTTcc                                              <  1:36587097/68‑1 (MQ=255)
                       aGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGAATTTTGATTTGCGATGCCGTTcc                                              <  1:9626434/68‑1 (MQ=255)
                       aGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGAATTTTGATTTGCGATGCCGTTc                                               >  1:32712008/1‑67 (MQ=255)
                       aGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGAATTTTGATTTGCGATGCCGTTc                                               >  1:34299291/1‑67 (MQ=255)
                        ggTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTccc                                             >  1:26980133/1‑68 (MQ=255)
                        ggTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTccc                                             >  1:25340246/1‑68 (MQ=255)
                        ggTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTccc                                             >  1:18263779/1‑68 (MQ=255)
                        ggTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGAATTTTGATTTGCGATGCCGTTccc                                             >  1:12704792/1‑68 (MQ=255)
                        ggTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGAATTTTGATTTGCGATGCCGTTccc                                             >  1:7856859/1‑68 (MQ=255)
                          tCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCat                                           >  1:18484385/1‑68 (MQ=255)
                          tCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCa                                            <  1:14526947/67‑1 (MQ=255)
                          tCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGAATTTTGATTTGCGATGCCGTTCCCat                                           >  1:31599917/1‑68 (MQ=255)
                          tCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGAATTTTGATTTGCGATGCCGTTCCCa                                            <  1:33061662/67‑1 (MQ=255)
                          tCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGAATTTTGATTTGCGATGCCGTTCCCa                                            <  1:21850810/67‑1 (MQ=255)
                           cAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCat                                           <  1:9713103/67‑1 (MQ=255)
                               aTGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAg                                       >  1:36255092/1‑67 (MQ=255)
                               aTGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAg                                       >  1:20675432/1‑67 (MQ=255)
                               aTGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAg                                       >  1:30360016/1‑67 (MQ=255)
                               aTGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGAATTTTGATTTGCGATGCCGTTCCCATATAg                                       >  1:24060992/1‑67 (MQ=255)
                               aTGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGAATTTTGATTTGCGATGCCGTTCCCATATAg                                       >  1:26029838/1‑67 (MQ=255)
                               aTGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGAATTTTGATTTGCGATGCCGTTCCCATATAg                                       <  1:8549500/67‑1 (MQ=255)
                                tGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCTATGCCGTTCCCATATAGc                                      >  1:17772947/1‑67 (MQ=255)
                                tGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAGcc                                     <  1:4330718/68‑1 (MQ=255)
                                tGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAGcc                                     <  1:23813367/68‑1 (MQ=255)
                                tGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAGc                                      <  1:20087442/67‑1 (MQ=255)
                                   tCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGAATTTTGATTTGCGATGCCGTTCCCATATAGCCCg                                   <  1:31928835/67‑1 (MQ=255)
                                    cTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAGCCtgtg                                 <  1:35969828/68‑1 (MQ=255)
                                    cTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTg                                 >  1:6045350/1‑68 (MQ=255)
                                    cTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTg                                 <  1:6317711/68‑1 (MQ=255)
                                    cTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTg                                 <  1:22188126/68‑1 (MQ=255)
                                         cGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTAc                            >  1:36823094/1‑68 (MQ=255)
                                         cGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTAc                            >  1:20858847/1‑68 (MQ=255)
                                         cGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTAc                            >  1:21367377/1‑68 (MQ=255)
                                         cGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTAc                            >  1:4968603/1‑68 (MQ=255)
                                         cGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTAc                            >  1:5496776/1‑68 (MQ=255)
                                         cGTAGCGCAGTCAGCCACTCGGCAAGAATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTAc                            >  1:6984151/1‑68 (MQ=255)
                                          gTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTACa                           >  1:1530725/1‑68 (MQ=255)
                                          gTAGCGCAGTCAGCCACTCGGCAAGAATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTACa                           <  1:25055947/68‑1 (MQ=255)
                                               gcAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCa                        <  1:16649294/66‑1 (MQ=255)
                                               gcAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCACg                      >  1:32278551/1‑68 (MQ=255)
                                               gcAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCACg                      >  1:2724195/1‑68 (MQ=255)
                                                 aGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCACGaa                    >  1:28699293/1‑68 (MQ=255)
                                                   tCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGt                  >  1:17249603/1‑68 (MQ=255)
                                                    cAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAg                   <  1:34074973/66‑1 (MQ=255)
                                                    cAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAg                   <  1:8288861/66‑1 (MQ=255)
                                                    cAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGtt                 >  1:13167942/1‑68 (MQ=255)
                                                        cACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGgc              <  1:32599304/67‑1 (MQ=255)
                                                         aCTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGgcgc            >  1:15011461/1‑68 (MQ=255)
                                                         aCTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGgcgc            >  1:34985315/1‑68 (MQ=255)
                                                            cgccAAGAATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGgcgcagc         <  1:8900024/65‑1 (MQ=255)
                                                            cGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGgcgcagc         <  1:19213413/68‑1 (MQ=255)
                                                            cGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGgcgcagc         <  1:12176668/68‑1 (MQ=255)
                                                            cGGCAAGAATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGgcgcagc         <  1:10001077/68‑1 (MQ=255)
                                                               cAAGGATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGgcgcagcgc       <  1:15985842/67‑1 (MQ=255)
                                                                   |                                                                   
TCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCAC  >  CP009273/2530215‑2530349

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