Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,939,673 T→C N174S (AAC→AGC)  lpxM ← myristoyl‑acyl carrier protein (ACP)‑dependent acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,939,6730TC96.1% 167.1 / ‑3.3 51N174S (AAC→AGC) lpxMmyristoyl‑acyl carrier protein (ACP)‑dependent acyltransferase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (24/25);  total (25/26)
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.94e-01

CTGTGCTCTGGGCCATGATCCTGATCGGGTAAATAATATCCCCAGTACCCCTGACGTACCGACTGGATGAATGGTTTAATACCGTCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGTTCCAGACATAATCAAAAACCGGGTTGCCCTGATTATGGAACATCGCTGCCATTTTCTGCCCTTGCGAGGCCATCAGCATGGCAGGAATATCGACGGCCCAACCGTGCGGCACCAGAAAGA  >  NC_000913/1939544‑1939793
                                                                                                                                 |                                                                                                                        
cTGTGCTCTGGGCCATGATCCTGATCGGGTAAATAATATCCCCAGTACCCCTGACGTACCGACTGGATGAATGGTTTAATACCGTCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACAt                                                                                                                 <  1:435305/139‑1 (MQ=255)
     ctctGGGCCATGATCCTGATCGGGTAAATAATATCCCCAGTACCCCTGACGTACCGACTGGATGAATGGTTTAATACCGTCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCa                                                                                                            >  1:406238/1‑139 (MQ=255)
               tGATCCTGATCGGGTAAATAATATCCCCAGTACCCCTGACGTACCGACTGTATGAATGGTTTAATACCGTCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCggg                                                                                                   <  2:1025462/138‑1 (MQ=255)
               tGATCCTGATCGGGTAAATAATATCCCCAGTACCCCTGACGTACCGACTGGATGAATGGTTTAATACCGTCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGTTCCAGACATAATCAAAAACCGGGt                                                                                                  >  1:940882/1‑139 (MQ=255)
               tGATCCTGATCGGGTAAATAATATCCCCAGTACCCCTGACGTACCGACTGGATGAATGGTTTAATACCGTCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCggg                                                                                                   >  1:1025462/1‑138 (MQ=255)
               tGATCCTGATCGGGTAAATAATATCCCCAGTACCCCTGACGTACCGACTGGATGAATGGTTTAATACCGTCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGGGt                                                                                                  <  2:406238/139‑1 (MQ=255)
                   ccTGATCGGGTAAATAATATCCCCAGTACCCCTGACGTACCGACTGGATGAATGGTTTAATACCGTCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGGGTTGcc                                                                                              <  2:902163/139‑1 (MQ=255)
                    cTGATCGGGTAAATAATATCACCAGTACCCCTGACGTACCGACTGGATGAATGGTTTAATACCGTCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGAGTTAccc                                                                                             >  1:861154/1‑139 (MQ=255)
                      gATCGGGTAAATAATATCCCCAGTACCCCTGACGTACCGACTGGATGAATGGTTTAATACCGTCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGGGTTGCCCTg                                                                                           <  1:1119688/139‑1 (MQ=255)
                           ggTAAATAATATCCCCAGTACCCCTGACGTACCGACTGGATGAATGGTTTAATACCGTCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGGGTTGcc                                                                                              >  1:1053623/1‑131 (MQ=255)
                           ggTAAATAATATCCCCAGTACCCCTGACGTACCGACTGGATGAATGGTTTAATACCGTCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGGGTTGcc                                                                                              <  2:1053623/131‑1 (MQ=255)
                                    tatCCCCAGTACCCCTGACGTACCGACTGGATGAATGGTTTAATACCGTCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGTTCCAGACATAATCAAAAACCGGGTTGCCCTGATTATGGAACATCg                                                                             <  2:940882/139‑1 (MQ=255)
                                    tatCCCCAGTACCCCTGACGTACCGACTGGATGAATGGTTTAATACCGTCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGGGTTGCCCTGATTATGGAACATCg                                                                             <  2:1152738/139‑1 (MQ=255)
                                          cAGTACCCCTGACGTACCGACTGGATGAATGGTTTAATACCGTCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGGGTTGCCCTGATTATGGAACATCGCTGCCa                                                                       >  1:876887/1‑139 (MQ=255)
                                           aGTACCCCTGACGTACCGACTGGATGAATGGTTTAATACCGTCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGGGTTGCCCTGATTATGGAACATCGCTGCCAt                                                                      >  2:1265182/1‑139 (MQ=255)
                                            gTACCCCTGACGTACCGACTGGATGAATGGTTTAATACCGTCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGCTACAGACATAATCAAAAACCGGGTTGCCCTGATTATGGAACATCGCTGCCAtt                                                                     >  1:201900/1‑139 (MQ=255)
                                                  cTGACGTACCGACTGGATGAATGGTTTAATACCGTCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACTAATCAAAAACCGGGTTGCCCTGATTATGGAACATCGCTGCCATTTTCTGcc                                                              <  1:369748/139‑1 (MQ=255)
                                                   tGACGTACCGACTGGATGAATGGTTTAATACCGTCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGGGTTGCCCTGATTATGGAACATCGCTGCCATTTTCTGcc                                                              <  1:792187/139‑1 (MQ=255)
                                                     aCGTACCGACTGGATGAATGGTTTAATACCGTCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATc                                                                                                             >  2:1219319/1‑90 (MQ=255)
                                                     aCGTACCGACTGGATGAATGGTTTAATACCGTCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATc                                                                                                             <  1:1219319/90‑1 (MQ=255)
                                                         aCCGACTGGATGAATGGTTTAATACCGTCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGa                                                                                                                    <  1:757456/79‑1 (MQ=255)
                                                         aCCGACTGGATGAATGGTTTAATACCGTCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGa                                                                                                                    >  2:757456/1‑79 (MQ=255)
                                                         aCCGACTGGATGAATGGTTTAATACCGTCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGGGTTGCCCt                                                                                            >  1:1580083/1‑103 (MQ=255)
                                                         aCCGACTGGATGAATGGTTTAATACCGTCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGGGTTGCCCt                                                                                            <  2:1580083/103‑1 (MQ=255)
                                                         aCCGACTGGATGAATGGTTTAATACCGTCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGGGTTGCCCTGATTATGGAACATCGCTGCCATTTTCTGCCCTTGCg                                                        <  1:1265182/139‑1 (MQ=255)
                                                           cGACTGGATGAATGGTTTAATACCGTCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGGGTTGCCCTGATTATGGAACATCGCTGCCATTTTCTGCCCTTGCGAg                                                      <  2:876887/139‑1 (MQ=255)
                                                                  aTGAATGGTTTAATACCGTCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGGGTTGcc                                                                                              <  2:678904/92‑1 (MQ=255)
                                                                  aTGAATGGTTTAATACCGTCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGGGTTGcc                                                                                              >  1:678904/1‑92 (MQ=255)
                                                                    gAATGGTTTAATACCGTCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGGGTTGCCCTGATTATGGAACATCGCTGCCATTTTCTGCCCTTGCGAGGCCATCAGc                                             <  2:861154/139‑1 (MQ=255)
                                                                            tAATACCGTCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGGGTTGCCCTGATTATGGAACATCGCTGCCATTTTCTg                                                                >  2:1086811/1‑112 (MQ=255)
                                                                            tAATACCGTCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGGGTTGCCCTGATTATGGAACATCGCTGCCATTTTCTg                                                                <  1:1086811/112‑1 (MQ=255)
                                                                                    tCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGGGTTGCCCTGGTTATGGAACATCGCTGCCATTTTCTGCCCTTGCGAGGCCATCAGCATGGCAGGAATATCGa                             >  2:1447921/1‑139 (MQ=255)
                                                                                    tCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGGGTTGCCCTGATTATGGAACATCGCTGc                                                                         >  1:1150722/1‑95 (MQ=255)
                                                                                    tCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGGGTTGCCCTGATTATGGAACATCGCTGc                                                                         <  2:1150722/95‑1 (MQ=255)
                                                                                    tCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGGGTTGCCCTGATTATGGAACATCGCTGCCATTTTCTGCCCTTGCGAGGCCATCAGCATGGCAGGAATATCGa                             >  2:1309514/1‑139 (MQ=255)
                                                                                    tCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGGGTTGCCCTGATTATGGAACATCGCTGCCATTTTCTGCCCTTGCGAGGCCATCAGCATGGCAGGAATATCGa                             <  2:201900/139‑1 (MQ=255)
                                                                                                tGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCggg                                                                                                   >  1:1334279/1‑57 (MQ=255)
                                                                                                tGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCggg                                                                                                   <  2:1334279/57‑1 (MQ=255)
                                                                                                tGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGGGTTGCCCTGATTATGGAACATCGCTGCCATTTTCTGCCCTTGCg                                                        >  2:127357/1‑100 (MQ=255)
                                                                                                tGCAGACGACCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGGGTTGCCCTGATTATGGAACATCGCTGCCATTTTCTGCCCTTGCg                                                        <  1:127357/100‑1 (MQ=255)
                                                                                                      cgacCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGGGTTGCCCTGATTATGGAACATCGCTGCCATTTTCTGCCCTTGCGAGGCCATCAGCATGGCAGGAATATCGACGGCCCAACCGTGCGGCa           >  2:1569027/1‑139 (MQ=255)
                                                                                                       gacCGCCAACGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGGGTTGCCCTGATTATGGAACATCGCTGCCATTTTCTGCCCTTGCGAGGCCATCAGCATGGCAGGAATATCGACGGCCCAACCGTGCGGCAc          >  2:285447/1‑139 (MQ=255)
                                                                                                        acCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGGGTTGcc                                                                                              >  2:652822/1‑54 (MQ=255)
                                                                                                        acCGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGGGTTGcc                                                                                              <  1:652822/54‑1 (MQ=255)
                                                                                                         ccGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCggg                                                                                                   <  1:912179/48‑1 (MQ=255)
                                                                                                         ccGCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCggg                                                                                                   >  2:912179/1‑48 (MQ=255)
                                                                                                           gCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGGGt                                                                                                  <  2:378231/47‑1 (MQ=255)
                                                                                                           gCCAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGGGt                                                                                                  >  1:378231/1‑47 (MQ=255)
                                                                                                            ccAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGGGTTGCCCTGATTATGGAACATCGCTGCCATTTTCTGCCCTTGCGAg                                                      <  2:1264399/90‑1 (MQ=255)
                                                                                                            ccAAAGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGGGTTGCCCTGATTATGGAACATCGCTGCCATTTTCTGCCCTTGCGAg                                                      >  1:1264399/1‑90 (MQ=255)
                                                                                                               aaGCGACGACGCACCGTGCTCCAGACATAATCAAAAACCGGGTTGCCCTGATTATGGAACATCGCTGCCATTTTCTGCCCTTGCGAGGCCATCAGCATGGCAGGAATATCGACGGCCCAACCGTGCGGCACCAGAAAGa  >  1:106958/1‑139 (MQ=255)
                                                                                                                                 |                                                                                                                        
CTGTGCTCTGGGCCATGATCCTGATCGGGTAAATAATATCCCCAGTACCCCTGACGTACCGACTGGATGAATGGTTTAATACCGTCATTTCTCGCATGCAGACGACCGCCAAAGCGACGACGCACCGTGTTCCAGACATAATCAAAAACCGGGTTGCCCTGATTATGGAACATCGCTGCCATTTTCTGCCCTTGCGAGGCCATCAGCATGGCAGGAATATCGACGGCCCAACCGTGCGGCACCAGAAAGA  >  NC_000913/1939544‑1939793

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