Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 657,102 G→T M182I (ATG→ATT pagP → Lipid IVA palmitoyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913657,1020GT97.9% 150.6 / ‑2.8 47M182I (ATG→ATTpagPLipid IVA palmitoyltransferase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (26/20);  total (27/20)
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.94e-01

ACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGT  >  NC_000913/656978‑657237
                                                                                                                            |                                                                                                                                       
aCTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATTCGCTTTCAGTTTTga                                                                                                                          <  1:938560/140‑1 (MQ=255)
aCTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATTCGCTTTCAGTTTTga                                                                                                                          >  1:327181/1‑140 (MQ=255)
  tGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTTCCTACAACAATGGCAATGTGTACTTTGCCTGGATTCGCTTTCAGTTTTgaga                                                                                                                        <  2:575528/140‑1 (MQ=255)
        tACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATTCGCTTTCAGTTTTGAGACAAATg                                                                                                                  <  1:12450/140‑1 (MQ=255)
            tCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATTCGCTTTCAGTTTTGAGACAAATGAAGt                                                                                                              <  1:926485/140‑1 (MQ=255)
               ctctCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATTCGCTTTCAGTTTTGAGACAAATGAAGtttt                                                                                                           >  2:885968/1‑140 (MQ=255)
               ctctCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATTCGCTTTCAGTTTTGAGACAAATGAAGtttt                                                                                                           >  2:572781/1‑140 (MQ=255)
               ctctCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATTCGCTTTCAGTTTTGAGACAAATGAAGtttt                                                                                                           >  2:146268/1‑140 (MQ=255)
                  tcCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGt                                                                                                        <  2:957538/140‑1 (MQ=255)
                    ccGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTaa                                                                                                      >  1:812537/1‑140 (MQ=255)
                                  aTTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCa                                                                                        >  1:399259/1‑140 (MQ=255)
                                  aTTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCa                                                                                        <  2:877450/140‑1 (MQ=255)
                                      gCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAg                                                                                    <  1:146268/140‑1 (MQ=255)
                                      gCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAg                                                                                    <  1:230569/140‑1 (MQ=255)
                                           cGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTaaa                                                                               >  2:94924/1‑140 (MQ=255)
                                           cGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTaaa                                                                               >  2:260394/1‑140 (MQ=255)
                                           cGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGAACTTTGGCTGGATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTaaa                                                                               >  2:948951/1‑140 (MQ=255)
                                                gTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATTCGCTTTCAGTTTCGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAg                                                                          >  2:659245/1‑140 (MQ=255)
                                                 ttATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTTTTTAAAATCAATAGCTAAAATAAGt                                                                         <  2:688494/140‑1 (MQ=255)
                                                  tATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTa                                                                        <  1:408014/140‑1 (MQ=255)
                                                     ggCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACa                                                                     >  1:146567/1‑140 (MQ=255)
                                                          aGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCaaa                                                                >  2:702955/1‑140 (MQ=255)
                                                           gTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCaaaa                                                               >  2:29460/1‑140 (MQ=255)
                                                                 tttCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAAcg                                                         <  1:29460/140‑1 (MQ=255)
                                                                                    cGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTa                                      >  1:797133/1‑140 (MQ=255)
                                                                                     gggTACCTACAACAATGGCAATGTGTACTTTGCCTGGATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAg                                     <  2:522833/140‑1 (MQ=255)
                                                                                      ggTACCTACAACAATGGCAATGTGTACTTTGCCTGGATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGt                                    <  2:903118/140‑1 (MQ=255)
                                                                                      ggTACCTACAACAATGGCAATGTGTACTTTGCCTGGATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGt                                    >  2:311701/1‑140 (MQ=255)
                                                                                          ccTACAACAATGGCAATGTGTACTTTGCCTGGATTCGCTTTCAGTTTTGATACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAg                                >  2:213768/1‑140 (MQ=255)
                                                                                              caacaaTGGCAATGTGTACTTTGCCTGGATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTac                            >  1:57023/1‑140 (MQ=255)
                                                                                               aacaaTGGCAATGTGTACTTTGCCTGGATTCGCTTTCAGTTTTGATACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTaca                           <  1:213768/140‑1 (MQ=255)
                                                                                                     ggCAATGTGTACTTTGCCTGGATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGa                      <  2:640618/139‑1 (MQ=255)
                                                                                                     ggCAATGTGTACTTTGCCTGGATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGa                      >  1:640618/1‑139 (MQ=255)
                                                                                                     ggCAATGTGTACTTTGCCTGGATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGAc                     >  2:859240/1‑140 (MQ=255)
                                                                                                      gCAATGTGTACTTTGCCTGGATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACa                    <  1:311701/140‑1 (MQ=255)
                                                                                                        aaTGTGTACTTTGCCTGGATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACaaa                  <  2:57023/140‑1 (MQ=255)
                                                                                                         aTGTGTACTTTGCCTGGATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAg                 <  1:859240/140‑1 (MQ=255)
                                                                                                            tgtACTTTGCCTGGATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGt              >  1:507655/1‑140 (MQ=255)
                                                                                                               aCTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTAcc           >  1:470850/1‑140 (MQ=255)
                                                                                                                cTTTGCCTGGATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAGAATCAATAGCTAAAACAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGCCAAAg                 <  2:850005/133‑1 (MQ=255)
                                                                                                                cTTTGCCTGGATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAGAATCAATAGCTAAAACAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAg                 >  1:850005/1‑133 (MQ=255)
                                                                                                                    gCCTGGATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAAt      <  2:797133/140‑1 (MQ=255)
                                                                                                                     ccTGGATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAAtt     >  1:331962/1‑140 (MQ=255)
                                                                                                                     ccTGGATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAAtt     >  2:815098/1‑140 (MQ=255)
                                                                                                                     ccTGGATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAAtt     >  1:664216/1‑140 (MQ=255)
                                                                                                                     ccTGGATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAAtt     >  1:20736/1‑140 (MQ=255)
                                                                                                                        ggATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGgt  <  1:527521/140‑1 (MQ=255)
                                                                                                                        ggATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGgt  >  1:310123/1‑140 (MQ=255)
                                                                                                                            |                                                                                                                                       
ACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGT  >  NC_000913/656978‑657237

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