Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,796,969 G→A 10.0% R39H (CGT→CAT)  fabZ → 3‑hydroxyacyl‑ACP dehydratase FabZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,796,9690GA10.0% 146.1 / 3.7 55R39H (CGT→CAT) fabZ3‑hydroxyacyl‑ACP dehydratase FabZ
Reads supporting (aligned to +/- strand):  ref base G (24/25);  new base A (4/2);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGCTGGCTCCTCTTTGGATTGCAAAAGGAGCGTGTGGTCAGCACCTGCGCTCCCTATTCTTATTACAGGCTTCCCCACCGAAATGATGGACATCAACGAGATTCGCGAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTCAGCATCAATGAGCCGTTTTTCAATGGCCACTTCCCGGCACACCCGATCATGCCGGGTGTCCTGATCATCGAAGCCATGGCCCAGGCGGCCGGCATC  >  NC_002947/1796822‑1797117
                                                                                                                                                   |                                                                                                                                                    
tAGCTGGCTCCTCTTTGGATTGCAAAAGGAGCGTGTGGTCAGCACCTGCGCTCCCTATTCTTATTACAGGCTTCCCCACCGAAATGATGGACATCAACGAGATTCGCGAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCgtg                                                                                                                                                    <  1:1035360/150‑1 (MQ=255)
    tGGCTCCTCTTTGGATTGCAAAAGGAGCGTGTGGTCAGCACCTGCGCTCCCTATTCTTATTACAGGCTTCCCCACCGCCCTGATGGACATCAACGACATACTCGGAGATCTGCCGCGCCGCTACCCGTTCCTGCTGGTAGACCGTGTGAc                                                                                                                                                <  1:387789/150‑1 (MQ=255)
     ggCTCCTCTTTGGATTGCAAAAGGAGCGGGTGGTCAGCACCTGCGCTCCCTATTCTTATTACAGGCTTCCCCACCGAAATGATGGACATCAACGAGATTCGCGAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACg                                                                                                                                               >  1:1320423/1‑150 (MQ=255)
           tctTTGGATTGCAAAAGGAGCGTGTGGTCAGCACCTGCGCTCCCTATTCTTATTACAGGCTTCCCCACCGAAATGATGGACATCAACGAGATTCGCGAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTg                                                                                                                                         <  1:1151964/150‑1 (MQ=255)
             tttGGATTGCAAAAGGAGCGTGTGGTCAGCACCTGCGCTCCCTATTCTTATTACAGGCTTCCCCACCGAAATGATGGACATCAGCGAGATTCGCGAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGa                                                                                                                                       <  2:1136676/150‑1 (MQ=255)
             tttGGATTGCAAAAGGAGCGTGTGGTCAGCACCTGCGCTCCCTATTCTTATTACAGGCTTCCCCACCGAAATGATGGACATCAACGAGATTCGCGAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGa                                                                                                                                       <  2:68581/150‑1 (MQ=255)
                 gATTGCAAAAGGAGCGTGTGGTCAGCACCTGCGCTCCCTATTCTTATTACAGGCTTCCCCACCGAAATGATGGACATCAACGAGATTCGCGAATACCTGCCTCACCGCTCGAAACGCCTGCTGGTAGACCGTGTGACGGATCTGGACTTc                                                                                                                                   >  1:784909/1‑150 (MQ=255)
                     gCAAAAGGAGCGTGTGGTCAGCACCTGCGCTCCCTATTCTTATTACAGGCTTCCCCACCGAAATGATGGACATCAACGAGATTCGCGAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAgg                                                                                                                               <  1:13537/150‑1 (MQ=255)
                              gCGTGTGGTCAGCACCTGCGCTCCCTATTCTTATTACAGGCTTCCCCACCGAAATGATGGACATCAACGAGATTCGCGAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGGCCTCGAAGTCCagatcc                                                                                                                      <  1:1260969/150‑4 (MQ=255)
                                   tgGTCAGCACCTGCGCTCCCTATTCTTATTACAGGCTTCCCCACCGAAATGATGGACATCAACGAGATTCGCGAATACCTGCCTCACCGCTACCAGTTCCTGATGGTAGACCAGGTGACGGACGTGGACTGCGATGCCCAGAGCGTTCGt                                                                                                                 >  2:867225/1‑150 (MQ=255)
                                                  cTCCCTATTCTTATTACAGGCTTCCCCACCGAAATGATGGACATCAACGAGATTCGCGAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTc                                                                                                  <  2:901374/150‑1 (MQ=255)
                                                   tCCCTATTCTTATTACAGGCTTCCCCACCGAAATGATGGACATCAACGAGATTCGCGAATACCTGCCGCACCGCTACCCGTTCCTGCTGGTAGACCATATGACGGATCTGGCCTTCGAGGCCCAGAGCATTCGTGCCTACCAGAATGTcc                                                                                                 >  2:149694/1‑149 (MQ=255)
                                                     ccTATTCTTATTACAGGCTTCCCCACCGAAATGATGGACATCAACGAGATTCGCGAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTCAGc                                                                                               >  2:604264/1‑150 (MQ=255)
                                                           cTTATTACAGGCTTCCCCACCGAAATGATGGACATCAACGAGATTCGCGAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTCAGCATCAAt                                                                                         >  1:710549/1‑150 (MQ=255)
                                                                     gCTTCCCCACCGAAATGATGGACATCAACGAGATTCGCGAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTCAGCATCAATGAGCCGtttt                                                                               >  2:1240048/1‑150 (MQ=255)
                                                                       ttCCCCACCGAAATGATGGACATCAACGAGATTCGCGAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTCAGCATCAATGAGCCGTTTTTc                                                                             >  2:333884/1‑150 (MQ=255)
                                                                          cccACCGAAATGATGGACATCAACGAGATTCGCGAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTCAGCATCAATGAGCCGTTTTTCAAt                                                                          >  1:298284/1‑150 (MQ=255)
                                                                           ccACCGAAATGATGGACATCAACGAGATTCGCGAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAAt                                                                                                     >  2:1064990/1‑122 (MQ=255)
                                                                           ccACCGAAATGATGGACATCAACGAGATTCGCGAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAAt                                                                                                     <  1:1064990/122‑1 (MQ=255)
                                                                           ccACCGAAATGATGGACATCAACGAGATTCGCGAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTCAGCATCAATGAGCCGTTTTTCAATg                                                                         <  1:867225/150‑1 (MQ=255)
                                                                           ccACCGAAATGATGGACATCAACGAGATTCGCGAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTCAGCATCAATGAGCCGTTTTTCAATg                                                                         >  2:1321465/1‑150 (MQ=255)
                                                                             aCCGAAATGATGGACATCAACGAGATTCGCGAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTCAGCATCAATGAGCCGTTTTTCAATGGc                                                                       <  1:604264/150‑1 (MQ=255)
                                                                             aCCGAAATGATGGACATCAACGAGATTCGCGAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTCAGCATCAATGAGCCGTTTTTCAATGGc                                                                       >  1:1127891/1‑150 (MQ=255)
                                                                              ccGAAATGATGGACATCAACGAGATTCGCGAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTCAGCATCAATGAGCCGTTTTTCAATGGcc                                                                      <  2:1205457/150‑1 (MQ=255)
                                                                                gAAATGATGGACATCAACGAGATTCGCGAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTg                                                                                                                <  1:1072802/106‑1 (MQ=255)
                                                                                gAAATGATGGACATCAACGAGATTCGCGAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTg                                                                                                                >  2:1072802/1‑106 (MQ=255)
                                                                                        ggACATCAACGAAATTCGCGAAGACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTCAGCATCAATGAGCCGTTTTTCAATGGCCACTTCCCGGc                                                            >  2:448772/1‑150 (MQ=255)
                                                                                               aaCGAGATTCGCGAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTCAGCATCAATGAGCCGTTTTTCAATGGCCACTTCCCGGCACACCCg                                                     <  2:710549/150‑1 (MQ=255)
                                                                                               aaCGAGATTCGCGAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTCAGCATCAATGAGCCGTTTTTCAATGGCCACTTCCCGGCACACCCg                                                     >  1:1377493/1‑150 (MQ=255)
                                                                                                aCGAGATTCGCGAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTCAGCATCAATGAGCCGTTTTTCAATGGCCACTTCCCGGCACACCCGa                                                    <  2:1377493/150‑1 (MQ=255)
                                                                                                aCGAGATTCGCGAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTCAGCATCAATGAGCCGTTTTTCAATGGCCACTTCCCGGCACACCCGa                                                    >  2:50161/1‑150 (MQ=255)
                                                                                                aCGAGATTCGCGAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTCAGCATCAATGAGCCGTTTTTCAATGGCCACTTCCCGGCACACCCGa                                                    <  1:333884/150‑1 (MQ=255)
                                                                                                     aTTCGCGAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCATGTGACGGATCTGGACTTCGAGGcc                                                                                                                             >  1:668861/1‑72 (MQ=255)
                                                                                                     aTTCGCGAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCATGTGACGGATCTGGACTTCGAGGcc                                                                                                                             >  1:669920/1‑72 (MQ=255)
                                                                                                     aTTCGCGAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCATGTGACGGATCTGGACTTCGAGGcc                                                                                                                             <  2:669920/72‑1 (MQ=255)
                                                                                                     aTTCGCGAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCATGTGACGGATCTGGACTTCGAGGcc                                                                                                                             <  2:668861/72‑1 (MQ=255)
                                                                                                      ttCGCGAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTCAGCATCAATGAGCCGTTTTTCAATGGCCACTTCCCGGCACACCCGATCATGc                                              <  1:149694/150‑1 (MQ=255)
                                                                                                           gAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTCAGCATCAATGAGCCGTTTTTCAATGGCCACTTCCCGGCACACCCGATCATGCCGGgt                                         <  1:1240048/150‑1 (MQ=255)
                                                                                                              tACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGt                                                                                                   >  2:964858/1‑89 (MQ=255)
                                                                                                              tACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGt                                                                                                   <  1:964858/89‑1 (MQ=255)
                                                                                                              tACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTCAGCATCAATGAGCCGTTTTTCAATGGCCACTTCCCGGCACACCCGATCATGCCGGGTGTc                                      <  2:956462/150‑1 (MQ=255)
                                                                                                                   gCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACaa                                                                                                         >  2:365309/1‑78 (MQ=255)
                                                                                                                   gCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACaa                                                                                                         <  1:365309/78‑1 (MQ=255)
                                                                                                                    ccTCACCGCTACCCGTTCCTGCTGGTAGGCCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTCAGCATCAATGAGCCGTTTTTCAATGGCCACTTCCCGGCACACCCGATCATGCCGGGTGTCCTGatc                                >  1:928803/1‑150 (MQ=255)
                                                                                                                     cTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTCAGCATCAATGAGCCGTTTTTCAATGGTCACTTCCCGGCacac                                                        <  2:278536/125‑1 (MQ=255)
                                                                                                                     cTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTCAGCATCAATGAGCCGTTTTTCAATGGTCACTTCCCGGCacac                                                        >  1:278536/1‑125 (MQ=255)
                                                                                                                         ccGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTCAGCATCAATGATCCGTTTTTCAATGGCCACTTCCCGGCACACCCGATCATGCCGGGTGTCCTGATCATCGa                           >  1:144737/1‑150 (MQ=255)
                                                                                                                          cGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTCAGCATCAATGa                                                                                       <  2:30674/89‑1 (MQ=255)
                                                                                                                          cGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTCAGCATCAATGa                                                                                       >  1:30674/1‑89 (MQ=255)
                                                                                                                          cGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTCAGCATCAATGAGCCGTTTTTCAATGGCCACTTCCCGGCACACCCGATCATGCCGGGTGTCCTGATCATCGaa                          >  2:1260916/1‑150 (MQ=255)
                                                                                                                             tACCCGTTCCTGCTGGTAGACCGTGTGACCGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTCAGCATCAATGAGCCGTTTTTCAATGGCCACTTCCCGGCACACCCGATCATGCCGGGTGTCCTGATCATCGAAGcc                       <  2:784909/150‑1 (MQ=255)
                                                                                                                              aCCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTCAGCATCAATGAGCCGTTTTTCAATGGCCACTTCCCGGCACACCCGATCATGCCGGGTGTCCTGATCATCGAAGCCa                      <  2:1320423/150‑1 (MQ=255)
                                                                                                                                ccGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTCAGCATCAATGAGCCGTTTTTCAATGGCCACTTCCCGGCACACCCGATCATGCCGGGTGTCCTGATCATCGAAGCCATg                    >  1:1234297/1‑150 (MQ=255)
                                                                                                                                ccGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTCAGCATCAATGAGCCGCCTTCCAAGGGCCGCGTCCCGGCACACCCTATCATGCCGGGTGTCCTGATCATCGAAGCCATg                    >  2:1403917/1‑150 (MQ=255)
                                                                                                                                                  cGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTCAGCATCAATGAGCCGTTTTTCAATGGCCACTTCCCGGCACACCCGATCATGCCGGGTGTCCTGATCATCGAAGCCATGGCCCAGGCGGCCGGCATc  <  1:1260916/150‑1 (MQ=255)
                                                                                                                                                  cGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTCAGCATCAATGAGCCGTTTTTCAATGGCCACTTCCCGGCACACCCGATCATGCCGGGTGTCCTGATCATCGAAGCCATGGCCCAGGCGGCCGGCATc  >  2:1208093/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TAGCTGGCTCCTCTTTGGATTGCAAAAGGAGCGTGTGGTCAGCACCTGCGCTCCCTATTCTTATTACAGGCTTCCCCACCGAAATGATGGACATCAACGAGATTCGCGAATACCTGCCTCACCGCTACCCGTTCCTGCTGGTAGACCGTGTGACGGATCTGGACTTCGAGGCCCAGAGCATTCGTGCCTACAAGAATGTCAGCATCAATGAGCCGTTTTTCAATGGCCACTTCCCGGCACACCCGATCATGCCGGGTGTCCTGATCATCGAAGCCATGGCCCAGGCGGCCGGCATC  >  NC_002947/1796822‑1797117

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