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,1200TC95.2% 67.9 / ‑3.0 21D673D (GAT→GACacnBbifunctional aconitate hydratase 2/2‑methylisocitrate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (9/11);  minor base A (1/0);  total (10/11)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.80e-01

GTGGATGATCGCGGAAGGTTACGGCGATCGTCGTACCCTGGAACGTCGTATTCAGGGCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGATCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTACAGGGTGAGAAGATCGACGAAGTGTTTATCGGTTCCTGCATGACCAACATC  >  NZ_CP009273/129985‑130243
                                                                                                                                       |                                                                                                                           
gTGGATGATCGCGGAAGGTTACGGCGATCGTCGTACCCTGGAACGTCGTATTCAGGGCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTg                                                                                                                          >  2:394571/1‑139 (MQ=255)
gTGGATGATCGCGGAAGGTTACGGCGATCGTCGTACCCTGGAACGTCGTATTCAGGGCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTg                                                                                                                          <  2:158786/139‑1 (MQ=255)
gTGGATGATCGCGGAAGGTTACGGCGATAGTCGTACCCTGGAACGTCGTATTCAGGGCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACCTCGAACTg                                                                                                                          >  2:173068/1‑139 (MQ=255)
     tgatCGCGGAAGGTTACGGCGATCGTCGTACCCTGGAACGTCGTATTCAGGGCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGa                                                                                                                     >  2:221357/1‑139 (MQ=255)
           cgGAAGGTTACGGCGATCGTCGTACCCTGGAACGTCGTATTCAGGGCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTaa                                                                                                               <  1:394571/139‑1 (MQ=255)
           cgGAAGGTTACGGCGATCGTCGTACCCTGGAACGTCGTATTCAGGGCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTaa                                                                                                               <  2:10285/139‑1 (MQ=255)
                   tACGGCGATCGTCGTACCCTGGAACGTCGTATTCAGGGCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCaa                                                                                                       >  1:318621/1‑139 (MQ=255)
                         gATCGTCGTACCCTGGAACGTCGTATTCAGGGCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCtgt                                                                                                 <  1:102005/139‑1 (MQ=255)
                                        gAACGTCGTATTCAGGGCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGAcc                                                                                  >  2:203662/1‑139 (MQ=255)
                                                          aTGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCTGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGc                                                                >  1:166964/1‑139 (MQ=255)
                                                                 aaTGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGc                                                         <  1:2177/139‑1 (MQ=255)
                                                                   tggctggcGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCgg                                                       >  2:110659/1‑139 (MQ=255)
                                                                                   aGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTACAGGGTgag                                            >  2:389643/1‑134 (MQ=255)
                                                                                                  aTGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTACAGGGTGAGAAGATCGACGAAGTGTTTAt                        <  2:281161/139‑1 (MQ=255)
                                                                                                  aTGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTACAGGGTGAGAAGATCGACGAAGTGTTTAt                        <  1:400441/139‑1 (MQ=255)
                                                                                                           cGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTACAGGGTGAGAAGATCGACGAAGTGTTTATCGGTTCCTg               >  2:226521/1‑139 (MQ=255)
                                                                                                            ggAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTACAGGGTGAGAAGATCGACGAAGTGTTTATCGGTTCCTGc              <  2:431831/139‑1 (MQ=255)
                                                                                                            ggAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTACAGGGTGAGAAGATCGACGAAGTGTTTATCGGTTCCTGc              <  2:240032/139‑1 (MQ=255)
                                                                                                            ggAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTACAGGGTGAGAAGATCGACGAAGTGTTTATCGGTTCCTGc              <  1:110659/139‑1 (MQ=255)
                                                                                                                      gCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTACAGGGTGAGAAGATCGACGAAGTGTTTATCGGTTCCTGCATGACCAACa    <  1:232582/139‑1 (MQ=255)
                                                                                                                        aGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTACAGGGTGAGAAGATCGACGAAGTGTTTATCGGTTCCTGCATGACCAACATc  >  1:226758/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                           
GTGGATGATCGCGGAAGGTTACGGCGATCGTCGTACCCTGGAACGTCGTATTCAGGGCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGATCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTACAGGGTGAGAAGATCGACGAAGTGTTTATCGGTTCCTGCATGACCAACATC  >  NZ_CP009273/129985‑130243

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