Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,121,319 G→A V159V (GTC→GTT gmd ← GDP‑mannose 4,6‑dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,121,3190GA98.2% 214.5 / ‑4.2 57V159V (GTC→GTTgmdGDP‑mannose 4,6‑dehydratase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (25/31);  total (25/32)
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.61e-01

TTGCGCGGGTGATTTTGCGGGTAACGAAGGTTTCGCCGCGGCGCGGGGATTCATGGTTGAAGAGAATTCCGTTACAGGCGTACATGCCGTAGGATTCACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCGACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATACAGTTCAGAGGTGGAAGCCTGATAGAAACGAGTTTTCTTTTCCAGACCGAGGAAGCGGATCGC  >  NZ_CP009273/2121183‑2121450
                                                                                                                                        |                                                                                                                                   
ttGCGCGGGTGATTTTGCGGGTAACGAAGGTTTCGCCGCGGCGCGGGGATTCATGGTTGAAGAGAATTCCGTTACAGGCGTACATGCCGTAGGATTCACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAAc                                                                                                                                   <  1:1148903/139‑1 (MQ=255)
 tGCGCGGGTGATTTTGCGGGTAACGAAGGTTTCGCCGCGGCGCGGGGATTCATGGTTGAAGAGAATTCCGTTACAGGCGTACATGCCGTAGGATTCACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAAc                                                                                                                                   <  2:242361/138‑1 (MQ=255)
 tGCGCGGGTGATTTTGCGGGTAACGAAGGTTTCGCCGCGGCGCGGGGATTCATGGTTGAAGAGAATTCCGTTACAGGCGTACATGCCGTAGGATTCACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAAc                                                                                                                                   >  1:242361/1‑138 (MQ=255)
        gTGATTTTGCGGGTAACGAAGGTTTCGCCGCGGCGCGGGGATTCATGGTTGAAGAGAATTCCGTTACAGGCGTACATGCCGTAGGATTCACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACg                                                                                                                           <  1:30921/139‑1 (MQ=255)
         tGATTTTGCGGGTAACGAAGGTTTCGCCGCGGCGCGGGGATTCATGGTTGAAGAGAATTCCGTTACAGGCGTACATGCCGTAGGATTCACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACgg                                                                                                                          <  1:863695/139‑1 (MQ=255)
              ttGCGGGTAACGAAGGTTTCGCCGCGGCGCGGGATTCATGGTTGAAGAGAATTCCGTTACAGGCGTACATGCCGTAGGATTCACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATc                                                                                                                     <  1:474953/138‑1 (MQ=255)
               tGCGGGTAACGAAGGTTTCGCCGCGGCGCGGGGATTCATGGTTGAAGAGAATTCCGTTACAGGCGTACATGCCGTAGGATTCACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATc                                                                                                                     <  2:297126/138‑1 (MQ=255)
               tGCGGGTAACGAAGGTTTCGCCGCGGCGCGGGGATTCATGGTTGAAGAGAATTCCGTTACAGGCGTACATGCCGTAGGATTCACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATc                                                                                                                     <  1:160297/138‑1 (MQ=255)
                  gggTAACGAAGGTTTCGCCGCGGCGCGGGGATTCATGGTTGAAGAGAATTCCGTTACAGGCGTACATGCCGTAGGATTCACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATggc                                                                                                                   >  2:86597/1‑137 (MQ=255)
                  gggTAACGAAGGTTTCGCCGCGGCGCGGGGATTCATGGTTGAAGAGAATTCCGTTACAGGCGTACATGCCGTAGGATTCACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATcgc                                                                                                                   <  1:86597/137‑1 (MQ=255)
                             gTTTCGCCGCGGCGCGGGGATTCATGGTTGAAGAGAATTCCGTTACAGGCGTACATGCCGTAGGATTCACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCGACCGCATACGGAGATCGCGGGTAGAACGGCg                                                                                                      <  1:301775/139‑1 (MQ=255)
                             gTTTCGCCGCGGCGCGGGGATTCATGGTTGAAGAGAATTCCGTTACAGGCGTACATGCCGTAGGATTCACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCg                                                                                                      <  1:1160148/139‑1 (MQ=255)
                               ttCGCCGCGGCGCGGGGATTCATGGTTGAAGAGAATTCCGTTACAGGCGTACATGCCGTAGGATTCACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCGTg                                                                                                    >  2:340637/1‑139 (MQ=255)
                                tCGCCGCGGCGCGGGGATTCATGGTTGAAGAGAATTCCGTTACAGGCGTACATGCCGTAGGATTCACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCGTgg                                                                                                   >  1:1149954/1‑139 (MQ=255)
                                      cggcgcggGGATTCATGGTTGAAGAGAATTCCGTTACAGGCGTACATGCCGTAGGATTCACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCtt                                                                                             >  2:345107/1‑139 (MQ=255)
                                              ggATTCATGGTTGAAGAGAATTCCGTTACAGGCGTACATGCCGTAGGATTCACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGa                                                                                     <  2:128692/139‑1 (MQ=255)
                                                    aTGGTTGAAGAGAATTCCGTTACAGGCGTACATGCCGTAGGATTCACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTcc                                                                               <  1:1079756/139‑1 (MQ=255)
                                                    aTGGTTGAAGAGAATTCCGTTACAGGCGTACATGCCGTAGGATTCACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTcc                                                                               <  2:414028/139‑1 (MQ=255)
                                                    aTGGTTGAAGAGAATTCCGTTACAGGCGTACATGCCGTAGGATTCACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTc                                                                                <  1:1179208/138‑1 (MQ=255)
                                                    aTGGTTGAAGAGAATTCCGTTACAGGCGTACATGCCGTAGGATTCACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTc                                                                                <  2:597838/138‑1 (MQ=255)
                                                    aTGGTTGAAGAGAATTCCGTTACAGGCGTACATGCCGTAGGATTCACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTc                                                                                <  2:552088/138‑1 (MQ=255)
                                                        ttGAAGAGAATTCCGTGACAGGCGTTCATGCCGTAGGATTCACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAGCGGCGTGGTCTCTTTCTGCGGAATTTACTgtg                                                                           <  2:1052326/139‑3 (MQ=255)
                                                          gAAGAGAATTCCGTTACAGGCGTACATGCCGTAGGATTCACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCAc                                                                          >  2:151100/1‑138 (MQ=255)
                                                                  ttCCGTTACAGGCGTACATGCCGTAGGATTCACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATa                                                                 <  2:1024591/139‑1 (MQ=255)
                                                                  ttCCGTTACAGGCGTACATGCCGTAGGATTCACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATa                                                                 >  1:916506/1‑139 (MQ=255)
                                                                  ttCCGTTACAGGCGTACATGCCGTAGGATTCACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATa                                                                 <  1:664820/139‑1 (MQ=255)
                                                                  ttCCGTTACAGGCGTACATGCCGTAGGATTCACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATa                                                                 <  1:340637/139‑1 (MQ=255)
                                                                  ttCCGTTACAGGCGTACATGCCGTAGGATTCACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATa                                                                 >  1:300368/1‑139 (MQ=255)
                                                                  ttCCGTTACAGGCGTACATGCCGTAGGATTCACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATa                                                                 <  1:151100/139‑1 (MQ=255)
                                                                            ggCGTACATGCCGTAGGATTCACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATACAGTTCagag                                                       >  1:648162/1‑139 (MQ=255)
                                                                                         tAGGATTCACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCg                                                                                                      <  2:1040513/79‑1 (MQ=255)
                                                                                         tAGGATTCACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCg                                                                                                      >  1:1040513/1‑79 (MQ=255)
                                                                                                cACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATACAGTTCAGAGGTGGAAGCCTGATAGAAAc                                    >  2:822403/1‑138 (MQ=255)
                                                                                                 aCGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATc                                                                                                                     >  1:87090/1‑56 (MQ=255)
                                                                                                 aCGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATc                                                                                                                     <  2:87090/56‑1 (MQ=255)
                                                                                                 aCGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATa                                                                 >  2:107717/1‑108 (MQ=255)
                                                                                                 aCGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATa                                                                 <  1:107717/108‑1 (MQ=255)
                                                                                                 aCGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATACAGTTCAGAGGTGGAAGCCTGATAGAAAc                                    >  1:853271/1‑137 (MQ=255)
                                                                                                  cGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATACAGTTCAGAGGTGGAAGCCTGATAGAAACGAg                                 <  2:1035814/139‑1 (MQ=255)
                                                                                                        ttAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCg                                                                                                      >  1:728239/1‑64 (MQ=255)
                                                                                                        ttAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCg                                                                                                      <  2:1286751/64‑1 (MQ=255)
                                                                                                        ttAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCg                                                                                                      >  1:1286751/1‑64 (MQ=255)
                                                                                                        ttAACGGTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCg                                                                                                      <  2:728239/64‑1 (MQ=255)
                                                                                                             ggTGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATACAGTTCAGAGGTGGAAGCCTGATAGAAACGAGttttcttttc                       >  1:536426/1‑138 (MQ=255)
                                                                                                               tGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCGTGGtc                                                                                                 <  1:826524/62‑1 (MQ=255)
                                                                                                               tGATCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCGTGGtc                                                                                                 >  2:826524/1‑62 (MQ=255)
                                                                                                                 aTCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCa                                                                           >  2:528811/1‑82 (MQ=255)
                                                                                                                 aTCCAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCa                                                                           <  1:528811/82‑1 (MQ=255)
                                                                                                                    cAGTAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATACAGTTCAGAGGTGGAAGCCTGATAGAAACGAGTTTTCTTTTCCAGACCg                <  2:811091/138‑1 (MQ=255)
                                                                                                                       tAGGCGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATACAGTTCAGAGGTGGAAGCCTGATAGAAACGAGTTTTCTTTTCCAGACCGAgg             >  1:114333/1‑138 (MQ=255)
                                                                                                                           cGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATACAGTTCag                                                         <  1:1064964/90‑1 (MQ=255)
                                                                                                                           cGTACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATACAGTTCag                                                         >  2:1064964/1‑90 (MQ=255)
                                                                                                                             tACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATACAGTTCAGAGGTGGAAGCCTGATAGAAACGAGTTTTCTTTTCCAGACCGAGGAAGCGGa      >  1:1085336/1‑139 (MQ=255)
                                                                                                                             tACAGTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATACAGTTCAGAGGTGGAAGCCTGATAGAAACGAGTTTTCTTTTCCAGACCGAGGAAGCGGa      >  2:916921/1‑139 (MQ=255)
                                                                                                                                 gTTTGGCAACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATACAGTTCAGAGGTGGAAGCCTGATAGAAACGAGTTTTCTTTTCCAGACCGAGGAAGCGGATCGc  >  2:95856/1‑139 (MQ=255)
                                                                                                                                     ggCAACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATACAGTTc                                                           <  2:1195723/78‑1 (MQ=255)
                                                                                                                                     ggCAACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATACAGTTc                                                           >  1:1195723/1‑78 (MQ=255)
                                                                                                                                        |                                                                                                                                   
TTGCGCGGGTGATTTTGCGGGTAACGAAGGTTTCGCCGCGGCGCGGGGATTCATGGTTGAAGAGAATTCCGTTACAGGCGTACATGCCGTAGGATTCACGGTAGTTAACGGTGATCCAGTAGGCGTACAGTTTGGCGACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATACAGTTCAGAGGTGGAAGCCTGATAGAAACGAGTTTTCTTTTCCAGACCGAGGAAGCGGATCGC  >  NZ_CP009273/2121183‑2121450

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