Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 130,120 T→C D673D (GAT→GAC acnB → bifunctional aconitate hydratase 2/2‑methylisocitrate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273130,1200TC97.9% 167.5 / ‑3.8 48D673D (GAT→GACacnBbifunctional aconitate hydratase 2/2‑methylisocitrate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (25/22);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.81e-01

GATGATCGCGGAAGGTTACGGCGATCGTCGTACCCTGGAACGTCGTATTCAGGGCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGATCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTACAGGGTGAGAAGATCGACGAAGTGTTTATCGGTTCCTGCATGACCAACATCGGTCACTTCCGT  >  NZ_CP009273/129988‑130255
                                                                                                                                    |                                                                                                                                       
gatgatCGCGGAAGGTTACGGCGATCGTCGTACCCTGGAACGTCGTATTCAGGGCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCg                                                                                                                                   >  2:647502/1‑139 (MQ=255)
               ttACGGCGATCGTCGTACCCTGGAACGTCGTATTCAGGGCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCa                                                                                                                    <  2:352687/139‑1 (MQ=255)
                  cGGCGATCGTCGTACCCTGGAACGTCGTATTCAGGGCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGATCTGGCGGATATTAAAGAGCCAAt                                                                                                                  <  1:345531/138‑1 (MQ=255)
                  cGGCGATCGTCGTACCCTGGAACGTCGTATTCAGGGCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAAt                                                                                                                  <  1:337083/138‑1 (MQ=255)
                      gATCGTCGTACCCTGGAACGTCGTATTCAGGGCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCtgt                                                                                                             <  1:190013/139‑1 (MQ=255)
                        tcgtcgTACCCTGGAACGTCGTATTCAGGGCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCtgtgt                                                                                                           >  1:700585/1‑139 (MQ=255)
                                   tGGAACGTCGTATTCAGGGCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGa                                                                                                >  1:657812/1‑139 (MQ=255)
                                     gAACGTCGTATTCAGGGCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGAcc                                                                                              <  2:469217/139‑1 (MQ=255)
                                     gAACGTCGTATTCAGGGCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGAcc                                                                                              >  2:550457/1‑139 (MQ=255)
                                     gAACGTCGTATTCAGGGCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGAcc                                                                                              >  2:162668/1‑139 (MQ=255)
                                      aaCGTCGTATTCAGGGCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGAccc                                                                                             >  1:295183/1‑139 (MQ=255)
                                        cgtcgtATTCAGGGCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCg                                                                                            >  1:666702/1‑138 (MQ=255)
                                          tcgtATTCAGGGCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGAt                                                                                         >  2:22217/1‑139 (MQ=255)
                                           cgtATTCAGGGCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGAt                                                                                         >  1:492276/1‑138 (MQ=255)
                                             tATTCAGGGCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGAc                                                                                      >  2:113350/1‑139 (MQ=255)
                                                    ggCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTc                                                                               >  2:112740/1‑139 (MQ=255)
                                                    ggCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTc                                                                               >  2:27480/1‑139 (MQ=255)
                                                    ggCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTc                                                                               >  2:636328/1‑139 (MQ=255)
                                                    ggCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAAAGACCCGGATGACGCGCGTc                                                                               >  2:467354/1‑139 (MQ=255)
                                                      cATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCg                                                                             <  1:658679/139‑1 (MQ=255)
                                                       aTGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGc                                                                            >  2:716777/1‑139 (MQ=255)
                                                        tGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCt                                                                           <  2:657812/139‑1 (MQ=255)
                                                            aaaaTGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCt                                                                       <  1:22217/139‑1 (MQ=255)
                                                             aaaTGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTg                                                                      <  1:570677/139‑1 (MQ=255)
                                                             aaaTGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTg                                                                      <  1:162668/139‑1 (MQ=255)
                                                             aaaTGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTg                                                                      <  2:306704/139‑1 (MQ=255)
                                                                  gctggcGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTa                                                                 <  2:492276/139‑1 (MQ=255)
                                                                   ctggcGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTAc                                                                <  1:778385/139‑1 (MQ=255)
                                                                   ctggcGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTAc                                                                <  1:311805/139‑1 (MQ=255)
                                                                       cGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCt                                                                           <  2:138691/124‑1 (MQ=255)
                                                                       cGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCt                                                                           >  1:138691/1‑124 (MQ=255)
                                                                       cGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGc                                                                     <  2:809740/130‑1 (MQ=255)
                                                                       cGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGc                                                                     >  1:809740/1‑130 (MQ=255)
                                                                                     ctgGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTACAGGGTGAGAAGATcgacg                                              <  1:432559/139‑1 (MQ=255)
                                                                                     ctgGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGCACAGGGTGAGAAGATcgac                                               >  2:179871/1‑138 (MQ=255)
                                                                                       ggAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTACAGGGTGAGAAGATcgacg                                              >  2:740640/1‑137 (MQ=255)
                                                                                               aTGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTACAGGGTGAGAAGATCGACGAAGTGTTTAt                                    <  2:295183/139‑1 (MQ=255)
                                                                                               aTGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTACAGGGTGAGAAGATCGACGAAGTGTTTAt                                    >  2:466179/1‑139 (MQ=255)
                                                                                               aTGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTACAGGGTGAGAAGATCGACGAAGTGTTTAt                                    <  1:636328/139‑1 (MQ=255)
                                                                                               aTGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTACAGGGTGAGAAGATCGACGAAGTGTTTAt                                    <  2:608359/139‑1 (MQ=255)
                                                                                                         ggAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTACAGGGTGAGAAGATCGACGAAGTGTTTATCGGTTCCTGc                          <  1:121178/139‑1 (MQ=255)
                                                                                                                   gCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTACAGGGTGAGAAGATCGACGAAGTGTTTATCGGTTCCTGCATGACCAACa                <  1:550457/139‑1 (MQ=255)
                                                                                                                       tGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTACAGGGTGAGAAGATCGACGAAGTGTTTATCGGTTCCTGCATGACCAACATCgg            >  1:410153/1‑139 (MQ=255)
                                                                                                                       tGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTACAGGGTGAGAAGATCGACGAAGTGTTTATCGGTTCCTGCATGACCAACATCgg            >  1:414627/1‑139 (MQ=255)
                                                                                                                         atcgacatcgacCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTACAGGGTGAGAAGATCGACGAAGTGTTTATCGGTTCCTGCATGACCAACATCGGt           <  1:112740/138‑1 (MQ=255)
                                                                                                                         atcgacatcgacCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTACAGGGTGAGAAGATCGACGAAGTGTTTATCGGTTCCTGCATGACCAACATCGGt           <  2:511058/138‑1 (MQ=255)
                                                                                                                                 cGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTACAGGGTGAGAAGATCGACGAAGTGTTTAt                                    >  1:494623/1‑105 (MQ=255)
                                                                                                                                 cGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTACAGGGTGAGAAGATCGACGAAGTGTTTAt                                    <  2:494623/105‑1 (MQ=255)
                                                                                                                                 cGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTACAGGGTGAGAAGATCGACGAAGTGTTTATCGGTTCCTGCATGACCAACATCGGTCACTTCCGt  >  1:540581/1‑139 (MQ=255)
                                                                                                                                 cGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTACAGGGTGAGAAGATCGACGAAGTGTTTATCGGTTCCTGCATGACCAACATCGGTCACTTCCGt  >  2:370690/1‑139 (MQ=255)
                                                                                                                                    |                                                                                                                                       
GATGATCGCGGAAGGTTACGGCGATCGTCGTACCCTGGAACGTCGTATTCAGGGCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGATCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTACAGGGTGAGAAGATCGACGAAGTGTTTATCGGTTCCTGCATGACCAACATCGGTCACTTCCGT  >  NZ_CP009273/129988‑130255

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