Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,158,088 G→A 7.6% R132H (CGT→CAT)  fabG → 3‑oxoacyl‑ACP reductase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,158,0880GA 7.6% 138.8 / 2.8 53R132H (CGT→CAT) fabG3‑oxoacyl‑ACP reductase subunit
Reads supporting (aligned to +/- strand):  ref base G (27/21);  new base A (2/2);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

GGCCATCTTGGTCAACAATGCCGGTATCACCCGCGACAACCTCATGCTGCGCATGAAAGACGACGAGTGGTTTGATGTCATCGACACCAACCTGAACAGCCTCTACCGTCTGTCCAAGGGCGTGCTGCGTGGCATGACCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATGGGTAACGCAGGTCAGGCCAACTACGCGGCTGCCAAGGCCGGTCTGGAAGGTTTCAGCCGCGCCCTGGCGCGTGAAGTGGGTTCGCGTGGTATCACCGTCAACTCGG  >  NC_002947/2157942‑2158237
                                                                                                                                                  |                                                                                                                                                     
ggCCATCTTGGTCAACAATGCCGGTATCACCCGCGACAACCTCATGCTGCGCATGAAAGACGACGAGTGGTTTGATGTCATCGACACCAACCTGAACAGCCTCTACCGTCTGTCCAAGGGCGTGCTGCGTGGCATGACCAAGGCGCGTTg                                                                                                                                                    <  2:481028/150‑1 (MQ=255)
   cATCTTGGTCAACAATGCCGGTATCACCCGCGACAACCTCATGCTGCGCATGAAAGACGACGAGTGGTTTGATGTCATCGACACCGACCTGAACAGCCTCTATCGTCTGTCCAGAGTCGTGCCGCGTGGCATGACCAAGACGCGGCgggg                                                                                                                                                 >  2:280309/1‑150 (MQ=255)
   cATCTTGGTCAACAATGCCGGTATCACCCGCGACAACCTCATGCTGCGCATGAAAGACGACGAGTGGTTTGATGTCATCGACACCAACCTGAACAGCCTCTACCGTCTGTCCAAGGGCGTGCTGCGTGGCATGACCAAGGCGCGTTgggg                                                                                                                                                 >  1:64235/1‑150 (MQ=255)
       ttGGTCAACAATGCCGGTATCACCCGCGACAACCTCATGCTGCGCATGAAAGACGACGAGTGGTTTGATGTCATCGACACCAACCTGAACAGCCTCTACCGTCTGTCCAAGGGCGTGCTGCGTGGCATGACCAAGGCGCGTTGGGGTCGt                                                                                                                                             >  2:170327/1‑150 (MQ=255)
        tGGTCAACAATGCCGGTATCACCCGCGACAACCTCATGCTGCGCATGAAAGACGACGAGTGGTTTGATGTCATCGACACCAACCTGAACAGCCTCTACCGTCTGTCCAAGGGCGTGCTGCGTGGCATGACCAAGGCGCGTTGGGGTCGTa                                                                                                                                            <  1:916837/150‑1 (MQ=255)
           tCAACAATGCCGGTATCACCCGCGACAACCTCATGCTGCGCATGAAAGACGACGAGTGGTTTGATGTCATCGACACCAACCTGAACAGCCTCTACCGTCTGTCCAAGGGCGTGCTGCGTGGCATGACCAAGGCGCGTTGGGGTCGTatca                                                                                                                                         <  2:666438/150‑1 (MQ=255)
             aacaaTGCCGGTATCACCCGCGACAACCTCATGCTGCGCATGAAAGACGACGAGTGGTTTGATGTCATCGACACCAACCTGAACAGCCTCTACCGTCTGTCCAAGGGCGTGCTGCGTGGCATGACCAAGGCGCGTTGGGGTCGTatcatc                                                                                                                                       >  1:1118607/1‑150 (MQ=255)
               caaTGCCGGTATCACCCGCGACAACCTCATGCTGCGCATGCAAGACGACGAGTGGTTTGATGTCATCGACACCAACCTGAACAGCCTCTACCGTCTGTCCAAGGGTGTGCTGCGTGGCATGACCAAGGCGCGTTGGGGACGTATCATCAg                                                                                                                                     >  1:1225573/1‑150 (MQ=255)
                        tATCACCCGCGACAACCTCATGCTGCGCATGAAAGACGACGAGTGGTTTGATGTCATCGACACCAACCTGAACAGCCTCTACCGTCTGTCCAAGGGCGTGCTGCGTGGCATGACCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTc                                                                                                                            <  1:280309/150‑1 (MQ=255)
                             cccGCGACAACCTCATGCTGCGCATGAAAGACGACGAGTGGTTTGATGTCATCGACACCAACCTGAACAGCCTCTACCGTCTGTCCAAGGGCGTGCTGCGTGGCATGACCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCg                                                                                                                       >  2:401013/1‑150 (MQ=255)
                              ccGCGACAACCTCATGCTGCGCATGAAAGACGACGAGTGGTTTGATGTCATCGACACCAACCTGAACAGCCTCTACCGTCTGTCCAAGGGCGTGCTGCGTGGCATGACCAAGGCGCGTTggg                                                                                                                                                  <  2:868936/122‑1 (MQ=255)
                              ccGCGACAACCTCATGCTGCGCATGAAAGACGACGAGTGGTTTGATGTCATCGACACCAACCTGAACAGCCTCTACCGTCTGTCCAAGGGCGTGCTGCGTGGCATGACCAAGGCGCGTTggg                                                                                                                                                  >  1:868936/1‑122 (MQ=255)
                              ccGCGACAACCTCATGCTGCGCATGAAAGACGACGAGTGGTTTGATGTCATCGACACCAACCTGAACAGCCTCTACCGTCTGTCCAAGGGCGTGCTGCGTGGCATGACCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGt                                                                                                                      >  1:240704/1‑150 (MQ=255)
                                gcgACAACCTCATGCTGCGCATGAAAGACGACGAGTGGTTTGATGTCATCGACACCAACCTGAACAGCCTCTACCGTCTGTCCAAGGGCGTGCTGCGTGGCATGACCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTg                                                                                                                    >  2:551488/1‑150 (MQ=255)
                                       ccTCATGCTGCGCATGAAAGACGACGAGTGGTTTGATGTCATCGACACCAACCTGAACAGCCTCTACCGTCTGTCCAAGGGCGTGCTGCGTGGCATGACCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCAt                                                                                                             >  1:924983/1‑150 (MQ=255)
                                       ccTCATGCTGCGCATGAAAGACGACGAGTGGTTTGATGTCATCGACACCAACCTGAACAGCCTCTACCGTCTGTCCAAGGGCGTGCTGCGTGGCATGACCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCAt                                                                                                             >  2:894838/1‑150 (MQ=255)
                                         tCATGCTGCGCATGAAAGACGACGAGTGGTTTGATGTCATCGACACCAACCTGAACAGCCTCTACCGTCTGTCCAAGGGCGTGCTGCGTGGCATGACCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATgg                                                                                                           <  2:9732/150‑1 (MQ=255)
                                         tCATGCTGCGCATGAAAGACGACGAGTGGTTTGATGTCATCGACACCAACCTGAACAGCCTCTACCGTCTGTCCAAGGGCGTGCTGCGTGGCATGACCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATgg                                                                                                           <  1:170327/150‑1 (MQ=255)
                                                             gacgaGTGGTTTGATGTCATCGACACCAACCTGAACAGCCTCTACCGTCTGTCCAAGGGCGTGCTGCGTGGCATGACCAAGGCGCTTTGGGGTCGTATCATCAGCATCGTCTCGGTCGTTGTTGCCATGGTTAACTCAGGTCAGGCCAAc                                                                                       >  2:269973/1‑150 (MQ=255)
                                                                   tGGTTTGATGTCATCGACACCAACCTGAACAGCCTCTACCGTCTGTCCAAGGGCGTGCTGCGTGGCATGACCAAGGCGCGTTGGGGTCGTatcatc                                                                                                                                       <  1:271502/96‑1 (MQ=255)
                                                                   tGGTTTGATGTCATCGACACCAACCTGAACAGCCTCTACCGTCTGTCCAAGGGCGTGCTGCGTGGCATGACCAAGGCGCGTTGGGGTCGTatcatc                                                                                                                                       >  2:271502/1‑96 (MQ=255)
                                                                      tttGATGTCATCGACACCAACCTGAACAGCCTCTACCGTCTGTCCAAGGGCGTGCTGCGTGGCATGACCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATGGGTAACGCAGGTCAGGCCAACTACGCGGCt                                                                              >  1:132709/1‑150 (MQ=255)
                                                                             tCATCGACACCAACCTGAACAGCCTCTACCGTCTGTCCAAGGGCGTGCTGCGTGGCATGACCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATGGGTAACGCAGGTCAGGCCAACTACGCGGCTGCCAAgg                                                                       >  1:741807/1‑150 (MQ=255)
                                                                                tCGACACCAACCTGAACAGCCTCTACCGTCTGTCCAAGGGTGTGCTGCGTGGCATGACCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATGGGTAACGCAGGTCAGGCCAACTACGCGGCTGCCAAGGCCg                                                                    <  2:1225573/150‑1 (MQ=255)
                                                                                  gACACCAACCTGAACAGCCTCTACCGTCTGTCCAAGGGCGTGCTGCGTGGCATGACCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATGGGTAACGCAGGTCAGGCCAACTACGCGGCTGCCAAGGCCGGt                                                                  <  2:924983/150‑1 (MQ=255)
                                                                                  gACACCAACCTGAACAGCCTCTACCGTCTGTCCAAGGGCGTGCTGCGTGGCATGACCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATGGGTAACGCAGGTCAGGCCAACTACGCGGCTGCCAAGGCCGGt                                                                  <  2:741807/150‑1 (MQ=255)
                                                                                     acCAACCTGAACAGCCTCTACCGTCTGTCCAAGGGCGTGCTGCGTGGCATGACCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATGGGTAACGCAGGTCAGGCCAACTACGCGGCTGCCAAGGCCGGTCTg                                                               >  2:800330/1‑150 (MQ=255)
                                                                                        aaCCTGAACAGCCTCTACCGTCTGTCCAAGGGCGTGCTGCGTGGCATGACCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATGGGTAACGCAGGTCAGGCCAACTACGCGGCTGCCAAGGCCGGTCTGGaa                                                            >  1:1082956/1‑150 (MQ=255)
                                                                                         aCCTGAACAGCCTCTACCGTCTGTCCAAGGGCGTGCTGCGTGGCATGACCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATGGGTAACGCAGGTCAGGCCAACTACGCGGCTGCCAAGGCCGGTCTGGAAg                                                           <  2:1118607/150‑1 (MQ=255)
                                                                                         aCCTGAACAGCCTCAACCGTCTGCCCAAGCGCGTGCTGAGTGGCATGACGAAGCCGCGTTGCGGTCGTATGATCAGCCGCGGTCAGGTCGGTGGTGCAACGGGTAACGCAGGTCAAGCCAACTACGCGGCTGCCAAGGCCAGTCTGGAAg                                                           <  2:240704/150‑1 (MQ=255)
                                                                                          ccTGAACAGCCTCTACCGTCTGTCCAAGGGCGTGCTGCGTGGCATGACCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATGGGTAACGCAGGTCAGGCCAACTACGCGGCTGCCAAGGCCGGTCTGGAAgg                                                          <  1:405859/150‑1 (MQ=255)
                                                                                                   ccTCTACCGTCTGTCCAAGGGCGTGCTGCGTGGCATGACCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATGGGTAACGCAGGTCAGGCCAACTACGCGGCTGCCAAGGCCGGTCTGGGCGGTTGCAGCcg                                                 >  2:732719/1‑150 (MQ=255)
                                                                                                     tctACCGTCTGTCCAAGGGCGTGCTGCGTGGCATGACCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATGGGTAACGCAGGTCAGGCCAACTACGCGGCTGCCAAGGCCGGTCTGGAAGGTTTCAGCcgcg                                               >  1:712615/1‑150 (MQ=255)
                                                                                                       tACCGTCTGTCCAAGGGCGTGCTGCGTGGCATGACCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATGGGTAACGCAGGTCAGGCCAACTACGCGGCTGCCAAGGCCGGTCTGGAAGGTTTCAGCCGCGcc                                             >  2:1256095/1‑150 (MQ=255)
                                                                                                                   aaGGGCGTGCTGCGTGGCATGACCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATGGGTAACGCAGGTCAGGCCAACTACGCGGCTGc                                                                            >  2:1199489/1‑107 (MQ=255)
                                                                                                                   aaGGGCGTGCTGCGTGGCATGACCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATGGGTAACGCAGGTCAGGCCAACTACGCGGCTGc                                                                            <  1:1199489/107‑1 (MQ=255)
                                                                                                                   aaGGGCGTGCTGCGTGGCATGACCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATGGGTAACGCAGGTCAGGCCAACTACGCGGCTGCCAAGGCCGGTCTGGAAGGTTTCAGCCGCGCCCTGGCGCGTGaa                                 >  2:151985/1‑150 (MQ=255)
                                                                                                                       gCGTGCTGCGTGGCATGACCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATGGGTAACGCAGGTCAGGCCAACTACGCGGCTGCCAAGGCCGGTCTGGAAGGTTTCAGCCGCGCCCTGGCGCGTGAAGTgg                             <  1:401013/150‑1 (MQ=255)
                                                                                                                       gCGTGCTGCGTGGCATGACCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATGGGTAACGCAGGTCAGGCCAACTACGCGGCTGCCAAGGCCGGTCTGGAAGGTTTCAGCCGCGCCCTGGCGCGTGAAGTgg                             <  2:712615/150‑1 (MQ=255)
                                                                                                                         gTGCTGCGTGGCATGACCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATGGGTAACGCAGGTCAGGCCAACTACGCGGCTGCCAAGGCCGGTCTGGAAGGTTTCAGCCGCGCCCTGGCGCGTGAAGTGGGt                           >  2:569192/1‑150 (MQ=255)
                                                                                                                             tgcGTGGCATGACCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTg                                                                                                                    >  1:680639/1‑57 (MQ=255)
                                                                                                                             tgcGTGGCATGACCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTg                                                                                                                    <  2:680639/57‑1 (MQ=255)
                                                                                                                             tgcGTGGCATGACCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATGGGTAACGCAGGTCAGGCCAACTACGCGGCTGCCAAGGCCGGTCTGGAAGGTTTCAGCCGCGCCCTGGCGCGTGAAGTGGGTTcgc                       >  1:340036/1‑150 (MQ=255)
                                                                                                                                gTGGCATGACCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATGGGTAACGCAGGTCAGGCCAACTACGCGGCTGCCAAGGCCGGTCTGGAAGGTTTCAGCCGCGCCCTGGCGCGTGAAGTGGGTTCGCGTg                    <  2:1082956/150‑1 (MQ=255)
                                                                                                                                     aTGACCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATGGGTAACGCAGGTCAGGCCAACTACGCGGCTGCCAAGGCCGGTCTGGAAGGTTTCAGCCGCGCCCTGGCGCGTGAAGTGGGTTCGCGTGGTATc               >  2:1275464/1‑150 (MQ=255)
                                                                                                                                     aTGACCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATGGGTAACGCAGGTCAGGCCAACTACGCGGCTGCCAAGGCCGGTCTGGAAGGTTTCAGCCGCGCCCTGGCGCGTGAAGTGGGTTCGCGTGGTATc               >  1:934854/1‑150 (MQ=255)
                                                                                                                                     aTGACCAAGGCGCATTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATGGGTAACGCAGGTCAGGCCAACTACGCGGCt                                                                              >  2:1274050/1‑87 (MQ=255)
                                                                                                                                     aTGACCAAGGCGCATTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATGGGTAACGCAGGTCAGGCCAACTACGCGGCt                                                                              <  1:1274050/87‑1 (MQ=255)
                                                                                                                                     aTGACCAAGGCGCATTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATGGGTAACGCAGGTCAGGCCAACTACGCGGCt                                                                              <  1:1274357/87‑1 (MQ=255)
                                                                                                                                     aTGACCAAGGCGCATTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATGGGTAACGCAGGTCAGGCCAACTACGCGGCt                                                                              >  2:1274357/1‑87 (MQ=255)
                                                                                                                                        aCCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATGGGTAACGCAGGTCAGGCCAACTACGCGGCTGCCAAGGCCGGTCTGGAAGGTTTCAGCCGCGCCCTGGCGCGTGAAGTGGGTTCGCGTGGTATCAcc            <  1:118923/150‑1 (MQ=255)
                                                                                                                                            aGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATGGGTAACGCAGGTCAGGCCAACTACGCGGCTGCCAAGGCCGGTCTGGAAGGTTTCAGCCGCGCCCTGGCGCGTGAAGTGGGTTCGCGt                     >  1:482664/1‑137 (MQ=255)
                                                                                                                                            aGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATGGGTAACGCAGGTCAGGCCAACTACGCGGCTGCCAAGGCCGGTCTGGAAGGTTTCAGCCGCGCCCTGGCGCGTGAAGTGGGTTCGCGt                     <  2:482664/137‑1 (MQ=255)
                                                                                                                                                 cgTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATGGGTCACGCAGGTCAGGCCAACTACGCGGCTGCCAAGGCCGGTCTGGAAGGTTTCAGCCGCGCCCTGGCGCGTGAAGTGGGTTCGCGTGGTATCACCCTCGACTCg   <  2:196370/150‑1 (MQ=255)
                                                                                                                                                  gTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATTGGTAACGCAGGTCAGGCCAACTACGCGGCTGCCAAGGCCGGTCTGGAAGGTTTCAGCCGCGCCCTGGCGCGTGAAGTGGGTTCGCGTGGTATCACCGTCAACTCgg  <  2:64235/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GGCCATCTTGGTCAACAATGCCGGTATCACCCGCGACAACCTCATGCTGCGCATGAAAGACGACGAGTGGTTTGATGTCATCGACACCAACCTGAACAGCCTCTACCGTCTGTCCAAGGGCGTGCTGCGTGGCATGACCAAGGCGCGTTGGGGTCGTATCATCAGCATCGGCTCGGTCGTTGGTGCCATGGGTAACGCAGGTCAGGCCAACTACGCGGCTGCCAAGGCCGGTCTGGAAGGTTTCAGCCGCGCCCTGGCGCGTGAAGTGGGTTCGCGTGGTATCACCGTCAACTCGG  >  NC_002947/2157942‑2158237

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